R help - Jan 2008 - R object as a function

I want to use a function as an argument to ingtegrate it twice.
See the following (senseless) example of a double integration:

test<-function(sf,lo,up,rest) {
  innerFkn<-function(sf,lo) {
    inte=integrate(f=sf,lower=lo,upper=4)
    return( inte$value )
  }
  integral=integrate(f=innerFkn,lower=1,upper=2,sf=sf,lo=lo,up=up)
  return( integral$vlaue+rest )
}
test(sf=stepfun(c(0,1),c(2,-1,3)),lo=0,up=2,rest=12)

Why isn't it possible to define the "innerFkn" inside
"test"?
"sf" is a stepfun, but it should possibly be any function.
How can I define some R object like a stepfun (depending on variables)
which can be evaluated like a function at some "lo"?
Thanks for help,
Thomas

Okay, let me try to better say what I meant:

myfun1<-function(x=5,pa) {
  return(pa[1]*x^2+pa[2]*x+pa[3])
}
myfun2<-function(x=5,param,myfunk1) {
  return(param[1]*myfunk1(x)+param[2]*myfunk1(x))
}
test<-function(pars1,pars2,lo,up){
  mf=myfun1(x=2,pa=8*pars1)
  integ=integrate(f=myfun2,lower=lo,upper=up,param=pars2,myfunk1=mf)
  return( 2*integ$value )
}
test(pars1=c(1,2,3),pars2=c(-1,1),lo=2,up=7)

So the function myfun2 will be integrated and has a function (myfunk1)
itself as an argument.
Why is this not possible?
Thanks,
Thomas

PS: the "x=5" and "pa=8*" has no influence...

On 22/01/2008 5:30 AM, Thomas Steiner wrote:
> I want to use a function as an argument to ingtegrate it twice.
> See the following (senseless) example of a double integration:
> 
> test<-function(sf,lo,up,rest) {
>   innerFkn<-function(sf,lo) {
>     inte=integrate(f=sf,lower=lo,upper=4)
>     return( inte$value )
>   }
>   integral=integrate(f=innerFkn,lower=1,upper=2,sf=sf,lo=lo,up=up)
>   return( integral$vlaue+rest )
> }
> test(sf=stepfun(c(0,1),c(2,-1,3)),lo=0,up=2,rest=12)
> 
> Why isn't it possible to define the "innerFkn" inside
"test"?
It is possible, but it needs to take arguments in the order that 
integrate expects, i.e. the first argument needs to be the values at 
which it will be evaluated.  You don't specify any "x" value.

The other problem is that integrate is not vectorized, it can only take 
scalar values for lower and upper, so you'll need a loop within innerFkn 
to do the integral over the values being passed in.

Duncan Murdoch
> "sf" is a stepfun, but it should possibly be any function.
> How can I define some R object like a stepfun (depending on variables)
> which can be evaluated like a function at some "lo"?
> Thanks for help,
> Thomas
> 
> ______________________________________________
> 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.