tests/testthat/test-PlotEmpiricalOperatingCharacteristics.R

In dpc10ster/rjafroc: Artificial Intelligence Systems and Observer Performance

context("PlotEmpiricalOperatingCharacteristics: ROC & FROC & vectors & lists")
test_that("ROC & FROC & vectors & lists", {

  fn <- paste0(test_path(), "/goodValues361/Plots/ds04-ROC-vectors", ".rds")
  if (!file.exists(fn)) {
    warning(paste0("File not found - generating new ",fn))
    p1 <- PlotEmpiricalOperatingCharacteristics(dataset04, trts = seq(1,5), rdrs = seq(1,4), opChType = "wAFROC")
    saveRDS(p1, file = fn)
  }
  p1 <- readRDS(fn)
  p2 <- PlotEmpiricalOperatingCharacteristics(dataset04, trts = seq(1,5), rdrs = seq(1,4), opChType = "wAFROC")
  expect_is(p2$Plot, "ggplot")
  expect_equal(p2$Points, p1$Points, check.environment = FALSE)
  
  plotT <- list(1, 2, c(1:2), c(1:2))
  plotR <- list(2, c(2:3), c(1:3), 1)
  fn <- paste0(test_path(), "/goodValues361/Plots/ds04-ROC-lists", ".rds")
  if (!file.exists(fn)) {
    warning(paste0("File not found - generating new ",fn))
    p1 <- PlotEmpiricalOperatingCharacteristics(dataset04, trts = plotT, rdrs = plotR, opChType = "ROC")
    saveRDS(p1, file = fn)
  }
  p1 <- readRDS(fn)
  p2 <- PlotEmpiricalOperatingCharacteristics(dataset04, trts = plotT, rdrs = plotR, opChType = "ROC")
  expect_is(p2$Plot, "ggplot")
  expect_equal(p2$Points, p1$Points, check.environment = FALSE)

  fn <- paste0(test_path(), "/goodValues361/Plots/ds04-FROC", ".rds")
  if (!file.exists(fn)) {
    warning(paste0("File not found - generating new ",fn))
    p1 <- PlotEmpiricalOperatingCharacteristics(dataset04, trts = plotT, rdrs = plotR, opChType = "FROC")
    saveRDS(p1, file = fn)
  }
  p1 <- readRDS(fn)
  p2 <- PlotEmpiricalOperatingCharacteristics(dataset04, trts = plotT, rdrs = plotR, opChType = "FROC")
  expect_is(p2$Plot, "ggplot")
  expect_equal(p2$Points, p1$Points, check.environment = FALSE)

  fn <- paste0(test_path(), "/goodValues361/Plots/ds04-AFROC", ".rds")
  if (!file.exists(fn)) {
    warning(paste0("File not found - generating new ",fn))
    p1 <- PlotEmpiricalOperatingCharacteristics(dataset04, trts = plotT, rdrs = plotR, opChType = "AFROC")
    saveRDS(p1, file = fn)
  }
  p1 <- readRDS(fn)
  p2 <- PlotEmpiricalOperatingCharacteristics(dataset04, trts = plotT, rdrs = plotR, opChType = "AFROC")
  expect_is(p2$Plot, "ggplot")
  expect_equal(p2$Points, p1$Points, check.environment = FALSE)

  fn <- paste0(test_path(), "/goodValues361/Plots/ds04-wAFROC", ".rds")
  if (!file.exists(fn)) {
    warning(paste0("File not found - generating new ",fn))
    p1 <- PlotEmpiricalOperatingCharacteristics(dataset04, trts = plotT, rdrs = plotR, opChType = "wAFROC")
    saveRDS(p1, file = fn)
  }
  p1 <- readRDS(fn)
  p2 <- PlotEmpiricalOperatingCharacteristics(dataset04, trts = plotT, rdrs = plotR, opChType = "wAFROC")
  expect_is(p2$Plot, "ggplot")
  expect_equal(p2$Points, p1$Points, check.environment = FALSE)

  fn <- paste0(test_path(), "/goodValues361/Plots/ds04-AFROC1", ".rds")
  if (!file.exists(fn)) {
    warning(paste0("File not found - generating new ",fn))
    p1 <- PlotEmpiricalOperatingCharacteristics(dataset04, trts = plotT, rdrs = plotR, opChType = "AFROC1")
    saveRDS(p1, file = fn)
  }
  p1 <- readRDS(fn)
  p2 <- PlotEmpiricalOperatingCharacteristics(dataset04, trts = plotT, rdrs = plotR, opChType = "AFROC1")
  expect_is(p2$Plot, "ggplot")
  expect_equal(p2$Points, p1$Points, check.environment = FALSE)

  fn <- paste0(test_path(), "/goodValues361/Plots/ds04-wAFROC1", ".rds")
  if (!file.exists(fn)) {
    warning(paste0("File not found - generating new ",fn))
    p1 <- PlotEmpiricalOperatingCharacteristics(dataset04, trts = plotT, rdrs = plotR, opChType = "wAFROC1")
    saveRDS(p1, file = fn)
  }
  p1 <- readRDS(fn)
  p2 <- PlotEmpiricalOperatingCharacteristics(dataset04, trts = plotT, rdrs = plotR, opChType = "wAFROC1")
  expect_is(p2$Plot, "ggplot")
  expect_equal(p2$Points, p1$Points, check.environment = FALSE)

})  




context("PlotEmpiricalOperatingCharacteristics: LROC")
test_that("PlotOperatingCharacteristics-LROC", {
  
  K1 <- 10;K2 <- 10;mu <- 1;lambda <- 1;nu <- 0.8;zeta1 <- -3;perCase <- rep(1, K2)
  frocData <- SimulateFrocDataset(mu, lambda, nu, zeta1, I = 2, J = 5, K1, K2, perCase, seed = 5)
  lrocData <- DfFroc2Lroc(frocData)
  fn <- paste0(test_path(), "/goodValues361/Plots/lrocData-ROC", ".rds")
  if (!file.exists(fn)) {
    warning(paste0("File not found - generating new ",fn))
    p <- PlotEmpiricalOperatingCharacteristics(lrocData, opChType = "ROC" )
    saveRDS(p, file = fn)
  }
  p1 <- readRDS(fn)
  p2 <- PlotEmpiricalOperatingCharacteristics(lrocData, opChType = "ROC" )
  expect_is(p2$Plot, "ggplot")
  expect_equal(p2$Points, p1$Points,check.environment = FALSE)

  fn <- paste0(test_path(), "/goodValues361/Plots/lrocData-LROC", ".rds")
  if (!file.exists(fn)) {
    warning(paste0("File not found - generating new ",fn))
    p <- PlotEmpiricalOperatingCharacteristics(lrocData, opChType = "LROC" )
    saveRDS(p, file = fn)
  }
  p1 <- readRDS(fn)
  p2<-PlotEmpiricalOperatingCharacteristics(lrocData, opChType = "LROC" )
  expect_is(p2$Plot, "ggplot")
  expect_equal(p2$Points, p1$Points, check.environment = FALSE)
  
  plotT <- list(1, 2)
  plotR <- list(seq(1,5), seq(1,5)) # 5 readers
  fn <- paste0(test_path(), "/goodValues361/Plots/lrocData-ROC-lists", ".rds")
  if (!file.exists(fn)) {
    warning(paste0("File not found - generating new ",fn))
    p <- PlotEmpiricalOperatingCharacteristics(lrocData, trts = plotT, rdrs = plotR,  opChType = "ROC")
    saveRDS(p, file = fn)
  }
  p1 <- readRDS(fn)
  p2 <- PlotEmpiricalOperatingCharacteristics(lrocData, trts = plotT, rdrs = plotR,  opChType = "ROC")
  expect_is(p2$Plot, "ggplot")
  expect_equal(p2$Points, p1$Points, check.environment = FALSE)
  
  
  fn <- paste0(test_path(), "/goodValues361/Plots/lrocData-LROC-lists", ".rds")
  if (!file.exists(fn)) {
    warning(paste0("File not found - generating new ",fn))
    p <- PlotEmpiricalOperatingCharacteristics(lrocData, trts = plotT, rdrs = plotR,  opChType = "LROC")
    saveRDS(p, file = fn)
  }
  p1 <- readRDS(fn)
  p2<-PlotEmpiricalOperatingCharacteristics(lrocData, trts = plotT, rdrs = plotR,  opChType = "LROC")
  expect_is(p2$Plot, "ggplot")
  expect_equal(p2$Points, p1$Points, check.environment = FALSE)
  
  fn <- paste0(test_path(), "/goodValues361/Plots/lrocData-LROC-vectors", ".rds")
  if (!file.exists(fn)) {
    warning(paste0("File not found - generating new ",fn))
    p <- PlotEmpiricalOperatingCharacteristics(lrocData, trts = c(1,2), rdrs = seq(1,5),  opChType = "LROC")
    saveRDS(p, file = fn)
  }
  p1 <- readRDS(fn)
  p2<- PlotEmpiricalOperatingCharacteristics(lrocData, trts = c(1,2), rdrs = seq(1,5),  opChType = "LROC")
  expect_is(p2$Plot, "ggplot")
  expect_equal(p2$Points, p1$Points, check.environment = FALSE)
  
})