Question: Error in tmp_haplo[hap2_index, ] <- hap2
0
gravatar for kas
4 weeks ago by
kas0
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",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.

ihs rehh fastphaseformat R • 103 views
ADD COMMENTlink modified 4 weeks ago by zx87548.3k • written 4 weeks ago by kas0

Please use edit to add more details.

ADD REPLYlink written 4 weeks ago by zx87548.3k
Please log in to add an answer.

Help
Access

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