2016-10-16 6 views
-1

私はStataのplogit adoを使って、回帰係数から回帰係数のベクトルを取得しました。 カラム名が存在する回帰係数の1 x 60ベクトルです。 私は、各係数を各被験者の変数値で乗算することによって線形予測子を得る必要があります。スタータの回帰係数の1つのベクトルを線形予測子に変換する

それは係数ベクトルを乗じてこれはかなりあいまいなユーザーが記述したコマンドであり、そしてそれがために参考にされているだろう被験者ごとに

+0

'predict'を使用できませんか? – chan1142

答えて

0

を線形予測を提供しますベクターに列名を抽出することが可能ですあなたはits murky originsを明らかにする。あなたはすぐに答えを得ているかもしれません!

それはそう言いませんが

plogitを使用すると、インデックス機能を取得するために predictでそれをフォローアップすることができ、および予測確率を取得すると、その周りに invlogit()機能を包むだけの問題ではありません。私も楽しいために手でそれをやった。

sysuse auto, clear 
plogit foreign weight mpg, lasso lambda(2) 

/* A. With Predict */ 
predict double phat, xb 
replace phat = invlogit(phat) 

/* B. By Hand */ 
/* (1) get a list of coefficients */ 
local coefs: colnames e(b) 

/* (2) Build the index function expression */ 
foreach x of local coefs { 
    if "`x'" != "_cons" { 
     local index "`index' `=_b[`x']'*`x' +" 
    } 
    else { 
     local index "`index' `=_b[`x']'" 
    } 
} 

/* (3) predict phat */ 
gen double phat2 = invlogit(`index') 
関連する問題