Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft of selective_sweep API #828

Merged
merged 1 commit into from Nov 16, 2021
Merged

Draft of selective_sweep API #828

merged 1 commit into from Nov 16, 2021

Conversation

molpopgen
Copy link
Owner

No description provided.

@molpopgen
Copy link
Owner Author

cc @apragsdale

Copy link
Collaborator

@apragsdale apragsdale left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I didn't run these vignettes locally, but they look very helpful! Just a few comments (which I may be confused about)

pdict = {
"recregions": [fwdpy11.PoissonInterval(0, 1, 1e-1)],
# Here, 2 means that fitness is multiplicative
# over 1, 1+hs, 1+2s.
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Should this be 1, 1+hs, 1+s or 1, 1+2hs, 1+2s, with h=0.5 meaning additive selection? (maybe I'm getting confused with how it interacts with Multiplicative(2.0)...)

Copy link
Owner Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

It is always 1+hs, so h changes range as scaling changes. yes, that is likely a sticking point...

pdict = {
"recregions": [fwdpy11.PoissonInterval(0, 1, 1e-1)],
# Here, 2 means that fitness is multiplicative
# over 1, 1+hs, 1+2s.
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Same here?

# ancestry, which we can get either
# from a burn-in or from msprime.
initial_ts = msprime.sim_ancestry(
samples=500,
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Do you want samples to equal population_size?

Copy link
Owner Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I need to double-check this.

# ancestry, which we can get either
# from a burn-in or from msprime.
initial_ts = msprime.sim_ancestry(
samples=500,
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Same here?

Copy link
Owner Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

again, will double-check. Some of this is confusion due to the change in defaults to diploidy for sim_ancestry.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

2 participants