next up previous
Next: Playing with Substitution Matrices Up: Calculating Evolutionary Distance, Part Previous: Calculating % Change

Using the Poisson Correction

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