RSS 2.0 Feed

» Welcome Guest Log In :: Register

Pages: (14) < [1] 2 3 4 5 6 ... >   
  Topic: Evolutionary Computation, Stuff that drives AEs nuts< Next Oldest | Next Newest >  
Wesley R. Elsberry



Posts: 4937
Joined: May 2002

(Permalink) Posted: Feb. 11 2017,16:30   

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 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

    
  399 replies since Mar. 17 2009,11:00 < Next Oldest | Next Newest >  

Pages: (14) < [1] 2 3 4 5 6 ... >   


Track this topic Email this topic Print this topic

[ Read the Board Rules ] | [Useful Links] | [Evolving Designs]