>
> Andrew
>
> The problem is with your model specification. Tr|1/ide(dam) is now written
> as at(Tr,1).ide(dam).
>
> Regards
> Neil Cullen
I doubt this would fix your problem. The old syntax is still permitted
in simple cases as here.
>
> > -----Original Message-----
> > From: Andrew Swan [SMTP:andrew@chiswick.anprod.csiro.au]
> > Sent: Wednesday, 3 February 1999 11:28
> > To: asreml@ram.chiswick.anprod.csiro.au
> > Subject: two trait maternal model
> >
> > hi - can anyone tell me what's wrong with the following syntax? i am
> > trying to fit a bivariate model including direct and maternal genetic
> > effects for both traits (and a p.e. effect for gfw10):
> >
> > gfw10 cfw10 ~ Trait Tr.line Tr.damage Tr.brtype,
> > !r Tr.animal Tr.dam Tr|1/ide(dam) 0.0858 !f Tr.flock
> >
> > 1 2 1
> > 0
> > 2 0 US !+3
> > 0.0510783
> > 0.037 0.0420101
> >
> > Tr.animal 2
> > 4 0 US !+10
> > 0.0367198
> > -0.00801379 0.00681885
> > 0.0267 -0.00576 0.0304376
> > -0.0085 0.0064 -0.0101964 0.0123053
> > 9870 0 0
> >
> > problem is i get a "Structure / Factor mismatch " message
> >
> > (i think i have also ordered the starting values in the g-structure
> > incorrectly - should be gfw10_d, gfw10_m, cfw10_d, cfw10_m ??)
> >
> > thanks
> > --
> > Andrew Swan
I can't see an obvious error. The 'last line read' should giv a clue
in conjunction with the .asr file as to how far through ASREML has got.
[The last line read will not include the initial values covered by !+3 and !+10
as these are read directly.
[ASREML reports each G structure as it reads it in]
I suggest replacing
9870 0 0
with
animal 0 -1
Also, in the factor definitions make sure both animal and dam have !P
specified. Otherwise, dam will not be coded as per animal and will not
agree in size so that the structure size for Tr.ani Tr.dam will not
be 4 times 9870 but something else.
I trust this helps
Arthur
[It reports each G structure as it reads it in]
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
Arthur Gilmour PhD email: Arthur.Gilmour@agric.nsw.gov.au
Senior 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 6362 0046
ASREML is still free 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
to asreml-request@chiswick.anprod.CSIRO.au
The address for messages to the list is asreml@chiswick.anprod.CSIRO.au
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml
<> <> <> <> <> <> <>
"You (Father) have given Him (Your Son) authority over all flesh,
that He should give eternal life to as many as You have given Him.
And this is eternal life,
that they may know You, the only true God,
and Jesus Christ whom You have sent." John 17: 2-3.
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml