[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