|Wesley R. Elsberry
Joined: May 2002
|Quote (WebHopper @ Feb. 11 2017,05:27)|
|Wesley, if I understand you well the expected correct bases after mutation is|
E = C + (u * (L - C) / K) - (u * C * (K - 1) / K)
C = number of correct bases
K = alphabet length
L = number of base pairs
u = mutation rate
I guess with "expected" E you refer to the mean. So if C=L we have
E = L (K + u - K u)/K
This is a linear law with respect to L if the other parameters are held constant. However, this does not match a Monte Carlo simulation, which shows that E is exponential with respect to L. So there must be some error in your calculation or I don't understand you well.
DiEb: Thanks for the link, I will look into this. Seems to be complicated however...
Edit: I'll look this over again. I derived it with Monte Carlo ground-truthing, so I'm not sure where we are diverging in expectation.
Edited by Wesley R. Elsberry on Feb. 11 2017,16:35
"You can't teach an old dogma new tricks." - Dorothy Parker