[R-sig-ME] Unrealistic coefficient values from an MCMCglmm mixed model

Mon Apr 1 18:45:43 CEST 2019

Hi Ben,

In my model, I have set the intercepts for GeneticType so that B isn't
relative to A (using "-1"). There are only 2 levels to GeneticType.
Apologies for the confusion.

The correct model should read:
model<- MCMCglmm(LRS~GeneticType*NAO+
                             GeneticType *Temp-1,
                          random = ~Year_of_Spawning,
                          family = "poisson",
                          data = data,
                          verbose = TRUE,
                          nitt = 1010000, burnin = 1000, thin = 1000)

Regards Temp, it is mean-centred whereas NAO wasn't. I re-ran the model
with a mean-centred NAO variable and it made no difference.

Cheers,
Ronan

-- 
Ronan O'Sullivan | Ph.D student | School of Biological, Earth and
Environmental Sciences, University College Cork, Ireland |
http://fisheye.ucc.ie/toms-team/  <http://fisheye.ucc.ie/toms-team/>

Irish Ecological Association - Ordinary Committee Member

	[[alternative HTML version deleted]]