Logit multinomial em R: mlogit versus nnet
Quero executar um logit multinomial em R e usei duas bibliotecas,nnet
emlogit
, que produzem resultados diferentes e relatam diferentes tipos de estatísticas. Minhas perguntas são:
Qual é a fonte de discrepância entre os coeficientes e os erros padrão relatados pornnet
e aqueles relatados pormlogit
?
Gostaria de relatar meus resultados a umLatex
arquivo usandostargazer
. Ao fazer isso, há uma troca problemática:
Se eu usar os resultados demlogit
então recebo as estatísticas que desejo, como psuedo R ao quadrado, no entanto, a saída está no formato longo (veja o exemplo abaixo).
Se eu usar os resultados dennet
então o formato é o esperado, mas relata estatísticas nas quais não estou interessado, como o AIC, mas não inclui, por exemplo, psuedo R ao quadrado.
Eu gostaria de ter as estatísticas relatadas pormlogit
na formatação dennet
quando eu usostargazer
.
Aqui está um exemplo reproduzível, com três alternativas de escolha:
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)
A saída tex, quando compilada, é do que me refiro como "formato longo", que considero indesejável:
Agora, usandonnet
:
library(nnet)
mlogit.model2 = multinom(y ~ 1 + col1+col2, data=mydata)
stargazer(mlogit.model2)
Dá a saída tex:
que é do formato "amplo" que eu desejo. Observe os diferentes coeficientes e erros padrão.