Difference between revisions of "METAL Quick Start"
Line 33: | Line 33: | ||
|} | |} | ||
− | + | Although some effort was made to harmonize analysis strategies (for example, by excluding individuals with a diagnosis of diabetes as well as other individuals with elevated fasting glucose levels), you will notice that the three files are formatted somewhat differently. In order to combine results across studies, one critical piece of information that METAL will need are details of the formating used for each file. You probably also noticed that the FUSION input files are a bit smaller than those for the other studies. This is because they have been compressed with [www.gzip.org gzip]. This is not a problem, because METAL can transparently handle gzip-compressed files. | |
+ | |||
+ | == The METAL Script == | ||
+ | |||
+ | We will now walk through the [[METAL Glucose Example Script]]. The first thing to know is that METAL scripts can include comments and that these are indicating by using a hash sign '''#''' as the first character in a line. Thus: | ||
− | + | # This is a comment. |
Revision as of 07:05, 30 March 2010
Learning by Example
Many people find that it easiest to learn by example. So, if reading through our extensive documentation is not for you, you might find it more appealing to walk through a simple METAL analysis.
The Glucose Data
This example examines evidence for association between fasting glucose levels and genetic markers in the G6PC2 (chr. 2), GCK (chr. 7) and MTNR1B (chr. 11) regions. It uses results from 3 genomewide association studies: FUSION, SardiNIA and DGI. Genetic variants in the three loci impact fasting glucose levels and, in the case of MTRN1B, also impact the risk of type 2 diabetes.
You can download a copy of the input files used in this example analysis from the METAL Download Page.
Input Files
Initially, each of the three studies analyzed association between fasting glucose levels and genotyped and imputed variants in each of the three loci. Study specific results are stored in the following files:
'STUDY | INPUT FILE | FILE SIZE |
FUSION | MAGIC_FUSION_Results.txt.gz | 46 kb |
SardiNIA | magic_SARDINIA.tbl | 188 kb |
DGI | DGI_three_regions.txt | 188 kb |
Although some effort was made to harmonize analysis strategies (for example, by excluding individuals with a diagnosis of diabetes as well as other individuals with elevated fasting glucose levels), you will notice that the three files are formatted somewhat differently. In order to combine results across studies, one critical piece of information that METAL will need are details of the formating used for each file. You probably also noticed that the FUSION input files are a bit smaller than those for the other studies. This is because they have been compressed with [www.gzip.org gzip]. This is not a problem, because METAL can transparently handle gzip-compressed files.
The METAL Script
We will now walk through the METAL Glucose Example Script. The first thing to know is that METAL scripts can include comments and that these are indicating by using a hash sign # as the first character in a line. Thus:
# This is a comment.