[R] Simultaneously summarizing many models

ONKELINX, Thierry Thierry.ONKELINX at inbo.be
Fri Feb 22 13:00:22 CET 2008


Store the models into a list and the use an apply-function.

x=1:10 #this creates some example data
y=rnorm(10)
m1=lm(x~y)
m2=lm(x~1)
Models <- list(m1, m2)
lapply(Models, summary)

HTH,

Thierry  


----------------------------------------------------------------------------
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest
Cel biometrie, methodologie en kwaliteitszorg / Section biometrics, methodology and quality assurance
Gaverstraat 4
9500 Geraardsbergen
Belgium 
tel. + 32 54/436 185
Thierry.Onkelinx op inbo.be 
www.inbo.be 

Do not put your faith in what statistics say until you have carefully considered what they do not say.  ~William W. Watt
A statistical analysis, properly conducted, is a delicate dissection of uncertainties, a surgery of suppositions. ~M.J.Moroney

-----Oorspronkelijk bericht-----
Van: r-help-bounces op r-project.org [mailto:r-help-bounces op r-project.org] Namens Christoph Scherber
Verzonden: vrijdag 22 februari 2008 12:46
Aan: r-help op stat.math.ethz.ch
Onderwerp: [R] Simultaneously summarizing many models

Dear R users,

Let´s say I have 10 models, each named m1,m2,m3..., and I would like to summarize them automatically 
and simultaneously  - e.g., to extract parameter estimates later on from all models; how can I do that?

I have tried:

x=1:10 #this creates some example data
y=rnorm(10)

m1=lm(x~y)
m2=lm(x~1)

sum.lms=function(x)summary(paste("m",x,sep=""))
sum.lms(1:2)

but the paste() function always returns "m1" "m2" in quotation marks instead of m1 m2

Do you know a way around this? I guess I would need to use expression() or substitute() but I don´t 
know how to construct m1, m2... without the quotation marks.

Thanks very much for your help!

Best wishes,

Christoph


(using R 2.5.1 on Windows XP SP2)



-- 
Dr. Christoph Scherber
DNPW, Agroecology
University of Goettingen
Waldweg 26
D-37073 Goettingen
Germany

phone +49(0)551 39 8807
fax   +49(0)551 39 8806
homepage www.gwdg.de/~cscherb1

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



More information about the R-help mailing list