模擬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