deriv3 example on Venables/Ripley page 263 - R-help

Michael Wang · 2001-05-06T13:17:46Z

What does it require to run the example on page 263 of Venables/Ripley book using R 1.3.0? I got the following error and I have no clue how to correct it. Thank you very much. (If you reply to the list, please send a copy to me.) > lmix2 <- deriv3( + ~ -log(p*dnorm((x-u1)/s1)/s1 + (1-p)*dnorm((x-u2)/s2)/s2), + c("p", "u1", "s1", "u2", "s2"), + function(x, p, u1, s1, u2, s2) NULL) Error in deriv3.formula(~-log(p * dnorm((x - u1)/s1)/s1 + (1 - p) * dnorm((x - : Function H is not in

Brian Ripley

Sun, May 6, 2001 7:03 AM #

On Sun, 6 May 2001, Michael Wang wrote:

Time travel and a large dose of optimism!

1) R 1.3.0 is not released, and it is certainly not finished yet.

2) Its deriv3 function is not compatible with that in the S-PLUS MASS
library, in particular it is not currently user-extensible as the V&R
version is.

3) The S-PLUS library MASS contains extensions to the engine behind
S's deriv and MASS's deriv3 to cover this example.  If by release time
R's equivalents are user-extensible it would need the extensions to be
written.

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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._