Tutorial: RAREMETAL

From Genome Analysis Wiki
Jump to navigationJump to search

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