[R-sig-ME] fit fixed-effect GLMs within each level of the random factors

Ben Bolker bolker at ufl.edu
Thu Jun 3 18:33:55 CEST 2010


  Not quite as convenient as lmList, but something like

  L <- lapply(split(yourdata,yourdata$random_factor),glm,
              family=...)
  sapply(L,coef)

  should work fine.

Shujuan Feng wrote:
> Hi all,
> 
> 
> 
> Is there a function similar to"lmList" for GLM?  I'd like to fit 
> fixed-effect GLMs within each level of the random factors to take a look at 
> the distribution of estimated parameters.
> 
> 
> 
> Thanks
> 
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models


-- 
Ben Bolker
Associate professor, Biology Dep't, Univ. of Florida
*** NEW E-MAIL ADDRESSES:
***   bbolker at gmail.com , bolker at math.mcmaster.ca
bolker at ufl.edu / people.biology.ufl.edu/bolker
GPG key: people.biology.ufl.edu/bolker/benbolker-publickey.asc




More information about the R-sig-mixed-models mailing list