# test script for quietTCRgenes.R - testcases are NOT comprehensive!
test_that("quietTCRgenes works", {
library(SeuratObject)
data("trex_example")
features <- rownames(trex_example@assays$RNA$counts)
expect_equal(
quietTCRgenes(features),
getdata("quietTCRgenes", "quietTCRgenes_feature.vector")
)
VariableFeatures(trex_example) <- features
trex_example <- quietTCRgenes(trex_example)
expect_equal(
quietTCRgenes(features),
VariableFeatures(trex_example)
)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.