I've got a bunch of similar datasets, all of which I've fit to linear
models.  I'd like to easily create arrays of a specific parameter from each
linear model (e.g., all of the intercepts in one array).  I figured I'd put
the model objects into an array, and then (somehow) I could easily create
corresponding arrays of intercepts or residuals or whatever, but I can't the
parameters back out.

 

Right now I've stored the model objects in a 2-D array:
> lms.ASP <- array(list(), c(3,4))

 

Then I fill the array element-by-element:
> surf105.lm. ASP <- lm(ASP ~ time)
> lms.ASP[1,1] <- list(surf105.lm.ASP)

 

Something is successfully being stored in the array:
> test <- lms.tx.ASP[1,1]
> test
[[1]]
Call:
lm(formula = ASP ~ time)
Coefficients:
     (Intercept)  elapsed.time..hr  
     0.430732          0.004073  

 

But I can't seem to call extraction functions on the linear models:
> fitted(lms.ASP[1,1])
NULL

 

It seems like something less than the actual linear model object is being
stored in the array, but I don't understand what's happening, or how to
easily batch-extract parameters of linear models.  Any advice?

 

 

____________________________________

Andrew D. Steen, Ph.D.

Center for Geomicrobiology, Aarhus University

Ny Munkegade 114
DK-8000 Aarhus C
Denmark
Tel: +45 8942 3241
Fax: +45 8942 2722

andrew.steen@biology.au.dk

 


	[[alternative HTML version deleted]]

