Defines functions congruence

Documented in congruence

congruence <- function(Z, R, from=NA, to=NA, steps=10000) {
  # Congruence
  # are arguments provided, or do I have to infer them from the input strings?
  if (is.na(from))  {from <- 0}
  if (is.na(to))    {to   <- length(Z)}
  Zf <- ecdf(Z) # empirical CDF
  Rf <- ecdf(R) # empirical CDF
  r <- to - from
  z <- seq(from, to, by = r/steps)        # values at which to evaluate ecdf
  a <- sum(abs(Zf(z) - Rf(z)))/length(z)  # area between the two curves (approximated)

Try the polrep package in your browser

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

polrep documentation built on Feb. 24, 2018, 3:01 a.m.