Re: [問題] stata 模擬
看板Statistics作者bcs (= ="frailty..gggg XD)時間17年前 (2009/04/27 18:36)推噓1(1推 0噓 0→)留言1則, 1人參與討論串2/2 (看更多)
: 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
推
04/27 18:58, , 1F
04/27 18:58, 1F
討論串 (同標題文章)
完整討論串 (本文為第 2 之 2 篇)：