Next: Playing with Substitution Matrices
Up: Calculating Evolutionary Distance, Part
Previous: Calculating % Change
In this exercise, you'll use your % Change values to compute corrected
evolutionary distances. You'll need to find someway to compute the
natural log. You can use a command called perl to do that, if you
don't have a calculator. The formula for the Poisson
corrected-distance D is D = -log (1 - %ID/100).
To compute this in perl, say for a % Change of 15%, use the following
command:
% perl -e 'print (- log(.85)),"\n"'
Write down your corrected values in your table. Do the same for the
real data. Do the values fit evolutionary distances for the
simulated data better than the % Change values?
David Ardell
2005-01-26