From Genome Analysis Wiki
Jump to navigationJump to search
23 bytes removed
, 10:52, 29 January 2014
Line 39: |
Line 39: |
| | | |
| raremetalworker --ped example1.ped --dat example1.dat --vcf example1.vcf.gz --traitName QT1 | | raremetalworker --ped example1.ped --dat example1.dat --vcf example1.vcf.gz --traitName QT1 |
− | --inverseNormal --makeResiduals --kinSave --kinGeno --prefix STUDY1 --labelHits | + | --inverseNormal --makeResiduals --kinSave --kinGeno --prefix STUDY1 |
| | | |
| This command transforms phenotypes to normality (--inverseNormal), calculates trait residuals after adjusting for covariates (--makeResiduals), estimates relatedness between individuals and saves this for later use (--kinSave and --kinGeno), and even generates some annotated PDF files summarizing results (--labelHits). | | This command transforms phenotypes to normality (--inverseNormal), calculates trait residuals after adjusting for covariates (--makeResiduals), estimates relatedness between individuals and saves this for later use (--kinSave and --kinGeno), and even generates some annotated PDF files summarizing results (--labelHits). |
Line 54: |
Line 54: |
| | | |
| raremetalworker --ped example2.ped --dat example2.dat --vcf example2.vcf.gz --traitName QT1 | | raremetalworker --ped example2.ped --dat example2.dat --vcf example2.vcf.gz --traitName QT1 |
− | --inverseNormal --makeResiduals --kinSave --kinGeno --prefix STUDY2 --labelHits | + | --inverseNormal --makeResiduals --kinSave --kinGeno --prefix STUDY2 |
| | | |
| ==STEP 3: Run RAREMETAL for Meta-Analysis== | | ==STEP 3: Run RAREMETAL for Meta-Analysis== |