[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