[R-sig-ME] "Animal Model" plus logistic regression
juliet.hannah at gmail.com
Thu Jun 10 23:57:07 CEST 2010
For a continuous trait, the animal model can be fit using the kinship
package. I don't think it works with discrete traits yet. My
understanding is that you need to be able to specify a custom
correlation matrix. Can glmer do that?
On Mon, Jun 7, 2010 at 11:41 AM, Bertrand Servin
<bertrand.servin at toulouse.inra.fr> wrote:
> Dear list,
> I tried searching in the archive, but cannot seem to find a relevant
> thread on my problem. I would like to fit a generalized linear
> model (logit link) where the random effects are individuals and
> I estimated (from genetic data) the expected correlation between
> these random effects. This is sometimes called "Animal model" because
> it is quite used in livestock genetics.
> This could be a simple dataset:
> The model is just logit(y) ~ mu + e
> with Var(e)=V.sigma^2_a+I.sigma^2_e
> So V is not to be estimated within the model fit, I just want to
> get an estimate of sigma^2_a, and the model residuals.
> I was thinking among these lines:
> and then modify something within this model.
> Related question, how do I actually fit the model once
> I made the modifications needed ?
> Is this the right way to go ?
> thanks for your help
> Bertrand Servin
> R-sig-mixed-models at r-project.org mailing list
More information about the R-sig-mixed-models