Difference between revisions of "GREGOR"
Line 8: | Line 8: | ||
=== Download from webpage === | === Download from webpage === | ||
− | Through this link [http://www.sph.umich.edu/csg/jich/FallInBed/ FallInBed Download], you can download a copy of | + | Through this link [http://www.sph.umich.edu/csg/jich/FallInBed/ FallInBed Download], you can download a copy of GREGOR. |
=== Download from GitHub with Git === | === Download from GitHub with Git === | ||
You can create your own git clone(copy) using: | You can create your own git clone(copy) using: | ||
− | git clone https://github.com/jinchen-umich/ | + | git clone https://github.com/jinchen-umich/GREGOR.git |
or | or | ||
− | git clone git://github.com/jinchen-umich/ | + | git clone git://github.com/jinchen-umich/GREGOR.git |
− | Either of these two commands creates a directory called | + | Either of these two commands creates a directory called GREGOR in the current directory. |
=== Update your copy === | === Update your copy === | ||
If you have already downloaded your copy, use the following commands to update: | If you have already downloaded your copy, use the following commands to update: | ||
− | 1. cd pathToYourCopy/ | + | 1. cd pathToYourCopy/GREGOR |
2. git pull | 2. git pull | ||
Line 28: | Line 28: | ||
# Latest Code (master branch) | # Latest Code (master branch) | ||
#: via Website | #: via Website | ||
− | #:# Go to : https://github.com/jinchen-umich/ | + | #:# Go to : https://github.com/jinchen-umich/GREGOR |
#:# Click on the <code>Download ZIP</code> button on the right side panel. | #:# Click on the <code>Download ZIP</code> button on the right side panel. | ||
#: via Command Line | #: via Command Line | ||
− | #:: <code>wget https://github.com/jinchen-umich/ | + | #:: <code>wget https://github.com/jinchen-umich/GREGOR/archive/master.zip</code> |
After downloading the file, uncompress (unzip/untar) it. The directory created will be named <code>FallInBed</code>. | After downloading the file, uncompress (unzip/untar) it. The directory created will be named <code>FallInBed</code>. | ||
Line 39: | Line 39: | ||
To build GREGOR, copy the GREGOR package to the directory you want, and then run the following command: | To build GREGOR, copy the GREGOR package to the directory you want, and then run the following command: | ||
− | tar xzvf | + | tar xzvf GREGOR.tar.gz |
− | After you unzip, you can find 3 directories in " | + | After you unzip, you can find 3 directories in "GREGOR" (./example ./lib ./script). |
== Download reference files == | == Download reference files == | ||
Line 47: | Line 47: | ||
Download the reference files from this link [http://www.sph.umich.edu/csg/jich/FallInBed/ FallInBed Download], then un-package the file | Download the reference files from this link [http://www.sph.umich.edu/csg/jich/FallInBed/ FallInBed Download], then un-package the file | ||
− | tar xzvf | + | tar xzvf GREGOR.ref.tar.gz |
After unzip, you will get 47 reference files in the directory "~/ref" | After unzip, you will get 47 reference files in the directory "~/ref" | ||
Line 55: | Line 55: | ||
When you run | When you run | ||
− | perl | + | perl GREGOR.pl |
− | you will get some information about | + | you will get some information about GREGOR |
---------------------------------------------------------------------------------- | ---------------------------------------------------------------------------------- | ||
− | + | GREGOR.pl : Functional annotation of trait-associated variants | |
---------------------------------------------------------------------------------- | ---------------------------------------------------------------------------------- | ||
This program tests for enrichment of an input list of trait-associated index | This program tests for enrichment of an input list of trait-associated index | ||
Line 73: | Line 73: | ||
Report Bug(s) : jich[at]umich[dot]edu | Report Bug(s) : jich[at]umich[dot]edu | ||
---------------------------------------------------------------------------------- | ---------------------------------------------------------------------------------- | ||
− | Usage : perl | + | Usage : perl GREGOR.pl --conf [conf.file] |
---------------------------------------------------------------------------------- | ---------------------------------------------------------------------------------- | ||
Line 79: | Line 79: | ||
The following command is a typical command line: | The following command is a typical command line: | ||
− | perl | + | perl GREGOR.pl --conf [conf.file] |
Example configuration file can be found in example directory. Users have to modify the configurations before running. | Example configuration file can be found in example directory. Users have to modify the configurations before running. | ||
== Configuration File == | == Configuration File == | ||
− | The example configuration file below illustrate how to configure the | + | The example configuration file below illustrate how to configure the GREGOR configuration file. |
############################################################################### | ############################################################################### | ||
Line 131: | Line 131: | ||
== Testing GREGOR == | == Testing GREGOR == | ||
− | There is an example directory in ~/ | + | There is an example directory in ~/GREGOR. You can find index SNP file, 3 bed files, bed file index and example config file. |
After change your config file, you can run a test. | After change your config file, you can run a test. | ||
− | perl ~/ | + | perl ~/GREGOR/script/GREGOR.pl --conf ~/GREGOR/example/example.conf |
After running 2 minutes +/- 1 minutes. You will get result file "StatisticSummaryFile.txt" in your defined output directory. | After running 2 minutes +/- 1 minutes. You will get result file "StatisticSummaryFile.txt" in your defined output directory. |
Revision as of 06:59, 24 September 2014
GREGOR
GREGOR (Genomic Regulatory Elements and Gwas Overlap algoRithm) is a tool built to evaluate global enrichment of trait-associated variants in experimentally annotated epigenomic regulatory features.
Because all reference data are version hg19, please make sure that your index SNP list and BED files are also version hg19.
Get GREGOR Source Codes
Download from webpage
Through this link FallInBed Download, you can download a copy of GREGOR.
Download from GitHub with Git
You can create your own git clone(copy) using:
git clone https://github.com/jinchen-umich/GREGOR.git
or
git clone git://github.com/jinchen-umich/GREGOR.git
Either of these two commands creates a directory called GREGOR in the current directory.
Update your copy
If you have already downloaded your copy, use the following commands to update:
1. cd pathToYourCopy/GREGOR 2. git pull
Download From GitHub without Git
If there is no git in your system, you can still download from GitHub:
- Latest Code (master branch)
- via Website
- Go to : https://github.com/jinchen-umich/GREGOR
- Click on the
Download ZIP
button on the right side panel.
- via Command Line
- via Website
After downloading the file, uncompress (unzip/untar) it. The directory created will be named FallInBed
.
Build GREGOR
To build GREGOR, copy the GREGOR package to the directory you want, and then run the following command:
tar xzvf GREGOR.tar.gz
After you unzip, you can find 3 directories in "GREGOR" (./example ./lib ./script).
Download reference files
Download the reference files from this link FallInBed Download, then un-package the file
tar xzvf GREGOR.ref.tar.gz
After unzip, you will get 47 reference files in the directory "~/ref"
Basic Usage Example
When you run
perl GREGOR.pl
you will get some information about GREGOR
GREGOR.pl : Functional annotation of trait-associated variants
This program tests for enrichment of an input list of trait-associated index SNPs ([chr:pos] format or rsID, hg19) in experimentally annotated regulatory domains (BED files).
Note: the index SNPs should be hg19 version. All maf and LD data are from 1000G EUR samples! (Release date : May 21, 2011)
Version : 1.1.0
Report Bug(s) : jich[at]umich[dot]edu
Usage : perl GREGOR.pl --conf [conf.file]
The following command is a typical command line:
perl GREGOR.pl --conf [conf.file]
Example configuration file can be found in example directory. Users have to modify the configurations before running.
Configuration File
The example configuration file below illustrate how to configure the GREGOR configuration file.
############################################################################### # CHIPSEQ ENRICHMENT CONFIGURATION FILE # This configuration file contains run-time configuration of # CHIP_SEQ ENRICHMENT ############################################################################### ## KEY ELEMENTS TO CONFIGURE : NEED TO MODIFY ############################################################################### INDEX_SNP_FILE = /net/dumbo/home/jchen/prj/chiseq/codes/FallInBed_Binomial/example/example.index.snps.rsid.list.txt ## e.g. /home/myid/data/FallInBed/example.snp.txt BED_FILE_INDEX = /net/dumbo/home/jchen/prj/chiseq/codes/FallInBed_Binomial/example/example.bed.file.index ## e.g. /home/myid/data/FallInBed/bedfiles.index REF_DIR = /net/dumbo/home/jchen/prj/chiseq/codes/FallInBed_Binomial/ref/ ## reference directory R2THRESHOLD = 0.7 LDWINDOWSIZE = 1000000 OUT_DIR = /net/dumbo/home/jchen/prj/chiseq/codes/FallInBed_Binomial/example/example.rsid.20130808/ ## e.g. /home/myid/data/FallInBed/result/ MIN_NEIGHBOR_NUM = 500 BEDFILE_IS_SORTED = True MOSRUN = mosbatch -E/tmp -i -m2000 -j20,43,122,135,137,138,149,151,153,154,155,156,162,163 sh -c
In the config file, there are several parameters to adjust:
INDEX_SNP_FILE: This file contains a single column of trait-associated input SNPs, without a header. Variants can be listed in rsid or hg19 chr:pos format.
BED_FILE_INDEX: This file lists the datasets (e.g. BED files) to be used for enrichment analysis. Use complete paths to file locations and make sure positions are in hg19 format.
REF_DIR: Define reference file directory which you download at here.
R2THRESHOLD and LDWINDOWSIZE: These two parameters define the index SNP LD proxies by R2 threshold and LD window size.
OUT_DIR: All result files are saved to this folder, where the script will create multiple sub-directories. Index SNPs are in the folder "index_SNP"; Random SNPs are in the folder "random_SNP".
MIN_NEIGHBOR_NUM: Define neighbor number around index SNP. Script will find no less than this number around every index SNP.
BEDFILE_IS_SORTED: True or false, depending on whether the BED files listed in the index file are sorted.
Results Output
The file StatisticSummaryFile.txt in the output directory contains enrichment results with the following information:
Bed_File: The individual datasets used in the enrichment analysis
InBed_Index_SNP: Number of index SNPs or their LD proxies that overlaps regulatory regions in each dataset
Pvalue: P-value calculated assuming a sum of binomial distributions to represent the number of index SNPs (or LD proxies) that overlap a dataset compared to the expectation observed in the matched control sets
- Note: SNPs that cannot be converted from rsID to chr:pos format are listed in the output file rsid.index.snp.txt. SNPs for which there are no LD proxies or no MAF data available are listed in the output file nonannoted.index.snp.txt.
Testing GREGOR
There is an example directory in ~/GREGOR. You can find index SNP file, 3 bed files, bed file index and example config file. After change your config file, you can run a test.
perl ~/GREGOR/script/GREGOR.pl --conf ~/GREGOR/example/example.conf
After running 2 minutes +/- 1 minutes. You will get result file "StatisticSummaryFile.txt" in your defined output directory.
Acknowledgements
GREGOR is the result of collaborative efforts by Cristen Willer, Jin Chen, Wei Zhou, Ellen Schmidt, He Zhang, and Goncalo Abecasis. Please email Cristen Willer [cristen@umich.edu] with any questions.