[R] problem with deriv3?

[Prof Brian D Ripley]

On Sun, 10 Mar 2002, B. D. Mccullough wrote:

> Using S+ 6 for Linux and
> R 1.4.1
>
> Taking verbatim some lines from S-Plus that work
> perfectly, and running them on R, I get an
> unexpected error message:

This is uncredited (by you) work of Bill Venables in extending S's deriv
which we contributed to S-PLUS.  For several years this extension was only
in the MASS library, and may well still be for your S-PLUS (you didn't
mention the version).

> R : Copyright 2002, The R Development Core Team
> Version 1.4.1  (2002-01-30)
>
>
> > library(MASS)
> > data <- read.table("/usr/local/data/mcmanus.dat",header=T)
> > mcmanus.frm <-
> deriv3(~-log(pnorm((2*D1-1)*(b0+b1*TT+b2*Y+b3*LF+b4*NW+b5*D2))),c("b0","b1","b
> 2","b3","b4","b5"),function(b0,b1,b2,b3,b4,b5,D1,TT,Y,LF,NW,D2) NULL)
> Error in deriv3.formula(~-log(pnorm((2 * D1 - 1) * (b0 + b1 * TT + b2 *  :
> 	Function   ^A is not in the derivatives table
> >
>
>
> I did a search on the R website and at the e-mail archive,
> and found nothing of help.

>From the NEWS file of R-devel

    o   deriv/deriv3() can now handle expressions involving pnorm and
        dnorm (with a single argument), as in S-PLUS.

> How can I decipher "^A" and make the appropriate corrections
> to the derivative table in R?

Just download and build R-devel.  Your example works there.


> PS: I believe that in R I should delete the "log" and instead use
> the "log.p=T" option; presumably it is numerically more efficient.
> If this is incorrect, please let me know.

Not if you want to use deriv3: see the NEWS item.

-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272860 (secr)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._