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

Hello, I am trying to get the XP-EHH values by using HP-EXX. (https://github.com/joepickrell/xpehh) . 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!!

Regards

Vic

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.

Help
Access

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