R/interval_posterior_nominator2.r

In bpp: Computations Around Bayesian Predictive Power

interval_posterior_nominator2 <- function(x, prior = "normal", IntEffBoundary, IntFutBoundary, interimSE, priormean, ...){
  # list input arguments
  inp <- list(...)

  res1 <- rep(NA, length(x))
  for (i in 1:length(x)){
    covmat2 <- rbind(c(interimSE[1], interimSE[2]), 
                     c(interimSE[2], interimSE[2])) ^ 2
    t2 <- as.numeric(pmvnorm(lower = IntEffBoundary, upper = IntFutBoundary, mean = c(x[i], x[i]), sigma = covmat2))
    if (prior == "normal"){res1[i] <- t2 * dnorm(x[i], mean = priormean, sd = inp$priorsigma)}
  }
  
  return(res1)
}