13 views

### Paul May

May 12, 2022, 8:34:51 AMMay 12
to R-inla discussion group
Hello,

Is there any way for R-inla to provide an identity link for Gamma GLMMs? I understand the issues with this (need to insure linear predictor \eta > 0), but that log transformation is awful strong.

In my particular application, I have repeat observation for each realization of the linear predictor

y_ij ~ Gamma(g^{-1}(\eta_j), scale);   i =1 , ... , n_j

so I can get a decent approximation of \mu_j = g^{-1}(\eta_j) with the sample mean \bar{y}_j. These look plausibly normally distributed, and a log transformation introduces a significant left tail.

Thanks,
Paul