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

From: Moshe Olshansky <m_olshansky_at_yahoo.com>
Date: Tue, 03 Jun 2008 21:24:17 -0700 (PDT)


Try something like this:

fits <- list(500)
for (i in 1:500)
{
if (sum(table1$id == i) == 0) fits[[i]] <- NA else fits[[i]] <- lm(y~t,data=table1,subset=(id==i)) }

> From: Manli Yan <manliyanrhelp_at_gmail.com>
> Subject: [R] linear model in the repeated data type~
> To: r-help@r-project.org
> Received: Wednesday, 4 June, 2008, 2:06 PM
> 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 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.



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 Wed 04 Jun 2008 - 05:14:27 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 Wed 04 Jun 2008 - 06:30:38 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