Dear R community,
I am trying to code a fairly complex equation for optim(). My current
approach 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 evaluation or do I need to
reformulate my approach?
Thanks in advance for your help,
Scott C. Stark
University of Arizona
[[alternative HTML version deleted]]
______________________________________________
R-help@stat.math.ethz.ch 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.
