test/test_m6A_cluster_profiler.R

In ZhenWei10/m6ALogisticModel: Linear model analysis on RNA modification data.

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
  )
})