Re: heterogeneous error structure

From: Arthur Gilmour <arthur.gilmour_at_AGRIC.NSW.GOV.AU>
Date: Wed, 9 Mar 2005 12:26:50 +1100

Dear Alistair,

I enclose a random regression example for your consideration.

Following is an amended version of your job.

TRAIT ~ mu SEX AGEF !r pol(sAGE,2).ANIMAL pol(sAGE,2).ide(ANIMAL) MOTHER
 !f mv
0 0 2
#1 1 3

#0 0 ID 3 # single homogeneous error variance

pol(sAGE,2).ANIMAL 2
pol(sAGE,2) 0 US !GP
0.96 7.801
0.01 0.01 0.01

pol(sAGE,2).ide(ANIMAL) 2
pol(sAGE,2) 0 US !GP
0.96 7.801
0.01 0.01 0.01
ide(ANIMAL) 0 ID

#MOTHER 0 IDV 0.226

This assumes an animal model with repeat observations in time
for the animals. This has a common independent residual variance plus
the animal level random regression covariance structure at the residual
and the genetic level.

You could extend it to have a random regression for the maternal

It is sometimes difficult to get solutions for the US variance structures
as it is difficult to identify the good starting values. !GP should

Best wishes in the Name of Jesus Christ, my Saviour.
Check out

Arthur Gilmour PhD
Principal Research Scientist (Biometrics)
NSW Department of Primary Industries
Orange Agricultural Institute, Forest Rd, ORANGE, 2800, AUSTRALIA
fax: 02 6391 3899; 02 6391 3922 Australia +61
telephone work: 02 6391 3815; home: 02 6364 3288;
The ASReml discussion group
   has reactivated. To join,

Ellen is doing well after her Hip replacement on 24 Feb but
she will need crutches for a while yet.

This message is intended for the addressee named and may contain
confidential information. If you are not the intended recipient or
received it in error, please delete the message and notify sender. Views
expressed are those of the individual sender and are not necessarily the
views of their organisation.

Received on Tue Mar 09 2005 - 12:26:50 EST

This webpage is part of the ASReml-l discussion list archives 2004-2010. More information on ASReml can be found at the VSN website. This discussion list is now deprecated - please use the VSN forum for discussion on ASReml. (These online archives were generated using the hypermail package.)