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@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.
