There is also clm() (for cumulative link models) from package ordinal
that has much the same interface that polr() has, but it does give you
p-values for the regression parameters. A simple example from
examples(clm):
library(ordinal)
data(wine)
fm1 <- clm(rating ~ contact + temp, data=wine)
(summ <- summary(fm1))
giving the output:
formula: rating ~ contact + temp
data: wine
link threshold nobs logLik AIC niter max.grad cond.H
logit flexible 72 -86.49 184.98 6(0) 4.01e-12 2.7e+01
Coefficients:
Estimate Std. Error z value Pr(>|z|)
contactyes 1.5278 0.4766 3.205 0.00135 **
tempwarm 2.5031 0.5287 4.735 2.19e-06 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Threshold coefficients:
Estimate Std. Error z value
1|2 -1.3444 0.5171 -2.600
2|3 1.2508 0.4379 2.857
3|4 3.4669 0.5978 5.800
4|5 5.0064 0.7309 6.850
The p-values are extracted with:
> coef(summ)[5:6, 4]
contactyes tempwarm
1.348440e-03 2.194605e-06
Cheers,
Rune
On 30 August 2011 03:30, Frank Harrell <f.harr...@vanderbilt.edu> wrote:
> Use the rms package to replace Design. Run anova(fit object from lrm) which
> produces a matrix from which you can extract P-values. This also handles
> the case of multiple degrees of freedom per predictor.
> Frank
>
> Debs Majumdar wrote:
>>
>> Hi,
>>
>> Are there any packages which prints out p-values for OLR's (like
>> `ologit' from Stata)? I want to run a bunch of OLRs and print the p-value
>> for the first coefficient from each of them.
>>
>>
>> I checked polr() under MASS and it doesn't.
>>
>> There's a lrm() function under Design which does print out p-values but I
>> couldn't extract p-values from the output.
>>
>>
>> Thanks,
>>
>> Debs
>>
>>
>> ______________________________________________
>> R-help@r-project.org 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.
>>
>
>
> -----
> Frank Harrell
> Department of Biostatistics, Vanderbilt University
> --
> View this message in context:
> http://r.789695.n4.nabble.com/Ordinal-logistic-regression-p-values-tp3777674p3777768.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help@r-project.org 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.
>
--
Rune Haubo Bojesen Christensen
Ph.D. Student, M.Sc. Eng.
Phone: (+45) 45 25 33 63
Mobile: (+45) 30 26 45 54
DTU Informatics, Section for Statistics
Technical University of Denmark, Build. 305, Room 122,
DK-2800 Kgs. Lyngby, Denmark
______________________________________________
R-help@r-project.org 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.
