R/MAW1.R

MAW1 <- function (U = U, V = V, k = 1000){
  Empirical <- MW1(U = U, V = V)
  random <- lapply(1:k, function(x) MW1(U = sample(U), V = sample(V)))
  EV <- mean(unlist(random))
  return((Empirical - EV)/(1 - EV))
}

Try the mri package in your browser

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

mri documentation built on May 1, 2019, 10:21 p.m.