Hello, I am using brms for the first time since some people mentioned that it is able to account for multiple response variables and continuous random term (in the predictors). My main goal is to estimate genetic correlation among some traits given a certain kinship matrix.
My data looks like this:
Family Sex Diet humerus ulna femur
67441 fam_1 M H 12.255 14.260 15.725
67442 fam_1 F L 11.870 12.960 15.140
67443 fam_1 F H 11.080 12.745 14.190
While my kinship matrix looks like this:
67441 67442 67443
67441 1.0106853 0.4224381 0.5476577
67442 0.4224381 0.9045946 0.4349701
67443 0.5476577 0.4349701 0.9937757
I read in the Vignettes that to include multiple response variables I should use the mvbind() function, but I could not find a way to include my kinship matrix as a random term in my model. If anyone has any idea on how to do this, it would be really helpful.
Thank you