Wednesday 13 December 2017

Thinking about the bootstrap


  1. Bootstrap samples experimental units but in phylogenetics you sample the VARIABLES - sites.
  2. How should we treat sites?
    1. Remove totall variant?
    2. Remove sites where a row is missing?
  3. You cannot say that parametric and non-parametric are the same thing. They are correlated but not directly comparable.
    1. Carry out FastTree with H5N8, then H5 then N8
    2. Use the parametric and non-parametric bootstraps
    3. Use the CONSEL measures as well.
  4. Having more bootstraps than 100 makes NO difference to the bootstrap values. They converge quickly empirically.
    1. This is far below the theoretical numbers needed by Efron says that this is usual.
    2. Suggests that sites are linked and so there is less independent variability than it appears.
    3. Need to experiment with conserved sites.
    4. Need to experiment with the substitution models to look at sensitivity and also gamma.
  5. There is a lack of independence between sites in the evolutionary models but this is IGNORED in the bootstrap calculations. You should bootstrap codons and not individual bases.

Need to create synthetic data where the true tree is known. This can be used to test:
  1. Effects of sampling by censoring the data.
  2. Evaluate modeltest.
  3. Check trees from bad evolutionary models against the best models (probably the same!!!)

No comments: