stopifnot(require("testthat"), require("lcc"))
context("running sampled CCC")
test_that("runing sampled CCC",{
require(MASS)
S<-matrix(c(2.0, 0.3,
0.3, 1.8), byrow = TRUE, ncol=2)
data<-mvrnorm(n=100, mu=c(10,9.5), Sigma=S)
ccc_lin<-2*cov(data[,1],data[,2])/(var(data[,1])+var(data[,2])+(mean(data[,1])-mean(data[,2]))^2)
ccc_lcc <- lcc:::CCC(data[,1],data[,2])
expect_equal(ccc_lin, ccc_lcc)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.