ChunkChromosome: Difference between revisions
From Genome Analysis Wiki
Jump to navigationJump to search
| Line 14: | Line 14: | ||
<source lang="bash"> | <source lang="bash"> | ||
mach1 -d chr1.dat -p chr1.ped --rounds 20 --states 200 --phase --interim 5 --sample 5 | mach1 -d chr1.dat -p chr1.ped --rounds 20 --states 200 --phase --interim 5 --sample 5 --prefix chr1 | ||
minimac --refHaps | minimac --refHaps 1000genomes.chr1.hap.gz --refSnps 1000genomes.chr1.snps --haps chr1.hap.gz --snps chr1.snps --rounds 5 --states 200 --prefix imputation-results | ||
</source> | </source> | ||
== Download == | == Download == | ||
Revision as of 07:26, 5 August 2011
ChunkChromosome is a helper utility for minimac and MaCH. It can be used to facilitate analyses of very large datasets in overlapping slices.
Parameters
ChunkChromosome expects three parameters:
- A data file (specified with the -d command line option), listing all the markers along one chromosome. The data file can optionally include phenotype and other information, which is safely ignored.
- A desired core chunk size, in markers (specified with the -n command line option and defaulting to 5000 markers).
- A desired overlap between chunks, also in markers (specified with the -o command line option and defaulting to 500 markers).
Usage
Suppose you plan to run 1000 Genomes Imputation using MaCH and minimac. Typically, you'd accomplish this by running the two commands:
mach1 -d chr1.dat -p chr1.ped --rounds 20 --states 200 --phase --interim 5 --sample 5 --prefix chr1
minimac --refHaps 1000genomes.chr1.hap.gz --refSnps 1000genomes.chr1.snps --haps chr1.hap.gz --snps chr1.snps --rounds 5 --states 200 --prefix imputation-results
