作者Lucas5566 (5566)
看板Statistics
标题[问题] stata 模拟
时间Mon Apr 27 10:30:31 2009
我想要模拟某个分配的中位数的抽样分配
以常态(0,1)来说
我写下面的程式
program define normsim,rclass
syntax[,obs(integer 1) mu(real 0) sigma(real 1)]
drop_all
set obs `obs'
tempvar z
gen `z'=`mu'+`sigma'*invnorm(uniform())
问题在这边，我不知道怎麽求出中位数
我本来是用
gen median=median(`z')
return scalar M=median
end
然後
simulate M=median,reps(1000):normsim,obs(1000) mu(0) sigma(1)
可是它执行是说我的syntax有问题
我不知道我错出在哪
很有可能是我return中位数的地方有问题(我觉得应该是= =+)
但是我不知道要怎麽算中位数(codebook指令可以有50%的percentile，但
是不知道怎麽return)
请高手帮忙
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 218.168.201.7
※ 编辑: Lucas5566 来自: 218.168.201.7 (04/27 16:47)