Skip to content
Prev 17290 / 20628 Next

binary trait correlation across environments (experimental trials) using MCMCglmm?

Hi,

Your model bi_model_trial is the correct one. However you must fix the residual variance at one in the prior. You have estimated the genetic (co)variance matrix for the two trials, from which you can obtain the genetic correlation.  Alternatively, you could fit animal+animal:Trial which assumes the genetic variances in the two trials are the same and the correlation is positive. The correlation in this latter model is obtained as VAR(animal)/(VAR(animal)+VAR(animal:Trial)).  Also, there seems to be a problem with your Surv data as it has three levels rather than 2 and so a cutpoint is being estimated.

Cheers,

Jarrod
On 2 Jan 2019, at 16:33, Plough, Louis <lplough at umces.edu<mailto:lplough at umces.edu>> wrote:
bi_model_trial

The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336.