library(testthat)
test_that("meth k means cluster", {
library(m6ALogisticModel)
SE_CQN <- readRDS("SE_CQN_filtered.rds")
set.seed(1)
Y <- plot_row_joint(SE_CQN,RETURN_INDX = T)
expect_that(Y, is_a("integer"))
expect_that(length(Y) == nrow(SE_CQN), is_true())
})
test_that("multinomial model", {
glm_multinomial(Y,mcols(SE_CQN),"GLM_bg")
})
test_that("go multinomial", {
library(TxDb.Hsapiens.UCSC.hg19.knownGene)
library(org.Hs.eg.db)
go_multinomial(Y,
row_ranges = rowRanges(SE_CQN),
txdb = TxDb.Hsapiens.UCSC.hg19.knownGene,
orgDb = org.Hs.eg.db,
HDER = "GO_test",
GO_Slim = T
)
})
test_that("multinomial general", {
glm_regular(Y,
row_ranges = mcols(rowgr),
HDER = "glm",
family = "poisson",
CUT_OFF = 5
)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.