SKATplus | R Documentation |
SKATlus
provides enhanced power over SKAT by properly estimating the null distribution of SKAT.This version uses only subjects with lower phenotypic values for estimating the null distribution. That is, the "controls" are those of lower phenotypic values. When "controls" are of higher phenotypic values, change the sign of the phenotypic values in order to use this function.
SKATplus(
y,
G,
X = NULL,
out_type = "D",
tau = NULL,
permutation = FALSE,
B = 1000
)
y |
a vector of phenotype on |
G |
an |
X |
a matrix of covariates. It has |
out_type |
an indicator of the outcome type. |
tau |
proportion of selected subjects used for SKATplus. |
permutation |
an indicator. Use permutation for p-value or not. |
B |
the number of permutations if |
A list with class "htest
" containing the following components:
* statistic the value of the test statistic, which is the same as SKAT statistic.
* parameter sample size and the number of SNPs
* p.value the p-value for SKATplus computed using Davies' method.
* method a character string indicting the test performed.
* data.name a character string giving the name of the data.
Kai Wang <kai-wang@uiowa.edu>
Wang, K. (2016) Boosting the power of the sequence kernel association test (SKAT) almost surely by properly estimating its null distribution. A J Hum Genet. 99 (1), 104-114.
n=1000
y = c(rep(1, n/2), rep(0, n/2))
maf = seq(0.05, 0.5, 0.05)
g = NULL
for (j in 1:10){
geno.freq = c(maf[j]^2, 2*maf[j]*(1-maf[j]), (1-maf[j])^2)
g = cbind(g, sample(c(0,1,2), n, replace=TRUE, prob=geno.freq))
}
SKATplus(y, g, X=NULL, out_type="D", tau=NULL, permutation=FALSE, B=1000)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.