R help - Oct 2011 - cuhre usage ?? multidimensional integration

my=function(x){
len=1
for(i in 1:len){
y[i]=x[i]
}
g=1
w=NULL
t=NULL
for(i in 1:len)w[i]=x[i+len]
for(i in 1:len)t[i]=x[i+2*len]
for(i in 1:len)g=g*dnorm(y[i])*dnorm(w[i])*dnorm(z[i])
return(g)
}
cuhre(6,1,my,rep(-100,6),rep(100,6))

Error in crff(match.call(), integrand, "cuhre", libargs, ...) : 
  Additional argument  not expected in the integrand function

function change to my=function(x,g,i,j)
result is not right. it should be 1, but it turns out to be 0.039...

How can I make this work?

Thank you!

--
View this message in context:
http://r.789695.n4.nabble.com/cuhre-usage-multidimensional-integration-tp3873478p3873478.html
Sent from the R help mailing list archive at Nabble.com.

Perhaps you should start by writing vectorized code; as it stands,
your code suggests you don't understand what simple operations like y
<- x actually do.

More to your question: what are cuhre  & crff ? They are not in base R
nor in any packages I have current loaded.

Michael


On Wed, Oct 5, 2011 at 1:26 AM, sevenfrost <linshuang11 at gmail.com>
wrote:
> my=function(x){
> len=1
> for(i in 1:len){
> y[i]=x[i]
> }
> g=1
> w=NULL
> t=NULL
> for(i in 1:len)w[i]=x[i+len]
> for(i in 1:len)t[i]=x[i+2*len]
> for(i in 1:len)g=g*dnorm(y[i])*dnorm(w[i])*dnorm(z[i])
> return(g)
> }
> cuhre(6,1,my,rep(-100,6),rep(100,6))
>
> Error in crff(match.call(), integrand, "cuhre", libargs, ...) :
> ?Additional argument ?not expected in the integrand function
>
> function change to my=function(x,g,i,j)
> result is not right. it should be 1, but it turns out to be 0.039...
>
> How can I make this work?
>
> Thank you!
>
> --
> View this message in context:
http://r.789695.n4.nabble.com/cuhre-usage-multidimensional-integration-tp3873478p3873478.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.
>

cuhre is a function in package R2Cuba for multidimensional integration.

There are three kind of variables to integrate. I just use y,w,t to
distinguish them from each other. Is there any problem ?

--
View this message in context:
http://r.789695.n4.nabble.com/cuhre-usage-multidimensional-integration-tp3873478p3876182.html
Sent from the R help mailing list archive at Nabble.com.