vartools
offers nine varitan association tests:
R CMD INSTALL vartools
library(vartools)
pgdata <- as.matrix(read.table(system.file("extdata","phengen.dat",package="vartools"), as.is=T, skip = 1))
cmc.pvalue <- cmc(table=pgdata)
print(cmc.pvalue)
library(vartools)
alpha <- 0.05
num.permutation <- 3000
quiet <- 1
alternative <- 1
maf.upper.bound <- 0.05
casectrl.dat <- read.table(system.file("extdata","phengen.dat",package="vartools"), skip = 1)
kbac.pvalue <- kbac(table=casectrl.dat, alpha, num.permutation, quiet, maf.upper.bound, alternative)
print(kbac.pvalue)
library(vartools)
casectrl.dat <- read.table(system.file("extdata","phengen.dat",package="vartools"), skip = 1)
calpha.stat <- calpha(casectrl.dat)
calpha.stat
library(vartools)
?vt
casectrl.dat <- read.table(system.file("extdata","phengen2.dat",package="vartools"), skip = 1)
vt.stat <- vt(casectrl.dat)
vt.stat
library(vartools)
?rarecover
casectrl.dat <- read.table(system.file("extdata","phengen.dat",package="vartools"), skip = 1)
rarecover.stat <- rarecover(casectrl.dat)
rarecover.stat
library(vartools)
?rbt
casectrl.dat <- read.table(system.file("extdata","phengen.dat",package="vartools"), skip = 1)
rbt.stat <- rbt(casectrl.dat)
rbt.stat
In additional, vartools
offers SKAT, WSS, aSum tests.Please check the user manual to use them.
Ilya Y. Zhbannikov, i.zhbannikov@mail.ru, ilyaz@uidaho.edu
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.