Re: AS-REML for stillbirth
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: AS-REML for stillbirth



Dear Mike,


From the information you supply, ASREML can handle the problem for a moderate
(say < 40,000) calves.  From your description, I presume you will prepare a sire
pedigree file (which includes all sires and maternal grandsires)

i.e.  sire.ped  contains 3 fields
sire sire_of_sire  dam or mgs  

Assuming I understand the model you want (as there are many possabilities)
the ASREML job might be something like

Still birth study
 calf
 sire !P
 mgs !P
 sex 2     
 hys  2222    #  where 2222 is number of HYS
 stillbirth
sire.ped   !skip 1    # !mgs
calf.dat
 
stillbirth ~ sex !r sire mgs !f hys    # Normal trait
  or
   
stillbirth !bin ~ sex !r sire mgs !f hys # Binomial trait



I would fit this first, and when you have values for the
sire and mgs variance components,  extend the model to estimate their 
covariance by appending

0 0 1
sire 2
2 0 US   s-s   s-mgs   mgs-mgs
sire 0 AINV

where  s-s is replaced by the sire variance component,
       s-mgs  by the covariance
       mgs-mgs by the maternal variance component.
       
       

> X-Authentication-Warning: lamb.chiswick.anprod.csiro.au: petidomo set sender 
to asreml-owner@ram.chiswick.anprod.csiro.au using -f
> From: "Mingfeng Luo" <MLUO@aps.uoguelph.ca>
> To: asreml@ram.chiswick.anprod.csiro.au
> Date: Wed, 9 Jun 1999 10:33:15 EDT
> Subject: AS-REML for stillbirth
> 
>  I try to estimate the genetic parameters of stillbirth of dairy 
>  cattle cattle. The model includes the sex effect of calf, HYS as 
>  random effects, sire and maternal grandsire as random genetic 
>  effects.
>  
>  Can AS-REML this problem?
>  
>  Thanks in advance.
>  
>  
>  Mike Luo
> **********************************
> *  Mike F. Luo                   *
> *  Animal & Poultry Science      * 
> *  Unversity of Guelph           *
> *  Guelph, ON, N1G 2W1           * 
> *  Tel: (519)824-4120 X 8354     *         
> *  Fax: (519)767-0573            * 
> **********************************                                   
>               
> 
> --
> 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 mailinglist 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