模拟Stata 8集群引导回归

Question

我试图沿着引导回归模型的系数存储一系列标量。下面的代码看起来像来自Stata [P]程序手册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模型。我对这种方法有些狡猾，但没有任何东西阻止我复制分析。

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