Question: the XP-EHH instruction
0
gravatar for Victor
2.4 years 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 • 1.2k views
ADD COMMENTlink modified 5 months ago by shanshuwen900 • written 2.4 years ago by Victor10

Hi, I have the same question, do you have a solution now?

ADD REPLYlink written 5 months ago by shanshuwen900
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: 1073 users visited in the last hour