SKAT test

vtools show test SKAT

Name:          SKAT
Description:   SKAT (Wu et al 2011) wrapper of its original R implementation
usage: vtools associate --method SKAT [-h] [--name NAME]
                                      [-k {linear,linear.weighted,IBS,IBS.weighted,quadratic,2wayIX}]
                                      [--beta_param BETA_PARAM BETA_PARAM]
                                      [-m {davies,liu,liu.mod,optimal}]
                                      [-i {fixed,random}]
                                      [--logistic_weights PARAM PARAM]
                                      [-r [CORR [CORR ...]]]
                                      [--missing_cutoff MISSING_CUTOFF]
                                      [--resampling N] [--small_sample]
                                      [--resampling_kurtosis N]
                                      {quantitative,disease}

SNP-set (Sequence) Kernel Association Test (Wu et al 2011). This is a wrapper
for the R package "SKAT" implemented & maintained by Dr. Seunggeun Lee, with a
similar interface and minimal descriptions based on the SKAT package
documentation (May 11, 2012). Please refer to
http://http://cran.r-project.org/web/packages/SKAT/ for details of usage. To
use this test you should have R installed with SKAT v0.75 or higher. The SKAT
commands applied to the data will be recorded and saved in the project log
file.

positional arguments:
  {quantitative,disease}
                        Phenotype is quantitative trait or disease trait (0/1
                        coding). Default set to "quantitative"

optional arguments:
  -h, --help            show this help message and exit
  --name NAME           Name of the test that will be appended to names of
                        output fields, usually used to differentiate output of
                        different tests, or the same test with different
                        parameters.
  -k {linear,linear.weighted,IBS,IBS.weighted,quadratic,2wayIX}, --kernel {linear,linear.weighted,IBS,IBS.weighted,quadratic,2wayIX}
                        A type of kernel. Default set to "linear.weighted".
                        Please refer to SKAT documentation for details.
  --beta_param BETA_PARAM BETA_PARAM
                        Parameters for beta weights. It is only used with
                        weighted kernels. Default set to (1,25). Please refer
                        to SKAT documentation for details.
  -m {davies,liu,liu.mod,optimal}, --method {davies,liu,liu.mod,optimal}
                        A method to compute the p-value. Default set to
                        "davies". Please refer to SKAT documentation for
                        details.
  -i {fixed,random}, --impute {fixed,random}
                        A method to impute missing genotypes. Default set to
                        "fixed". Please refer to SKAT documentation for
                        details.
  --logistic_weights PARAM PARAM
                        This option, if specified, will get the logistic
                        weights from genotype matrix Z and apply this weight
                        to SKAT. It requires two input parameters par1 and
                        par2. To use the SKAT default setting, type
                        `--logistic_weights 0.07 150'. Please refer to SKAT
                        documentation for details.
  -r [CORR [CORR ...]], --corr [CORR [CORR ...]]
                        The pho parameter of SKAT test. Default is 0. Please
                        refer to SKAT documentation for details.
  --missing_cutoff MISSING_CUTOFF
                        a cutoff of the missing rates of SNPs. Any SNPs with
                        missing rates higher than cutoff will be excluded from
                        the analysis. Default set to 0.15
  --resampling N        Number of resampling using bootstrap method. Set it to
                        '0' if you do not want to apply resampling.
  --small_sample        This option, if evoked, will apply small sample
                        adjustment "SKAT_Null_Model_MomentAdjust" for small
                        sample size and binary trait. Please refer to SKAT
                        documentation for details.
  --resampling_kurtosis N
                        Number of resampling to estimate kurtosis, for small
                        sample size adjustment. Set it to '0' if you do not
                        wnat to apply the adjustment. The SKAT default setting
                        is 10000. Please refer to SKAT documentation for
                        details.

vtools associate rare status -m "SKAT --name skat quantitative" --group_by refGene.name2 --\
to_db skat -j8 > skat.txt

INFO: 3180 samples are found
INFO: 2632 groups are found
INFO: Starting 8 processes to load genotypes
Loading genotypes: 100% [=========================================================================================================================================] 3,180 32.8/s in 00:01:36
Testing for association: 100% [================================================================================================================================] 2,632/147 8.9/s in 00:04:56
INFO: Association tests on 2632 groups have completed. 147 failed.
INFO: Using annotation DB skat in project test.
INFO: Annotation database used to record results of association tests. Created on Wed, 30 Jan 2013 21:34:23


vtools show fields | grep skat

skat.refGene_name2           refGene_name2
skat.sample_size_skat        Sample size
skat.Q_stats_skat            Test statistic for SKAT, "Q"
skat.pvalue_skat             p-value


head skat.txt

refGene_name2	sample_size_skat	Q_stats_skat	pvalue_skat
AADACL4	3180	33707.7	0.379148
ABCD3	3180	1178.25	0.961708
AAMP	3180	5905.71	0.612598
ABCB10	3180	55121.9	0.109206
ABCB6	3180	16500.2	0.812062
ABCG5	3180	9829.17	0.76832
ABI2	3180	42491.9	0.0100467
ABHD1	3180	1315.49	0.880286
ABL2	3180	794.385	0.963097