Nothing
context("Testing downstream analysis functions")
gn1 <- paste0("N_", 1:100)
gn2 <- paste0("N_", seq(2, 200, length=100))
sn1 <- paste0("S_", 1:50)
sn2 <- paste0("S_", 1:75)
coexpSets <- list(
a=matrix(rnorm(100*100), 100, dimnames=list(gn1, gn1)),
b=matrix(rnorm(100*100), 100, dimnames=list(gn2, gn2))
)
adjSets <- coexpSets
exprSets <- list(
a=matrix(rnorm(50*100), 50, dimnames=list(sn1, gn1)),
b=matrix(rnorm(75*100), 75, dimnames=list(sn2, gn2))
)
moduleAssignments <- list(a=sample(1:7, 100, replace=TRUE), b=NULL)
names(moduleAssignments[[1]]) <- gn1
# Only analyse modules with > 2 genes
modules <- moduleAssignments[[1]][intersect(names(moduleAssignments[[1]]),
colnames(adjSets[[2]]))]
modules <- table(modules)
modules <- names(modules[modules > 2])
nModules <- length(modules)
test_that("'networkProperties' function runs without error", {
expect_is(
networkProperties(
adjSets, exprSets, coexpSets, moduleAssignments, modules=modules[1],
verbose=FALSE
), "list"
)
props <- networkProperties(
adjSets[[1]][1:10, 1:10], exprSets[[1]][,1:10], coexpSets[[1]][1:10, 1:10],
verbose=FALSE
)
expect_is(props, "list")
})
test_that("'nodeOrder' function runs without error", {
n <- nodeOrder(
adjSets, exprSets, coexpSets, moduleAssignments, modules=modules[1],
verbose=FALSE
)
expect_is(n, "character")
n <- nodeOrder(
adjSets, NULL, coexpSets, moduleAssignments, modules=modules[1:2],
orderModules=FALSE, verbose=FALSE
)
expect_is(n, "character")
})
test_that("'sampleOrder' function runs without error", {
s <- sampleOrder(
adjSets, exprSets, coexpSets, moduleAssignments, modules=modules[1],
verbose=FALSE
)
expect_is(s, "character")
expect_error(
sampleOrder(
adjSets, NULL, coexpSets, moduleAssignments, modules=modules[1:2],
simplify=FALSE, verbose=FALSE
)
)
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.