Question: the XP-EHH instruction
gravatar for Victor
22 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 • 969 views
ADD COMMENTlink modified 22 months ago • written 22 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: 1618 users visited in the last hour