Integrate over an integral in R

Question

I want to solve the following in R:

∫₀^H [π(t) ∫_t^H A(x) dx ] dt <

Answer 1

In your integrand, lower = t is not vectorised, so the call to integrate is not doing what you expected*. Vectorising over t fixes this issue,

expected_loss <- function(H){
  integrand <- function(t) prior(t) * integrate(A, lower = t, upper = H)$value
  vint <- Vectorize(integrand, "t")
  loss <- integrate(vint, lower = 0, upper = H)$value
  return(loss)
} 

expected_loss(.5)
# [1] 0.7946429
expected_loss(1)
# [1] 0.8571429

*: a closer look at integrate revealed that passing vectors to lower and/or upper was silently allowed, but only the first value was taken into account. When integrating over a wider interval the quadrature scheme picked a first point further from the origin, resulting in the unintuitive decrease that you observed.

After reporting this behaviour to r-devel, this user-error will now be caught by integrate thanks to Martin Maechler (R-devel).