inst/unitTests/test_agg.R

In hansenlab/mpra: Analyze massively parallel reporter assays

test_agg <- function() {
    data(mpraSetExample)
    data(mpraSetAggExample)
    agg_manual <- getRNA(mpraSetExample, aggregate = TRUE)
    agg_pre <- getRNA(mpraSetAggExample)
    stopifnot(identical(dim(agg_manual), dim(agg_pre)))
    agg_pre <- agg_pre[rownames(agg_manual),]
    colnames(agg_pre) <- gsub("cond_sample_", "", colnames(agg_pre))
    checkEquals(agg_manual, agg_pre)
}

test_data <- function() {
    data(mpraSetExample)
    checkTrue(validObject(mpraSetExample))
}