Choosing projection values for a 2D SFS

196 views
Skip to first unread message

Aidan Short

unread,
Aug 3, 2021, 4:53:11 PM8/3/21
to dadi-user
Hi Prof. Gutenkunst,

This isn't really an issue, but kind of a general question about choosing projection values. In the past I have tried to maximize the number of segregating sites based on the segregating SNPs in the two populations.
For example:
fs = dadi.Spectrum.from_data_dict(dd, ['population 1', 'population 2'], [100,50])
print(fs.S())

However, I noticed that tools such as easySFS suggest selecting the projection values based on the segregating SNPs in each of the single populations. 
For example: 
fs = dadi.Spectrum.from_data_dict(dd, ['population 1'], [100])
print(fs.S())
fs = dadi.Spectrum.from_data_dict(dd, ['population 2'], [50])
print(fs.S())

I guess I am wondering if it would be better for me to maximize the number of segregating sites within each population or between the two populations.

Sorry if this confusing, I am still new to using dadi.

Sincerely, 

Aidan Short

Ryan Gutenkunst

unread,
Aug 3, 2021, 6:37:07 PM8/3/21
to dadi-user
Hello Aidan,

In our own analyses, we tend to maximize the segregating sites across the whole spectrum. But to my knowledge no one has rigorously tested this rule of thumb. So do whichever seems more logical to you; I doubt it makes a big difference in any reasonable case.

Best,
Ryan
> --
> 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 on the web visit https://groups.google.com/d/msgid/dadi-user/ef2741ab-71c9-4d10-b236-f377a267ba2dn%40googlegroups.com.

Aidan Short

unread,
Aug 4, 2021, 9:30:31 AM8/4/21
to dadi...@googlegroups.com
Thanks Ryan! I appreciate your advice!

Sincerely,
Aidan

You received this message because you are subscribed to a topic in the Google Groups "dadi-user" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/dadi-user/mfgoorieqTE/unsubscribe.
To unsubscribe from this group and all its topics, send an email to dadi-user+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/dadi-user/1685EC2C-F55A-4AEE-AB58-3996C8742962%40gmail.com.
Reply all
Reply to author
Forward
0 new messages