Dear respected members,
Anybody using "rehh" package, kindly help me. While calculating iHS values, its showing like this,
hh_map <- data2haplohh(hap_file = "bta1__hapguess_switch.out", map_file = "12breedsmap.inp", allele_coding = "map",chr.name = 1,popsel = 12, verbose = FALSE,) * Reading input file(s) * Map info: 1770 markers declared for chromosome 1. Haplotype input file in fastPHASE format assumed. Haplotypes in the fastPHASE output file originate from 12 populations. Error in tmp_haplo[hap2_index, ] <- hap2 : number of items to replace is not a multiple of replacement length
If anybody has used the iHS method for selection signature analysis, please provide those scripts to do those.