ywh19860616 发表于 2014-6-7 10:10
上面是十等分分组统计均值。
谢谢您的解答,我根据上面的程序和数据情况,修改了obs以及recode里面的数据,但是结果不太对,能否帮忙看看问题出在哪里?clear
set obs 330820
gen lntfp = uniform()
gen eng = rnormal()
cumul eng,gen(index) eq
recode index (min/5=1)(5/5.5=2)(5.5/6=3)(6/6.5=4)(6.5/7=5)(7/7.5=6)(7.5/8=7)(8/8.5=8)(8.5/9=9)(9/9.5=10)(9.5/10=11)(10/10.5=12)(10.5/11=13)(11/11.5=14)(11.5/12=15)(12/12.5=16)(12.5/13=17)(13/13.5=18)(13.5/14=19)(14/14.5=20)(14.5/15=21)(15/15.5=22)(15.5/16=23)(16/16.5=24)(16.5/17=25)(17/17.5=26)(17.5/18=27)(18/max=28)
statsby mean=r(mean), clear by(index): summarize lntfp