[R] Two functions as parametrs of a function.

Duncan Murdoch murdoch.duncan at gmail.com
Thu Apr 7 14:04:02 CEST 2011


On 07/04/2011 7:45 AM, Kenneth Roy Cabrera Torres wrote:
> Hi R users:
>
> I'm trying to make a function where two of the parameters are
> functions, but I don't know how to put each set of parameters for
> each function.
>
> What am I missing?
>
> I try this code:
>
> f2<-function(n=2,nsim=100,fun1=rnorm,par1=list(),fun2=rnorm,par2=list()){
>      force(fun1)
>      force(fun2)
>      force(n)
>      p1<-unlist(par1)
>      p2<-unlist(par2)

You don't want to unlist those.  Leave them as lists, and use the 
do.call() function for the calls.  For example,

localfun1 <- function() do.call(fun1, c(list(n), par1))


>      force(p1)
>      force(p2)

Calling force() is only useful for arguments to your function:  it kicks 
the evaluator, so it's not so lazy.   Since p1 and p2 were calculated 
earlier, there's no need to force them.  You might need to force n, 
since it is only used in the localfun1 and localfun2 calls, and if the 
argument to f2 changes before you make those calls, the wrong value will 
be used.  That's pretty unlikely with your code, but since fun1 and fun2 
can do anything, it's safest to fix the value of n before you call 
them.  The same argument applies to fun1 and fun2 themselves.

Duncan Murdoch

>      localfun1<- function() fun1(n, p1)
>      localfun2<- function() fun2(n, p2)
>      vp<-replicate(nsim,t.test(localfun1(), localfun2())$p.value)
>      return(vp)
> }
>
> f2(fun1=rbeta,par1=list(shape1=2,shape2=2),fun2=rbeta,par2=list(shape1=1,shape2=2))
>
> Thank you for your help.
>
> Kenneth
>
> ______________________________________________
> 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