Nothing
test_that("SRD returns correct results",
{
cov.matrix.1 <- cov(matrix(rnorm(30*10), 30, 10))
cov.matrix.2 <- cov(matrix(rnorm(30*10), 30, 10))
srd.output <- SRD(cov.matrix.1, cov.matrix.2)
expect_that(srd.output, is_a("SRD"))
expect_that(names(srd.output), equals(c("output", "pc1", "model")))
expect_that(names(srd.output[[3]]), equals(c("quantiles", "interval", "code")))
expect_that(dim(srd.output[[1]]), equals(c(10, 6)))
expect_that(length(srd.output[[2]]), equals(10))
expect_that(length(srd.output[[3]]), equals(3))
expect_that(srd.output[[1]], is_a("matrix"))
expect_that(srd.output[[2]], is_a("numeric"))
expect_that(srd.output[[3]], is_a("list"))
srd.scores.bool <- sapply(srd.output[[1]][,1], function(x) isTRUE(x > -1 & x < 1))
srd.lower.bool <- sapply(srd.output[[1]][,2], function(x) isTRUE(x > -1 & x < 1))
srd.upper.bool <- sapply(srd.output[[1]][,3], function(x) isTRUE(x > -1 & x < 1))
srd.sd.bool <- sapply(srd.output[[1]][,4], function(x) isTRUE(x > 0 & x < 1))
expect_that(sum(srd.scores.bool), equals(length(srd.scores.bool)))
expect_that(sum(srd.lower.bool), equals(length(srd.scores.bool)))
expect_that(sum(srd.upper.bool), equals(length(srd.scores.bool)))
expect_that(mean(srd.output[[1]][,5]), equals(0))
expect_that(mean(srd.output[[1]][,6]), equals(0))
expect_that(mean(srd.output[[2]]), equals(0))
expect_that(as.numeric(srd.output[[2]] < srd.output[[3]][[2]]), equals(srd.output[[3]][[3]]))
}
)
test_that("SRD returns correct results on lists",
{
m.list <- RandomMatrix(10, 4)
suppressWarnings(RNGversion("3.5.0"))
set.seed(42)
srd.all <- SRD(m.list)
suppressWarnings(RNGversion("3.5.0"))
set.seed(42)
one.and.two <- srd.all[1, 2][[1]]
one.and.two.default <- SRD(m.list[[1]], m.list[[2]])
expect_that(one.and.two, is_equivalent_to(one.and.two.default))
expect_that(srd.all, is_a('matrix'))
expect_true(isSymmetric(srd.all))
}
)
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.