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 this particular case, you don't need to Vectorize since the integral of dbeta is already implemented in R through pbeta. Try this:

prior <- function(t) dbeta(t, 1, 24)
#define the integral of the A function instead
Aint     <- function(x,H) pbeta(H, 1, 4) - pbeta(x,1,4)
expected_loss <- function(H){
  integrand<-function(x) Aint(x,H)*prior(x)
  loss          <- integrate(integrand, lower = 0, upper = H)$value
  return(loss)
}
expected_loss(.5)
#[1] 0.7946429
expected_loss(1)
#[1] 0.8571429