Jeff Evans
2009-Feb-26 18:04 UTC
[R] generalized linear mixed models with a beta distribution
Has there been any follow up to this question? I have found myself wondering the same thing: How then does SAS fit a beta distributed GLMM? It also fits the negative binomial distribution. Both of these would be useful in glmer/lmer if they aren't 'illegal' as Brian suggested. Especially as SAS indicates a favorable delta BIC of over 1000 when I fit the beta to my data (could be the beginning of a great song..) versus my original binomial fit. Jeff Evans Michigan State University
Douglas Bates
2009-Feb-26 20:49 UTC
[R] generalized linear mixed models with a beta distribution
On Thu, Feb 26, 2009 at 12:04 PM, Jeff Evans <evansj18 at msu.edu> wrote:> Has there been any follow up to this question? I have found myself wondering > the same thing: How then does SAS fit a beta distributed GLMM? It also fits > the negative binomial distribution.When SAS decides to open-source their code we'll be able to find out.> Both of these would be useful in glmer/lmer if they aren't 'illegal' as > Brian suggested. Especially as SAS indicates a favorable delta BIC of over > 1000 when I fit the beta to my data (could be the beginning of a great > song..) versus my original binomial fit.Definitions of generalized linear mixed models are not entirely straightforward, at least for me. I'm making some progress but, as always, it is slower than one would like it to be.
dave fournier
2009-Feb-27 02:20 UTC
[R] generalized linear mixed models with a beta distribution
You can fit this kind of model (and negative binomial) and more difficult mixed models with AD Model Builder's random effects module which is now freely available at http://admb-project.org/ -- David A. Fournier P.O. Box 2040, Sidney, B.C. V8l 3S3 Canada Phone/FAX 250-655-3364 http://otter-rsch.com
Apparently Analagous Threads
- Nonlinear statistical modeling -- a comparison of R and AD Model Builder
- Robust estimation of variance components for a nested design
- Maximum likelihood estimation of parameters make no biological sense
- singular convergence in glmmPQL
- glmmpql and lmer keep failing