[R] Problem with adapt package

Martin Maechler maechler at stat.math.ethz.ch
Fri Jun 29 14:03:50 CEST 2001

>>>>> "MM" == Martin Maechler <maechler at stat.math.ethz.ch> writes:

>>>>> "JonR" == Jonathan Rougier <J.C.Rougier at durham.ac.uk> writes:
    JonR> Hi Anthony,
    JonR> Anthony Tate wrote:
    >>> I am having problems using the adapt add on package to integrate a
    >>> function over 2 dimensions, using the following code.
    >>> library(adapt)
    >>> alb1 <- function(p,X,th,alpha){
    >>> y <- p[1]; z <- p[2]
    >>> c <- sqrt((X + z*tan(th))^2 + z^2)
    >>> r <- sqrt(c^2 + y^2)
    >>> exp(-alpha*r) / r^3
    >>> }
    >>> x <- 0.05
    >>> y <- x
    >>> w <- 0.5
    >>> delta <- 0.1
    >>> czang <- pi/4
    >>> alpha <- 15.0
    >>> y[i] <- adapt(ndim=2, lower=c(0,0), upper=c(w,delta),
    >>> minpts=1000, maxpts=10^6, funct=alb1,  eps=0.005,
    >>> X=x, th=czang, alpha=alpha)$value
    >>> The code gives the following error
    >>> Error in function (p, X, th, alpha)  : Argument "X" is missing, with no default

    JonR> This is a bug that I thought we had squashed.  Crudely there is a
    JonR> "..."  missing in the line "ff<-functn".  I have a feeling that
    JonR> Martin fixed this with some tricky code using ".Alias", but the
    JonR> fix has not appeared on CRAN.  Martin?!

    MM> Yes, this definitely rings a bell!  Indeed I had fixed something
    MM> like that, submitted to Thomas, who had put it into the integrate
    MM> package then.

    MM> "adapt" has been packaged by Brian Ripley and I haven't checked
all wrong!  I apologize.
B.R. and Thomas Lumley and me were just discussing the topic, 
but Thomas, not Brian volunteered to re-package the adapt() function!
    MM>  ...
    MM> Okay, yes, apply (something like) the following patch

    MM>  <...>
 (see last message, the patch is fine).

Martin Maechler <maechler at stat.math.ethz.ch>	http://stat.ethz.ch/~maechler/
Seminar fuer Statistik, ETH-Zentrum  LEO D10	Leonhardstr. 27
ETH (Federal Inst. Technology)	8092 Zurich	SWITZERLAND
phone: x-41-1-632-3408		fax: ...-1228			<><
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