This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error in tmp_haplo[hap2_index, ] <- hap2

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.

rehh fastphaseformat r ihs

Please use edit to add more details.

1 answer

If anybody has used the iHS method for selection signature analysis, please provide those scripts to do those.

Log in to answer this question.