Displaying 1 result from an estimated 1 matches for "my_gamma".
2009 Sep 17
2
R functions with array arguments
...rking 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...