NEWS0.md

MSKAT

Testing the rare variant set association with multiple continuous phenotypes

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

Testing the rare variant set association across multiple traits.



baolinwu/MSKAT documentation built on May 28, 2019, 6:37 p.m.