Minimac: Tutorial - Revision history

Goncalo: /* Imputation into Phased Haplotypes - minimac(2) */

2017-01-25T16:34:40Z

Imputation into Phased Haplotypes - minimac(2)

← Older revision		Revision as of 12:34, 25 January 2017
Line 49:		Line 49:
	=== Imputation into Phased Haplotypes - minimac(2)===		=== Imputation into Phased Haplotypes - minimac(2)===

	Imputing genotypes using '''minimac(2)''' is a straightforward process: after selecting a set of reference haplotypes, plugging-in the target haplotypes from the previous step and setting the number of rounds to use for estimating model parameters (which describe the length and conservation of haplotype stretches shared between the reference panel and your study samples), imputation should proceed rapidly. ~~Because marker names can change between dbSNP versions~~, it ~~is usually a good idea to include~~ '~~'aliases'' file that provides mappings between earlier marker names and~~ the ~~current preferred name for each polymorphism~~.		Imputing genotypes using '''minimac(2)''' is a straightforward process: after selecting a set of reference haplotypes, plugging-in the target haplotypes from the previous step and setting the number of rounds to use for estimating model parameters (which describe the length and conservation of haplotype stretches shared between the reference panel and your study samples), imputation should proceed rapidly.

			Minimac needs a file listing the variants in your sample. If your directory already includes a "sample.snps" file, no worries. If it doesn't, you can generate one using "sample.dat" as input with the following command:

			cut -f 2 -d " " sample.dat > sample.snps

	The minimac command line would look like this:		The minimac command line would look like this:

Cfuchsb at 07:38, 28 January 2015

2015-01-28T07:38:44Z

← Older revision		Revision as of 03:38, 28 January 2015
Line 3:		Line 3:
	== Getting Started ==		== Getting Started ==

	Download [http://~~www~~.sph.umich.edu~~/csg~~/abecasis/MaCH/download/ MaCH] and [http://genome.sph.umich.edu/wiki/Minimac#Download Minimac] or [http://genome.sph.umich.edu/wiki/Minimac2#Download Minimac2]. Furthermore, example data used in this tutorial can be found [http://~~www~~.sph.umich.edu~~/csg~~/cfuchsb/minimac2_example.tgz here]		Download [http://csg.sph.umich.edu/abecasis/MaCH/download/ MaCH] and [http://genome.sph.umich.edu/wiki/Minimac#Download Minimac] or [http://genome.sph.umich.edu/wiki/Minimac2#Download Minimac2]. Furthermore, example data used in this tutorial can be found [http://csg.sph.umich.edu/cfuchsb/minimac2_example.tgz here]

	== Minimac and Minimac2 Imputation ==		== Minimac and Minimac2 Imputation ==
Line 15:		Line 15:
	./mach1 -d sample.dat -p sample.ped --rounds 20 --states 50 --phase --interim 5 --sample 5 --prefix sample.pp \| tee mach.log		./mach1 -d sample.dat -p sample.ped --rounds 20 --states 50 --phase --interim 5 --sample 5 --prefix sample.pp \| tee mach.log

	This will request that MaCH estimate haplotypes for your sample, using 20 iterations of its Markov sampler and conditioning each update on up to 50 haplotypes. A summary description of these parameters follows (but for a more complete description, you should go to the [http://~~www~~.sph.umich.edu~~/csg~~/abecasis/MaCH/ MaCH website]):		This will request that MaCH estimate haplotypes for your sample, using 20 iterations of its Markov sampler and conditioning each update on up to 50 haplotypes. A summary description of these parameters follows (but for a more complete description, you should go to the [http://csg.sph.umich.edu/abecasis/MaCH/ MaCH website]):

	{\| class="wikitable" border="1" cellpadding="2"		{\| class="wikitable" border="1" cellpadding="2"
Line 23:		Line 23:
	\|-		\|-
	\|style=white-space:nowrap\|<code>-d sample.dat</code>		\|style=white-space:nowrap\|<code>-d sample.dat</code>
	\| Data file in [http://~~www~~.sph.umich.edu~~/csg~~/abecasis/Merlin/tour/input_files.html Merlin format]. Markers should be listed according to their order along the chromosome.		\| Data file in [http://csg.sph.umich.edu/abecasis/Merlin/tour/input_files.html Merlin format]. Markers should be listed according to their order along the chromosome.
	\|-		\|-
	\| <code>-p sample.ped</code>		\| <code>-p sample.ped</code>
	\| Pedigree file in [http://~~www~~.sph.umich.edu~~/csg~~/abecasis/Merlin/tour/input_files.html Merlin format]. Alleles should be labeled on the forward strand.		\| Pedigree file in [http://csg.sph.umich.edu/abecasis/Merlin/tour/input_files.html Merlin format]. Alleles should be labeled on the forward strand.
	\|-		\|-
	\| <code>--states 200</code>		\| <code>--states 200</code>

Cfuchsb: /* Imputation into Phased Haplotypes - minimac(2) */

2014-09-16T07:24:05Z

Imputation into Phased Haplotypes - minimac(2)

← Older revision		Revision as of 03:24, 16 September 2014
Line 57:		Line 57:
	or		or

	./minimac2 --refHaps hapmap.hap --refSnps hapmap.snps --haps sample.pp.gz --snps sample.snps --prefix ~~sample~~.imp \| tee minimac2.log		./minimac2 --refHaps hapmap.hap --refSnps hapmap.snps --haps sample.pp.gz --snps sample.snps --prefix sample2.imp \| tee minimac2.log


Line 101:		Line 101:
	or		or

	./minimac2-omp --cpus 4 --refHaps hapmap.hap --refSnps hapmap.snps --haps sample.pp.gz --snps sample.snps --prefix ~~sample~~.imp \| tee minimac-omp.log		./minimac2-omp --cpus 4 --refHaps hapmap.hap --refSnps hapmap.snps --haps sample.pp.gz --snps sample.snps --prefix sample2.imp \| tee minimac-omp.log

	== Imputation quality evaluation ==		== Imputation quality evaluation ==

Cfuchsb: /* Getting Started */

2014-09-16T07:23:41Z

Getting Started

← Older revision		Revision as of 03:23, 16 September 2014
Line 3:		Line 3:
	== Getting Started ==		== Getting Started ==

	Download [http://www.sph.umich.edu/csg/abecasis/MaCH/download/ MaCH] and [http://genome.sph.umich.edu/wiki/Minimac#Download Minimac] or [http://genome.sph.umich.edu/wiki/Minimac2#Download Minimac2]. Furthermore, example data used in this tutorial can be found [http://www.sph.umich.edu/csg/~~abecasis/downloads~~/minimac2_example.tgz here]		Download [http://www.sph.umich.edu/csg/abecasis/MaCH/download/ MaCH] and [http://genome.sph.umich.edu/wiki/Minimac#Download Minimac] or [http://genome.sph.umich.edu/wiki/Minimac2#Download Minimac2]. Furthermore, example data used in this tutorial can be found [http://www.sph.umich.edu/csg/cfuchsb/minimac2_example.tgz here]

	== Minimac and Minimac2 Imputation ==		== Minimac and Minimac2 Imputation ==

Cfuchsb: /* Getting Started */

2014-09-16T07:23:19Z

Getting Started

← Older revision		Revision as of 03:23, 16 September 2014
Line 3:		Line 3:
	== Getting Started ==		== Getting Started ==

	Download [http://www.sph.umich.edu/csg/abecasis/MaCH/download/ MaCH] and [http://genome.sph.umich.edu/wiki/Minimac#Download Minimac] or [http://genome.sph.umich.edu/wiki/Minimac2#Download Minimac2]. Furthermore, example data used in this tutorial can be found [http://www.sph.umich.edu/csg/abecasis/downloads/~~minimac_example~~.~~tar.gz~~ here]		Download [http://www.sph.umich.edu/csg/abecasis/MaCH/download/ MaCH] and [http://genome.sph.umich.edu/wiki/Minimac#Download Minimac] or [http://genome.sph.umich.edu/wiki/Minimac2#Download Minimac2]. Furthermore, example data used in this tutorial can be found [http://www.sph.umich.edu/csg/abecasis/downloads/minimac2_example.tgz here]

	== Minimac and Minimac2 Imputation ==		== Minimac and Minimac2 Imputation ==

Cfuchsb: /* Minimac Imputation */

2014-09-16T06:41:42Z

Minimac Imputation

← Older revision		Revision as of 02:41, 16 September 2014
Line 5:		Line 5:
	Download [http://www.sph.umich.edu/csg/abecasis/MaCH/download/ MaCH] and [http://genome.sph.umich.edu/wiki/Minimac#Download Minimac] or [http://genome.sph.umich.edu/wiki/Minimac2#Download Minimac2]. Furthermore, example data used in this tutorial can be found [http://www.sph.umich.edu/csg/abecasis/downloads/minimac_example.tar.gz here]		Download [http://www.sph.umich.edu/csg/abecasis/MaCH/download/ MaCH] and [http://genome.sph.umich.edu/wiki/Minimac#Download Minimac] or [http://genome.sph.umich.edu/wiki/Minimac2#Download Minimac2]. Furthermore, example data used in this tutorial can be found [http://www.sph.umich.edu/csg/abecasis/downloads/minimac_example.tar.gz here]

	== Minimac Imputation ==		== Minimac and Minimac2 Imputation ==

	[[Minimac]] and [[Minimac2]] relies on a two step approach. First, the samples that are to be analyzed must be phased into a series of estimated haplotypes. Second, imputation is carried out directly into these phased haplotypes. As newer reference panels become available, only the second step must be repeated.		[[Minimac]] and [[Minimac2]] relies on a two step approach. First, the samples that are to be analyzed must be phased into a series of estimated haplotypes. Second, imputation is carried out directly into these phased haplotypes. As newer reference panels become available, only the second step must be repeated.

Cfuchsb: /* Reference */

2014-09-16T06:41:22Z

Reference

← Older revision		Revision as of 02:41, 16 September 2014
Line 115:		Line 115:
	= Reference =		= Reference =

	If you use minimac, please cite:		If you use minimac or minimac2, please cite:

	Howie B, Fuchsberger C, Stephens M, Marchini J, and Abecasis GR.		Howie B, Fuchsberger C, Stephens M, Marchini J, and Abecasis GR.

Cfuchsb: /* Minimac Imputation */

2014-09-16T06:40:51Z

Minimac Imputation

← Older revision		Revision as of 02:40, 16 September 2014
Line 7:		Line 7:
	== Minimac Imputation ==		== Minimac Imputation ==

	[[Minimac]] relies on a two step approach. First, the samples that are to be analyzed must be phased into a series of estimated haplotypes. Second, imputation is carried out directly into these phased haplotypes. As newer reference panels become available, only the second step must be repeated.		[[Minimac]] and [[Minimac2]] relies on a two step approach. First, the samples that are to be analyzed must be phased into a series of estimated haplotypes. Second, imputation is carried out directly into these phased haplotypes. As newer reference panels become available, only the second step must be repeated.

	=== Pre-phasing - MaCH ===		=== Pre-phasing - MaCH ===
Line 47:		Line 47:
	\|}		\|}

	=== Imputation into Phased Haplotypes - minimac ===		=== Imputation into Phased Haplotypes - minimac(2)===

	Imputing genotypes using '''minimac''' is a straightforward process: after selecting a set of reference haplotypes, plugging-in the target haplotypes from the previous step and setting the number of rounds to use for estimating model parameters (which describe the length and conservation of haplotype stretches shared between the reference panel and your study samples), imputation should proceed rapidly. Because marker names can change between dbSNP versions, it is usually a good idea to include ''aliases'' file that provides mappings between earlier marker names and the current preferred name for each polymorphism.		Imputing genotypes using '''minimac(2)''' is a straightforward process: after selecting a set of reference haplotypes, plugging-in the target haplotypes from the previous step and setting the number of rounds to use for estimating model parameters (which describe the length and conservation of haplotype stretches shared between the reference panel and your study samples), imputation should proceed rapidly. Because marker names can change between dbSNP versions, it is usually a good idea to include ''aliases'' file that provides mappings between earlier marker names and the current preferred name for each polymorphism.

	The minimac command line would look like this:		The minimac command line would look like this:
Line 55:		Line 55:
	./minimac --refHaps hapmap.hap --refSnps hapmap.snps --haps sample.pp.gz --snps sample.snps --prefix sample.imp \| tee minimac.log		./minimac --refHaps hapmap.hap --refSnps hapmap.snps --haps sample.pp.gz --snps sample.snps --prefix sample.imp \| tee minimac.log

	A detailed description of all minimac options is available [[Minimac Command Reference\|elsewhere]]. Here is a brief description of the above parameters:		or

			./minimac2 --refHaps hapmap.hap --refSnps hapmap.snps --haps sample.pp.gz --snps sample.snps --prefix sample.imp \| tee minimac2.log


			A detailed description of all minimac(2) options is available [[Minimac Command Reference\|elsewhere]]. Here is a brief description of the above parameters:

	{\| class="wikitable" border="1" cellpadding="2"		{\| class="wikitable" border="1" cellpadding="2"
Line 90:		Line 95:


	You can speed-up things by running minimac in parallel by launching the [http://genome.sph.umich.edu/wiki/~~Minimac~~#Multiprocessor_Version ~~minimac~~-omp] version. On our cluster 4 cpus per minimac is optimal (--cpus 4).		You can speed-up things by running minimac in parallel by launching the [http://genome.sph.umich.edu/wiki/Minimac2#Multiprocessor_Version minimac2-omp] version. On our cluster 4 cpus per minimac(2) is optimal (--cpus 4).

	./minimac-omp --cpus 4 --refHaps hapmap.hap --refSnps hapmap.snps --haps sample.pp.gz --snps sample.snps --prefix sample.imp \| tee minimac-omp.log		./minimac-omp --cpus 4 --refHaps hapmap.hap --refSnps hapmap.snps --haps sample.pp.gz --snps sample.snps --prefix sample.imp \| tee minimac-omp.log

			or

			./minimac2-omp --cpus 4 --refHaps hapmap.hap --refSnps hapmap.snps --haps sample.pp.gz --snps sample.snps --prefix sample.imp \| tee minimac-omp.log

	== Imputation quality evaluation ==		== Imputation quality evaluation ==

Cfuchsb: /* Getting Started */

2014-09-16T06:39:02Z

Getting Started

← Older revision		Revision as of 02:39, 16 September 2014
Line 3:		Line 3:
	== Getting Started ==		== Getting Started ==

	Download [http://www.sph.umich.edu/csg/abecasis/MaCH/download/ MaCH] and [http://genome.sph.umich.edu/wiki/Minimac#Download Minimac]. Furthermore, example data used in this tutorial can be found [http://www.sph.umich.edu/csg/abecasis/downloads/minimac_example.tar.gz here]		Download [http://www.sph.umich.edu/csg/abecasis/MaCH/download/ MaCH] and [http://genome.sph.umich.edu/wiki/Minimac#Download Minimac] or [http://genome.sph.umich.edu/wiki/Minimac2#Download Minimac2]. Furthermore, example data used in this tutorial can be found [http://www.sph.umich.edu/csg/abecasis/downloads/minimac_example.tar.gz here]

	== Minimac Imputation ==		== Minimac Imputation ==

Cfuchsb at 06:38, 16 September 2014

2014-09-16T06:38:32Z

← Older revision		Revision as of 02:38, 16 September 2014
Line 1:		Line 1:
	Before reading this tutorial, you might find it useful to spend a few minutes reading through the main [[Minimac]] documentation.		Before reading this tutorial, you might find it useful to spend a few minutes reading through the main [[Minimac]] and [[Minimac2]] documentation.

	== Getting Started ==		== Getting Started ==