R で多項ロジットを実行したいのですが、2 つのライブラリ と を使用しましたnnet。mlogitこれらは異なる結果を生成し、異なるタイプの統計を報告します。私の質問は次のとおりです。
nnetによって報告された係数と標準誤差と によって報告された係数と標準誤差の不一致の原因は何mlogitですか?Latexを使用して結果をファイルに報告したいと思いますstargazer。そうする場合、問題のあるトレードオフがあります。その結果を使用すると、
mlogit疑似 R 二乗などの必要な統計が得られますが、出力は長い形式になります (以下の例を参照)。結果を使用すると
nnet、形式は期待どおりですが、AIC などの興味のない統計が報告されますが、たとえば、疑似 R 二乗は含まれません。
を使用する場合
mlogit、 のフォーマットで統計を報告してもらいたいです。nnetstargazer
以下は再現可能な例で、3 つの選択肢があります。
library(mlogit)
df = data.frame(c(0,1,1,2,0,1,0), c(1,6,7,4,2,2,1), c(683,276,756,487,776,100,982))
colnames(df) <- c('y', 'col1', 'col2')
mydata = df
mldata <- mlogit.data(mydata, choice="y", shape="wide")
mlogit.model1 <- mlogit(y ~ 1| col1+col2, data=mldata)
コンパイル時の tex 出力は、望ましくないと思われる「長い形式」と呼ばれるものです。
今、使用してnnet：
library(nnet)
mlogit.model2 = multinom(y ~ 1 + col1+col2, data=mydata)
stargazer(mlogit.model2)
tex出力を与えます：
これは私が望む「ワイド」フォーマットです。異なる係数と標準誤差に注意してください。