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
1.64
0.96 7.801
0.01 0.01 0.01
ANIMAL 0 AINV
pol(sAGE,2).ide(ANIMAL) 2
pol(sAGE,2) 0 US !GP
1.64
0.96 7.801
0.01 0.01 0.01
ide(ANIMAL) 0 ID
#MOTHER 1
#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
component.
It is sometimes difficult to get solutions for the US variance structures
as it is difficult to identify the good starting values. !GP should
help.
Best wishes in the Name of Jesus Christ, my Saviour.
Check out http://members.dodo.net.au/~cargovale
Arthur Gilmour PhD
Principal Research Scientist (Biometrics)
NSW Department of Primary Industries
Orange Agricultural Institute, Forest Rd, ORANGE, 2800, AUSTRALIA
mailto:Arthur.Gilmour_at_agric.nsw.gov.au
fax: 02 6391 3899; 02 6391 3922 Australia +61
telephone work: 02 6391 3815; home: 02 6364 3288;
http://www.agric.nsw.gov.au/reader/17263
http://www.asreml.co.uk/
The ASReml discussion group ASREML-L_at_agric.nsw.gov.au
has reactivated. To join, mailto:arthur.gilmour_at_agric.nsw.gov.au
Cookbook: http://uncronopio.org/luis/asreml_cookbook.html
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.
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.)