Re: Spline predict points]

David,

Thanks for the tip. Moving 'predictpoints' to the asreml call did the trick.

I also tried your other suggestion and downloaded the latest version from:
http://www.vsni.co.uk/downloads/asreml/as2windows.php?data=upgrade

But 'predictpoints' in the 'predict' call didn't work. See the
version info below.

Did I find the latest version, or is it somewhere else?

Best,

Kevin

R> asreml.About()
Mixed Models for S Language Environments
About asreml(): Mon Oct 29 09:17:32 2007

asreml() functions visible in: c:/R/R-2.6.0/library/asreml/
asreml() version: 2.00
   Build date: 28 August 2007
ASReml object library: 2.00bb
   Build date: 28 August 2007

On 10/28/07, Butler, David <David.Butler_at_dpi.qld.gov.au> wrote:
> Hi Kevin,
>
> predict() uses update() and at some point there was a bug in update().
> either trying the latest version or shifting the predictpoints argument
> to the asreml() call from predict() should work.
>
> Dave.
>
> On Sat, 2007-10-27 at 04:39, Kevin Wright wrote:
> > I'm exploring the use of splines with a mixed model in asreml-r
> > following the example of the orange tree growth for Tree 1.
> >
> > Page 117 of ASREML-R manual says:
> > "The predictpoints argument adds the nominated points to the design
> > matrix for prediction purposes (Figure 8.13)."
> >
> > Does this mean that I will obtain predictions at the design points AND
>
> > the predictpoints, or only at predictpoints, or something else?
> >
> > Here's a simple example:
> >
> > o1 <- asreml(circ ~ 1 + lin(x),
> > random= ~spl(x, 7),
> > splinepoints=list(x=c(118, 484, 664, 1004, 1231, 1372,
> 1582)),
> > data=orange, subset=Tree==1) o.pred <- predict(o1,
> > classify = "x", predictpoints=list(x=200)) o.pred$predictions$x$pvals
> >
> > Notes:
> > - The cells of the hypertable are calculated from all model terms
> constructed
> > solely from factors in the averaging and classify sets.
> >
> > x predicted.value standard.error est.status
> > 1 118.0 30.76677 6.044060 Estimable
> > 2 191.2 37.54277 5.304859 Estimable
> > 3 264.4 44.31327 4.738176 Estimable
> > 4 337.6 51.07276 4.357754 Estimable
> > 5 410.8 57.81574 4.126452 Estimable
> > 6 484.0 64.53670 3.978459 Estimable
> > 7 557.2 71.22233 3.880751 Estimable
> > 8 630.4 77.82801 3.849271 Estimable
> > 9 664.0 80.81897 3.846828 Estimable
> > 10 737.2 87.22137 3.829022 Estimable
> > etc
> >
> > Since I specified a prediction at x=200, I was expecting to see one
> > row for "x=200" in the table.
> >
> > Would appreciate any clarifications of this.
> >
> > Thanks, K Wright
> >
>
> ********************************DISCLAIMER****
> **********
> The information contained in the above e-mail message or messages (which
> includes any attachments) is confidential and may be legally privileged.
> It is intended only for the use of the person or entity to which it is
> addressed. If you are not the addressee any form of disclosure,
> copying, modification, distribution or any action taken or omitted in
> reliance on the information is unauthorised. Opinions contained in the
> message(s) do not necessarily reflect the opinions of the Queensland
> Government and its authorities. If you received this communication in
> error, please notify the sender immediately and delete it from your
> computer system network.
>
>
> ********************************DISCLAIMER****************************
> The information contained in the above e-mail message or messages
> (which includes any attachments) is confidential and may be legally
> privileged. It is intended only for the use of the person or entity
> to which it is addressed. If you are not the addressee any form of
> disclosure, copying, modification, distribution or any action taken
> or omitted in reliance on the information is unauthorised. Opinions
> contained in the message(s) do not necessarily reflect the opinions
> of the Queensland Government and its authorities. If you received
> this communication in error, please notify the sender immediately and
> delete it from your computer system network.
>
Received on Fri Oct 29 2007 - 09:48:46 EST