Re: Interpreting asreml output

Hi,

The random effect structure is not multivariate hence only a single
variance component for animal and dam effects are fitted. Try
something like:

random=~us(trait):ped(Animal)

for the random structures and it should make more sense.

Cheers,

Jarrod

On 11 Oct 2007, at 23:09, Derek Roff wrote:

> Please excuse this question,the answer to which should be
> obvious, but is baffling me. Using the asreml splus add on
> I wish to determine the variances and covariances for two
> traits. The command issued is
> HS.asr <- asreml(cbind(TRAIT1,TRAIT2)~trait,
>
> random=~ped(Animal,var=T)+Dam,
>
> rcov=~units:us(trait),
>
> ginverse=list(Animal=AINV), data=Data)
> summary(HS.asr)$varcomp
> which gives
> gamma component std.error
> z.ratio
> ped(Animal, var = T) 0.04206207 0.04206207 0.014288383
> 2.9437953
> Dam 0.00377920 0.00377920 0.006370836
> 0.5932031
> R!variance 1.00000000 1.00000000
> NA NA
> R!trait.TRAIT1:TRAIT1 1.75769701 1.75769701 0.054551955
> 32.2206051
> R!trait.TRAIT2:TRAIT1 0.22529479 0.22529479 0.015193119
> 14.8287390
> R!trait.TRAIT2:TRAIT2 0.14791726 0.14791726 0.008980163
> 16.4715568
>
> constraint
> ped(Animal, var = T) Positive
> Dam Positive
> R!variance Fixed
> R!trait.TRAIT1:TRAIT1 Unconstrained
> R!trait.TRAIT2:TRAIT1 Unconstrained
> R!trait.TRAIT2:TRAIT2 Unconstrained
>
> Which are the phenotypic and genetic (co)variance components
> for the two traits?
> Sincerely,
> Derek
> Derek Roff,
> Department of Biology,
> University of California,
> Riverside,
> CA 92521
> USA
>
Received on Tue Oct 12 2007 - 10:31:30 EST