[R] Noobie question, regression across levels
Ben Bolker
bolker at ufl.edu
Thu Dec 18 14:56:58 CET 2008
AllenL wrote:
>
> This is what I ended up using:
>
> Data.subset<-data.frame(PlotFinal,YearFinal,BioFinal,ritFinal)
> ###Subset of main dataset
> Data.length<-sapply(split(Data.subset,PlotFinal),nrow) ### The
> number of data points in each plot
>
> Data.sub<-split(Data.subset,PlotFinal) ##Split Data.subset by
> PlotFinal
> Good.levels<-Data.sub[Data.length>10] ##Take only those Plots that
> have >10 years of observation
>
> lmfun<-function(Good.levels)
> {lm(ritFinal~BioFinal+YearFinal,data=Good.levels)} ##Set up the
> ###regression function
>
Potentially slightly confusing to have the same name within the lmfun
function ... I probably would
have written this as
lmfun<-function(dat) {lm(ritFinal~BioFinal+YearFinal,data=dat)}
but the results would be identical.
> lmList<-lapply(Good.levels,lmfun) ###Apply above
> function to all "good levels"
>
> coef.List<-lapply(lmList,coef) ###List of the coefficients of
> EACH above regression
>
> This seems to have worked beautifully.
> So, my current problem is extracting the coefficients from "coef.List"
> (and I think this is getting to the roof of my general confusion).
> coef.List seems to be a list of a list of lm objects, right? How come them
> I cannot call them by their position directly (ie. coef.List[1] gives me
> the first lm object. I cannot seem to call the first coefficient in this
> object (the intercept).)
> My goal is to have a bunch of new lists, each one listing a subset of
> coefficients from each regression (ie. one will be a list of all the
> intercepts, another the mean of all the "YearFinal" effects, etc).
>
lmList (not coef.List) should be a list of lm objects, coef.List is a list
of vectors.
I actually think your confusion is simpler/more fundamental. If you want
the
first element within the first element of coef.List you actually need
coef.List[[1]][1],
not coef.List[1][1] ...
You should also consider sapply instead of lapply above, that will give you
a matrix of coefficients (which may be easier to deal with)
--
View this message in context: http://www.nabble.com/Noobie-question%2C-regression-across-levels-tp21020222p21073282.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list