求救STATA做Bootstrap方法-需要用到两个回归方程的系数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

正在学习用STATA做Bootstrap的方法，，如果我想要通过Bootstrap得到一个统计值的BC置信区间，而这个统计值是通过两个回归方程得到的，请问应该如何做？

比如：第一个回归： PM=a0+a1SL+a2SR
         第二个回归：OC=b0+b1SR+b2PM+b3SL+b4SR*SL+b5SR*PM

我想要求a1*b5的BC置信区间，应该如何做？

万分感谢！！！

扫码加我 拉你入群

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

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

分享0 收藏11 回帖

关键词：Bootstrap方法 Bootstrap Bootstra Stata boots 置信区间 如何 统计

请各位高手踊跃帮忙啊~~

.....

等待回答。。。。

帮忙啊～～～

求救！！

I guess similar to this one

http://www.ats.ucla.edu/stat/stata/faq/mulmediation.htm

capture drop program bootmm
program bootmm, rclass
  syntax [if] [in]
  sureg (read math)(write math)(science read write math) `if' `in'
  return scalar indread  = [read]_b[math]*[science]_b[read]
  return scalar indwrite = [write]_b[math]*[science]_b[write]
  return scalar indtotal = [read]_b[math]*[science]_b[read]+ ///
                           [write]_b[math]*[science]_b[write]
end

bootstrap r(indread) r(indwrite) r(indtotal), bca reps(5000) nodots: bootmm

in your case

capture drop program bootaaa
program bootaaa, rclass
sureg ( PM SL SR) (OC SR PM SL SR*SL SR*PM) `if' `in'
return scalar aaa  = [PM]_b[SL]*[OC]_b[SR*PM]
end
bootstrap r(aaa) , bca reps(5000) nodots: bootaaa

sophiafinn 发表于 2013-2-7 08:25
in your case

capture drop program bootaaa

十分感谢！！我自己做了一下，应该就是这个样子了。不过在sureg语句上面需要再加上syntax [if][in] 否则运行结果就会显示insufficient observations to compute jackknife standard errors no results will be saved。

再次感谢！！