Dear Tracey,
leg(DIM,1).Sire.Lactation
6 0 US !+21
0.0328252
-0.00346642 0.0425384
0.0 0.0 0.0634254
0.0 0.0 0.0106747 0.0704660
0.0 0.0 0.0 0.0 0.0634254
0.0 0.0 0.0 0.0 0.0106747 0.0704660
Sire 0 AINV
leg(DIM,1).Cow.Lactation
6 0 US !+21
0.489681
0.00261273 0.497089
0.0 0.0 0.603377
0.0 0.0 0.00561644 0.617541
0.0 0.0 0.0 0.0 0.603377
0.0 0.0 0.0 0.0 0.00561644 0.617541
Cow 0 0
Now leg(DIM,1).Sire.Lactation
is odered 2 by Nsires by 3
but you have written the structure presumably
as 2 by 3 by Nsires (i.e. 6 by Nsires)
hence it fails.
You need to declare the term as
leg(DIM,1).Lactation.Sire
in the model
and in the G structure header, and then the order will match
Similarly leg(DIM,1).Cow.Lactation
needs to be reordered to leg(DIM,1).Lactation.Cow
The .asr output points to this error, albeit obtusely
Structure for leg(DIM,1).Sire.Lactation has 966054 levels defined
Warning: leg(DIM,1).Sire.Lactation is a 2 x483027 interaction.
Are you sure the G structure matches it? 6 x161009
Sire.Lactation is a161009 x 3 interaction.
------------------------
Arthur Gilmour
Retired Principal Research Scientist (Biometrics)
-------------------- m2f --------------------
Sent using Mail2Forum (http://www.mail2forum.com).
Read this topic online here:
http://www.vsni.co.uk/forum/viewtopic.php?p=878#878
-------------------- m2f --------------------
Received on Sat Sep 25 2009 - 00:31:46 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.)