tests/testthat/test-renderer3-roc.R

acontext("TestROC")

data(TestROC, package="animint2")

##RColorBrewer::display.brewer.all()
## x <- RColorBrewer::brewer.pal(12, "Paired")
## cat(deparse(x))

algo.colors <-
  c(hmcan="#A6CEE3",
    hmcan.broad="#1F78B4",
    PeakSegJoint="grey40",
    macs="#FDBF6F",
    macs.broad="#FF7F00")

path.before.params <- list(
  title="Test ROC curves, predicted peaks and errors",
  roc=ggplot()+
  geom_path(aes(FPR, TPR, group=Model, key=Model, color=Model),
            data=TestROC$roc)+
  scale_color_manual(values=algo.colors, breaks=names(algo.colors))+
  geom_point(aes(FPR, TPR, color=Model, key=paste(model, parameter),
                 size=parameter,
                 fill=parameter),
             shape=21,
             data=subset(TestROC$parameters, parameter=="learned"))+
  geom_point(aes(FPR, TPR, color=Model, key=paste(model, parameter),
                 size=parameter,
                 fill=parameter),
             shape=21,
             data=subset(TestROC$parameters, parameter=="default"))+
  geom_point(aes(FPR, TPR, color=Model,
                 key=ModelParam),
             clickSelects="ModelParam",
             size=4,
             alpha=0.9,
             data=subset(TestROC$roc, same.as.prev==FALSE))+
  scale_fill_manual(values=c(default="black", learned="white"))+
  scale_size_manual(values=c(default=7, learned=10))+
  scale_shape_manual(values=c(default=20, learned=1))+
  coord_equal()+
  theme_grey()+
  scale_x_continuous(paste(
    "False positive rate in test labels",
    "= Probability(peak | no peak)"),
                     breaks=seq(0, 1, by=0.2))+
  scale_y_continuous(paste(
    "True positive rate in test labels",
    "= Probability(peak | peak)"),
                     breaks=seq(0, 1, by=0.2))+
  theme(panel.margin=grid::unit(0, "lines")),
  first=list(ModelParam="PeakSegJoint -0.8"))

test_that("path before params, 5 paths rendered", {
  info <- animint2HTML(path.before.params)
  entry.list <- info$plots$roc$legend$Model$entries
  linetype.computed <- sapply(entry.list, "[[", "pathlinetype")
  linetype.expected <- rep(1, 5)
  expect_identical(linetype.computed, linetype.expected)
  path.list <- getNodeSet(info$html, '//g[@class="geom1_path_roc"]//path')
  expect_equal(length(path.list), 5)
  stroke.vec <- getStyleValue(
    info$html, '//td[@class="roc_legend"]//line', "stroke")
  expect_color(stroke.vec, algo.colors)
})

path.after.params <- list(
  title="Test ROC curves, predicted peaks and errors",
  roc=ggplot()+
  geom_point(aes(FPR, TPR, color=Model, key=paste(model, parameter),
                 size=parameter,
                 fill=parameter),
             shape=21,
             data=subset(TestROC$parameters, parameter=="learned"))+
  geom_point(aes(FPR, TPR, color=Model, key=paste(model, parameter),
                 size=parameter,
                 fill=parameter),
             shape=21,
             data=subset(TestROC$parameters, parameter=="default"))+
  scale_color_manual(values=algo.colors, breaks=names(algo.colors))+
  geom_path(aes(FPR, TPR, group=Model, key=Model, color=Model),
            data=TestROC$roc)+
  geom_point(aes(FPR, TPR, color=Model,
                 key=ModelParam),
             clickSelects="ModelParam",
             size=4,
             alpha=0.9,
             data=subset(TestROC$roc, same.as.prev==FALSE))+
  scale_fill_manual(values=c(default="black", learned="white"))+
  scale_size_manual(values=c(default=7, learned=10))+
  scale_shape_manual(values=c(default=20, learned=1))+
  coord_equal()+
  theme_grey()+
  scale_x_continuous(paste(
    "False positive rate in test labels",
    "= Probability(peak | no peak)"),
                     breaks=seq(0, 1, by=0.2))+
  scale_y_continuous(paste(
    "True positive rate in test labels",
    "= Probability(peak | peak)"),
                     breaks=seq(0, 1, by=0.2))+
  theme(panel.margin=grid::unit(0, "lines")),
  first=list(ModelParam="PeakSegJoint -0.8"))

test_that("path after params, 5 paths rendered", {
  info <- animint2HTML(path.after.params)
  entry.list <- info$plots$roc$legend$Model$entries
  linetype.computed <- sapply(entry.list, "[[", "pathlinetype")
  linetype.expected <- rep(1, 5)
  expect_identical(linetype.computed, linetype.expected)
  path.list <- getNodeSet(info$html, '//g[@class="geom3_path_roc"]//path')
  expect_equal(length(path.list), 5)
  stroke.vec <- getStyleValue(
    info$html, '//td[@class="roc_legend"]//line', "stroke")
  expect_color(stroke.vec, algo.colors)
})

Try the animint2 package in your browser

Any scripts or data that you put into this service are public.

animint2 documentation built on Nov. 22, 2023, 1:07 a.m.