RE: Bivariate anaysis
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
RE: Bivariate anaysis
It looks to me that the !ASUV is inappropriate. It messes-up the default
sizes (ASREML reports to expect double the size of what is actually
present).
> Luc
> Dr. L.L.G. Janss
> ID-Lelystad Institute for Animal Science and Health
> P.O. Box 65, 8200 AB Lelystad, The Netherlands
> phone: + 31 (0) 320 238 259, fax: + 31 (0) 320 238 050
> e-mail: l.l.g.janss@id.wag-ur.nl, internet: www.id.wageningen-ur.nl
> visiting address: Edelhertweg 15, Lelystad
> This message is intended exclusively for the addressee. It may contain
> information that is confidential. Any use or publication of this e-mail
> message without permission of ID-Lelystad is not allowed. If you are not
> the intended recipient please notify us and destroy this message.
>
>
> -----Original Message-----
> From: Abbas.Pakdel@alg.vf.wau.nl [SMTP:Abbas.Pakdel@alg.vf.wau.nl]
> Sent: Thursday, November 29, 2001 3:14 PM
> To: asreml@chiswick.anprod.csiro.au
> Subject: 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
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml