tests/testthat/test-createDistMat.R

## Tests for createDistMat
## library(Rtpca); library(testthat)
context("createDistMat")

library(dplyr)
library(Biobase)

m1 <- matrix(1:12, ncol = 4)
m2 <- matrix(2:13, ncol = 4)
m3 <- matrix(c(2:10, 1:7), ncol = 4)

rownames(m1) <- 1:3
rownames(m2) <- 2:4
rownames(m3) <- 2:5

colnames(m1) <- paste0("X", 1:4)
colnames(m2) <- paste0("X", 1:4)
colnames(m3) <- paste0("X", 1:4)

mat_list <- list(
    m1, m2, m3
)

ctrl_dist_mat <- matrix(c(0, 2 ,2, 0),
                        ncol = 2)
rownames(ctrl_dist_mat) <- 
    colnames(ctrl_dist_mat) <- c("2", "3")

expect_identical(
    createDistMat(mat_list),
    ctrl_dist_mat
)

expr1 <- ExpressionSet(m1)
expr2 <- ExpressionSet(m2)
expr3 <- ExpressionSet(m3)

exprSet_list <- list(
    expr1, expr2, expr3
)

expect_identical(
    createDistMat(exprSet_list),
    ctrl_dist_mat
)

df_list <- list(
    data.frame(m1), 
    data.frame(m2), 
    data.frame(m3)
)

expect_identical(
    createDistMat(exprSet_list), 
    ctrl_dist_mat
)

tbl_list <- list(
    data.frame(m1) %>% 
        as_tibble() %>% 
        mutate(gene_name = 1:3), 
    data.frame(m2) %>% 
        as_tibble() %>% 
        mutate(gene_name = 2:4), 
    data.frame(m3) %>% 
        as_tibble() %>% 
        mutate(gene_name = 2:5)
)

expect_identical(
    createDistMat(tbl_list, 
                      rownameCol = "gene_name"), 
    ctrl_dist_mat
)

Try the Rtpca package in your browser

Any scripts or data that you put into this service are public.

Rtpca documentation built on Nov. 8, 2020, 7:44 p.m.