Question: Imputation Without A Reference Panel
3
gravatar for 14134125465346445
6.0 years ago by
United Kingdom
141341254653464453.4k wrote:

I would like to know if there are there any tools to impute missing SNPs for a resequencing dataset of hundreds of individuals for which there is no reference panel?

The closest I got is:

(1) in using a subset of the most complete individuals with the "include reference" MACH (http://www.sph.umich.edu/csg/abecasis/MaCH) option:

INCLUDE REFERENCE (e.g. HAPMAP) GENOTYPES TO YOUR DATASET:

If you select this option, you should simply create one large pooled dataset. Some individuals will have missing data and others will have much more complete genotyping information.

In addition to estimating the most likely genotype for each individual, you can use the command line options --dosage and --quality options to request additional information about each inferred genotype.

Other tools I looked at:

(2) in using Impute (https://mathgen.stats.ox.ac.uk/impute/impute_v2.html):

We have proposed a simple and universal solution to this problem: we provide all available reference haplotypes to IMPUTE2, then let the software choose a "custom" reference panel for each individual to be imputed. There are several advantages to this approach

Beagle (http://faculty.washington.edu/browning/beagle/beagle.html) doesn't seem to mention options for including the reference in the input.

snps imputation • 4.6k views
ADD COMMENTlink modified 5.1 years ago by kumar.vinod81280 • written 6.0 years ago by 141341254653464453.4k
1

Mach and Beagle do not require a reference panel. You will known that if you try them once.

ADD REPLYlink written 6.0 years ago by lh331k

duplicate of Does Imputation Always Need Reference Panel ?

ADD REPLYlink written 6.0 years ago by Michael Dondrup45k
0
gravatar for pd3
6.0 years ago by
pd3340
pd3340 wrote:

How many sites do you need to fill? Both IMPUTE2 and BEAGLE can be used without reference genotypes. Of course, the more "holes" (missing data) there are, the less accurate the method get.

ADD COMMENTlink modified 6.0 years ago • written 6.0 years ago by pd3340
0
gravatar for kumar.vinod81
5.1 years ago by
kumar.vinod81280
New Delhi
kumar.vinod81280 wrote:

What command line to use for missing genotype imputation in MACH without a reference panel. Is this?

To estimate model parameters ...

prompt> mach1 -d sample.dat -p sample.ped -s hapmap.snps -h hapmap.haplos --greedy --rounds 10 --prefix round1

To fill in missing genotypes ...

prompt> mach1 -d sample.dat -p sample.ped -s hapmap.snps -h hapmap.haplos --greedy --errormap round1.erate --cross round1.rec --mle --mldetails

Or we have to left some commands......

Thanks

ADD COMMENTlink written 5.1 years ago by kumar.vinod81280
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2234 users visited in the last hour