Assertion `RSDDataset->numberOfSamples>=1'

[ChrisK1988]

Hello,

I have been using RAiSD on a modern human dataset with ~4100 individuals to assess positive selection segments across the genome. I have separated my VCFs by population and am getting an error with only one population. It errors out on the following command with the below issue and I am not sure how to fix it.

Command: ./RAiSD -n FIN -I Path/to/VCF/FIN_only.ALLChroms.vcf -R

RAiSD: sources/RAiSD_Dataset.c:1806: RSDDataset_getNumberOfSamples_vcf: Assertion `RSDDataset->numberOfSamples>=1' failed.
Aborted

The VCF files were created at different times but from the same original data source. They followed the same construction commands so I am not sure why this specific population is throwing an error. Based on the error it is suggesting that I have some issues with the samples, but that is not the case as the samples appear to be correct and have corresponding phased genotypes (like the rest of the data did).

Thank you for your help!

Chris

[alachins]

Hi Chris, This error might be caused when the vcf contains no samples or no SNPs. There has to be at least one valid SNP in your data to prevent this assertion from failing. Nikos

…

On Fri, Apr 12, 2024 at 6:37 PM ChrisK1988 ***@***.***> wrote: Hello, I have been using RAiSD on a modern human dataset with ~4100 individuals to assess positive selection segments across the genome. I have separated my VCFs by population and am getting an error with only one population. It errors out on the following command with the below issue and I am not sure how to fix it. Command: ./RAiSD -n FIN -I Path/to/VCF/FIN_only.ALLChroms.vcf -R *RAiSD: sources/RAiSD_Dataset.c:1806: RSDDataset_getNumberOfSamples_vcf: Assertion `RSDDataset->numberOfSamples>=1' failed. Aborted* The VCF files were created at different times but from the same original data source. They followed the same construction commands so I am not sure why this specific population is throwing an error. Thank you for your help! Chris — Reply to this email directly, view it on GitHub <#51>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AALKWCTQ6TOBAYV6F7ZANDLY5AES5AVCNFSM6AAAAABGEOEAA2VHI2DSMVQWIX3LMV43ASLTON2WKOZSGI2DANJQG42TSNI> . You are receiving this because you are subscribed to this thread.Message ID: ***@***.***>

-- Nikolaos Alachiotis

[ChrisK1988]

Hi. Sorry for the very late comment back. The VCF file does contain samples and SNPs. I am not sure the underlying mechanics behind it, but it seems like the issue was resolved by updating a version of some of the dependencies, which is interesting because it worked before without the update.

At any rate, it is resolved and I appreciate the response.