## ----include = FALSE---------------------------------------------------------- knitr::opts_chunk$set( collapse = TRUE, comment = "#>" ) ## ----setup-------------------------------------------------------------------- library(waddR) spec.output<-c("pval","d.wass^2","perc.loc","perc.size","perc.shape") ## ----diff_loc----------------------------------------------------------------- set.seed(24) ctrl<-rnorm(300,0,1) dd1<-rnorm(300,1,1) set.seed(32) wasserstein.test(ctrl,dd1,method="SP",permnum=10000)[spec.output] wasserstein.test(ctrl,dd1,method="ASY")[spec.output] ## ----diff_size---------------------------------------------------------------- set.seed(24) ctrl<-rnorm(300,0,1) dd2<-rnorm(300,0,2) set.seed(32) wasserstein.test(ctrl,dd2)[spec.output] wasserstein.test(ctrl,dd2,method="ASY")[spec.output] ## ----diff_shape--------------------------------------------------------------- set.seed(24) ctrl<-rnorm(300,6.5,sqrt(13.25)) sam1<-rnorm(300,3,1) sam2<-rnorm(300,10,1) dd3<-sapply(1:300, function(n) { sample(c(sam1[n],sam2[n]),1,prob=c(0.5, 0.5))}) set.seed(32) wasserstein.test(ctrl,dd3)[spec.output] wasserstein.test(ctrl,dd3,method="ASY")[spec.output] ## ----no_diff------------------------------------------------------------------ set.seed(24) ctrl<-rnorm(300,0,1) nodd<-rnorm(300,0,1) set.seed(32) wasserstein.test(ctrl,nodd)[spec.output] wasserstein.test(ctrl,nodd,method="ASY")[spec.output] ## ----session-info------------------------------------------------------------- sessionInfo()