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