[R] individual intercept and slope
Phil Spector
spector at stat.berkeley.edu
Tue Nov 2 16:53:19 CET 2010
You didn't say what form you wanted the output in, but
here's one way:
> sapply(split(dat,dat$individual),function(s)lm(height~time,data=s)$coef)
1 2
(Intercept) 8.466667 19.866667
time 2.485714 -2.057143
- Phil Spector
Statistical Computing Facility
Department of Statistics
UC Berkeley
spector at stat.berkeley.edu
On Tue, 2 Nov 2010, Rosario Garcia Gil wrote:
> Hello
>
> I would like to extract the estimates for the intercept and slope by individual for growth from a lm fit.
> Any advice?
>
> Individual Time point Height
> 1 1 10
> 1 2 11
> 1 3 23
> 1 4 15
> 1 5 21
> 1 6 23
> 2 1 24
> 2 2 12
> 2 3 9
> 2 4 10
> 2 5 11
> 2 6 10
> ...
>
> Thanks
>
> ______________________________________________
> 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.
>
More information about the R-help
mailing list