I am trying to fit a multivariate sire maternal grandsire model. The solutions
I get seem to be OK. However I would like to estimate the covariance between
the sire and maternal component of the same trait. Is this possible? I have
used the following .as file:
calveas 6
weight
gest
race 6
sexparity 4
yearseason 100 !I
m_p 90 !P
f_gp 129 !I
../pest/mgs.ped !MGS
../pest/mgsall.dat !skip 1
calveas weight gest ~ Trait Tr.race Tr.sexparity !r Tr.m_p Tr.f_gp !f
Tr.yearseason
1 2 2
0 0 0
3 0 US 0.4 0.9 16 0.5 6.5 17
Tr.m_p
3 0 US 0.02 0.2 2.2 0.08 0.5 3.7
m_p
Tr.f_p
3 0 US 0.02 0.05 1.6 0.03 0.4 1.3
f_p
For the covariances I get the following estimates:
Sire
Maternal
Trait 1 (Calveas) 0.015 0.16 0.085
Trait 2 (Weight) 2.0 0.52
Trait 3 (Gest) 3.2
Trait 1
0.019 0.059 0.035
Trait 2
1.4 0.37
Trait 3
1.1
But I would like to have the following estimates as well:
Sire
Maternal
Trait 1 (Calveas) 0.015 0.16 0.085 covtrait1
Trait 2 (Weight) 2.0 0.52
covtrait2
Trait 3 (Gest) 3.2
covtrait3
Trait 1
0.019 0.059 0.035
Trait 2
1.4 0.37
Trait 3
1.1
In fact it might be possible to estimate the (non-symmetric) 3x3 matrix with
covariances between the maternal and sire components.
Could anyone help me with this problem?
Alex Hulzebosch
Department of Animal Breeding and Genetics
Wageningen, the Netherlands
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml