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]]

