[R] p value for coefficients in multinomial model

From: Xingwang Ye <xwye_at_sibs.ac.cn>
Date: Thu 05 Apr 2007 - 13:00:58 GMT

Dear all,    

1)how can I easily get p value for the coefficients of factors in a multinomial model? 2)why the p values for "type III" test with Anova are not identical to that from SAS?  

for example:

 A,B and C are categorical variables,but the proportions of each level in each categorical variables are not balance. Y is a nominal variables(>=3 categories);       

To do multinomial logistic regression in SAS, I tried proc logistic data=a;

    class A B C Y/param = ref ref = first;     model Y=A B C/link=logit;

In R, I tried this:

    b <- as.data.frame(lapply(a[,c("A","B","C","Y")], as.factor))     mod.mul<- (Y~A+B+C,data=b)

I can get the identical coefficients and std.errors in R as in SAS. how can I easily get the p value for each coeefficients? However,the "LR Chisq" for each variables are similar but not identical to that from SAS, why?          

BTW,please forward your kind answers to my email. Thank you.      

Best wishes
yours, sincerely
Xingwang Ye

PhD candidate     
Research Group of Nutrition Related Cancers and Other Chronic Diseases      
Institute for Nutritional Sciences,  
Shanghai Institutes of Biological Sciences,     
Chinese Academy of Sciences     
P.O.Box 32     
294 Taiyuan Road     
Shanghai 200031     


R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Thu Apr 05 23:03:36 2007

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Thu 05 Apr 2007 - 15:30:47 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.