Description Usage Arguments Value Author(s) Examples
Calculate HSIC values
1 |
A |
is one matrix. |
B |
is one matrix. |
P |
is "phenoype", a vector. |
... |
not used. |
a matrix
Beibei Jiang beibei_jiang@psych.mpg.de
1 2 3 4 5 6 7 8 9 10 11 |
# simulate some data
set.seed(123)
geno1 <- matrix(sample(0:2, size = 1000, replace = TRUE, prob = c(0.5, 0.3, 0.2)), ncol = 10)
geno2 <- matrix(sample(0:2, size = 2000, replace = TRUE, prob = c(0.4, 0.3, 0.3)), ncol = 20)
dimnames(geno1) <- list(row = paste0("IND", 1:nrow(geno1)), col = paste0("rs", 1:ncol(geno1)))
dimnames(geno2) <- list(row = paste0("IND", 1:nrow(geno2)), col = paste0("exm", 1:ncol(geno2)))
epiHSIC(A = scale(geno1),
B = scale(geno2),
P = rnorm(100))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.