Ich führte ein multinomiales Logit-Modell in JMP durch und erhielt Ergebnisse, die den AIC sowie Chi-Quadrat-p-Werte für jede Parameterschätzung enthielten. Das Modell hat ein kategoriales Ergebnis und 7 kategoriale erklärende Varianten.
Ich passe dann an, was ich dachte, würde dasselbe Modell in R bauen, indem ich die multinom
Funktion im nnet- Paket verwende.
Der Code war im Grunde:
fit1 <- multinom(y ~ x1+x2+...xn,data=mydata);
summary(fit1);
Die beiden geben jedoch unterschiedliche Ergebnisse. Mit JMP ist der AIC 2923,21 und mit nnet::multinom
dem AIC 3116,588.
Meine erste Frage lautet also: Ist eines der Modelle falsch?
Die zweite Sache ist, dass JMP Chi-Quadrat-p-Werte für jede Parameterschätzung angibt, die ich benötige. Zusammenfassung auf dem Multinom fit1
läuft nicht - es gibt nur die Schätzungen, AIC und Abweichung.
Meine zweite Frage lautet daher: Gibt es eine Möglichkeit, die p-Werte für das Modell und die Schätzungen bei der Verwendung zu ermitteln nnet::multinom
?
Ich weiß, dass mlogit ein weiteres R-Paket dafür ist, und es sieht so aus, als ob seine Ausgabe die p-Werte enthält. Ich konnte jedoch nicht mlogit
mit meinen Daten ausgeführt werden. Ich glaube, ich hatte die Daten richtig formatiert, aber es hieß, ich hätte eine ungültige Formel. Ich habe die gleiche Formel verwendet, die ich verwendet habe multinom
, aber es scheint, dass ein anderes Format für die Verwendung einer Pipe erforderlich ist, und ich verstehe nicht, wie das funktioniert.
Vielen Dank.