[R] R functions with array arguments

[baptiste auguie]

Try this,

sapply(Mabslim , my_gamma, alpha=-1, xstar = -21, xmax = -27)

or wrap it with ?Vectorize,

vmy_gamma = Vectorize(my_gamma, vectorize.args = "xlim")
vmy_gamma(alpha=-1, xstar = -21, xlim= Mabslim, xmax = -27)


HTH,

baptiste


2009/9/17 Maurizio Paolillo <paolillo at na.infn.it>:
> Dear R users,
> I'm trying to implement a self-defined function with multiple
> arguments, one of which is an array, but I find that the result is a
> single value instead of an array.
> This is the example I'm working on:
>
>
> # define integration limit vector
>  > Mabslim <- c(-17.95, -16.65, -17.27, -17.62, -16.76, -17.07, -17.02)
>
> # Define function
>  > schech<-function(x,alpha,xstar) (10^(0.4*(alpha+1)*(xstar-
> x))*exp(-10^(0.4*(xstar-x))))
>
> # Define new function integrating the previous one
>  > my_gamma<-function(alpha,xstar,xlim,xmax)
> integrate(schech,xmax,xlim,alpha,xstar)$value
>
>  > my_gamma(-1,-21,Mabslim,-27)
> [1] 2.487746
>
> Note that 'Mabslim' is used as an upper integration limit within the
> function. I tried to use "sapply(....)" but it looks to me as if this
> can be used only if the array is the first argument of the function.
> I'm a beginner to R, apologies if I'm just using the wrong approach.
> Any suggestion on how to solve this?
> Thanks,
> Maurizio
>
>
>        [[alternative HTML version deleted]]
>
> ______________________________________________
> 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.
>