Question: the XP-EHH instruction
gravatar for Victor
9 months ago by
Victor10 wrote:

Hello, I am trying to get the XP-EHH values by using HP-EXX. ( . I only have map and ped plink files in two populations. So 1. I used plink function --recode vcf, to make vcf files (preQC had been done, plus --snp-only -just acgt). 2. Eeagle was applied to phase my dataset. (it should have no-missing data in my output files, and all alleles become 0 and 1 only) 3 Then I wrote a R-script to format a map files for XP-EHH.

I have some questions about XP-EHH. 1. How many map file I should make for XP-EHH command? The example ./xpehh -m example.ihsmap.................

But I have 2 populations. Should I merge the 2 map files into one? Or I should use the overlap snp? Why I ask this it is because I try to use the big map file merged from 2 map files and always get the error message like:"borken the pipeline or out of range........" 2. I check the ihshap file.I saw no row name or column name in the files. I am wondering how do the script process the files. The output of example shown the chr ID. Sould I make the same order for ihsmap and ihshap?

Are there anyone has experience in process XP-EHH

Thank you for your time!!



phased plink haplotype xp-ehh vcf • 551 views
ADD COMMENTlink modified 9 months ago • written 9 months ago by Victor10
Please log in to add an answer.


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