Re: [R] linear model in the repeated data type~

From: Austin, Matt <maustin_at_amgen.com>
Date: Thu, 05 Jun 2008 15:03:08 -0700

allFits <- lmList(y ~ t|id, data=table1, pool=FALSE)

allCoefs <- sapply(allFits, coef) ## preferred by me

or

allCoefs <- list(length(allFits))
for(i in 1:length(allFits)) allCoef[[i]] <- coef(allFits[[i]])

--Matt



From: Manli Yan [mailto:manliyanrhelp_at_gmail.com] Sent: Wednesday, June 04, 2008 9:12 PM
To: Austin, Matt
Cc: r-help_at_r-project.org
Subject: Re: [R] linear model in the repeated data type~

hi:lot thanks,how to use list to extract,I type allFit$coefficents,it came to nothing, such as I need to extract the estimates,how to do it by using list

2008/6/3 Austin, Matt <maustin_at_amgen.com<mailto:maustin_at_amgen.com>>: How about

library(nlme)
allFits <- lmList(y ~ t|id, data=table1, pool=FALSE)

or

allFits <- by(table1, table1$id, function(x) lm(y ~ t, data=x))

Both ways store the results as a list, so you can access individual results using list extraction.

--Matt

-----Original Message-----

From: r-help-bounces_at_r-project.org<mailto:r-help-bounces_at_r-project.org> [mailto:r-help-bounces_at_r-project.org<mailto:r-help-bounces_at_r-project.org>] On Behalf Of Manli Yan Sent: Tuesday, June 03, 2008 9:07 PM
To: r-help_at_r-project.org<mailto:r-help_at_r-project.org> Subject: [R] linear model in the repeated data type~

 here is the data:

 y<-c(5,2,3,7,9,0,1,4,5)
id<-c(1,1,6,6,7,8,15,15,19)
t<-c(50,56,50,56,50,50,50,60,50)

table1<-data.frame(y,id,t)//longitudinal data

what I want to do is to use the linear model for each id ,then get the estimate value,like:

fit1<-lm(y~t,data=table1,subset=(id==1))

but ,you can see the variable "id" is quite irregular,they are not arranaged in order and many number missing,if I write a loop by using "for",it will give me a lot "NA", and for sure ,I dont want to type id=## for about 500 times,any one know how to deal with it?

       [[alternative HTML version deleted]]



R-help_at_r-project.org<mailto: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<http://www.r-project.org/posting-guide.html> and provide commented, minimal, self-contained, reproducible code.

        [[alternative HTML version deleted]]



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. Received on Thu 05 Jun 2008 - 23:06:28 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 06 Jun 2008 - 00:30:41 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive