Error in tmp_haplo[hap2_index, ] <- hap2
1
0
Entering edit mode
4.6 years ago
kas • 0

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 • 722 views
ADD COMMENT
0
Entering edit mode

Please use edit to add more details.

ADD REPLY

Login before adding your answer.

Traffic: 1899 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6