[R] cubature

Jeff Newmiller jdnewmil at dcn.davis.ca.us
Sun Mar 25 09:46:16 CEST 2012

Who knows? You mention different machines, yet fail to hint what they are much less provide sessionInfo() output, and x is not defined in your global environment in your snippet. Please read the posting guide.

I am not familiar with this library, but I would be surprised if it were capable of handling variable limits. I didn't see any hint that such would be supported in the help page. Are you sure the same code is being used on different machines?
Jeff Newmiller                        The     .....       .....  Go Live...
DCN:<jdnewmil at dcn.davis.ca.us>        Basics: ##.#.       ##.#.  Live Go...
                                      Live:   OO#.. Dead: OO#..  Playing
Research Engineer (Solar/Batteries            O.O#.       #.O#.  with
/Software/Embedded Controllers)               .OO#.       .OO#.  rocks...1k
Sent from my phone. Please excuse my brevity.

JMDS <jessica.dunn at rmit.edu.au> wrote:

>I am using adaptIntegrate from Cubature to do numerical integration on
>double integral with a 1 x 2 vector x.
>Say the function is something simple to start like f(x)=x1*x2 and I
>wish to
>integrate x1 over (0,365-x2) and x2 over (0,365)
>f <- function(x) {(x[2])*(x[1])} # "x" is vector
>int1<-adaptIntegrate(f, lowerLimit = c(0, 0), upperLimit = c(365-x[2],
>I recieve the following error:
>Error in adaptIntegrate(f, lowerLimit = c(0, 0), upperLimit = c(365 - 
>  object 'x' not found
>The problem is that this code works fine on other machines and I wonder
>If I
>am missing some referenced package somewhere.
>Any help would be greatly appreciated.
>View this message in context:
>Sent from the R help mailing list archive at Nabble.com.
>R-help at r-project.org mailing list
>PLEASE do read the posting guide
>and provide commented, minimal, self-contained, reproducible code.

More information about the R-help mailing list