Nothing
nTax <- 6L
dataset <- matrix(seq_len(nTax), nrow = nTax, ncol = 6,
dimnames = list(LETTERS[1:nTax], 1:6))
dataset[4, 3:5] <- "-"
dataset[5:6, 4:6] <- "?"
dataset[6, 3] <- "?"
test_that("Errors are handled", {
expect_error(ArtificialExtinction(dataset, "A", "E", replaceAmbiguous = "INVALID"))
expect_error(ArtificialExtinction(dataset, "A", "E", replaceCoded = "INVALID"))
})
test_that("Replacements ok", {
expectation <- matrix(rep("?", 6), 2, 3, dimnames = list(c("A", "B"), 4:6))
expect_equal(expectation,
ArtificialExtinction(dataset, subject = c("A", "B"),
template = "E",
replaceAmbiguous = "ambig")[1:2, 4:6])
expectation[] <- 1
colnames(expectation) <- NULL
expect_equal(PhyDatToMatrix(ArtificialExtinction(
MatrixToPhyDat(dataset),
c("A", "B"), "F", "unif", sampleFrom = "A"
))[1:2, 4:6],
expectation
)
expect_true(all(!"?" == ArtificialExtinction(
dataset[-6, ], 1:2, 5, "freq"
)[1:2, 4:6]))
expect_true(all(ArtificialExtinction(dataset, subject = "E", template = "F",
replaceAmbiguous = "binary")["E", 3:6] %in% 0:1))
expect_equal(setNames(rep("?", 3), 4:6),
ArtificialExtinction(dataset, subject = "E", template = "F",
replaceAmbiguous = "binary",
replaceAll = FALSE)["E", 4:6])
nChar <- 100
dataset <- rbind(subj = c(rep(2, nChar)),
templ = c(rep(3, nChar / 2), rep("?", nChar / 2)))
suppressWarnings(expect_gt(chisq.test(as.integer(
ArtEx(dataset, 1, 2, "binary", "binary")[1, ]))$p.value, 0.05))
})
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.