gamma distribution/PHYLIP

Kent E. Holsinger kent at darwin.darwin.eeb.uconn.edu
Thu Feb 29 09:30:08 EST 1996


>>>>> "Joe" == Joe Felsenstein <joe at evolution.genetics.washington.edu>
writes: In article <4giqf9$m06 at nntp3.u.washington.edu>
joe at evolution.genetics.washington.edu (Joe Felsenstein) writes:


    Joe> In likelihood methods (such as Ziheng Yang's methods that use
    Joe> gamma-distributed rates in his PAML package, or my Hidden
    Joe> Markov Model approach in DNAML) one can change the rate
    Joe> variation parameters until the overall likelihood of the tree
    Joe> is maximized.  With distance methods, I am not as sure how to
    Joe> do this.  Perhaps one could try various values until the
    Joe> goodness of fit of the tree (say, the sum of squares) is
    Joe> optimized, but doing this may have biases in it.

I don't think comparing goodness-of-fit values would work, unless they
were normalized in some way. The distance between sequences A and B
increases as the coefficient of variation increases, so the total tree
length would also increase as the coefficient of variation
increases. It's not obvious to me that there is a scale-invariant way
of comparing the goodness-of-fit statistics, but maybe I'm missing
something. (It wouldn't be the first time.)

-- Kent


-- 
Kent E. Holsinger                Kent at Darwin.EEB.UConn.Edu
-- Department of Ecology & Evolutionary Biology          
-- University of Connecticut, U-43                                       
-- Storrs, CT   06269-3043                                               



More information about the Mol-evol mailing list