Entering edit mode

2.1 years ago

taniamahmood38
▴
50

Dear all,

I am using rehh package on R software to compute site frequency spectrum (sfs) statistics on a single chromosome of my pooled-seq data of 12 individuals. First, I converted my unphased vcf file to an object of class haplohh using the command

```
hp_vcf <- data2haplohh(hap_file = "sample.vcf", polarize_vcf = FALSE, vcf_reader = "vcfR", chr.name = "CM004588.1", verbose = TRUE)
```

Then, I used the function calc_sfs_tests. My command was

```
Sample <- calc_sfs_tests(hp_vcf, polarized = FALSE, window_size = 150000, overlap = 75000, right = TRUE, verbose = TRUE)
```

Below is the image of output generated.

Now, I am stuck at this point. Kindly help me out that why there are zeros in last three columns and same values in theta's columns? where I am getting wrong? and what is the solution now? Any help would be highly appreciated.