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.