replace Mean with Difference

Author: qddyy

NAMESPACE

@@ -6,12 +6,12 @@ export(CDF)
 export(ChiSquare)
 export(ContingencyTableTest)
 export(Correlation)
+export(Difference)
 export(Friedman)
 export(JonckheereTerpstra)
 export(KSampleTest)
 export(KolmogorovSmirnov)
 export(KruskalWallis)
-export(Mean)
 export(MultiCompT)
 export(MultipleComparison)
 export(OneSampleTest)

R/Difference.R

@@ -0,0 +1,41 @@
+#' @title `r Difference$private_fields$.name`
+#' 
+#' @description Performs mean based two sample permutation test on data vectors. 
+#' 
+#' @aliases twosample.difference
+#' 
+#' @export
+#' 
+#' @importFrom R6 R6Class
+
+
+Difference <- R6Class(
+    classname = "Difference",
+    inherit = TwoSampleTest,
+    cloneable = FALSE,
+    public = list(
+        #' @description Create a new `Difference` object. 
+        #' 
+        #' @template init_params
+        #' 
+        #' @return A `Difference` object. 
+        initialize = function(
+            method = c("mean", "median"),
+            alternative = c("two_sided", "less", "greater"), n_permu = NULL
+        ) {
+            private$.method <- match.arg(method)
+
+            super$initialize(alternative = match.arg(alternative), n_permu = n_permu)
+        }
+    ),
+    private = list(
+        .name = "Two Sample Test Based on Mean or Median",
+
+        .define_statistic = function() {
+            private$.statistic_func <- switch(private$.method,
+                mean = function(x, y) mean(x) - mean(y),
+                median = function(x, y) median(x) - median(y)
+            )
+        }
+    )
+)

R/Mean.R

@@ -9,7 +9,7 @@ tests <- list(
     onesample.quantile = Quantile,
     onesample.cdf = CDF,
 
-    twosample.mean = Mean,
+    twosample.difference = Difference,
     twosample.wilcoxon = Wilcoxon,
     twosample.scoresum = ScoreSum,
     twosample.ansari = AnsariBradley,

R/pmt.R

@@ -77,9 +77,9 @@ t$ci
 ## Table 2.1.1
 
 ```{r}
-# See ?Mean or ?twosample.mean
+# See ?Difference or ?twosample.difference
 t <- pmt(
-    "twosample.mean", alternative = "greater"
+    "twosample.difference", alternative = "greater", method = "mean"
 )
 
 # See ?Table2.1.1
@@ -107,9 +107,9 @@ t$p_value
 ## Example 2.3.1
 
 ```{r}
-# See ?Mean or ?twosample.mean
+# See ?Difference or ?twosample.difference
 t <- pmt(
-    "twosample.mean", alternative = "greater"
+    "twosample.difference", alternative = "greater", method = "mean"
 )
 
 # See ?Table2.3.1
@@ -194,7 +194,7 @@ t$p_value
 ## Example 2.8.2
 
 ```{r}
-# See ?RatioMeanDeviance or ?twosample.rmd
+# See ?RatioDifferenceDeviance or ?twosample.rmd
 t <- pmt(
     "twosample.rmd", alternative = "greater"
 )