Open main menu

Genome Analysis Wiki β

Changes

RareMETALS

4,051 bytes added, 10:18, 18 May 2017
Where to download
rareMETALS is an R-package for performing single or gene-level tests for detecting rare variant associations. For questions regarding the use of this package, please contact Dajiang Liu (dajiang .liu at umich outlook dot educom) or Gonçalo Abecasis (goncalo at umich dot edu). The same methodology is also implemented in command line tools. Please see [http://genome.sph.umich.edu/wiki/Rare-Metal]
== Change Log ==
* 05/18/2017 Version 6.8 incorporates a number of new features and bug fixes. We included support for multi-allelic variants, the support for a new conditional analysis method, the support for cohort level genomic controls, and the bug fixes for calculating heterogeneity statistics such as Q and I2.
* 04/09/2016 Version 6.3 is released. Minor bug fix: Due to different level of missingness of variants in the gene, the single variant association statistics calculated using the covariance matrices of score statistics can be different than single variant association statistics calculated using vstat. This has lead to confusions. It has been fixed in version 6.3. The primary results from version 6.2 should be correct.
* 09/25/2015 Version 6.2 is released. Minor bug fix: Removed the incorrect warning information in version 6.1 when quantitative traits are meta-analyzed. The software incorrectly consider it as binary trait and suggested the use of rareMETALS2.
* 07/23/2015 Version 6.1 is released. Minor feature changes include output for VT the sites where the statistics are maximized; fixed a bug for determining monomorphic sites. Issue warnings when rareMETALS is used to analyze binary trait for meta-analysis.
* 05/19/2015 Version 6.0 is released. Minor feature addition: rareMETALS can now output of the set of variants that are analyzed in VT (i.e. the set of variants with MAF < the threshold where the VT statistic is maximized).
* 04/01/2015 Version 5.9 is released (not a April's fool joke)! A bug in calculating Cochran-Q statistic is fixed. A bug in conditional.rareMETALS.range.group is also fixed. No other analyses are affected.
* 01/24/2015 Version 5.8 is released, which fixed a serious bug for single variant unconditional association tests with group file. If you happen to run the analyses using rareMETALS.single.group() in version 5.7, the results are likely to be incorrect. Please rerun using version 5.8. Please note only rareMETALS.single.group function is affected. All other functions should not be affected by this error.
== Where to download ==
The R package can be downloaded from [[Media:rareMETALS_5rareMETALS_6.98.tar.gz | rareMETALS_5rareMETALS_6.98.tar.gz]]. It will be eventually released on the Comprehensive R-archive Network. If you want to perform gene-level association test using automatically generated annotations, you will also need [[Media:refFlat_hg19.txt.gz | refFlat_hg19.txt.gz]], which is the gene definition modified from refFlat.
== Documentation ==
library(rareMETALS)
res <- rareMETALS.single(score.stat.file,cov.file=NULL,range="19:11200093-11201275",alternative="two.sided",ix.gold=1,callrate.cutoff=0,hwe.cutoff=0)
###result can be explored as below###
> names(res)
[1] "p.value" "ref" "alt" "integratedData" "raw.data"
[6] "clean.data" "statistic" "direction.by.study" "anno" "maf"
[11] "QC.by.study" "no.sample" "beta1.est" "beta1.sd" "hsq.est"
[16] "nearby" "pos"
> print(res$pos)
[1] "19:11200093" "19:11200213" "19:11200235" "19:11200272" "19:11200282" "19:11200309" "19:11200412" "19:11200419"
[9] "19:11200431" "19:11200442" "19:11200475" "19:11200508" "19:11200514" "19:11200557" "19:11200579" "19:11200728"
[17] "19:11200753" "19:11200754" "19:11200806" "19:11200839" "19:11200840" "19:11200896" "19:11201124" "19:11201259"
[25] "19:11201274" "19:11201275"
> print(res$p.value)
[1] 0.551263675 0.056308558 0.172481571 0.734935815 0.922326732 0.053804524 0.886985353 0.903835162 0.005280228 0.266575301
[11] 0.196122312 0.157114376 0.951477852 0.840523624 0.759482777 0.112743041 0.414147263 0.825877149 0.006090142 0.096474975
[21] 0.096474975 0.956407850 0.038234190 0.253512486 0.550935361 0.482315038
== A Simple Tutorial for Using the rareMETALS.single.group function ==
 Dataset used to get the refaltList [[Media:groupFile.txt.gz]]
res.site<-read.table("groupFile.txt",header=T)
alternative = c("two.sided"), callrate.cutoff = 0,
hwe.cutoff = 0, correctFlip = TRUE, analyzeRefAltListOnly = TRUE)
  ###result can be explored as below### > names(res31) [1] "p.value" "ref" "alt" "integratedData" "raw.data" [6] "clean.data" "statistic" "direction.by.study" "anno" "maf" [11] "maf.byStudy" "maf.maxdiff.vec" "ix.maf.maxdiff.vec" "maf.sd.vec" "no.sample.mat" [16] "no.sample" "beta1.est" "beta1.sd" "QC.by.study" "hsq.est" [21] "nearby" "cochranQ.stat" "cochranQ.df" "cochranQ.pVal" "I2" [26] "log.mat" "pos" > print(res31$respos) [1] "19:11200093" "19:11200213" "19:11200235" "19:11200272" "19:11200282" "19:11200309" "19:11200412" "19:11200419" [9] "19:11200431" "19:11200442" "19:11200475" "19:11200508" "19:11200514" "19:11200557" "19:11200579" "19:11200728" [17] "19:11200753" "19:11200754" "19:11200806" "19:11200839" "19:11200840" "19:11200896" "19:11201124" "19:11201259" [25] "19:11201274" "19:11201275" > print(res31$p.outvalue) [1] NA NA NA NA 0.9223267 NA NA NA NA NA NA NA [13] NA NA NA NA NA NA NA NA NA NA NA NA [25] NA NA
== A Simple Tutorial for Using the rareMETALS.range function ==
131
edits