作者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