Re: Bivariate anaysis + maternal genetic effect
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Bivariate anaysis + maternal genetic effect
On Thu, 6 Dec 2001 Abbas.Pakdel@alg.vf.wau.nl wrote:
>
> 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
The preceeding line should read
1 2 1
> 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
>
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
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