蒙特卡罗模拟

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

最近在看蒙特卡罗模拟，对照书本练习，出现以下结果，不知问题在哪里（红色部分），求大神指点！
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

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：蒙特卡罗模拟 蒙特卡罗 蒙特卡 simulate function 蒙特卡罗 function occurred version double

resch2应该是Simulation卡方的一个外部命令，你加载下。

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不是外部命令，不知我的理解对不对，还麻烦您给予解答。（我对照其他书本练习也是这种情况，不知何故？？）

大念 发表于 2014-7-15 10:00
非常感谢您的解答，但书上（周广肃，Stata统计分析与应用）说，一开始定义了程序reschi2(program reschi2 ...

如果是自己定义的程序的话，可能是是这个程序定义的问题了。你看看程序的定义是不是出问题了，Stata我记得定义程序在10以后有一些区别。

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

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);

crystal8832 发表于 2014-7-15 13:15
如果是自己定义的程序的话，可能是是这个程序定义的问题了。你看看程序的定义是不是出问题了，Stata我记得 ...

谢谢您的回复！我仔细看了一下，程序定义没问题，用其他例子做还是这种情况，很是困惑！

大念 发表于 2014-7-15 14:47
谢谢您的回复！我仔细看了一下，程序定义没问题，用其他例子做还是这种情况，很是困惑！

我帮你跑了一下，我这边是可以做的哦。

crystal8832 发表于 2014-7-15 23:03
我帮你跑了一下，我这边是可以做的哦。

我重新做还是那种情况，唉！不过还是要感谢您的热心！

大念 发表于 2014-7-15 14:43
感谢您的及时回复，我仔细看了书本，程序定义没问题，我再对照另一本书（陈强，高级计量经济学及stata应用 ...

你好，我之前也碰到程序无法运行的问题，你可以把程序写到do文档里，再运行，就能成功了。我试过可以。