Bug not catched

[VincentGardeux]

Hello,

We recently realized that the VCF we were using was not correct (i.e. not same samples than in the count matrices).
However the phm program did not find a bug and ran and even produced results.

Could you implement a check (maybe on the column names of the VCF / count matrix) to be sure that the correct genotypes are taken into account? And raise an error, if this is not the case?

Thanks in advance

[natsuhiko]

Hi Vincent,

Sorry for the late reply.

Because we implemented the input count data as a binary file, it is not easy to perform the sample consistency check at this moment. However, I will consider this point in the next update.

Thanks for the useful feedback,
Natsuhiko