Reconstructing phylogeny by quadratically approximated maximum likelihood.

M D Woodhams, M D Hendy
Author Information
  1. M D Woodhams: Allan Wilson Centre for Molecular Ecology and Evolution, Massey University, Palmerston North, New Zealand. m.d.woodhams@massey.ac.nz

Abstract

Maximum likelihood (ML) for phylogenetic inference from sequence data remains a method of choice, but has computational limitations. In particular, it cannot be applied for a global search through all potential trees when the number of taxa is large, and hence a heuristic restriction in the search space is required. In this paper, we derive a quadratic approximation, QAML, to the likelihood function whose maximum is easily determined for a given tree. The derivation depends on Hadamard conjugation, and hence is limited to the simple symmetric models of Kimura and of Jukes and Cantor. Preliminary testing has demonstrated the accuracy of QAML is close to that of ML.

MeSH Term

Algorithms
Chromosome Mapping
Computer Simulation
Evolution, Molecular
Likelihood Functions
Models, Genetic
Models, Statistical
Phylogeny
Sequence Analysis, DNA