[R] help with polytomous logistic regression

John Fox jfox at mcmaster.ca
Fri Jan 7 15:56:55 CET 2005


Dear Anne,

There's an anova() method for multinom objects that will allow you to get
likelihood-ratio tests. Likewise, the Anova() function in the car package
has a method for multinom objects.

As to plotting, I have a paper, coauthored with Bob Andersen, on producing
"effect" plots (partial plots) for polytomous logistic regression, including
multinomial and proportional-odds logistic regression. The paper, at
<http://socserv.socsci.mcmaster.ca/jfox/logit-effect-displays.pdf>, has an R
function for computing such effects, including standard errors. There are
examples in the paper as well. Eventually, I intend to incorporate this into
the effects package to make the graphs more or less automatic.

Regards,
 John

--------------------------------
John Fox
Department of Sociology
McMaster University
Hamilton, Ontario
Canada L8S 4M4
905-525-9140x23604
http://socserv.mcmaster.ca/jfox 
-------------------------------- 

> -----Original Message-----
> From: r-help-bounces at stat.math.ethz.ch 
> [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Anne
> Sent: Friday, January 07, 2005 8:53 AM
> To: R list
> Subject: [R] help with polytomous logistic regression
> 
> Hi!
> 
> I'm trying to do some ploytomous logistic regression using 
> multinom() in the nnet package, but am a bit confused about 
> interpretation of the results Is it possible to get the 
> following quantities: 
> 
> I:  maximum likelihood estimates to test for fit of model and 
> significance of each predictor
> 
> (I would like to produce a table of the following type)
> Analysis of Variance: MLE    (values are non sensical, I know!)
>                df      chi2           p
> intercept   2       57            .003
> v1            4       89            .876
> v2            2         7            .05
> LR        110     450            0.93
> 
> II:  chi square values and corresponding p-values for each 
> level of the predictor variates (I'm OK for the estimates and 
> se from summary,)
> 
>  
>                level estimate     se        chi2         p   
> intercept   1        5              .9         32          .002
>                 2       6.3            .8         31          .03
> v1             3       89            .876       43          .001
>                 4       65            .05         67          
> 0.001      
>  v2...
> 
> 
> Is there a convenient way to plot the results? (I'd like to 
> display visually the effects of the predictors: any sugestion?)
> 
> OK, I probably miss something here (no experience with (non 
> ordinal) polytomous logistic regression yet!)
> 
> Thank  very much to all
> 
> Anne
> 
> ----------------------------------------------------
> Anne Piotet
> Tel: +41 79 359 83 32 (mobile)
> Email: anne.piotet at m-td.com
> ---------------------------------------------------
> M-TD Modelling and Technology Development PSE-C
> CH-1015 Lausanne
> Switzerland
> Tel: +41 21 693 83 98
> Fax: +41 21 646 41 33
> --------------------------------------------------
>  
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at 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




More information about the R-help mailing list