Re: Spline predict points]

Dave,

Perhaps you could ask VSN to announce to the list when the new library
has been made available.

scott

-----Original Message-----
From: ASReml users discussion group [mailto:ASREML-L_at_AGRIC.NSW.GOV.AU]
On Behalf Of Butler, David
Sent: Wednesday, 31 October 2007 8:47 AM
To: ASREML-L_at_AGRIC.NSW.GOV.AU
Subject: Re: Spline predict points]

Hi Kevin,

Looks like VSN don't have the latest version on their site - it should
be library 2.00bi. Probably not their fault as I went on leave and asked
someone send it to them on my behalf. I'll check it out.

Dave.

-----Original Message-----
From: ASReml users discussion group [mailto:ASREML-L_at_AGRIC.NSW.GOV.AU]
On Behalf Of Kevin Wright
Sent: Tuesday, 30 October 2007 12:49 AM
To: ASREML-L_at_AGRIC.NSW.GOV.AU
Subject: 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
   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*
ion 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 Sun Oct 31 2007 - 09:35:58 EST