consensusHap: num2base failed #1
Description
Dir professor Li :
The R vresion : 3.5.3 . I run the command in Rstudio.
I have tested the built-in data in Hapi package , whilch can perform good , and get the same consensusHap result .
while I use 5 gametes data ,which is my own data and as same format as Example gmt.rda , to find the crossover pisition . Use this command to input 5 gametes data :
gmt <- read.table(file="test-5-gmt.txt",rownames=1).Then I use the same command to run ,but I get the different result when I transed the consensusHap 0/1 to A/T/C/G . The result is the same as this：
The result shows :
1. the hap 0/1 trans to the 1/2/3/4 but not the A/T/C/G ,
2. at some locations it shows the hap1 same to the hap2 .
The file I have uploaded .
so I'm confused. Would you please explain it to me? Thank you very much !