Hello,

 

I am trying to integrate a 2-dimensional function that already calls the
function adapt. More precisely, I am calling 

 

adapt(2,lower=c(-100,-100),upper=c(100,100),functn=function(s){1-exp(-50*Unc
enteredGauss(c(-10,10,-10,10),60,s)})

 

where UncenteredGauss is given by the following code in R:

 

UncenteredGauss=function(limsamp,sigma2,s) 

{

int=adapt(2,lower=c(limsamp[1],limsamp[3]),upper=c(limsamp[2],limsamp[4]),fu
nctn=function(u){h(sigma2,s,u)})

value=int$value

return(value)

}

h=function(sigma2,s,u){1/(2*pi*sigma2)*exp(-((u[1]-s[1])^2+(u[2]-s[2])^2)/(2
*sigma2))}

 

This does not return an error message but a value that is not the good one.

UncenteredGauss gives the good result, but I believe the problem comes from
calling in adapt a function that already calls adapt

Any idea on how to integrate a function that already contains an integral???

 

Thank you very much

 

Helene Morlon

Green Lab

Center for Ecology and Evolutionary Biology

University of Oregon, Eugene

morlon.helene@gmail.com


	[[alternative HTML version deleted]]

