epiHSIC {episcan}

Calculate HSIC values

Description

Calculate HSIC values

Usage

epiHSIC(A = NULL, B = NULL, P = NULL, ...)

Arguments

A

is one matrix.

B

is one matrix.

P

is "phenoype", a vector.

...

not used.

Value

a matrix

Author(s)

Beibei Jiang [e-mail]

Examples

 
# 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))