I seem to be able to extract everything but the t value and Std.error
out of the summary function for polr.
Category <- factor(rep(c(-1,0,1),2))
Freq<-c(2,4,4,3,3,4)
Sample<-factor(rep(c("x","y"), rep(3,2)))
plr<-polr(Category~Sample,weights=Freq)
coefficients(summary(plr)) ???
Can anyone help
Cheers
Daniel Yearsley
dyea002 at math.auckland.ac.nz
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at
stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
I seem to be able to extract everything but the t-value and std. error out
of the summary function for polr.
Category <- factor(rep(c(-1,0,1),2))
Freq <- c(2,4,4,3,3,4)
Sample <- factor(rep(c("x","y"), rep(3,2)))
plr <- polr(Category~Sample, weights = Freq)
coefficients(summary(plr))
Can anyone help??
Cheers
Daniel Yearsley
dyea002 at math.auckland.ac.nz
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at
stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
On Mon, 12 Feb 2001, Daniel Yearsley wrote:> I seem to be able to extract everything but the t value and Std.error > out of the summary function for polr. > > Category <- factor(rep(c(-1,0,1),2)) > Freq<-c(2,4,4,3,3,4) > Sample<-factor(rep(c("x","y"), rep(3,2))) > plr<-polr(Category~Sample,weights=Freq) > coefficients(summary(plr)) ???You do it exactly as for lm. The t-values are summary(plr)$coefficients[, "t value"] and the standard errors are summary(plr)$coefficients[, "Std. Error"] These work for me in your example. (R 1.2.1, MASS from VR_6.2-3) You can also use vcov(plr) to get the estimated covaraince matrix of the estimates. -- Brian D. Ripley, ripley at stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272860 (secr) Oxford OX1 3TG, UK Fax: +44 1865 272595 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
What version of R and what platform. I just pasted your example code
directly into R (1.2.1/Windows) and it ran without an error. I get the
t-values and the Std Errors with no problem.
It does precede the output by noting:
Re-fitting to get Hessian
Value Std. Error t value
Sampley -0.2056671 0.8285269 -0.2482322
-1|0 -1.1998700 0.6613101 -1.8143835
0|1 0.3079109 0.6023141 0.5112132
Was there something else you expected to happen?
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at
stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._