Displaying 2 results from an estimated 2 matches for "boundsy".
Did you mean:
bounds
2007 Aug 23
1
Expedite scalar f(x) evaluation over vectors
...proach is too slow for optim().
I have a function that takes a double integral (hopefully correctly) across
two terms, e.g.,
doubleint <- function(c1,c2) {integrate(function(y) {
sapply(y, function(y) {
integrate(function(x) {a*(c1+x)+b*(c2+x)}, boundsx[1], boundsx[2])$value
})
}, boundsy[1],boundsy[2])
}
I would like to rapidly evaluate this function over vectors c1 and c2 of
equal length where the double integral is calculated for each (matching)
element c1_i & c2_i. At present I get length mismatch errors. Furthermore,
mapply() takes too long. Can I expedite this evalu...
2007 Aug 23
1
Clarification: Expedite scalar f(x) evaluation over vectors
...uble integral (hopefully correctly) across
two terms, e.g.,
<<note replacement of example function by hypothetical f(x,y)>>
doubleint <- function(c1,c2) {integrate(function(y) {
sapply(y, function(y) {
integrate(function(x) {f(x,y)}, boundsx[1], boundsx[2])$value
})
}, boundsy[1],boundsy[2])
}
<<I do not have a closed form for the real function>>
I would like to rapidly evaluate this function over vectors c1 and c2 of
equal length where the double integral is calculated for each (matching)
element c1_i & c2_i. At present I get length mismatch errors...