R/kfuncCOP_Pd.R

"kfuncCOP_Pd" <- function(z, d=2) {
  f <- sapply(z, function(f) {
         f + f*sum(sapply(1:(d-1), function(k) ((-log(f))^k) / factorial(k))) })
  f[is.nan(f)] <- 0
  return(f)
}

Try the copBasic package in your browser

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

copBasic documentation built on July 23, 2026, 1:07 a.m.