Difference between revisions of "Test EPACTS for DIAGRAM"

From Genome Analysis Wiki
Jump to navigationJump to search
Line 1: Line 1:
 
= Download  EPACTS =
 
= Download  EPACTS =
  
EPACTS is available for download [http://www.sph.umich.edu/csg/kang/epacts/download/index.html here].
+
EPACTS is available for download [http://www.sph.umich.edu/csg/kang/epacts/download/epacts_v2_01.noref_binary.2012_07_06.tar.gz here].
 
 
  
 
Requirements
 
Requirements
Line 9: Line 8:
 
*Perl vX
 
*Perl vX
 
*gcc vX
 
*gcc vX
 
  
 
= Install EPACTS =
 
= Install EPACTS =
  
 
Uncompress EPACTS package to the directory you would like to install   
 
Uncompress EPACTS package to the directory you would like to install   
   tar xzvf epacts2.noref_binary.2012_07_03.tar.gz
+
   tar xzvf epacts_v2_01.noref_binary.2012_07_06.tar.gz
  
 
= Example =
 
= Example =

Revision as of 17:17, 24 September 2012

Download EPACTS

EPACTS is available for download here.

Requirements

  • Linux (e.g. Ubuntu)
  • Perl vX
  • gcc vX

Install EPACTS

Uncompress EPACTS package to the directory you would like to install

 tar xzvf epacts_v2_01.noref_binary.2012_07_06.tar.gz

Example

Once installed, test out the software by running a quick example using the test data provided in the "example" directory. The example VCF and PED files are:

${EPACTS_DIRECTORY}/example/1000G_exome_chr20_example_softFiltered.calls.vcf.gz

${EPACTS_DIRECTORY}/example/1000G_dummy_pheno.ped


Run the single variant score test on the example data using this command:

${EPACTS_DIR}/epacts single 
--vcf {EPACTS_DIR}/example/1000G_exome_chr20_example_softFiltered.calls.vcf.gz 
--ped {EPACTS_DIR}/example/1000G_dummy_pheno.ped 
--min-maf 0.001 --chr 20 --pheno DISEASE --cov AGE --cov SEX --test b.score --anno 
--out {OUTPUT_DIR}/test --run 2 &

This command will run the single variant test on the input VCF and PED files, with a minimum MAF threshold of 0.001.  The phenotype is "DISEASE" and we are adjusting the analysis with covariates AGE and SEX.  The output file directory prefix is {OUTPUT_DIR}/test.  Finally, EPACTS will run the analysis in parallel on 2 CPUs.