Test EPACTS for DIAGRAM: Difference between revisions

From Genome Analysis Wiki
Jump to navigationJump to search
Cfuchsb (talk | contribs)
No edit summary
Cfuchsb (talk | contribs)
No edit summary
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/index.html here].
Line 11: Line 11:




# 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 epacts2.noref_binary.2012_07_03.tar.gz


# 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:  
= 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:  
<pre>${EPACTS_DIRECTORY}/example/1000G_exome_chr20_example_softFiltered.calls.vcf.gz
<pre>${EPACTS_DIRECTORY}/example/1000G_exome_chr20_example_softFiltered.calls.vcf.gz



Revision as of 17:16, 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 epacts2.noref_binary.2012_07_03.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.