1

我正在尝试沿自举回归模型的系数存储一系列标量。下面的代码看起来像 Stata [P]rogramming 手册中的示例，postfile显然是为了与此类程序一起使用。

问题出// commented在线路上，无法正常工作。更具体地说，问题似乎是下面的语法在 Stata 8 中有效，但在bootstrap过程中进行了一些更改后无法在 Stata 9+ 中工作。

cap pr drop bsreg
pr de bsreg
    reg mpg weight gear_ratio

    predict yhat
    qui sum yhat

//  sca mu = r(mean)
//  post sim (mu)
end

sysuse auto, clear

postfile sim mu using results , replace

    bootstrap, cluster(foreign) reps(5) seed(6112): bsreg

postclose sim

use results, clear

添加version 8到代码中并没有解决问题。有谁知道这个过程有什么问题，以及如何修复它以在 Stata 9+ 中执行？过去和最近都提出过这个问题，但没有找到答案。

很抱歉描述太长，这是一个很长的问题。

我把这个问题当作一个编程问题来提出，因为我正在使用这段代码来复制一些健康不平等研究。有必要引导整个过程，而不仅仅是reg模型。我对这种方法有一些疑问，但没有什么能阻止我复制分析。

4

1 回答 1

0

添加noisily到引导程序显示该predict命令存在问题。这是使用tempvar宏的修复。

cap pr drop bsreg
pr de bsreg
    reg mpg weight gear_ratio

    tempvar yhat
    predict `yhat'
    qui sum `yhat'

    sca mu = r(mean)
    post sim (mu)
end

sysuse auto, clear

postfile sim mu using results , replace

    bootstrap, cluster(foreign) reps(5) seed(6112): bsreg

postclose sim

use results, clear
于 2013-10-30T20:07:35.877 回答