RE: Bivariate anaysis
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
RE: Bivariate anaysis
You need to use !ASUV tif you want to fit the R structure as I cross
CORR but I do not think there is any advantage of using CORR
in this case. As Luc suggests, yu may need to specify
the number of records in the R structureArthur
Fri, 30 Nov 2001, Janss, L.L.G. wrote:
> 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
>
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
Arthur Gilmour PhD mailto:Arthur.Gilmour@agric.nsw.gov.au
Principal Research Scientist (Biometrics) fax: <61> 2 6391 3899
NSW Agriculture <61> 2 6391 3922
Orange Agricultural Institute telephone work: <61> 2 6391 3815
Forest Rd, ORANGE, 2800, AUSTRALIA home: <61> 2 6364 3288
mobile: <61> 4 1964 7675
ASREML is available by anonymous ftp from pub/aar on ftp.res.bbsrc.ac.uk
or point your web browser at ftp://ftp.res.bbsrc.ac.uk/pub/aar/
To join the asreml discussion list, send the message
subscribe
mailto:asreml-request@chiswick.anprod.CSIRO.au
To send messages to the list, mailto:asreml@chiswick.anprod.CSIRO.au
Asreml list archive: http://www.chiswick.anprod.csiro.au/lists/asreml
<> <> <> <> <> <> <>
"Blessed are the peacemakers,
for they shall be called sons of God" Jesus; Matthew 5:9
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml