millanek / Dsuite

Fast calculation of Patterson's D (ABBA-BABA) and the f4-ratio statistics across many populations/species
160 stars 26 forks source link

Dinvestigate error with vcf. Not sure why #88

Closed drew-olson-tmru closed 8 months ago

drew-olson-tmru commented 8 months ago

Hi,

I am trying to use Dinvestigate to identify introgression between two populations (F. fujikuroi and F. prolfieratum), however, I am having issues with the tool. I think it may have to do with the VCF file.

The text files are: https://github.com/drew-olson-tmru/grape_isolate_analysis/blob/main/dsuite.sets.txt and https://github.com/drew-olson-tmru/grape_isolate_analysis/blob/main/dsuite_trios_fol.txt

The VCF file is too large to include on Github.

Here is the error below.

There are 3 sets (populations/species) excluding the Outgroup /builddir/build/BUILD/gcc-11.4.1-20230605/obj-x86_64-redhat-linux/x86_64-redhat-linux/libstdc++-v3/include/bits/basic_string.h:1040: std::cxx11::basic_string<_CharT, _Traits, _Alloc>::const_reference std::cxx11::basic_string<_CharT, _Traits, _Alloc>::operator[](std::cxx11::basic_string<_CharT, _Traits, _Alloc>::size_type) const [with _CharT = char; _Traits = std::char_traits; _Alloc = std::allocator; std::cxx11::basic_string<_CharT, _Traits, _Alloc>::const_reference = const char&; std::__cxx11::basic_string<_CharT, _Traits, _Alloc>::size_type = long unsigned int]: Assertion '__pos <= size()' failed. Aborted (core dumped)

Please let me know if you can determine an issue based on this, otherwise I can provide more direct info from the vcf file.