library(MSKAT)
Y = matrix(rnorm(2000), 1000,2)
X = matrix(rnorm(2000), 1000,2)
G = matrix(rbinom(10000,2,0.02), 1000,10)
MSKAT(MSKAT.cnull(Y,X), G, W.beta=c(1,25))
## Time consuming
X.list = vector('list',2)
for(i in 1:2) X.list[[i]] = rbind(1, t(X))
MKMR(t(Y),X.list, Gm=G, W.beta=c(1,25))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.