全部版块 我的主页
› 论坛 › 计量经济学与统计论坛 五区 › 计量经济学与统计软件 › Stata专版
蒙特卡罗模拟
大念
21146 21
2014-07-14
最近在看蒙特卡罗模拟，对照书本练习，出现以下结果，不知问题在哪里（红色部分），求大神指点！
program reschi2,rclass
  1. version 10
  2. drop _all
  3. set obs 50
  4. gen double x=reschi2(1)
  5. gen y=3+2*x+rchi2(1)-3
  6. reg y x
  7. return scalar t2=(_b[x]-2)/_se[x]
  8. return scalar r2=abs(return(t2))>invttail(48,0.025)
  9. end

. simulate reject=r(r2),reps(1000) nodots seed(101):reschi2
Unknown function reschi2()
an error occurred when simulate executed reschi2


二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复
crystal8832
2014-7-14 23:54:46
resch2应该是Simulation卡方的一个外部命令，你加载下。
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

大念
2014-7-15 10:00:28
crystal8832 发表于 2014-7-14 23:54
resch2应该是Simulation卡方的一个外部命令，你加载下。
非常感谢您的解答，但书上（周广肃，Stata统计分析与应用）说，一开始定义了程序reschi2(program reschi2,rclass),后面simulate reject=r(r2),reps(1000) nodots seed(101):reschi2解释是（书上原文）：
这里，每次模拟的命令来自于程序"reschi2"，选项resps(1000)表明模拟次数为1000次，seed(101)表明设定种子为101，nodots表示不显示模拟过程的点。表达式“reject=r(r2)"表明将每次模拟的返回值r(r2)保存在reject中。因此，我认为它的意思是执行这个程序(即reschi2)1000次，reschi2不是外部命令，不知我的理解对不对，还麻烦您给予解答。（我对照其他书本练习也是这种情况，不知何故？？）
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

crystal8832
2014-7-15 13:15:16
大念 发表于 2014-7-15 10:00
非常感谢您的解答，但书上（周广肃，Stata统计分析与应用）说，一开始定义了程序reschi2(program reschi2 ...
如果是自己定义的程序的话，可能是是这个程序定义的问题了。你看看程序的定义是不是出问题了，Stata我记得定义程序在10以后有一些区别。
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

大念
2014-7-15 14:35:40
crystal8832 发表于 2014-7-15 13:15
如果是自己定义的程序的话，可能是是这个程序定义的问题了。你看看程序的定义是不是出问题了，Stata我记得 ...
感谢您的及时回复，我仔细看了一下，程序的定义没问题，我对照另一本书（陈强，高级计量经济学及stata应用）的例子练习，仍然出现这种情况，很是困惑！（是版本的原因吗？我用的是stata10，人家用stata9也没问题）
. program onesample,rclass
  1. drop _all
  2. set obs 30
  3. gen x=runiform()
  4. sum x
  5. return scalar mean_sample=r(mean)
  6. end

. simulate xbar=r(mean_sample),seed(101) reps(10000):onesample
Unknown function runiform()
an error occurred when simulate executed onesample
r(133);

end of do-file
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

大念
2014-7-15 14:43:55
crystal8832 发表于 2014-7-15 13:15
如果是自己定义的程序的话，可能是是这个程序定义的问题了。你看看程序的定义是不是出问题了，Stata我记得 ...
感谢您的及时回复，我仔细看了书本，程序定义没问题，我再对照另一本书（陈强，高级计量经济学及stata应用）的例子练习，仍然出现那种情况，很是困惑！
program onesample,rclass
  1. drop _all
  2. set obs 30
  3. gen x=runiform()
  4. sum x
  5. return scalar mean_sample=r(mean)
  6. end

. simulate xbar=r(mean_sample),seed(101) reps(10000):onesample
Unknown function runiform()
an error occurred when simulate executed onesample
r(133);

end of do-file

r(133);
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

微信QQ空间QQ微博
扫码加好友，拉您进群
各岗位、行业、专业交流群