[問題] stata 模擬
我想要模擬某個分配的中位數的抽樣分配
以常態(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)
討論串 (同標題文章)
完整討論串 (本文為第 1 之 2 篇):