1

ブートストラップ回帰モデルの係数に沿って一連のスカラーを格納しようとしています。以下のコードは、Stata [P]rogramming manual for の例のように見えますが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 に答える