R/KS.R

pKS <- function(q) c(.Call("pKS2", q, tol = 10^(-6)))

qKS <- function(p) {
  pKSm <- function(y) return(pKS(y)-p)
  quan <- uniroot(pKSm, lower = 0.2, upper = 3)
  res <- quan$root
  return(res)
}

Try the robts package in your browser

Any scripts or data that you put into this service are public.

robts documentation built on May 2, 2019, 4:55 p.m.