Hi everyone,
I am trying to use the package rehh to check selection signature. I've got some inappropriate results when I use the scan_hh function. So I tried to apply this function using the example dataset. You can see the results below. Any recommendation how to solve the problem?
> make.example.files()
> data(haplohh_cgu_bta12)
> res.scan<-scan_hh(haplohh_cgu_bta12)
> dim(res.scan)
[1] 1424 7
> head(res.scan)
CHR POSITION freq_A iHH_A iHH_D iES_Tang_et_al_2007 iES_Sabeti_et_al_2007
F1200140 12 79823 0.1500000 NA NA NA NA
F1200150 12 125974 0.4071429 NA NA NA NA
F1200170 12 175087 0.3571429 NA NA NA NA
F1200180 12 219152 0.2214286 NA NA NA NA
F1200190 12 256896 0.1750000 NA NA NA NA
F1200210 12 316254 0.3892857 NA NA NA NA