Hi all,
I'm trying to amend a previous model that included gene x sex
interactions, but now i'm looking at gene x sex x environment
interactions. I'm getting an error message appear that says 'V is
the wrong dimension for some prior$G/prior$R elements'
Here is my code -
prior.model.2<-list(R=list(V=diag(2)/252, nu=0.01),G=list (
G1=list (V=diag(2)/252, nu=2, alpha.mu=c(0,0),
alpha.V=diag(2)/252)))
model.2 <- MCMCglmm(relative_fec ~
sex*rep-1,random=~us(sex:temp):line, rcov=~idh(sex):units,
family="gaussian", nitt = 100000, burnin = 30000, thin=30, data =
all_treatments, prior = prior.model.2, verbose = FALSE)
I have relative fecundity measurements for males and females from 14
lines, that were measured in 3 temperatures. Each line, sex and temp
combination had 6 replicates.
Can anyone guide me on my prior here please?
Many thanks,
Tanya
[[alternative HTML version deleted]]