MaCH FAQ - Revision history

Ppwhite: /* How do I get imputation quality estimates? */

2017-02-02T15:22:49Z

How do I get imputation quality estimates?

← Older revision		Revision as of 15:22, 2 February 2017
Line 62:		Line 62:
	Estimated per allele error rate is 0.0293		Estimated per allele error rate is 0.0293

	A better approach is to mask a small proportion of SNPs (vs. genotypes in the above simple approach). One can generate a mask.dat from the original .dat file by simply changing the flag of a subset of markers from M to S2 without duplicating the .ped file. Post-imputation, one can use   [http://genome.sph.umich.edu/wiki/CalcMatch CalcMatch ]and [http://csg.sph.umich.edu//ylwtx/doseR2.tgz doseR2.pl ]to estimate genotypic/allelic error rate and correlation respectively. Both programs can be downloaded from [http://csg.sph.umich.edu//ylwtx/software.html http://~~www~~.sph.umich.edu/~~csg~~/ylwtx/software.html].		A better approach is to mask a small proportion of SNPs (vs. genotypes in the above simple approach). One can generate a mask.dat from the original .dat file by simply changing the flag of a subset of markers from M to S2 without duplicating the .ped file. Post-imputation, one can use   [http://genome.sph.umich.edu/wiki/CalcMatch CalcMatch ]and [http://csg.sph.umich.edu//ylwtx/doseR2.tgz doseR2.pl ]to estimate genotypic/allelic error rate and correlation respectively. Both programs can be downloaded from [http://csg.sph.umich.edu//ylwtx/software.html http://csg.sph.umich.edu//ylwtx/software.html].

	'''Warning''': Imputation involving masked datasets should be performed separately for imputation quality estimation. For production, one should use all available information.		'''Warning''': Imputation involving masked datasets should be performed separately for imputation quality estimation. For production, one should use all available information.

Ppwhite: /* How do I get imputation quality estimates? */

2017-02-02T15:22:19Z

How do I get imputation quality estimates?

← Older revision		Revision as of 15:22, 2 February 2017
Line 62:		Line 62:
	Estimated per allele error rate is 0.0293		Estimated per allele error rate is 0.0293

	A better approach is to mask a small proportion of SNPs (vs. genotypes in the above simple approach). One can generate a mask.dat from the original .dat file by simply changing the flag of a subset of markers from M to S2 without duplicating the .ped file. Post-imputation, one can use   [http://genome.sph.umich.edu/wiki/CalcMatch CalcMatch ]and [http://~~www~~.sph.umich.edu/~~csg~~/ylwtx/doseR2.tgz doseR2.pl ]to estimate genotypic/allelic error rate and correlation respectively. Both programs can be downloaded from [http://~~www~~.sph.umich.edu/~~csg~~/ylwtx/software.html http://www.sph.umich.edu/csg/ylwtx/software.html].		A better approach is to mask a small proportion of SNPs (vs. genotypes in the above simple approach). One can generate a mask.dat from the original .dat file by simply changing the flag of a subset of markers from M to S2 without duplicating the .ped file. Post-imputation, one can use   [http://genome.sph.umich.edu/wiki/CalcMatch CalcMatch ]and [http://csg.sph.umich.edu//ylwtx/doseR2.tgz doseR2.pl ]to estimate genotypic/allelic error rate and correlation respectively. Both programs can be downloaded from [http://csg.sph.umich.edu//ylwtx/software.html http://www.sph.umich.edu/csg/ylwtx/software.html].

	'''Warning''': Imputation involving masked datasets should be performed separately for imputation quality estimation. For production, one should use all available information.		'''Warning''': Imputation involving masked datasets should be performed separately for imputation quality estimation. For production, one should use all available information.

Ppwhite: /* Install MaCH */

2017-02-02T15:21:46Z

Install MaCH

← Older revision		Revision as of 15:21, 2 February 2017
Line 216:		Line 216:

	== Install MaCH ==		== Install MaCH ==
	We have source codes available through the MaCH download page: http://~~www~~.sph.umich.edu/~~csg~~/yli/mach/download/ <br>		We have source codes available through the MaCH download page: http://csg.sph.umich.edu//yli/mach/download/ <br>

	== More questions? ==		== More questions? ==

	Email [mailto:yunli@med.unc.edu Yun Li] or [mailto:goncalo@umich.edu Goncalo Abecasis].		Email [mailto:yunli@med.unc.edu Yun Li] or [mailto:goncalo@umich.edu Goncalo Abecasis].

Ppwhite: /* How do I get reference files for an region of interest? */

2017-02-02T15:21:15Z

How do I get reference files for an region of interest?

← Older revision		Revision as of 15:21, 2 February 2017
Line 79:		Line 79:
	== How do I get reference files for an region of interest? ==		== How do I get reference files for an region of interest? ==

	Note that you do not need to extract regional pedigree files for your own samples because SNPs in pedigree but not in reference will be automatically discarded. <br> 1. For HapMapII format, download haplotypes from http://~~www~~.sph.umich.edu/~~csg~~/ylwtx/HapMapForMach.tgz <br> 2. For MACH format, you can do the following:		Note that you do not need to extract regional pedigree files for your own samples because SNPs in pedigree but not in reference will be automatically discarded. <br> 1. For HapMapII format, download haplotypes from http://csg.sph.umich.edu//ylwtx/HapMapForMach.tgz <br> 2. For MACH format, you can do the following:

	*First, find the first and last SNP in the region you are interested in. Say "rsFIRST" and "rsLAST", defined according to position.		*First, find the first and last SNP in the region you are interested in. Say "rsFIRST" and "rsLAST", defined according to position.

Ppwhite: /* How do I interpret the imputation quality estimates? */

2017-02-02T15:20:54Z

How do I interpret the imputation quality estimates?

← Older revision		Revision as of 15:20, 2 February 2017
Line 69:		Line 69:
	In the simple approach, you will only get concordance/error estimates. There are two aspects to check. (1) the ratio between the genotypic error and allelic error. We expect that only a small proportion of errors where one homozygote is imputed as the other homozygote. Therefore, a ~2:1 ratio is expected. (2) the absolute error rate. There are several factors influencing imputation quality including the population to be imputed, the reference population and the genotyping panel used. Typically, we expect <2% allelic error rate among Caucasians and East Asians; 3-5% among Africans and African Americans. Figure below show imputation quality from the Human Genome Diversity Project (HGDP) for 52 populations across the world and by different HapMap reference panel.		In the simple approach, you will only get concordance/error estimates. There are two aspects to check. (1) the ratio between the genotypic error and allelic error. We expect that only a small proportion of errors where one homozygote is imputed as the other homozygote. Therefore, a ~2:1 ratio is expected. (2) the absolute error rate. There are several factors influencing imputation quality including the population to be imputed, the reference population and the genotyping panel used. Typically, we expect <2% allelic error rate among Caucasians and East Asians; 3-5% among Africans and African Americans. Figure below show imputation quality from the Human Genome Diversity Project (HGDP) for 52 populations across the world and by different HapMap reference panel.

	http://~~www~~.sph.umich.edu/~~csg~~/yli/figure3.gif		http://csg.sph.umich.edu//yli/figure3.gif

	Table 3 in the MaCH 1.0 paper tabulates imputation quality by commercial panel in CEU, YRI, and CHB+JPT.		Table 3 in the MaCH 1.0 paper tabulates imputation quality by commercial panel in CEU, YRI, and CHB+JPT.

Ppwhite: /* Where can I find HapMap III / 1000 Genomes reference files? */

2017-02-02T15:20:33Z

Where can I find HapMap III / 1000 Genomes reference files?

← Older revision		Revision as of 15:20, 2 February 2017
Line 48:		Line 48:
	== Where can I find HapMap III / 1000 Genomes reference files? ==		== Where can I find HapMap III / 1000 Genomes reference files? ==

	You can find these at the MaCH download page, which is at http://~~www~~.sph.umich.edu/~~csg~~/yli/mach/download/		You can find these at the MaCH download page, which is at http://csg.sph.umich.edu//yli/mach/download/

	== Does --mle overwrite input genotypes? ==		== Does --mle overwrite input genotypes? ==

Ppwhite: /* Where can I find combined HapMap reference files? */

2017-02-02T15:20:16Z

Where can I find combined HapMap reference files?

← Older revision		Revision as of 15:20, 2 February 2017
Line 44:		Line 44:
	== Where can I find combined HapMap reference files? ==		== Where can I find combined HapMap reference files? ==

	You can find them at http://~~www~~.sph.umich.edu/~~csg~~/yli/mach/download/HapMap-r21.html or on the HapMap Project website.		You can find them at http://csg.sph.umich.edu//yli/mach/download/HapMap-r21.html or on the HapMap Project website.

	== Where can I find HapMap III / 1000 Genomes reference files? ==		== Where can I find HapMap III / 1000 Genomes reference files? ==

Ppwhite: /* minimac */

2017-02-02T15:19:37Z

minimac

← Older revision		Revision as of 15:19, 2 February 2017
Line 4:		Line 4:

	This is the new 2-step procedure we are recommending, particularly considering people that are performing imputation multiple times (using HapMap as reference, or using updated releases of the 1000 Genomes data as reference). <br>		This is the new 2-step procedure we are recommending, particularly considering people that are performing imputation multiple times (using HapMap as reference, or using updated releases of the 1000 Genomes data as reference). <br>
	The first step is a pre-phasing step using MaCH. This step does not need external reference. This is a time-consuming step BUT is a one-time investment. For computational reason, we recommend breaking the genome into small overlapping segments ( [http://genome.sph.umich.edu/wiki/MaCH_FAQ#Divide_and_Conquer Divide-and-Conquer]) for this step. In general, we recommend >500Kb overlapping region on each side. For example, for Affymetrix 6.0 panel, if we use core region of 10Mb and flanking/overlapping region of 1Mb on each side, it will correspond to ~3500 SNps in the core region and ~350 SNPs on each side. For 2000 individuals, one job with ~4,200 SNPs running with --states 200 and -r 50, this would take ~40 hours. For other combinations, using the following link to estimate computing time [http://~~www~~.sph.umich.edu/~~csg~~/yli/MaCH-Admix/runtime.php#est runtime estimate]. <br>		The first step is a pre-phasing step using MaCH. This step does not need external reference. This is a time-consuming step BUT is a one-time investment. For computational reason, we recommend breaking the genome into small overlapping segments ( [http://genome.sph.umich.edu/wiki/MaCH_FAQ#Divide_and_Conquer Divide-and-Conquer]) for this step. In general, we recommend >500Kb overlapping region on each side. For example, for Affymetrix 6.0 panel, if we use core region of 10Mb and flanking/overlapping region of 1Mb on each side, it will correspond to ~3500 SNps in the core region and ~350 SNPs on each side. For 2000 individuals, one job with ~4,200 SNPs running with --states 200 and -r 50, this would take ~40 hours. For other combinations, using the following link to estimate computing time [http://csg.sph.umich.edu//yli/MaCH-Admix/runtime.php#est runtime estimate]. <br>
	The second step is the actual imputation step using minimac. This step can run on whole chromosomes. Regarding computing time, one million markers for 1000 individuals using 100 reference haplotypes takes ~ 1 hour; and computing time increases linearly with all the above three parameters. See [http://genome.sph.umich.edu/wiki/Minimac minimac] for details.		The second step is the actual imputation step using minimac. This step can run on whole chromosomes. Regarding computing time, one million markers for 1000 individuals using 100 reference haplotypes takes ~ 1 hour; and computing time increases linearly with all the above three parameters. See [http://genome.sph.umich.edu/wiki/Minimac minimac] for details.

Ylwtx: /* MaCH-Admix */

2012-04-24T20:13:10Z

MaCH-Admix

← Older revision		Revision as of 20:13, 24 April 2012
Line 9:		Line 9:
	=== MaCH-Admix ===		=== MaCH-Admix ===

	If you are doing imputation only ~~once~~ (think twice if this is really true) or under an immediate time pressure, you can use MaCH-Admix, which does not require pre-phased data and takes ~1/7 of the computing time of that typically needed for pre-phasing. For large dataset, we recommend breaking the genome into small overlapping segments ( [http://genome.sph.umich.edu/wiki/MaCH_FAQ#Divide_and_Conquer Divide-and-Conquer]). Details see [http://www.unc.edu/~yunmli/MaCH-Admix/ MaCH-Admix].		If you are doing imputation only a few (<5) times (think twice if this is really true) or under an immediate time pressure, you can use MaCH-Admix, which does not require pre-phased data and takes ~1/7 of the computing time of that typically needed for pre-phasing. For large dataset, we recommend breaking the genome into small overlapping segments ( [http://genome.sph.umich.edu/wiki/MaCH_FAQ#Divide_and_Conquer Divide-and-Conquer]). Details see [http://www.unc.edu/~yunmli/MaCH-Admix/ MaCH-Admix].

	=== Divide and Conquer ===		=== Divide and Conquer ===

Ylwtx: /* MaCH-Admix */

2012-04-24T20:12:33Z

MaCH-Admix

← Older revision		Revision as of 20:12, 24 April 2012
Line 9:		Line 9:
	=== MaCH-Admix ===		=== MaCH-Admix ===

	If you are doing imputation only once (think twice if this is really true) or under an immediate time pressure, you can use MaCH-Admix, which does not require pre-phased data and takes ~1/7 of the computing time of ~~the original MaCH~~. For large dataset, we recommend breaking the genome into small overlapping segments ( [http://genome.sph.umich.edu/wiki/MaCH_FAQ#Divide_and_Conquer Divide-and-Conquer]). Details see [http://www.unc.edu/~yunmli/MaCH-Admix/ MaCH-Admix].		If you are doing imputation only once (think twice if this is really true) or under an immediate time pressure, you can use MaCH-Admix, which does not require pre-phased data and takes ~1/7 of the computing time of that typically needed for pre-phasing. For large dataset, we recommend breaking the genome into small overlapping segments ( [http://genome.sph.umich.edu/wiki/MaCH_FAQ#Divide_and_Conquer Divide-and-Conquer]). Details see [http://www.unc.edu/~yunmli/MaCH-Admix/ MaCH-Admix].

	=== Divide and Conquer ===		=== Divide and Conquer ===