Difference between revisions of "Tutorial: RAREMETAL"

From Genome Analysis Wiki
Jump to navigationJump to search
Line 27: Line 27:
  
 
   $yourPath/bin/raremetalworker --ped $yourLocalPath/rmw_tutorial/inputfiles/example1.ped --dat $yourLocalPath/rmw_tutorial/inputfiles/example1.dat --vcf   
 
   $yourPath/bin/raremetalworker --ped $yourLocalPath/rmw_tutorial/inputfiles/example1.ped --dat $yourLocalPath/rmw_tutorial/inputfiles/example1.dat --vcf   
  $yourLocalPath/rmw_tutorial/inputfiles/example1.vcf.gz --kinGeno --kinSave --traitName LDL --inverseNormal --makeResiduals --useCovariates --prefix example1
+
        $yourLocalPath/rmw_tutorial/inputfiles/example1.vcf.gz --kinGeno --kinSave --traitName LDL --inverseNormal --makeResiduals --useCovariates --prefix example1
  
 
==STEP 3: Run RareMETAL to do Meta Analysis==
 
==STEP 3: Run RareMETAL to do Meta Analysis==

Revision as of 19:23, 7 March 2013

RareMETAL Tutorial

Introduction

In this tutorial, we will show how to run RareMetalWorker and RareMETAL to perform meta analysis using example data sets.

RareMetalWorkr is the pre-processing tool to analyze data from individual studies and generate summary statistics for further meta analysis using RareMETAL.

RareMETAL is the tool to do gene-level meta analysis.

STEP 1: Install Software and Download Example Data Sets

 tar xvzf raremetalworker.tutorial.tgz #extract
 cd rmw_tutorial 
 tar xvzf raremetal.tutorial.tgz #extract
 cd raremetal_tutorial

STEP 2: Run RareMetalWorker on Individual Studies

  • The first example has 743 individuals coded as unrelated according to PED file (each person belongs to an individual family).
  • there are ~1000 markers included in the VCF file.
  • To analyze this sample accounting for hidden relatedness, an empirical kinship should be calculated.
  • By using the following command, covariates are also adjusted and residuals inverse normalized on the fly.
 $yourPath/bin/raremetalworker --ped $yourLocalPath/rmw_tutorial/inputfiles/example1.ped --dat $yourLocalPath/rmw_tutorial/inputfiles/example1.dat --vcf  
       $yourLocalPath/rmw_tutorial/inputfiles/example1.vcf.gz --kinGeno --kinSave --traitName LDL --inverseNormal --makeResiduals --useCovariates --prefix example1

STEP 3: Run RareMETAL to do Meta Analysis