tests/testthat/test-rss.AUC.test.R

In generalRSS: Statistical Tools for Balanced and Unbalanced Ranked Set Sampling

library(testthat)

test_that("rss.AUC.test works as expected",{
  set.seed(123)
  rss.data1<-rss.simulation(H=3,nsamp=c(6,6,6),dist="normal", rho=0.8,delta=0)
  rss.data2<-rss.simulation(H=3,nsamp=c(8,8,8),dist="normal", rho=0.8,delta=0.2)
  result=rss.AUC.test(data1=rss.data1, data2=rss.data2, alpha=0.05, delta0=0.5)
  
  expect_equal(result$RSS_AUC, 0.4930556, tolerance = 1e-3)
  expect_equal(result$p.value, 0.9227096, tolerance = 1e-3)
  expect_equal(result$CI, c(0.3571189, 0.6331168), tolerance = 1e-3)
  expect_equal(result$`-2*log.LR`, 0.009413097, tolerance = 1e-3)
})

test_that("rss.AUC.test handles wrong input with error",{
  rss.data1<-rss.simulation(H=3,nsamp=c(6,6,6),dist="normal", rho=0.8,delta=0)
  rss.data2<-rss.simulation(H=3,nsamp=c(7,8,9),dist="normal", rho=0.8,delta=0.2)
  
  expect_error(rss.AUC.test(data1=rss.data1, data2=rss.data2, alpha=5, delta0=0.5),
               "alpha is out of bound.")
  expect_error(rss.AUC.test(data1=rss.data1[,1], data2=rss.data2, alpha=0.05, delta0=0.5),
               "The input data must contain 'rank' and 'y' variables.")
  expect_error(rss.AUC.test(data1=rss.data1[,2], data2=rss.data2, alpha=0.05, delta0=0.5),
               "The input data must contain 'rank' and 'y' variables.")
  expect_error(rss.AUC.test(data1=rss.data1, data2=rss.data2[,1], alpha=0.05, delta0=0.5),
               "The input data must contain 'rank' and 'y' variables.")
  expect_error(rss.AUC.test(data1=rss.data1, data2=rss.data2[,2], alpha=0.05, delta0=0.5),
               "The input data must contain 'rank' and 'y' variables.")
})