Random Regression

From: maryam <asremlforum_at_VSNI.CO.UK>
Date: Sat, 13 Feb 2010 09:29:33 +0000

Dear All,
i want to use ASReml to run a random regression model for mothly egg production for hens, i have 13 measurment for duration, i had
prepare a as file but every time that i run it i cant get results,
can anybody help me to do this and tell my mistakes?
with the best wishes,
Maryam
********* my as file******
 ANIM 24609
 SIRE !P
 DAM !P
 HATCH 3
 BDCLASS 3
 ASSESS 13
 RECORD !M -99

INDEXPED.TXT !REPEAT !MAKE
alllnew.txt !MAXIT 20
 RECORD ~ ASSESS HATCH BDCLASS !r pol(ASSESS,2).ANIM !f mv
1 2 1 !ASMV 13 !STEP .1
1893
ASSESS 0 US !+91 !GP
.09
0 0
0 0 0
0 0 0 0
0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 .005 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 .09
 

pol(ASSESS,2).ANIM 2
pol(ASSESS,2) 0 US !+91 !GP

.008
0 0
0 0 0
0 0 0 0
0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 .0004 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 .004 0 0 0 0 0 0 0 0 0 0 .006

ANIM
************************and asr file****************

  1 ANIM 24609 Factor 1 1 1399.8061 2762 0 0
   2 SIRE 2518 Direct 2 2.000 88.54 1996. 0 0
   3 DAM 2518 Direct 3 1.000 730.4 2516. 0 0
   4 HATCH 3 Factor 4 1 1.9519 3 0 0
   5 BDCLASS 3 Factor 5 1 1.9884 3 0 0
   6 ASSESS 13 Factor 6 1 7.0000 13 0 0
   7 RECORD 1 Variate 7 1.000 11.87 33.00 5631 564
   8 pol(ASSESS, 3 Polynoml 6 1.000 7.000 13.00 0 0
   9 pol(ASSESS,273827 Interaction 8 pol(ASSE: 3 1 ANIM :24609
  10 mv_estimates 564 Missing value
   1893 identity
 Warning: US matrix was not positive definite: MODIFIED
     13 UnStructure 0.0900 1.4730 38.1437 0.5240 1.8638 49.8811
    0.8501 0.9453 2.6358 61.9508 0.7847 0.7166 1.5797 4.9878
   59.2439 0.6499 0.7323 1.3377 3.8417 4.5062 55.9215 0.1038
    0.4558 0.1638 -0.3175 -0.2282 0.2861 31.7640 0.1427 0.4200
    0.2648 0.2352 0.2886 0.6830 2.0387 36.0837 0.3653 0.5287
    0.7054 1.3444 1.3873 1.5470 1.1976 1.9858 36.2577 0.1932
    0.2139 0.1714 0.2440 0.2650 0.5457 1.5276 2.0085 2.4250
   39.1734 0.2633 0.2212 0.3140 0.3616 0.3559 0.6158 1.1374
    1.6070 2.1700 2.7885 39.0620 0.1303 0.3740 0.1997 0.2358
    0.2223 0.2608 1.3695 1.5668 1.8962 2.4675 2.7158 45.0518
    0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
    0.0000 0.0000 0.0000 0.0000 0.0900
   24609 records assumed sorted 13 within 1893
 Warning: US matrix was not positive definite: MODIFIED
      3 UnStructure 0.0080 0.0581 0.7303 0.0643 0.1447 0.8944

   5679 identity
 Fault 1 Structure / Factor mismatch
  Last line read was: 0 0 0 0 0 0 0 0 0
 29 alllnew.txt
  Model specification: TERM LEVELS GAMMAS
 BDCLASS 3 0.000
 HATCH 3 0.000
 ASSESS 13 0.000
 pol(ASSESS,2).ANIM 73827 0.100
 mv_estimates 564 0.000
SECTIONS 24609 6 1
    TYPE 0 0 0 2003 0
  STRUCT 1893 0 0 0 0 1 0
              13 91 9 7 0 1 0
               3 6 9 98 0 1 0
            5679 0 0 0 0 1 0
   10 factors defined [max 500].
  103 variance parameters [max 600]. 2 special structures
 Final parameter values 0.10000 0.0000 1.0000
  0.90000E-01 1.4730 38.144 0.52404 1.8638 49.881 0.85006
  0.94531 2.6358 61.951 0.78472 0.71662 1.5797 4.9878
   59.244 0.64987 0.73233 1.3377 3.8417 4.5062 55.922
  0.10378 0.45584 0.16381 -.31748 -.22818 0.28608 31.764
  0.14271 0.42004 0.26482 0.23518 0.28860 0.68303 2.0387
   36.084 0.36527 0.52871 0.70542 1.3444 1.3873 1.5470
   1.1976 1.9858 36.258 0.19315 0.21388 0.17145 0.24397
  0.26500 0.54574 1.5276 2.0085 2.4250 39.173 0.26332
  0.22118 0.31401 0.36160 0.35595 0.61576 1.1374 1.6070
   2.1700 2.7885 39.062 0.13027 0.37397 0.19968 0.23582
  0.22229 0.26078 1.3695 1.5668 1.8962 2.4675 2.7158
   45.052 0.78588E-070.35505E-05-.46220E-05-.15366E-05-.25507E-05-.52937E-06
  -.57371E-05-.56535E-050.22441E-050.41586E-050.46457E-050.40795E-050.90000E-01
  0.80000E-020.58146E-010.73030 0.64349E-010.14475 0.89443
  Last line read was: 0 0 0 0 0 0 0 0 0
    10 1 24609 24609 8000
 Finished: 16 Sep 2002 12:52:52.103 Structure / Factor mismatch

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

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

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

-------------------- m2f --------------------
Received on Sat Feb 13 2010 - 09:29:33 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.)