math test ds.skewness

Author: Leire Abarrategui Martinez

tests/testthat/definition_tests/def-ds.skewness.R

@@ -28,7 +28,7 @@
   expect_equal(skewness.server.3, skewness.local.3, tolerance = ds.test_env$tolerance)
 }
 
-.test.skewness.mean.median.mode.combined<-function(variable.name,method)
+.test.skewness.mean.median.mode<-function(variable.name,method)
 {
   #There are not DataSHIELD functions to calculate the median and mode.
   #We need danger function
@@ -40,7 +40,52 @@
   server.variable<-c(server.variable$study1,server.variable$study2,server.variable$study3)
 
   server.mean<-ds.mean(x = variable.name, type = "combine")
+  server.mean<-server.mean[[1]][1]
   server.median<-median(server.variable)
   server.mode<-DescTools::Mode(server.variable)
 
+  if(skewness.server[[1]][1]>0)
+  {
+    if(!is.na(server.mode[1]))
+      {
+        expect_true(server.mode[1]<server.median)  
+        
+      }
+    expect_true(server.median<server.mean)
+    
+  }else{
+    if(!is.na(server.mode[1]))
+    {
+      expect_true(server.mode[1]>server.median) 
+    }
+    expect_true(server.median>server.mean)
+  }
+  
+}
+
+.test.skewness.Q1.Q2.Q3<-function(variable.name,method)
+{
+  #There are not DataSHIELD functions to calculate Q1,Q2 and Q3.
+  #We need danger function
+  library(dsDangerClient)
+  
+  skewness.server <- ds.skewness(x = variable.name, method = method, type = "combine")
+  
+  server.variable<-ds.DANGERprint(variable.name)
+  server.variable<-c(server.variable$study1,server.variable$study2,server.variable$study3)
+  
+  server.quantile<-quantile(server.variable)
+  Q1<-server.quantile[2]
+  Q2<-server.quantile[3]
+  Q3<-server.quantile[4]
+  
+  if(skewness.server[[1]][1]>0)
+  {
+    expect_true(Q3-Q2>Q2-Q1)  
+    
+  }else{
+    expect_true(Q3-Q2<Q2-Q1) 
+  }
+  
 }
+    

tests/testthat/setup.R

@@ -15,6 +15,7 @@ context("setup")
 
 library(lme4)
 library(e1071)
+library(DescTools)
 
 library(DSOpal)
 library(DSLite)

tests/testthat/test-math-ds.skewness.R

@@ -1,2 +1,34 @@
 source("connection_to_datasets/init_testing_datasets.R")
-source("definition_tests/def-ds.skewness.R")
+source("definition_tests/def-ds.skewness.R")
+
+context("ds.skewness()::math::mean.median.mode::multiple")
+test_that("skewness >0 & mode<median<mean or skewness<0 & mean<median<mode",
+{
+ connect.all.datasets()
+            
+  .test.skewness.mean.median.mode('D$INTEGER', method = 1)
+  .test.skewness.mean.median.mode('D$NON_NEGATIVE_INTEGER', method = 2)
+  .test.skewness.mean.median.mode('D$POSITIVE_INTEGER', method = 3)
+  #.test.skewness.mean.median.mode('D$NEGATIVE_INTEGER', method = 1) 
+  #.test.skewness.mean.median.mode('D$NUMERIC', method = 2)
+  .test.skewness.mean.median.mode('D$NON_NEGATIVE_NUMERIC', method = 3)
+  .test.skewness.mean.median.mode('D$POSITIVE_NUMERIC', method = 1)
+  .test.skewness.mean.median.mode('D$NEGATIVE_NUMERIC', method = 2) 
+
+})
+
+context("ds.skewness()::math::Q1.Q2.Q3::multiple")
+test_that("skewness >0 & Q3-Q2>Q2-Q1 or skewness<0 & Q3-Q2<Q2-Q1",
+{
+  connect.all.datasets()
+            
+  .test.skewness.Q1.Q2.Q3('D$INTEGER', method = 1)
+  .test.skewness.Q1.Q2.Q3('D$NON_NEGATIVE_INTEGER', method = 2)
+  .test.skewness.Q1.Q2.Q3('D$POSITIVE_INTEGER', method = 3)
+  .test.skewness.Q1.Q2.Q3('D$NEGATIVE_INTEGER', method = 1) 
+  #.test.skewness.Q1.Q2.Q3('D$NUMERIC', method = 2)
+  .test.skewness.Q1.Q2.Q3('D$NON_NEGATIVE_NUMERIC', method = 3)
+  .test.skewness.Q1.Q2.Q3('D$POSITIVE_NUMERIC', method = 1)
+  .test.skewness.Q1.Q2.Q3('D$NEGATIVE_NUMERIC', method = 2) 
+       
+})