Changes

From Genome Analysis Wiki
Jump to navigationJump to search
no edit summary
Line 1: Line 1:  +
[[Category:RAREMETAL]]
 +
[[Category:RAREMETALWORKER]]
 
== Useful Wiki Pages ==
 
== Useful Wiki Pages ==
   Line 53: Line 55:  
* We next analyse the second study using a similar command:
 
* We next analyse the second study using a similar command:
   −
   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
 
                     --inverseNormal --makeResiduals --kinSave --kinGeno --prefix STUDY2
   Line 74: Line 76:  
   tabix -c "#" -s 1 -b 2 -e 2 STUDY2.QT1.singlevar.cov.txt.gz
 
   tabix -c "#" -s 1 -b 2 -e 2 STUDY2.QT1.singlevar.cov.txt.gz
   −
* We next create two text files that will drive the meta-analysis. The first file lists the input files with summary statistics. Let's call it summaryfiles. In most Linux workstations, you can use the command pico or nano to create this file. These should be the contents of "summaryfiles":
+
* We next create two text files that will drive the meta-analysis. The first file lists the input files with summary statistics. Let's call it summaryfiles. In most Linux workstations, you can use the command pico or nano to create this file. These should be the contents of "summaryfiles":
    
   STUDY1.QT1.singlevar.score.txt.gz
 
   STUDY1.QT1.singlevar.score.txt.gz
32

edits

Navigation menu