Difference between revisions of "R Hackathon 1/Trait Evolution SG"

From Phyloinformatics
Jump to: navigation, search
Line 5: Line 5:
 
# Investigate statistical power of discriminating among various models
 
# Investigate statistical power of discriminating among various models
 
# Identify gaps in current implementation
 
# Identify gaps in current implementation
Accomplishments
+
== Accomplishments. ==
 
# Evaluated the results of continuous character analyses in different packages
 
# Evaluated the results of continuous character analyses in different packages
 
#* Packages are mostly consistent
 
#* Packages are mostly consistent

Revision as of 16:26, 24 January 2008

  • Participants: Harmon, Hipp, Hunt

Targets.

  1. Compare various implementations of the same methods (ape, geiger, OUCH, Mesquite)
  2. Improve functionality of character fitting in r
  3. Investigate statistical power of discriminating among various models
  4. Identify gaps in current implementation

Accomplishments.

  1. Evaluated the results of continuous character analyses in different packages
    • Packages are mostly consistent
    • Discrepancies come from two sources:
      • Different approaches (e.g. marginal versus joint likelihood)
      • Difficulties in finding the ML solution
  2. Improve functionality or, at least, interpretability of output
    • GEIGER was modified to give more reliable results by a more thorough search of the likelihood surface
    • Some ape functions seem (to us) unreliable for large trees

Test

    • Clarified why one sometimes gets different results from different programs
    • Modified GEIGER package so that results for fitting models of character evolution are more robust
    • Modified functions to rescale trees and variance-covariance matrices for non-ultrametric trees