Question: Error in tmp_haplo[hap2_index, ] <- hap2
gravatar for kas
14 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 • 228 views
ADD COMMENTlink modified 14 months ago by zx87549.7k • written 14 months ago by kas0

Please use edit to add more details.

ADD REPLYlink written 14 months ago by zx87549.7k
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: 1923 users visited in the last hour