tests/testthat/test_diff_RNA_ucsc.R

In huerqiang/GeoTcgaData: Processing Various Types of Data on GEO and TCGA

test_that("can parse example diff_RNA_ucsc", {
    skip_on_cran()
    df <- matrix(rnbinom(400, mu = 4, size = 10), 25, 16)
    df <- as.data.frame(df)
    rownames(df) <- paste0("gene", 1:25)
    colnames(df) <- paste0("sample", 1:16)
    df <- log2(df + 1)
    group <- sample(c("group1", "group2"), 16, replace = TRUE)
    df <- cbind(rownames(df), df)
    result <- differential_RNA(counts = df, group = group, 
        filte = FALSE, method = "limma", ucscData = TRUE)
    expect_true( "P.Value"    %in% colnames(result))
})