Further to my comment yesterday, I add the following:
>
> >
> > Hi,
> >
> > I want to get predictions for the subrace effects for four different
> > sites (traits), using the following code:
> >
> > !PART 21 #Full multivariate Four sites
> > !maxit 12
> > massflow bawdflow shalflow sso3flow ~ Trait Trait.subrace at(Trait,1).rep,
> > at(Trait,2).bigrep at(Trait,3).rep,
> > !r at(Trait,2).plot 38 -at(Trait,3).rep.block,
> > Trait.csiroID !f mv
> > 1 2 1 !STEP .1
> > 0 0 0
> > Trait 0 US !+10 !GUFUFFUFFFP
> > 734
> > 0 346
> > 0 0 192
> > 0 0 0 913
> >
> > Trait.csiroID 2
> > Trait 0 CORR !+10 !GUUUFFUUUUU
> > 1.08
> > 0.89 0.68
> > 0 0 1.11
> > 267 335 174 245
> > csiroID
> >
> > predict subrace
> >
> > The run converges but it does not produce a .pvs file with the
> > predictions for subrace, although there are solutions for Trait and
> > Trait.subrace in the .sln file.
> >
> > Unable to predict at(Trait,3) 4 0 0 1 -13 0 -1
> 3
> > LogL=-4157.71 S2= 1.0000 5063 df
> >
> > Not all subraces are in all sites. There are 11, 11, 18 and 8 in each
> > site respectively. What is the problem?
> >
> >
> > Luis
> >
>
> Dear Luis,
>
> The problem probably is that you will need to be more specific
> about the prediction.
>
> The fixed terms in the model include Trait, subrace, rep and bigrep
> and mv.
>
> predict subrace
> requests ASREML to average over the other 4 dimensions of a 5way
> table but I gather this table has many gaps so the average is undefined.
>
> If you made the rep and bigrep terms random and used the statement
>
> predict Trait subrace
>
> you will have more chance of getting subrace predictions for each Trait
> (those where the subrace occurs for that trait).
>
> If you want to keep the current model, try the predict command
>
> predict Trait subrace !present Trait subrace rep bigrep
>
> which should form a Trait x subrace table averaging over rep and bigrep.
>
> This however may not work because rep and bigrep only appear in an interaction
> as indeed subrace and ASREML does not yet know how to extract prescence
> information
> from an interaction.
>
> Arthur
On further reflection,
predict Trait subrace
should work (with the rep terms fixed or eandom) provided
if rep is fixed, it has a complete set of levels for traits 1 and 3
if bigrep is fixed, it has a complete set of levels.
Arthur
>
> <><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
> 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 6362 0046
> Cargo: <61> 2 6364 3288
> I havn't finished building my house at Cargo but am back at work.
>
> 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
> 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
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
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 6362 0046
Cargo: <61> 2 6364 3288
I havn't finished building my house at Cargo but am back at work.
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
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