Failed while ordering equations

Hi,

I wondered if anyone could suggest what I could do when I have the message 'Failed while ordering equations.' In the manual it suggests either increasing the workspace or simplifying the model, so I have tried various things but still get the message. Or can anyone see that there is anything obviously wrong in the .as or .asr file. There is a warning that for cow there are more levels than specified but I have checked this and the figure I have written in the .as file is correct.

Many thanks for your time,

Tracey

AS FILE

 Sire 1654 !P
 Cow 43592 !I
 Lactation
 HTD 109454 !I
 SeasonSubClass 99 !I
 DIM !-154.5 !/150
 CalvingYear 12 !I
 Het1
 Recomb1
 SCC1000 !m0
 lnSCC !=SCC1000 !^0 !M0
../MastitisPedigree.ped
../SCC_01092009_Lactation1.dat !CSV !MAXIT 100 !STEP 0.001 !WORKSPACE
lnSCC ~ mu leg(DIM,4).SeasonSubClass CalvingYear Het1 Recomb1,
!r Sire.leg(DIM,2) Cow.leg(DIM,2) HTD,
!f mv
4 1 2
26947 0 0 !s2=1.07909 #lactation 1 dim 4 to 24
34332 0 0 !s2=0.659710 #lactation 1 dim 25 to 49
282018 0 0 !s2=0.433130 #lactation 1 dim 50 to 249
72177 0 0 !s2=0.348163 #lactation 1 dim 250 to 305
Sire.leg 2
3 0 US !+6
0.0328252
-0.00346642 0.0425384
0.0 0.0 0.01
Sire 0 AINV
Cow.leg 2
3 0 US !+6
0.489681
0.00261273 0.497089
0.0 0.0 0.01
Cow

ASR FILE

  Sire 1654 !P
  Cow 43592 !I
  HTD 109454 !I
  SeasonSubClass 99 !I
  DIM !-154.5 !/150
  CalvingYear 12 !I
  SCC1000 !m0
  lnSCC !=SCC1000 !^0 !M0
 Reading pedigree file ../MastitisPedigree.ped : skipping 0 lines
 Using an adapted version of Meuwissen & Luo GSE 1992 305-313: Specify !METHOD 1 for column method.
 PEDIGREE [../MastitisPedigree.ped ] has 161009 identities, 524292 Non zero elements
 QUALIFIERS: !CSV !MAXIT 100 !STEP 0.001 !WORKSPACE
 Reading ../SCC_01092009_Lactation1.dat FREE FORMAT skipping 0 lines

 Univariate analysis of lnSCC
 Using 415474 records of 415474 read

 Notice: Model term leg(DIM,4) is not the same term as leg(DIM,2) but may be
          equivalent. Check the knot points reported in the .res file

  Model term Size #miss #zero MinNon0 Mean MaxNon0
   1 Sire !P 161009 0 0 4343. 0.9965E+05 0.1171E+06
  Warning: More levels found in Cow than specified
   2 Cow 43952 0 0 1 21727.3655 43952
   3 Lactation 0 0 1.000 1.000 1.000
   4 HTD 109454 0 0 1 44947.6498 109454
   5 SeasonSubClass 99 0 0 1 35.7476 99
   6 DIM 0 0 -1.003 -0.5955E-02 1.003
   7 CalvingYear 12 0 0 1 5.4059 12
   8 Het1 0221010 0.3000E-01 6.947 100.0
   9 Recomb1 0173583 0.3000E-01 6.694 50.02
  10 SCC1000 3069 0 1000. 0.1167E+06 0.1892E+08
  11 lnSCC Variate 3069 0 6.908 10.95 16.76
  12 mu 1
  13 leg(DIM,4) 5 0 0 -1.003 -.5955E-02 1.003
  14 leg(DIM,4).SeasonSubClass 495 13 leg(DIM,4): 5 5 SeasonSubClass : 99
  15 leg(DIM,2) 3 0 0 -1.003 -.5955E-02 1.003
  16 Sire.leg(DIM,2) 483027 1 Sire :***** 15 leg(DIM,2) : 3
  17 Cow.leg(DIM,2) 131856 2 Cow :43952 15 leg(DIM,2) : 3
  18 mv_estimates 3069
  26947 identity
   26947 units in section 1; initial variance 1.0791
  34332 identity
   34332 units in section 2; initial variance 0.65971
 282018 identity
  282018 units in section 3; initial variance 0.43313
  72177 identity
   72177 units in section 4; initial variance 0.34816
      3 UnStructure 0.0328 -0.0035 0.0425 0.0000 0.0000 0.0100

 161009 Ainverse
 Structure for Sire.leg(DIM,2) has 483027 levels defined
 Warning: Sire.leg(DIM,2) is a161009 x 3 interaction.
          Are you sure the G structure matches it? 3 x161009
      3 UnStructure 0.4897 0.0026 0.4971 0.0000 0.0000 0.0100

  43952 identity
 Structure for Cow.leg(DIM,2) has 131856 levels defined
 Warning: Cow.leg(DIM,2) is a 43952 x 3 interaction.
          Are you sure the G structure matches it? 3 x 43952
 Forming 727916 equations: 510 dense.
 Initial updates will be shrunk by factor 0.001
 Notice: LogL values are reported relative to a base of -320000.00
   1 LogL=-4351.60 S2= 0.0000 0 df
 Warning: - The estimation was ABORTED because a fault occured.
          The following results are probably partly wrong
               but are printed to help you locate possible causes
          The Standard Errors of variance components and solutions
               and derived quantities are probably wrong.
          Use the DL command line options to get details in the .asl file.
           Rerun using !MAXIT to stop the job earlier
               to get correct standard errors.

 Source Model terms Gamma Component Comp/SE % C
 Residual 415474 0
 HTD 109454 0 0.100000 0.100000 0.04 3 P
 Variance[ 1] 26947 0 1.07909 1.07909 0.00 0 ?
 Variance[ 2] 34332 0 0.659710 0.659710 0.00 0 P
 Variance[ 3] 282018 0 0.433130 0.433130 0.08 7 P
 Variance[ 4] 72177 0 0.348163 0.348163 0.00 0 P
 Sire.leg(DIM,2) UnStructured 1 1 0.328252E-01 0.328252E-01 0.01 0 U
 Sire.leg(DIM,2) UnStructured 2 1 -0.346642E-02 -0.346642E-02 -0.00 0 U
 Sire.leg(DIM,2) UnStructured 2 2 0.425384E-01 0.425384E-01 0.01 1 U
 Sire.leg(DIM,2) UnStructured 3 1 0.100000E-06 0.100000E-06 0.00 0 U
 Sire.leg(DIM,2) UnStructured 3 2 0.100000E-06 0.100000E-06 0.00 0 U
 Sire.leg(DIM,2) UnStructured 3 3 0.100000E-01 0.100000E-01 0.00 0 U
 Cow.leg(DIM,2) UnStructured 1 1 0.489681 0.489681 0.28 28 U
 Cow.leg(DIM,2) UnStructured 2 1 0.261273E-02 0.261273E-02 0.00 0 U
 Cow.leg(DIM,2) UnStructured 2 2 0.497089 0.497089 0.00 0 U
 Cow.leg(DIM,2) UnStructured 3 1 0.100000E-06 0.100000E-06 0.00 0 U
 Cow.leg(DIM,2) UnStructured 3 2 0.100000E-06 0.100000E-06 0.00 0 U
 Cow.leg(DIM,2) UnStructured 3 3 0.100000E-01 0.100000E-01 0.00 0 U
 Warning: Code B - fixed at a boundary (!GP) F - fixed by user
               ? - liable to change from P to B P - positive definite
               C - Constrained by user (!VCC) U - unbounded
               S - Singular Information matrix
 S means there is no information in the data for this parameter.
 Very small components with Comp/SE ratios of zero sometimes indicate poor
           scaling. Consider rescaling the design matrix in such cases.
 Covariance/Variance/Correlation Matrix UnStructured Sire.leg(DIM,2)
   0.000 -0.9277E-01 0.5519E-05
   0.000 0.000 0.4849E-05
   0.000 0.000 0.000
 Covariance/Variance/Correlation Matrix UnStructured Cow.leg(DIM,2)
   0.000 0.5296E-02 0.1429E-05
   0.000 0.000 0.1418E-05
   0.000 0.000 0.000
 Score Test 7 37.0 6 df

 Analysis of Variance NumDF F_inc
 Fault 7 Failed while ordering equations.
  Last line read was: Record Yhat Residual Hat
 SCC_01092009_Lact1_2 ../SCC_01092009_Lactation1.dat
  Model specification: TERM LEVELS GAMMAS
 Recomb1 1
 Het1 1
 CalvingYear 12
 leg(DIM,4).SeasonSubClass 495
 mu 1
 SeasonSubClass 99
 mv_estimates 3069
 Cow 43952
 HTD 109454 0.100 [ 9]
 Cow.leg(DIM,2) 131856
 Sire 161009
 Sire.leg(DIM,2) 483027
SECTIONS 26947 13 1 34332 14 2 282018 15
               3 72177 16 4
    TYPE 0 0 0 0 0 0 0 0
               9 2007 0 2005
  STRUCT 26947 0 0 0 0 1 0
           34332 0 0 0 0 1 0
          282018 0 0 0 0 1 0
           72177 0 0 0 0 1 0
               3 6 9 17 6 1 0
          161009 0 -2 0 0 1 0
               3 6 9 23 6 1 0
           43952 0 0 0 0 1 0
   18 factors defined [max 500].
   28 variance parameters [max1500]. 2 special structures
 Final parameter values [ 9: 28] 0.10000 0.0000 0.0000
   0.0000 1.0791 0.65971 0.43313 0.34816 0.32825E-01-.34664E-02
  0.42538E-010.10000E-060.10000E-060.10000E-010.48968 0.26127E-020.49709
  0.10000E-060.10000E-060.10000E-01
  Last line read was: Record Yhat Residual Hat
    18 7 415474415474 512000
 Finished: 16 Sep 2009 15:56:51.048 Failed while ordering equations.

-------------------- m2f --------------------

Sent using Mail2Forum (http://www.mail2forum.com).

Read this topic online here:
http://www.vsni.co.uk/forum/viewtopic.php?p=845#845

-------------------- m2f --------------------
Received on Fri Sep 17 2009 - 07:44:33 EST