R/comparison.R

Defines functions comparison

comparison <- function(x,y){
  Fx <- ecdf(x)
  Fy <- ecdf(y)
  qx <- quantile(Fx, seq(.1,.9,.1))
  qy <- quantile(Fy, seq(.1,.9,.1))
  dist <- sqrt(sum((qx - qy)^2)) / sqrt(sum(qx^2))
  if(is.nan(dist)){
    dist <- 0
  }
  #print(sum((qx - qy)^2))
  #print(sqrt(sum(qx^2)))
  return(dist)
}

Try the EEEA package in your browser

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

EEEA documentation built on June 10, 2025, 9:13 a.m.