[R-sig-ME] Parameter-expanded priors with MCMCglmm

W Robert Long longrob604 at gmail.com
Sun Oct 7 12:35:12 CEST 2012


I would be grateful if anyone can advise me on how to use parameter 
expanded priors in the following context.

With random=~idh(1+year):hospital I initially had:

prior3b <- list(R = list(V = 1, fix = 1), G=list(G1=list(V=diag(7), nu=1)))

which worked OK, but I would like to use parameter expanded prior. I 
have tried this:

prior3a <- list(R = list(V = 1, fix = 1), G = list(G1 = list(V=diag(7), 
nu = 1, alpha.mu = 0, alpha.V = diag(25^2,7,7)))

But this gives the error:
alpha.V is the wrong dimension for some prior$G/prior$R elements

Thanks
Robert Long
PhD Student
University of Leeds



More information about the R-sig-mixed-models mailing list