covariance_structure_of_correlated_random_effects_for_multiple_traits
This is an old revision of the document!
The documentation is bit confusing on how to declare several correlated effects for several traits. Whereas renumf90
does this automatically for maternal effects models, it is not so obvious if one does it by itself.
For instance terminal crosses models (e.g. |here) use:
- trait_pure: y = u + e (animal)
- trait_crossbred : y = s + e (sire of crossbred)
Imagine that we have two such traits. Data file has “animal” in column 4 and “sire” in column 5.
... 4 0 3084 cross 0 5 3084 cross ... RANDOM_GROUP 4 5 RANDOM_TYPE add_animal FILE ../SNPgr/ped.txt (CO)VARIANCES 10 0 0 5 0 0 0 0 5 0 0 12
Where we have a 4×4 matrix for covariances. In fact
Effect 4 Effect 5 Trait 1 Trait2 Trait1 Trait2
Effect 4 Trait 1 Var(u) 0 0 Cov(u,s)
Trait2 0 0 0 0
Effect 5 Trait 1 0 0 0 0
Trait2 Cov(s,u) Var(s)
covariance_structure_of_correlated_random_effects_for_multiple_traits.1522349099.txt.gz · Last modified: 2024/03/25 18:22 (external edit)