[R] Alternatives to integrate?

[. .]

Hi Michael,

This is the problem:

func <- Vectorize(function(x, a, sad, samp="pois", trunc=0, ...) {
  result <- function(x) {
    f1 <- function(n) {
			f <- function() {
        dcom <- paste("d", sad, sep="")
        dots <- c(as.name("n"), list(...))
        do.call(dcom, dots)
			}
      g <- function() {
        dcom <- paste("d", samp, sep="")
        lambda <- a * n
        dots <- c(as.name("x"), as.name("lambda"))
        do.call(dcom, dots)
      }
      f() * g()
    }
    integrate(f1,0,2000)$value
#     adaptIntegrate(f1,0,2000)$integral

#     n <- 0:2000
#     trapz(n,f1(n))

#     area(f1, 0, 2000, limit=10000, eps=1e-100)
  }
  return(result(x) / (1 - result(trunc)))
}, "x")
func(200, 0.05, "exp", rate=0.001)

If you could propose something I will be gratefull.

Thanks in advance.

On Thu, Sep 1, 2011 at 10:55 AM, R. Michael Weylandt
<michael.weylandt at gmail.com> wrote:
> Mr ". .",
>
> MASS::area comes to mind but it may be more helpful if you could say what
> you are looking for / why integrate is not appropriate it is for whatever
> you are doing.
>
> Strictly speaking, I suppose there are all sorts of "alternatives" to
> integrate() if you are willing to be really creative and build something
> from scratch: diff(), cumsum(), lm(), hist(), t(), c(), ....
>
> Michael Weylandt
>
> On Thu, Sep 1, 2011 at 9:53 AM, B77S <bps0002 at auburn.edu> wrote:
>>
>> package "caTools"
>> see ?trapz
>>
>>
>> . wrote:
>> >
>> > Hi all,
>> >
>> > is there any alternative to the function integrate?
>> >
>> > Any comments are welcome.
>> >
>> > Thanks in advance.
>> >
>> > ______________________________________________
>> > 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.
>> >
>>
>> --
>> View this message in context:
>> http://r.789695.n4.nabble.com/Alternatives-to-integrate-tp3783624p3783645.html
>> Sent from the R help mailing list archive at Nabble.com.
>>
>> ______________________________________________
>> 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.
>
>