R/vorob_threshold.R

In KrigInv: Kriging-Based Inversion for Deterministic and Noisy Computer Experiments

vorob_threshold <- function(pn){
  if(is.null(pn)) {return(0.5)}
  
  M <- length(pn)
  integ.pn <- mean(pn)
  pn.sort <- sort(pn)
  index <- (1-integ.pn)*M
  index.floor <- floor(index)
  index.cap <- min(M,index.floor+1)
  wbig <- index - index.floor
  wsmall <- 1- wbig
  alpha <- pn.sort[index.floor]*wsmall + pn.sort[index.cap]*wbig
  return(alpha)
}