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

Re: Constraining variance components



Dear Bruce,

I trust you are getting a few more ASREML jobs going now.

The analysis you propose should be straight forward.
The message indicates a syntax error so it is a pity you did not
enclose the .as file [or its tail].

The following syntax [as a simple model of what I believe you have]
should work.

  Different traits in different years.
    animal !P
    sire
    dam
    year 2
    y1
    y2
  data.dat
  data.dat
  y1 y2 ~ year !r Trait.animal
  1 2 1
  0
  Trait 0 US  20 0 30 
  
  Trait.animal 2
  Trait 0 US 5 2 8
  animal
  
  
  
Since there is no information on the error correlation, it
should remain at zero.

An alternative formulation would be as follows 
[ requires modifying the data  - assumes  missing observations are
coded as zero so that y1 and y2 can be added.  Requires
data sorted so that y1 ccurs first, then y2 and that you know
how many of each[say 1000 and 1100].

  Different traits in different years.
    animal !P
    sire
    dam
    year 2
    y      !+V6       #  put both traits in same column
    Trait 2 !>0 !+1   #  code 1 for trait 1, 2 for trait 2
  data.dat
  data.dat
  y ~ year !r Trait.animal
  2 1 1
  1000   !S2=20    # need to nominate correct numbers of each trait
  1100  !S2=30 
  
  Trait.animal 2
  Trait 0 US 5 2 8
  animal
  
  
  I hope this is what you need
  
  Arthur


> Hi,
> I want to estimate the genetic correlation between two traits measured on 
> different animals in different years (essentially the same structure as one 
> trait measured only on males and another measured only on females).  I have 
> different number of observations per 'trait'.
> 
> What is the best way to handle this?
> I get the error "Failed to parse R structure line" with trying to constrain 
> the residual covariance to zero.
> 
> Thanks in advance,
> Bruce Southey
> 
> 
> --
> 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

                        <> <> <> <> <> <> <>
"You (Father) have given Him (Your Son) authority over all flesh,
   that He should give eternal life to as many as You have given Him.
And this is eternal life,
   that they may know You, the only true God,
   and Jesus Christ whom You have sent."  John 17: 2-3.
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>


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