R/quantile.2skm.R

In pedrohcgs/kmte: An R Package for Treatment Effects with Censored Outcomes

# Compute Quantile for 2SKM distributions
quantile.2skm <- function(F.ecdf, probs = c(0, 0.25, 0.5, 0.75, 1)) {
  x <- knots(F.ecdf)
  x <- x[order(x)]                              #here
  quantiles <- vapply(probs, FUN = function(p) {x[which(p <= F.ecdf(x))[1]]}, FUN.VALUE = 1.0)
  #correct small amount at very top of distribution (deffective distributions).
  #quantiles[is.na(quantiles)==1] <- x[which(F.ecdf(x) == max(F.ecdf(x)))[1]]

  nams <- paste(format(round(probs * 100,
                if (length(probs) > 1) 2 - log10(diff(range(probs))) else 2)),
                "%", sep = "")
  names(quantiles) <- nams
  return(quantiles)
}