tests/testthat/test_plot.R

In catavallejos/BASiCS: Bayesian Analysis of Single-Cell Sequencing data

test_that("plot of chain works", {
  data(ChainSCReg)
  pdf(NULL)
  expect_error({
      p <- plot(ChainSCReg, Gene = 1)  
      p <- plot(ChainSCReg, Param = "nu", Cell = 1)
      p <- plot(ChainSCReg, Param = "beta", RegressionTerm = 1)
      p <- plot(ChainSCReg, Param = "sigma2")
  }, NA)
  dev.off()
})

test_that("plot of BASiCS_Summary works without spikes", {
  data(ChainSC)
  S <- Summary(ChainSC)
  pdf(NULL)
  expect_error(plot(S), NA)
  expect_error(plot(S, Param = "nu"), NA)
  dev.off()
})

test_that("plot works for BASiCS_Chain (non-regression, spikes)", {
  data(ChainSC)
  pdf(NULL)
  expect_error({
    p <- plot(ChainSC, Param = "mu", Gene = 1)
    p <- plot(ChainSC, Param = "delta", Gene = 1)
    p <- plot(ChainSC, Param = "phi", Cell = 1)
    p <- plot(ChainSC, Param = "s", Cell = 1)
    p <- plot(ChainSC, Param = "nu", Cell = 1)
    p <- plot(ChainSC, Param = "theta")
  }, NA)
  dev.off()
})

test_that("plot works for BASiCS_Chain (regression, no spikes)", {
  data(ChainSCReg)
  pdf(NULL)
  expect_error({
    plot(ChainSCReg, Param = "mu", Gene = 1)
    plot(ChainSCReg, Param = "delta", Gene = 1)
    plot(ChainSCReg, Param = "epsilon", Gene = 1)
    plot(ChainSCReg, Param = "s", Cell = 1)
    plot(ChainSCReg, Param = "nu", Cell = 1)
    plot(ChainSCReg, Param = "theta")
    plot(ChainSCReg, Param = "beta", RegressionTerm = 1)
  }, NA)
  dev.off()
})

test_that("plot works for BASiCS_Summary with all valid combinations", {
  data(ChainSCReg)
  SChain <- Summary(ChainSCReg)
  pdf(NULL)
  expect_error({
    plot(SChain, Param = "mu", Param2 = "delta")
    plot(SChain, Param = "mu", Param2 = "epsilon")
    plot(SChain, Param = "delta", Param2 = "epsilon")
    plot(SChain, Param = "s", Param2 = "phi")
    plot(SChain, Param = "s", Param2 = "nu")
    plot(SChain, Param = "phi", Param2 = "nu")
    plot(SChain, Param = "sigma2")
    plot(SChain, Param = "beta")
  }, NA)
  dev.off()
})

test_that("BASiCS_ShowFit works", {
  data(ChainSCReg)
  pdf(NULL)
  expect_warning(BASiCS_showFit(ChainSCReg, smooth = FALSE), "deprecated")
  expect_error({
    BASiCS_ShowFit(ChainSCReg, smooth = TRUE)
    BASiCS_ShowFit(ChainSCReg, smooth = FALSE)
  }, NA)
  dev.off()
})

test_that("Diagnostic plot works", {
  data(ChainSCReg)
  pdf(NULL)
  expect_warning(BASiCS_diagPlot(ChainSCReg), "deprecated")
  dev.off()
  g <- BASiCS_DiagPlot(ChainSCReg)
  expect_s3_class(g, "ggplot")
  g <- BASiCS_DiagPlot(ChainSCReg, Param = "delta")
  expect_s3_class(g, "ggplot")
  g <- BASiCS_DiagPlot(ChainSCReg, Param = "epsilon")
  expect_s3_class(g, "ggplot")
  g <- BASiCS_DiagPlot(ChainSCReg, Param = "delta", x = "mu", y = "epsilon")
  expect_s3_class(g, "ggplot")
  expect_error(
    BASiCS_DiagPlot(ChainSCReg, Param = "nu", x = "mu", y = "epsilon"),
    "Invalid combination of parameters"
  )
  expect_error(
    BASiCS_DiagPlot(ChainSCReg, Param = "delta", x = "delta", y = "nu"),
    "Invalid combination of parameters:"
  )
  g <- BASiCS_DiagPlot(ChainSCReg, Param = "delta", x = "mu", y = "epsilon")
  expect_s3_class(g, "ggplot")

  ChainSCReg@parameters[["epsilon"]][, 1] <- NA
  g <- BASiCS_DiagPlot(ChainSCReg, Param = "epsilon", x = "mu", y = "epsilon")
  expect_s3_class(g, "ggplot")
  g <- BASiCS_DiagPlot(ChainSCReg, Param = "epsilon")
  expect_s3_class(g, "ggplot")
})

test_that("Diagnostic hist works", {
  data(ChainSCReg)
  data(ChainSC)
  expect_warning(g <- BASiCS_diagHist(ChainSCReg, Param = "mu"), "deprecated")
  expect_s3_class(g, "ggplot")
  g <- BASiCS_DiagHist(ChainSC, Param = "delta")
  expect_s3_class(g, "ggplot")
  g <- BASiCS_DiagHist(ChainSC, Param = "nu")
  expect_s3_class(g, "ggplot")
  g <- BASiCS_DiagHist(ChainSC, Param = "nu")
  expect_s3_class(g, "ggplot")
})

test_that("DE plots work (non-regression)", {
  data(ChainSC)
  data(ChainRNA)
  Test <- BASiCS_TestDE(ChainSC, ChainRNA, Plot = FALSE, PlotOffset = FALSE)
  expect_s3_class(
    g <- BASiCS_PlotDE(Test),
    "gg"
  )
  expect_s3_class(
    g <- BASiCS_PlotDE(Test@Results[[1]]),
    "gg"
  )
  expect_s3_class(
    g <- BASiCS_PlotDE(Test@Results[[2]]),
    "gg"
  )

  for (type in c("MA", "Grid", "Volcano")) {
    expect_s3_class(
      g <- BASiCS_PlotDE(Test@Results[[1]], Plots = type),
      "gg"
    )
    expect_s3_class(
      g <- BASiCS_PlotDE(Test@Results[[2]], Plots = type),
      "gg"
    )
  }
})

test_that("DE plots work (regression)", {
  data(ChainSCReg)
  data(ChainRNAReg)
  Test <- BASiCS_TestDE(
    ChainSCReg,
    ChainRNAReg,
    Plot = FALSE,
    PlotOffset = FALSE
  )
  expect_s3_class(g <- BASiCS_PlotDE(Test), "gg")
  expect_s3_class(g <- BASiCS_PlotDE(Test@Results[[1]]), "gg")
  expect_s3_class(g <- BASiCS_PlotDE(Test@Results[[2]]), "gg")
  expect_s3_class(g <- BASiCS_PlotDE(Test@Results[[3]]), "gg")
  for (type in c("MA", "Grid", "Volcano")) {
    expect_s3_class(
      g <- BASiCS_PlotDE(Test@Results[[1]], Plots = type),
      "gg"
    )
    expect_s3_class(
      g <- BASiCS_PlotDE(Test@Results[[2]], Plots = type),
      "gg"
    )
    expect_s3_class(
      g <- BASiCS_PlotDE(Test@Results[[3]], Plots = type),
      "gg"
    )
  }
})

test_that("BASiCS_PlotOffset", {
  data(ChainSC)
  data(ChainRNA)
  expect_s3_class(g <- BASiCS_PlotOffset(ChainSC, ChainRNA), "gg")
})

test_that("DiagPlot works", {
  data(ChainSC)
  expect_s3_class(g <- BASiCS_DiagPlot(ChainSC), "gg")
  expect_s3_class(g <- BASiCS_DiagPlot(ChainSC, HLine = TRUE), "gg")
  expect_s3_class(g <- BASiCS_DiagPlot(ChainSC, HLine = 150), "gg")
  expect_s3_class(g <- BASiCS_DiagPlot(ChainSC, HLine = FALSE), "gg")
})

test_that("DiagHist works", {
  data(ChainSC)
  expect_s3_class(g <- BASiCS_DiagHist(ChainSC), "gg")
  expect_s3_class(g <- BASiCS_DiagHist(ChainSC, VLine = TRUE), "gg")
  expect_s3_class(g <- BASiCS_DiagHist(ChainSC, VLine = 150), "gg")
  expect_s3_class(g <- BASiCS_DiagHist(ChainSC, VLine = FALSE), "gg")
})