regression analysis(2)

Hi,
i want to run a random regression model, but really i cant get result any time,
i had prepared the data in the single trait design (i have 13 repeated measurment and i set then beneath each other as:
animal dam sare..... ASSESS RECORD
1 1 12
1 2 14
1 3 22
.....
1 13 10
,
and i had used !repeat qualifier in pedigree line,
but with this as file i cant reach the desire,
could any one help me, kindly? :)

as fiel:

 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 0 AINV

asr fil:

  SIRE !P
  DAM !P
  RECORD !M -99
  Reading pedigree file INDEXPED.TXT : skipping 0 lines
 PEDIGREE [INDEXPED.TXT ] has 2518 identities, 6896 Non zero elements
 QUALIFIERS: !MAXIT 20
  Reading alllnew.txt FREE FORMAT skipping 0 lines
 Univariate analysis of RECORD
 Using 24609 records of 24609 read
  Model term Size Type COL Minimum Mean Maximum #zero #miss
   1 ANIM 2518 Direct 1 3.000 1383. 2518. 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,2 7554 Interaction 8 pol(ASSE: 3 1 ANIM : 2518
  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

 Fault 1 G structure size does not match Ainverse
  Last line read was: 0 0 0 0 0 0 0 0 0
 30 alllnew.txt
  Model specification: TERM LEVELS GAMMAS
 BDCLASS 3 0.000
 HATCH 3 0.000
 ASSESS 13 0.000
 pol(ASSESS,2).ANIM 7554 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
             581 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 1000
 Finished: 16 Sep 2002 13:13:29.503 G structure size does not match Ainverse
**********************************8
thanks of you!
maryam

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

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

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

-------------------- m2f --------------------
Received on Sat Feb 13 2010 - 09:50:36 EST