[R] problem with deriv3?
B. D. Mccullough
bdm25 at drexel.edu
Mon Mar 11 03:10:45 CET 2002
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:
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.
How can I decipher "^A" and make the appropriate corrections
to the derivative table in R?
Please reply to me. I will post a summary to the list.
Bruce
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.
B. D. McCullough
bdmccullough at drexel.edu
www.pages.drexel.edu/~bdm25
w: 215-895-2134 f: 215-895-2907
Associate Professor
Department of Decision Sciences
Academic Building room 230
Drexel University
Philadelphia, PA 19104
``Disorganized technical documentation, albeit comprehensive, signals
a failure of the marketing department since it alone (if any) can
hire English majors.'' W. Kahan, ``Marketing versus Mathematics''
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list