I have a large vcf file from next gen sequencing. I am not able to segregate DNA sites using the seg.sites function.

I am able to load my vcf file using Discolor<-read.vcfR("file_name"), but when I try to find segregated sites in my DNA using seg.sites(Discolor) I get Error in latag2n(x) : no method for coercing this S4 class to a vector.

When I do class(Discolor), I get the following.

[1] "vcfR"
attr(,"package")
[1] "vcfR"
0

There are 0 best solutions below