Question: Reformat minimac imputated files to plink file
3
gravatar for Floris Brenk
5.4 years ago by
Floris Brenk890
USA
Floris Brenk890 wrote:

Hello all,

I am imputating my genotyped dataset, genotyped with illumina hapmap610 chip with minimac and mach  according to standard protocols: http://genome.sph.umich.edu/wiki/Minimac:_1000_Genomes_Imputation_Cookbook

I was wondering if anyone knows how to convert the output files of this imputation packages into standard plink format (ped and map files) with as side note that some filtering can be performed on the probability score that minimac gives as output.

Thanks in advance!

 

plink imputation mach minimac • 6.0k views
ADD COMMENTlink modified 3.1 years ago by ISB0 • written 5.4 years ago by Floris Brenk890

I have met the same question !! Convert the .dos, .info formate of mimimac to the plink binary formate !! Have you solved that ??

ADD REPLYlink written 3.2 years ago by xujiayuan3225740610
3
gravatar for Floris Brenk
5.4 years ago by
Floris Brenk890
USA
Floris Brenk890 wrote:

Found it:  http://www.complextraitgenomics.com/software/gcta/data_management.html

gcta64  --dosage-mach  test.mldose.gz  test.mlinfo.gz  --make-bed --out test

 

 

ADD COMMENTlink written 5.4 years ago by Floris Brenk890

I have tried gcta64, but the generated .bim formate didn't have the rsID, only the chr:bp position, still can't the following analysis? Do you know how to convert to the binary .bim formate?

ADD REPLYlink written 3.2 years ago by xujiayuan3225740610

Is it because the dose file does not have rsid? There is an option --rsid in Minimac3, which will use rsid instead of chr:bp in the output files including dose and info file.

ADD REPLYlink written 2.9 years ago by Xu19120
2
gravatar for Maxime Lamontagne
5.4 years ago by
Québec
Maxime Lamontagne2.1k wrote:

I think dose2geno is the software you need (http://www.mybiosoftware.com/population-genetics/4571).

ADD COMMENTlink written 5.4 years ago by Maxime Lamontagne2.1k

                

ADD REPLYlink modified 5.4 years ago • written 5.4 years ago by Floris Brenk890
1
gravatar for oyhel
3.5 years ago by
oyhel20
oyhel20 wrote:

I am wondering if you managed to figure out how to filter only high quality imputed genotypes from minimac3 output. I am struggling with the same problem you did. The dose2geno does provide best guess genotypes, but it does not seem to consider low quality scores. Using PLINKs vcf-min-gp flag is not working since plink expects 0-1 scale GP but minimac3 outputs phred scale. Any help would be highly appreciated if you managed to sort this!

ADD COMMENTlink written 3.5 years ago by oyhel20
1

gcta --imput-rsq 0.3 is OK

ADD REPLYlink written 3.2 years ago by xujiayuan3225740610

Hmm cant really remember how sorry, but I did filtered out all snps with a r2 score of lower than 0.3

ADD REPLYlink written 3.5 years ago by Floris Brenk890
0
gravatar for Floris Brenk
5.4 years ago by
Floris Brenk890
USA
Floris Brenk890 wrote:

Yes this looks indeed interesting. Although still a bit of reformatting needs to be done to get the actual .ped and .map files... One thing that still misses is that the probability score of the minimac output is not used now and I dont know what dose2geno does with imputed low probability scored genotypes...

ADD COMMENTlink written 5.4 years ago by Floris Brenk890
0
gravatar for xujiayuan32257406
3.2 years ago by
xujiayuan3225740610 wrote:

I have met the same question !! Convert the .dos, .info formate of mimimac to the plink binary formate !! Have you solved that ??

ADD COMMENTlink written 3.2 years ago by xujiayuan3225740610
0
gravatar for ISB
3.1 years ago by
ISB0
ISB0 wrote:

Use fcgene to convert minimac output to PED/MAP files. You can also use rsq as filter: fcgene --minimac-prob minimac_output.prob --minimac-infominimac_output.info --oformat plink --rsq 0.3 --out plink_files

ADD COMMENTlink written 3.1 years ago by ISB0

Thanks but fcgene requires pedinfo.txt and snpinfo.txt as input files. How can I prepare these 2 files using Minimac3 output files? And Minimac3 seems to not output .prob file.

ADD REPLYlink written 2.9 years ago by xliu425510
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: 1287 users visited in the last hour