2016-12-15 21 views
0

AERのivregは、機器変数によるロジット回帰をサポートしていますか?例えばRのロジスティック回帰における機器変数(AERのivreg)

Mortバイナリ変数 ICU_AdmImmedICU_LOSあり、

IV = ivreg (Mort ~ Age + Sex + APACHE + PART_SendImmed + ICU_AdmImmed + ICU_LOS | Age + Sex + APACHE + PART_SendImmed + NurseOCC_Adm + NurseOCC_Disch, 
    data = test) 

は、内因性の変数であり、 NurseOCC_AdmNurseOCC_Dischは、二つの操作変数です。

答えて

0

ロジットではなくプロビットを使用することをお勧めします。これは、解釈能力を向上させるだけでなく、ivprobを使用できるようにするものです。

data("mtcars") 

library(ivprobit) 

ivprob(I(mtcars$vs==0), mtcars$cyl, I(mtcars$disp>150), mtcars$hp) 
$coefficients 
      [,1] 
[1,] -62.48480 
[2,] 17.09827 
[3,] -74.66346 

$se 
[1] 93844628 28724155 134217728 

$tval 
       [,1] 
[1,] -6.658325e-07 
[2,] 5.952576e-07 
[3,] -5.562862e-07 

$pval 
numeric(0) 

$df 
numeric(0) 

$names 
[1] "(Intercept)" "X1"   "yhat" 

私は実際に私はこれがeconjobrumors.comの結論だったように、1つが存在しないかなり確信している、ロジットを使用してRでそれを行うには缶詰のパッケージを知りませんあなたがプロビットへのロジットを好む理由があれば、いつでもプロセスmanuallyを実行することができます。

+0

お返事ありがとうございます。私はivprobitですべての変数を手配しようとしていましたが、いくつかのエラーがあります。私は次のようにしました:iv2 = ivprobit(Mort_code〜Age + Sex + APACHE + PART_SendImmed、ICU_AdmImmed + ICU_LOS、NurseOCC_Adm + NurseOCC_Disch、 data = ds_test) – Ali

+0

@Ali私は 'data ='や '〜'そのパッケージで私の例のように、コンマを区切ります。これは次のようなものです: 'ivprob(y、x1、y2、x)'。ドキュメンテーションはさらに詳しく説明されています。この '?ivprob'のような文書を表示することができます。ああ、複数の変数を実行するには、この 'x = mtcars $ hp + mtcars $ qsec'のようなものです。この回答が役に立つ場合は、それを解決策としてマークしてください。 –

+0

通常の 'ivreg'は2段階の最小二乗を使います。この場合、「ivprob」の2つの段階は何ですか?最初の段階でまだOLSですか? – useR

関連する問題