You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
I was using msprime for the last weeks an i wondered:
How can I set the recombination rate in msprime to ensure a specific expected number of recombination events occurs in the entire Ancestral Recombination Graph (ARG)?
I thought i could calculate the number of expected recombination events (E) like this: E=0.0667×30=2.001.
However, the code above results in over 700 nodes with the flag 131072 which means 350 recombination events right?
Im very happy about some responses since i got quite confused xD Thanks in advance!
reacted with thumbs up emoji reacted with thumbs down emoji reacted with laugh emoji reacted with hooray emoji reacted with confused emoji reacted with heart emoji reacted with rocket emoji reacted with eyes emoji
-
Hi there!
I was using msprime for the last weeks an i wondered:
How can I set the recombination rate in msprime to ensure a specific expected number of recombination events occurs in the entire Ancestral Recombination Graph (ARG)?
Heres what i tried so far:
ts = msprime.sim_ancestry(3,ploidy = 1,population_size = 100, recombination_rate = 0.0667, sequence_length = 30,additional_nodes = msprime.NodeType.RECOMBINANT,coalescing_segments_only = False)
I thought i could calculate the number of expected recombination events (E) like this: E=0.0667×30=2.001.
However, the code above results in over 700 nodes with the flag 131072 which means 350 recombination events right?
Im very happy about some responses since i got quite confused xD Thanks in advance!
Best,
NewAtl
Beta Was this translation helpful? Give feedback.
All reactions