2017-06-06 11 views
0

私はロジスティック回帰を実行しており、その結果を失業者に報告しています。私はapply.coef = ORオプションを適用すると(オッズ比が報告されるように)、意義の星は間違って報告されることに気付きました - 場合によっては3つではなく星があり、時には星がありますしてはいけません。 たとえば、ロッジのオッズ比を報告しているStargazerは、間違った有意性を報告しています。

stargazer(basic.logit.model, 
     type="html", 
     apply.coef = OR, 
     column.labels = c("Base"), 
     dep.var.labels.include = FALSE, 
     digits=2, out=("basic_model_only.htm")) 

利回り

X 0.33 
(0.23) 
Constant 0.03 
(0.11) 
Observations 6,532 
Log Likelihood -552.64 
Akaike Inf. Crit. 1,109.28 
Note: *p<0.1; **p<0.05; ***p<0.01 

と適用-COEFオプションを指定せずに、結果は以下のとおりです。

X -1.10*** 
(0.23) 
Constant -3.68*** 
(0.11) 
Observations 6,532 
Log Likelihood -552.64 
Akaike Inf. Crit. 1,109.28 
Note: *p<0.1; **p<0.05; ***p<0.01 

私が何をしないのですか?

答えて

1

一般に、結果を再現できるように最小限の作業例を提供するよう努めなければなりません。この場合は、basic.logit.modelを生成するために使用するデータとコードがあることが良いでしょう。

あなたの質問には、apply.coefは係数だけを変換しますが、結果が示すように標準誤差は変換しません。これらの変換されていないSEを使用して重要度レベルを計算します。

これを避けるには、p引数を使用して、スタガジャーにカスタムp値(元のモデルのもの)を提供してください。

これは

p.values <- list(summary(basic.logit.model)$coefficients[,4] 

stargazer(basic.logit.model, 
    type="html", 
    apply.coef = OR, 
    p = p.values, 
    column.labels = c("Base"), 
    dep.var.labels.include = FALSE, 
    digits=2, out=("basic_model_only.htm")) 
のために働く必要があります
関連する問題