[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:
>Hi,
>
>I am using adaptIntegrate from Cubature to do numerical integration on
>a
>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],
>365))
>
>
>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.
>
>Thanks!
>
>
>
>
>--
>View this message in context:
>http://r.789695.n4.nabble.com/cubature-tp4502516p4502516.html
>Sent from the R help mailing list archive at Nabble.com.
>
>______________________________________________
>R-help at r-project.org mailing list
>https://stat.ethz.ch/mailman/listinfo/r-help
>PLEASE do read the posting guide
>http://www.R-project.org/posting-guide.html
>and provide commented, minimal, self-contained, reproducible code.
More information about the R-help
mailing list