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]]


