Question

我正在尝试沿着自举回归模型的系数存储一系列标量。下面的代码类似于postfile的Stata [P]编程手册中的示例，该手册显然是用于此类程序。

问题在于// 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模型。我对这种方法有一些狡辩，但没有什么可以阻止我复制分析。

Answer 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

模拟Stata 8群集自举回归

1 个答案: