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