R help - Mar 2011 - rms: getting adjusted R^2 from ols object

How can I extract the adjusted R^2 value from an ols object (using rms package)?

library(rms)
x <- rnorm(10)
y <- x + rnorm(10)
ols1 <- ols(y ~ x)

Typing "ols1" displays adjusted R^2 among other things, but how can I
assign it to a variable? I tried str(ols1) but couldn't see where to
go from there.

Thanks,
Mark Seeto

Hi,

If you look in the print.ols function you can see that the adjusted R^2
is not stored in the ols object but calculated on the fly, the relevant
code snippet from that function looks to be:

    x <- ols1

    stats <- x$stats
    resid <- x$residuals
    n <- length(resid)
    ndf <- stats["d.f."]
    df <- c(ndf, n - ndf - 1, ndf)
    r2 <- stats["R2"]
    rdf <- df[2]
    rsqa <- 1 - (1 - r2) * (n - 1)/rdf


Martyn


-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]
On Behalf Of Mark Seeto
Sent: 09 March 2011 09:10
To: r-help at r-project.org
Subject: [R] rms: getting adjusted R^2 from ols object

How can I extract the adjusted R^2 value from an ols object (using rms
package)?

library(rms)
x <- rnorm(10)
y <- x + rnorm(10)
ols1 <- ols(y ~ x)

Typing "ols1" displays adjusted R^2 among other things, but how can I
assign it to a variable? I tried str(ols1) but couldn't see where to
go from there.

Thanks,
Mark Seeto

______________________________________________
R-help at 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.

________________________________________________________________________
This e-mail has been scanned for all viruses by Star.\ _...{{dropped:12}}

On Mar 09, 2011; 11:09am Mark Seto wrote:
>> How can I extract the adjusted R^2 value from an ols object (using rms
>> package)? 
>> library(rms) 
>> x <- rnorm(10) 
>> y <- x + rnorm(10) 
>> ols1 <- ols(y ~ x)
##
ols1$stats
ols1$stats[4]

Regards, Mark.

--
View this message in context:
http://r.789695.n4.nabble.com/rms-getting-adjusted-R-2-from-ols-object-tp3343118p3343220.html
Sent from the R help mailing list archive at Nabble.com.