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

In datashield/dsBaseClient: DataSHIELD Client Functions

.test.skewness.combined <- function(variable.name,some.values,method)
{
 skewness.local <- e1071::skewness(x = some.values, na.rm = TRUE, type = method) 
 
 skewness.server <- ds.skewness(x = variable.name, method = method, type = "combine")
 skewness.server <- round(as.numeric(skewness.server[[1]][1]), digits = 8)
 
 expect_equal(skewness.server, skewness.local, tolerance = ds.test_env$tolerance)
}

.test.skewness.split <- function(variable.name,some.values.1,some.values.2,some.values.3,method)
{
  skewness.local.1 <- e1071::skewness(x = some.values.1, na.rm = TRUE, type = method)
  skewness.local.2 <- e1071::skewness(x = some.values.2, na.rm = TRUE, type = method)
  skewness.local.3 <- e1071::skewness(x = some.values.3, na.rm = TRUE, type = method)
  
  
  skewness.server <- ds.skewness(x = variable.name, method = method, type = "split")
  
  skewness.server.1 <- round(as.numeric(skewness.server[[1]][1]), digits = 8)
  skewness.server.2 <- round(as.numeric(skewness.server[[1]][2]), digits = 8)
  skewness.server.3 <- round(as.numeric(skewness.server[[1]][3]), digits = 8)
  
  expect_equal(skewness.server.1, skewness.local.1, tolerance = ds.test_env$low_tolerance)
  expect_equal(skewness.server.2, skewness.local.2, tolerance = ds.test_env$low_tolerance)
  expect_equal(skewness.server.3, skewness.local.3, tolerance = ds.test_env$low_tolerance)
}

.test.skewness.mean.median.mode<-function(variable.name,method)
{
  #There are not DataSHIELD functions to calculate the median and mode.
  #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.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(as.numeric(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(as.numeric(skewness.server[[1]][1])>0)
  {
    expect_true(Q3-Q2>Q2-Q1)  
    
  }else{
    expect_true(Q3-Q2<Q2-Q1) 
  }
  
}