Unable to form R inverse component
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Unable to form R inverse component
Hi,
While I analyse covariance between two trait, I encounter an error. Even
refering manual page 102-103, I still can't solve this problem. Can some
one help me ?
Thanks
PS.
.as and .asr files:
Animal model for FEC data from Cowra 95 and 96
animal !P #1
sire #2
dam !P #3
dambd 2 #4
sirebd 4 #5
mategp 2 #6
sex 2 #7
pad 18 #8
brrank 7 #9
brank 6 #10
rrank 3 #11
dob #12
age_sam #13
age_wean #14
fec210
.......
cowrawt.txt
cowrawt.txt !MAXIT 20
fec210 wwt ~Trait Tr.sirebd Tr.sex Tr.year Tr.brank Tr.rrank,
Tr|2.age_wean Tr|1.age_sam !r Tr.animal Tr.dam Tr.ide(dam),
!f Tr.dambd Tr.sirebd.dambd Tr.year.pad Tr.year.mategp
1 2 1 !STEP .1
0
2 0 US !+3 !GPUP
24.1094 1 8.67232
Tr.animal 2
4 0 US !+10 !GPUPUUPUUUP
17.0071
.0001 2.89729
.0001 .0001 .000143501E-03
.0001 .0001 .0001 2.30679
animal 0 0
ASREML [20 Feb 1998] Animal model for FEC data from Cowra 95 and 96
07/15/98 10:14:04.24 8.00 Mbyte fecwwt.asr
A-inverse retrieved from ainverse.bin
PEDIGREE [cowrawt.txt ] has 2757 animals, 7830 Non zero elements
QUALIFIERS: !MAXIT 20
Reading cowrawt.txt FREE FORMAT skipping 0 lines
Bivariate analysis of fec210 and wwt
1916 records obtained from 1920
lines of cowrawt.txt 1916 records
retained
Model term Size Type COL Minimum Mean Maximum #
zero #miss
1 animal 2757 Direct 1 3.0000 1496.1 2757.0
0. 0
2 sire 1 Covariate 2 3.0000 23.199 42.000
0. 0
3 dam 2757 Direct 3 1.0000 1018.7 2749.0
0. 0
4 dambd 2 Factor 4 1.0000 1.3314 2.0000
0. 0
...
31 Trait 2 Traits/Variat
32 Tr.sirebd 8 Interaction 31 Trait : 2 5 sirebd : 4
33 Tr.sex 4 Interaction 31 Trait : 2 7 sex : 2
34 Tr.year 4 Interaction 31 Trait : 2 19 year : 2
35 Tr.brank 12 Interaction 31 Trait : 2 10 brank : 6
36 Tr.rrank 6 Interaction 31 Trait : 2 11 rrank : 3
37 Tr_2.age_wea 1 Conditional 31 Trait : 2 14 age_wean : 1
38 Tr_1.age_sam 1 Conditional 31 Trait : 1 13 age_sam : 1
39 Tr.animal 5514 Interaction 31 Trait : 2 1 animal : 2757
40 Tr.dam 5514 Interaction 31 Trait : 2 3 dam : 2757
41 ide(dam) 2757 Factor 3 1.0000 1018.7 2749.0
0. 0
42 Tr.ide(dam) 5514 Interaction 31 Trait : 2 41 ide(dam) : 2757
43 Tr.dambd 4 Interaction 31 Trait : 2 4 dambd : 2
44 sirebd.dambd 8 Interaction 5 sirebd: 4 4 dambd : 2
45 Tr.sirebd.da 16 Interaction 31 Trait : 2 44 sirebd.dam: 8
46 year.pad 36 Interaction 19 year : 2 8 pad : 18
47 Tr.year.pad 72 Interaction 31 Trait : 2 46 year.pad : 36
48 year.mategp 4 Interaction 19 year : 2 6 mategp : 2
49 Tr.year.mate 8 Interaction 31 Trait : 2 48 year.mateg: 4
1916 identity
2 US=UnStr 24.11 1.00 8.67
4 US=UnStr 17.01 0.00 2.90 0.00 0.00 0.00 0.00
0.00
0.00 2.31
2757 Ainverse
Structure of Tr.animal has 11028 levels defined
Forming 16681 equations: 39 dense
Initial updates will be shrunk by factor 0.100
Fault 1009Unable to form R inverse component
Last line read was: animal 0 0 0 0 0 0 0
49 1009 1916 3832 1000
fec210 cowrawt.txt
Model specification: TERM LEVELS GAMMAS
Tr_1.age_sam 1 0.00 Tr_2.age_wea 1 0.00 Tr.rrank 6 0.00
Tr.brank 12 0.00 Tr.year 4 0.00 Tr.sex 4 0.00
Tr.sirebd 8 0.00 Trait 2 0.00 Tr.animal 5514 0.00
Tr.dam 5514 0.00 Tr.ide(dam) 5514 0.10 Tr.dambd 4 0.00
Tr.sirebd.da 16 0.00 Tr.year.pad 72 0.00 Tr.year.mate 8 0.00
SECTIONS 3832 16 1
TYPE 0 0 0 0 0 0 0 0
203 -1 11 0 0 0 0
STRUCT 1916 0 0 0 0 1 0 0
2 3 9 17 3 2 0 0
4 10 9 20 10 2 0 0
2757 0 -2 0 0 1 2 0
Final parameter values 0.10000 0.00000 0.00000
0.00000 0.00000 1.00000 24.10940 1.00000 8.67232 17.00710
0.00010 2.89729 0.00010 0.00010 0.00000 0.00010 0.00010
0.00010 2.30679
Finished: 10:16:01.78 Unable to form R inverse component
************************************************************************
Yuandan Zhang
Department of Animal Science
University of New England
Armidale, NSW 2351, Australia
Phone +61 + 2 + 6773 2756(Lab)
Fax +61 + 2 + 6773 3275
E-mail yzhang@metz.une.edu.au
homepage http://metz.une.edu.au/~yzhang/
************************************************************************