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