模拟 Stata 8 集群自举回归答案

【问题标题】：Emulate Stata 8 clustered bootstrapped regression模拟 Stata 8 集群自举回归
【发布时间】：2013-10-29 22:05:48
【问题描述】：

我正在尝试沿自举回归模型的系数存储一系列标量。下面的代码类似于 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+ 中执行？问题已提出in the past 和more recently，但没有找到答案。

对不起，描述太长了，这是一个很长的问题。

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

【问题讨论】：

请注意，对于不是您真正关心的最小二乘法，预测的平均值是样本平均值。

标签： stata

【解决方案1】：

将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

【讨论】：