Re: Interpreting asreml output

Dear Jarrod,
Thanks for the answer. Obvious when you pointed it out!
Sincerely,
Derek

---- Original message ----
>Date: Fri, 12 Oct 2007 10:31:30 +0100
>From: Jarrod Hadfield <J.Hadfield_at_ED.AC.UK>
>Subject: Re: Interpreting asreml output
>To: ASREML-L_at_AGRIC.NSW.GOV.AU
>
>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
>>
Derek Roff,
Department of Biology,
University of California,
Riverside,
CA 92521
USA
Received on Tue Oct 12 2007 - 10:43:55 EST