Folded Spectrum in python is not recognizing some of the individuals in my VCF file

38 views
Skip to first unread message

Morgan L. Lynn

unread,
Oct 29, 2024, 4:00:06 PM10/29/24
to dadi-user
Hello, with a few of my species in my pairwise analyses no data comes out after putting the vcf file into a folded spectrum. When I change the amount of called alleles from 10 (the true value) to 8 ( a random value I inputed) the folded spectrum is visible but only accounts for 2, or 3 of the alleles. This occurs with two other species where if I adjust the amount of called alleles in my populations the folded spectrum appears to have data, but when I use true values the folded spectrums are devoid of anything. Please help, I have attached examples of the issues I am seeing. dadi-issues-adjusted-alleles.pngdadi-problem-no-data.png

Ryan Gutenkunst

unread,
Oct 29, 2024, 4:50:44 PM10/29/24
to dadi-user
Hello Morgan,

This may occur if you often have missing data. For a site to enter the spectrum for 10 alleles, you must have at least 10 successfully called alleles at that site. Is it possible (for example, due to low coverage) that you often have missing calls?

Best,
Ryan

On Oct 29, 2024, at 11:19 AM, 'Morgan L. Lynn' via dadi-user <dadi...@googlegroups.com> wrote:

Hello, with a few of my species in my pairwise analyses no data comes out after putting the vcf file into a folded spectrum. When I change the amount of called alleles from 10 (the true value) to 8 ( a random value I inputed) the folded spectrum is visible but only accounts for 2, or 3 of the alleles. This occurs with two other species where if I adjust the amount of called alleles in my populations the folded spectrum appears to have data, but when I use true values the folded spectrums are devoid of anything. Please help, I have attached examples of the issues I am seeing. <dadi-issues-adjusted-alleles.png><dadi-problem-no-data.png>

--
You received this message because you are subscribed to the Google Groups "dadi-user" group.
To unsubscribe from this group and stop receiving emails from it, send an email to dadi-user+...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/dadi-user/bf38b17e-7b7a-447a-a8d0-0fbec03b4cdbn%40googlegroups.com.
<dadi-problem-no-data.png><dadi-issues-adjusted-alleles.png>

Reply all
Reply to author
Forward
0 new messages