Dear all,

I have found some difference of the results between multinom() function in
R and multinomial logistic regression in SPSS software.

The input data, model and parameters are below:

choles <- c(94, 158, 133, 164, 162, 182, 140, 157, 146, 182);
sbp <- c(105, 121, 128, 149, 132, 103, 97, 128, 114, 129);
case <- c(1, 3, 3, 2, 1, 2, 3, 1, 2, 2);

result <- multinom(case ~ choles + sbp + choles:sbp, abstol=1.0e-20,
reltol=1.0e-20, MaxNWts=10000);

However, the estimated coeffcients and standard errors of the coefficeints
are different from the SPSS.

For instance,

the estimated coefficients of the variable "choles" are 0.1946555 and
0.6244513 from the above result, but the SPSS result are  0.213120 and
0.662575.

Standard errors are much more different.

Why these kind of discrepancies occur?

I'll be very appreciated to whom explain the above matter.

Thanks.

Jinho Yoo.

	[[alternative HTML version deleted]]

