Wesley R. Elsberry
Posts: 4937 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)
where 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 groundtruthing, 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
