[Inla-commits] Heritability under a zero-inflated Poisson model

[Peter Chi]

Hello all,

I am trying to estimate heritability of a trait under the assumption of a zero-inflated Poisson distribution, using the AnimalINLA extension package.

I am a bit confused about the specification of lambda as described in the Holand 2011 paper. They suggest estimating lambda as \sum{y_i} / \sum{n_i} which I think for me amounts simply to the mean of my trait value, as each of my n_i's are 1. 

So, 2 questions:
1) Can anyone explain the justification for this choice of lambda? I've looked at one of the reference papers (Vazquez 2009) but I'm still unclear. Holand et al. state in their eqn (5) that heritability is estimated as: 
\sigma^2_u / (\sigma^2_u + \lambda^{-1})

So I guess \lambda^{-1} is supposed to represent the environmental variance, but I'm not sure why this would be the case.

2) As shown above, the denominator has \lambda^{-1} in it. But, I have suspicions that the estimated posterior distributions of heritability from animal.inla have been calculated with simply \lambda in the denominator, not \lambda^{-1}... Does anyone have experience with this?

Thanks!


Peter