preblom with repeated
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
preblom with repeated
Arthur,
I am doing a repeated measures analysis. There are problems.
the .as file read like thsis
Analysis of bivaraite on FEC and wwt of lambs from cowra
anim !P
sire !P
dam !P
sirebr 5 !A
dambr 3 !A
sex 2
cgp 40 !A
br 3 !A
fage1 #age at 1st sampling
fec1
fecln1
fecst1
fecct1
fage2
fec2
fecln2
fecst2
fecct2
wage
wwt
fecwwt.ped
asfecwwt_lamb.dat
fecct1 fecct2 ~ Trait at(Trait,1).fage1 at(Trait,2).fage2 !f Trait.sex
Trait.sirebr !r anim units
!f Trait.cgp Trait.mv
1 2 0 !ASUV
anim #no of animals in the data
Trait
The error is invalid model factor Trait.mv.
following is the asr file.
ASREML [28 Oct 1998] Analysis of bivaraite on FEC and wwt of lambs from
cowra
Fri Mar 26 11:50:53 1999 16.00 Mbyte feccbrt12
A-inverse retrieved from ainverse.bin
PEDIGREE [fecwwt.ped ] has 2189 identities, 6021 Non zero elements
Reading asfecwwt_lamb.dat FREE FORMAT skipping 0 lines
Bivariate analysis of fecct1 and fecct2
Using 1421 records [of 1421 read from 1421 lines of
asfecwwt_lamb.dat ]
Model term Size Type COL Minimum Mean Maximum #zero
#miss
1 anim 2189 Direct 1 3.000 1168. 2189. 0
0
2 sire 2189 Direct 2 2.000 553.3 1283. 3
0
3 dam 2189 Direct 3 1.000 850.5 2188. 2
0
4 sirebr 5 Factor 4 1 2.5550 4 0
3
5 dambr 3 Factor 5 1 1.3340 2 0
2
6 sex 2 Factor 6 1 1.6735 2 0
0
7 cgp 40 Factor 7 1 20.9993 40 0
0
8 br 3 Factor 8 1 2.0934 3 0
8
9 fage1 1 Covariat 9 64.00 109.0 158.0 0
71
10 fec1 1 Covariat 10 10.00 652.5 5500. 0
70
11 fecln1 1 Covariat 11 2.303 5.906 8.610 0
70
12 fecst1 1 Covariat 12 3.162 22.75 74.16 0
70
13 fecct1 1 Variate 13 2.154 7.756 17.60 0
70
14 fage2 1 Covariat 14 113.0 138.6 153.0 0
498
15 fec2 1 Covariat 15 10.00 168.3 1800. 0
490
16 fecln2 1 Covariat 16 2.303 4.150 7.500 0
490
17 fecst2 1 Covariat 17 3.162 10.56 42.43 0
490
18 fecct2 1 Variate 18 2.154 4.534 12.13 0
490
19 wage 1 Covariat 19 178.0 192.5 218.0 0
9
20 wwt 1 Covariat 20 9.500 25.41 44.00 0
9
21 Trait 2 Traits/Variat
22 at(Trait,1) 1 True If trait 1
23 at(Trait,1). 1 Interaction 22 at(Tra: 1 9 fage1 : 1
24 at(Trait,2) 1 True If trait 2
25 at(Trait,2). 1 Interaction 24 at(Tra: 1 14 fage2 : 1
26 Trait.sex 4 Interaction 21 Trait : 2 6 sex : 2
27 Trait.sirebr 10 Interaction 21 Trait : 2 4 sirebr : 5
28 units 1421 Units
29 Trait.cgp 80 Interaction 21 Trait : 2 7 cgp : 40
30 mv_estimates 0 Missing value
31 Trait.mv 0 Covariance 21 Trait : 2 30 mv_estimates : 0
ault 0 Invalid model factor Trait.mv
Last line read was: 963899 25 12713 2 2 2 2218 2 73 1700 7.44 41.23 11.91
122 700 6.55 26.46 8
.86 209 27.00
31 0 1421 1421 2000
ecct1 asfecwwt_lamb.dat
Model specification: TERM LEVELS GAMMAS
rait.sirebr 10 0.00 Trait.sex 4 0.00 at(Trait,2). 1 0.00
t(Trait,1). 1 0.00 Trait 2 0.00 anim 2189 0.10
nits 1421 0.10 Trait.cgp 80 0.00
SECTIONS 0 0 0
31 factors define [max 500].
0 variance parameters [max 550]. 2 special structures
inal parameter values
ARNING: Extra lines on the end of the input file are ignored from
nim
Finished: Fri Mar 26 11:50:58 1999 Invalid model factor Trait.mv
HOWEVER WHEN I CHANGE the missing value model factor to mv. The error
message change to this.
A-inverse retrieved from ainverse.bin
PEDIGREE [fecwwt.ped ] has 2189 identities, 6021 Non zero elements
Bivariate analysis of fecct1 and fecct2
139 records in the data
0 2189 139
R structures imply 0 + 2189 records: only 139 exist
Last line read was: anim 0 0 0 0 0
30 1 1421 139 2000
fecct1 asfecwwt_lamb.dat
Model specification: TERM LEVELS GAMMAS
Trait.sirebr 10 0.00 Trait.sex 4 0.00 at(Trait,2). 1 0.00
at(Trait,1). 1 0.00 Trait 2 0.00 anim 2189 0.10
units 1421 0.10 Trait.cgp 80 0.00 mv_estimates 70 0.00
SECTIONS 2189 10 1
STRUCT 2189 0 0 0 0 1 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 0
0 0 0 0 0 0 0 0
30 factors define [max 500].
10 variance parameters [max 550]. 2 special structures
Final parameter values 0.10000 0.10000 0.0000
0.0000 -.10000E-36
Fault 1 R structures imply 0 + 2189 records: only 139
exis
t
WARNING: Extra lines on the end of the input file are ignored from
Trait
Finished: Fri Mar 26 12:18:02 1999 R structures imply 0 + 2189
records: only 139 exis
t
the data file pattern if followed. Could you help me check what is wrong
there.
Regards,
Yuandan
962681 37 10373 3 1 2 2205 1 99 400 5.99 20 7.35 NA NA NA NA NA 183 17.00
962786 36 10122 3 1 1 2218 1 95 1400 7.24 37.42 11.16 NA NA NA NA NA 187 23.00
962809 37 10739 3 1 1 2217 3 94 200 5.3 14.14 5.84 NA NA NA NA NA 188 17.50
962824 28 10613 3 1 1 2205 1 94 1200 7.09 34.64 10.6 NA NA NA NA NA 188 24.50
962927 4 12648 2 2 1 2109 3 91 700 6.55 26.46 8.86 NA NA NA NA NA 191 11.00
962940 3 12717 4 2 2 2106 3 91 200 5.3 14.14 5.84 NA NA NA NA NA 191 33.50
...
963833 12 10674 3 1 2 2217 2 NA NA NA NA NA 129 600 6.4 24.49 8.42 202 19.50
963869 40 10110 1 1 1 2216 3 NA NA NA NA NA NA NA NA NA NA 204 21.00
963870 36 10679 3 1 2 2200 NA NA NA NA NA NA NA 200 5.3 14.14 5.84 NA NA
963871 36 10679 3 1 2 2218 2 NA NA NA NA NA 127 100 4.62 10 4.63 204 15.00
963889 25 12806 2 2 1 2218 3 NA NA NA NA NA 124 100 4.62 10 4.63 207 17.00
963848 12 10424 3 1 1 2217 2 79 400 5.99 20 7.35 128 10 2.3026 3.162 2.154
203 19.50
963853 41 10466 1 1 2 2216 2 79 2400 7.78 48.99 13.35 128 200 5.3 14.14
5.84 203 24.00
963854 12 10319 3 1 2 2205 3 79 1900 7.55 43.59 12.35 128 200 5.3 14.14
5.84 203 21.50
963855 12 10319 3 1 2 2205 3 79 100 4.62 10 4.63 128 10 2.3026 3.162 2.154
203 17.00
963856 42 10881 1 1 2 2205 1 79 300 5.71 17.32 6.68 128 100 4.62 10 4.63
203 22.00
963857 40 10809 1 1 1 2205 1 78 600 6.4 24.49 8.42 127 100 4.62 10 4.63 204
25.00
963858 25 10830 2 1 2 2205 1 78 100 4.62 10 4.63 127 10 2.3026 3.162 2.154
204 24.50
963859 41 10063 1 1 2 2205 1 78 300 5.71 17.32 6.68 127 500 6.22 22.36 7.92
204 24.00
A-inverse retrieved from ainverse.bin
PEDIGREE [fecwwt.ped ] has 2189 identities, 6021 Non zero elements
Bivariate analysis of fecct1 and fecct2
9 31 0.100000001490116 0
Invalid model factor Trait.mv
Last line read was: 963899 25 12713 2 2 2 2218 2 73 1700 7.44 41.23
11.91 122 700 6.55 26.46 8
.86 209 27.00
31 0 1421 1421 2000
fecct1 asfecwwt_lamb.dat
Model specification: TERM LEVELS GAMMAS
Trait.sirebr 10 0.00 Trait.sex 4 0.00 at(Trait,2). 1 0.00
at(Trait,1). 1 0.00 Trait 2 0.00 anim 2189 0.10
units 1421 0.10 Trait.cgp 80 0.00
SECTIONS 0 0 0
31 factors define [max 500].
0 variance parameters [max 550]. 2 special structures
Final parameter values
WARNING: Extra lines on the end of the input file are ignored from
anim
Finished: Fri Mar 26 11:50:58 1999 Invalid model factor Trait.mv
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml