Nothing
library(testthat)
library(iMRMC)
context("simRoeMetz")
# This flag should always be false except when the tests are first created.
flagSave <- FALSE
init.lecuyerRNG()
# Create an MRMC data frame
config.gRoeMetz <- sim.gRoeMetz.config()
# Simulate data
result <- sim.gRoeMetz(config.gRoeMetz)
saveResult <- result
# Split the data into four pieces, not including truth
df.MRMC.Apos <- droplevels(result[grepl("pos", result$caseID) & grepl("A", result$modalityID), ])
df.MRMC.Aneg <- droplevels(result[grepl("neg", result$caseID) & grepl("A", result$modalityID), ])
df.MRMC.Bpos <- droplevels(result[grepl("pos", result$caseID) & grepl("B", result$modalityID), ])
df.MRMC.Bneg <- droplevels(result[grepl("neg", result$caseID) & grepl("B", result$modalityID), ])
print("")
cat("modality A mean shift, neg = ", mean(df.MRMC.Aneg$score), "\n")
cat("modality A mean shift, pos = ", mean(df.MRMC.Apos$score), "\n")
cat("modality B mean shift, neg = ", mean(df.MRMC.Bneg$score), "\n")
cat("modality B mean shift, pos = ", mean(df.MRMC.Bpos$score), "\n")
str(result)
#### TEST ####
# Save the result to a file for future comparisons
fileName <- "test_simRoeMetz.Rdata"
if (flagSave) {
save(saveResult, file = file.path("tests", "testthat", fileName))
}
# Recover the expected results
saveResult <- 0
if (!file.exists(fileName)) {
fileName <- file.path("tests", "testthat", fileName)
}
load(fileName)
test_that(
"sim.gRoeMetz does not change", {
expect_equal(saveResult, result,tolerance=1e-5)
}
)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.