Re: Specifying variance components
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Specifying variance components



Hi Bruce
> 
> Hi,
> I need some help with specifying variance components under a multivariate 
> model with maternal effects.  The data set does not have sufficient 
> information to estimate the additive-maternal covariance.  So I am just 
> fitting dam as a second effect.
> 
> Below is a my current .as file (simplified as the random effects are the 
> important part at present), is this a correct specification?
> Is this the same as saying:
> wt0 wtw ~ Trait !r Tr.tag at(Tr,1).dam at(Tr,2).dam ?

YES  well same as
  wt0 wtw ~ Trait !r Tr.tag at(Tr,1).dam 0.2 at(Tr,2).dam 4.2
and omitting the last G structure
> 
> If I decide to fit only a maternal effect to wt0, the model changes to:
> wt0 wtw ~ Trait !r Tr.tag at(Tr,1).dam

  wt0 wtw ~ Trait !r Tr.tag at(Tr,1).dam 0.2
> 
> How would I specify starting variance components for the dam component?
> 
> 
> Thanks in advance for any comments,
> Bruce Southey
> 
> 
> 

Best wishes 
Arthur
> Bivariate with maternal
>  tag  !p  # 2525 animals in pedigree
>  sire
>  dam  !p
>  br 6 !I
>  yb 6 !I
>  sx   2 !I
>  aod   4 !I
>  bd    !-12260
>  aw
> ped.file
> wt.dat  !MAXIT 1 # !CONTINUE
> wt0 wtw ~ Trait !r Tr.tag Tr.dam
> 1 2 2   # !step 0.0001
> 0 0 0
> Trait 0 US 0.29202 0.63988 3.0156
> Tr.tag 2
> Trait 0 US 0.28609 0.85243 4.6901
> tag 0 0
> Tr.dam 2
> 2 0 DIAG 0.2  4.2
> dam 0 0
> 
> 
> 
> --
> Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml


<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
Arthur Gilmour PhD                  mailto:Arthur.Gilmour@agric.nsw.gov.au
Senior Research Scientist (Biometrics)               fax: <61> 2 6391 3899
NSW Agriculture                                           <61> 2 6391 3922
Orange Agricultural Institute             telephone work: <61> 2 6391 3815
Forest Rd, ORANGE, 2800, AUSTRALIA                  home: <61> 2 6362 0046

ASREML is still free by anonymous ftp from pub/aar on ftp.res.bbsrc.ac.uk
    or point your web browser at ftp://ftp.res.bbsrc.ac.uk/pub/aar/ 

To join the asreml discussion list, send the message  
     subscribe
mailto:asreml-request@chiswick.anprod.CSIRO.au

To send messages to the list, mailto:asreml@chiswick.anprod.CSIRO.au

Asreml list archive: http://www.chiswick.anprod.csiro.au/lists/asreml

                        <> <> <> <> <> <> <>
"Christ Jesus came into the world to save sinners"   I Timothy 1:15.
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>


--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml