作者mathematics3 (馬克勞林)
站內Statistics
標題[問題] 請問STATA 驗證大數法則的程式
時間Tue Dec 2 00:07:28 2008
forvalues i = 1(1)100{
set obs 3000
gen x`i' = (invnorm(uniform()))*sqrt(0.5)+3
qui su x`i'
display r(mean)
}
不知道為什麼
跑不出來...
是我迴圈的語法寫錯嗎?
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.114.208.22
1F：→ bcs:LLN樣本愈大，均數會愈接近參數,obs可設為不同樣本數. 12/02 12:14
2F：→ bcs:mata 裡面先設立100*1的matrix，J(100,1,0) 12/02 12:17
3F：→ bcs:再利用loop的方式塞入mean(rnormal(size,1,3,0.5)) 12/02 12:18
4F：→ bcs:最後把mata的matrix轉到stata，再用sum即可查驗mean值 12/02 12:19
5F：→ mathematics3:謝謝樓上..可惜我學的STATA很粗淺 matrix也還沒學到 12/02 12:49