Bivariate anaysis + maternal genetic effect
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Bivariate anaysis + maternal genetic effect
Dear asreml users
Hi
Thanks for last reply, But I have another problem.I would like to set up an
ASREML to do a bivariate analysis for traits which some of them, maternal
genetic effect have effect .
I set up the following . as file .
In this model GEW is body weight that maternal effect has effect on this trait
and for the second trait (HCT) maternal genetic effect dose not any effect.
But in . asr file every time I see this message : G structure header: Factor,
order
Any help greatly appreciated.
Best regards
abbas
The .as file
Bivariate Analysis for hct & gew
IND !P
SIRE
DAM !P
IHD 46 !I
SEX 2 !I
K 2 !I
AFDEL 10 !I
HCT !M -99
GEW !M -99
RV !M -99
TV !M -99
RATIO !M -99
PERRV !M -99
PERTV !M -99
UITV !M -99
VOCHT !M -99
BORST !M -99
LEVER !M -99
HART !M -99 # DEPENDENT VARIABLES
ali2.prn !MAKE #pedigree file
ali1.prn !MAXIT 250 #ascites data
GEW HCT ~ Trait Tr.SEX Tr.K Tr.AFDEL Tr.IHD !r Tr.IND at(Tr,1).DAM !f mv
1 2 2
0 0 0
2 0 US 4.14 -2.58 8.22 !GP # residual variance
Tr.IND 2
3 0 DIAG 1.24 6.86 0.22 !GP # Additive variance 1,2 and maternal effect
respectively.
5327 0 0
The .asr file:
ASREML [11 Jun 1999] Bivariate Analysis for hct & gew
Thu Dec 6 11:43:49 2001 8.00 Mbyte mhct1
Reading ali2.prn : skipping 0 lines
PEDIGREE [ali2.prn ] has 5327 identities, 16760 Non zero elements
QUALIFIERS: !MAXIT 250
Reading ali1.prn FREE FORMAT skipping 0 lines
Bivariate analysis of GEW and HCT
Using 4202 records [of 4202 read from 4202 lines of ali1.prn
]
Model term Size Type COL Minimum Mean Maximum #zero
#miss
1 IND 5327 Direct 1 895.0 3102. 5327. 0
0
2 SIRE 1 Covariat 2 0.9413E+08 0.9420E+08 0.9429E+08 0
0
3 DAM 5327 Direct 3 66.00 453.2 892.0 0
0
4 IHD 46 Factor 4 1 22.6675 46 0
0
5 SEX 2 Factor 5 1 1.5869 2 0
0
6 K 2 Factor 6 1 1.8377 2 0
0
7 AFDEL 10 Factor 7 1 4.9912 10 0
0
8 HCT 1 Variate 8 23.00 35.40 55.00 0
655
9 GEW 1 Variate 9 3.540 16.04 25.63 0
509
10 RV 1 Covariat 10 0.4000 1.955 5.300 0
542
11 TV 1 Covariat 11 2.300 6.968 11.80 0
544
12 RATIO 1 Covariat 12 8.000 27.94 60.00 0
544
13 PERRV 1 Covariat 13 2.747 12.47 51.50 0
556
14 PERTV 1 Covariat 14 2.238 4.393 9.002 0
558
15 UITV 1 Covariat 15 1.000 0.1600 1.000 2095
1708
16 VOCHT 1 Covariat 16 1.000 0.7763E-01 2.000 3542
505
17 BORST 1 Covariat 17 1.000 0.3084E-01 2.000 3585
505
18 LEVER 1 Covariat 18 1.000 0.6708E-01 2.000 3488
505
19 HART 1 Covariat 19 1.000 0.5890 2.000 1781
506
20 Trait 2 Traits/Variat
21 Tr.SEX 4 Interaction 20 Trait : 2 5 SEX :
2
22 Tr.K 4 Interaction 20 Trait : 2 6 K :
2
23 Tr.AFDEL 20 Interaction 20 Trait : 2 7 AFDEL :
10
24 Tr.IHD 92 Interaction 20 Trait : 2 4 IHD :
46
25 Tr.IND 10654 Interaction 20 Trait : 2 1 IND :
5327
26 at(Tr,1) 1 True If trait 1
27 at(Tr,1).DAM 5327 Interaction 26 at(Tr,1): 1 3 DAM :
5327
28 mv_estimates 1164 Missing value
4202 identity
2 US=UnStr 4.14 -2.58 8.22
8404 records assumed sorted 2 within 4202
3 DIAG=Dia 1.24 6.86 0.22
5327 Ainverse
Structure of Tr.IND has 15981 levels defined
Fault 1 G structure header: Factor, order
GEW ali1.prn
Model specification: TERM LEVELS GAMMAS
Tr.IHD 92 0.00 Tr.AFDEL 20 0.00 Tr.K 4 0.00
Tr.SEX 4 0.00 Trait 2 0.00 Tr.IND 10654 0.00
at(Tr,1).DAM 5327 0.00 mv_estimates 1164 0.00
SECTIONS 8404 9 1
TYPE 0 0 0 0 0 2003 -1 0
STRUCT 4202 0 0 0 0 1 0 0
2 3 9 10 0 1 0 0
3 3 8 13 0 1 0 0
5327 0 -2 0 0 1 2 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 0 0 0 0 0 0 0
28 factors define [max 500].
15 variance parameters [max 550]. 2 special structures
Final parameter values 0.0000 0.0000 0.0000
1.0000 4.1400 -2.5800 8.2200 1.2400 6.8600 0.22000
Last line read was:
28 1 4202 8404 1000
Finished: Thu Dec 6 11:44:01 2001 G structure header: Factor, order
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml