Re: Failed while ordering equations

From: Arthur <asremlforum_at_VSNI.CO.UK>
Date: Fri, 25 Sep 2009 00:31:46 +0100

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.)