Question: Error in tmp_haplo[hap2_index, ] <- hap2
gravatar for kas
9 months ago by
kas0 wrote:

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", = 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.

ihs rehh fastphaseformat R • 196 views
ADD COMMENTlink modified 9 months ago by zx87549.3k • written 9 months ago by kas0

Please use edit to add more details.

ADD REPLYlink written 9 months ago by zx87549.3k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 788 users visited in the last hour