tests/testthat/test-bayNorm.r

In WT215/bayNorm: Single-cell RNA sequencing data normalization

set.seed(12300)
data("EXAMPLE_DATA_list")

# bay_array_N<-array(rpois(1000*50*2,17),dim=c(1000,50,2))
# bay_array_C<-array(rpois(1000*50*2,58),dim=c(1000,50,2))
# \dontrun{
# noisy_output<-NOISY_FUN(bay_array_N,bay_array_C)

test_that("bayNorm runs properly",
          {
set.seed(1258484)
q1<- bayNorm(Data=EXAMPLE_DATA_list$inputdata[seq(1,10),seq(1,30)],BETA_vec = EXAMPLE_DATA_list$inputbeta[seq(1,30)],parallel = FALSE)
set.seed(1258484)
q2<- bayNorm(Data=EXAMPLE_DATA_list$inputdata[seq(1,10),seq(1,30)],BETA_vec = EXAMPLE_DATA_list$inputbeta[seq(1,30)],parallel = FALSE)
expect_identical(q1, q2)
}
)


test_that("bayNorm_sup runs properly",
          {
            set.seed(1258484)
            q1<- bayNorm(Data=EXAMPLE_DATA_list$inputdata[seq(1,10),seq(1,30)],BETA_vec = EXAMPLE_DATA_list$inputbeta[seq(1,30)],parallel = FALSE)

            set.seed(1258484)
            q2<- bayNorm_sup(Data=EXAMPLE_DATA_list$inputdata[seq(1,10),seq(1,30)],PRIORS=q1$PRIORS,input_params = q1$input_params,parallel = FALSE)
            expect_identical(q1$PRIORS, q2$PRIORS)
            }
)