Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
From | Marta García-Granero <mgarciagranero@gmail.com> |
To | statalist@hsphsun2.harvard.edu |
Subject | st: Adding "if group==1" option to a short program -- or ",by(group)" |
Date | Tue, 02 Jul 2013 09:51:05 +0200 |
Good morning. I have written this short program: program skr2 version 12.1 args x confirm variable `x' _pctile `x', nquantiles(8) local sk2=(r(r6)+r(r2)-2*r(r4))/(r(r6)-r(r2)) local kr2=(r(r7)-r(r5)+r(r3)-r(r1))/(r(r6)-r(r2)) - 1.2333 display "Bowley's skewness [-1; 1]: " %6.2f `sk2' display "Moor's kurtosis [-1.23; inf.]: " %6.2f `kr2' endAs it can be seen from its contents, it computes robust skewness and kurtosis coefficients for just one variable.
I promise I have tried to read the documentation on programming, but I couldn't find how to modify it to be able to accept the selection of subgroups based on a second variable, like:
skr2 iron if group==1 or maybe skr2 iron, by(group)Any hints would be very much appreciated. I don't need the complete solution (I'd like to face the task and learn more about Stata programming), just a gentle push in the correct direction.
Thanks a lot, Marta GG * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/faqs/resources/statalist-faq/ * http://www.ats.ucla.edu/stat/stata/