Re: [問題] stata 模擬

[bcs (= ="frailty..gggg XD)]

: program define normsim,rclass
: version 10.1
: syntax[,obs(integer 1) mu(real 0) sigma(real 1)]
: drop_all
: set obs `obs'
: tempvar z
: gen `z'=`mu'+`sigma'*invnorm(uniform())
:-> sum `z',d
:-> return scalar med=r(p50)
: end

:-> simulate x=r(med),reps(1000):normsim,obs(1000) mu(0) sigma(1)


改 '->' 即可以跑出來
 good luck :)

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.112.86.139
推 Lucas5566:謝謝妳!我會了!!