Bivariate anaysis
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Bivariate anaysis
Dear artour and asreml users
Hi
I would like to set up an ASREML to do a bivariate analysis . I have some
normal traits and some score traits. I run univariate analysis for all traits
atfirst , But I don't know why the following model did not run.
Any help greatly appreciated.
The .as file
Bivariate Analysis for hct & gew
IND !P
SIRE
DAM
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 50 !ASUV #ascites data
HCT RV ~ Trait Tr.SEX Tr.K Tr.AFDEL Tr.IHD !r Tr.IND !f mv
1 2 1
0
Trait 0 CORR 0.1 13.5 5.8
Tr.IND 2
Trait 0 CORR 0.4 4.12 1.24
IND
and the .asr file is :
ASREML [11 Jun 1999] Bivariate Analysis for hct & gew
Thu Nov 29 15:11:11 2001 8.00 Mbyte bi1
Reading ali2.prn : skipping 0 lines
PEDIGREE [ali2.prn ] has 5327 identities, 16760 Non zero elements
QUALIFIERS: !MAXIT 50 !ASUV
Reading ali1.prn FREE FORMAT skipping 0 lines
Bivariate analysis of HCT and RV
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 1 Covariat 3 0.9413E+08 0.9420E+08 0.9429E+08 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 Covariat 9 3.540 16.04 25.63 0
509
10 RV 1 Variate 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 mv_estimates 1197 Missing value
8404 identity
2 CORRelat 0.10 13.50 5.80
Fault 1 R structures imply 0 + 16808 records: only 8404 exist
HCT 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.10
mv_estimates 1197 0.00
SECTIONS 16808 8 1
TYPE 0 0 0 0 0 0 0
STRUCT 8404 0 0 0 0 1 0 0
2 3 5 9 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
26 factors define [max 500].
11 variance parameters [max 550]. 2 special structures
Final parameter values 0.10000 0.0000
-.10000E-36
0.10000 13.500 5.8000
Last line read was: Trait 0 CORR 0.1 13.5 5.8 0 0 0 0 0
26 1 4202 8404 1000
Finished: Thu Nov 29 15:11:23 2001 R structures imply 0 + 16808
records
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml