1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
module Probability.Distribution.Tree.Moves where
import MCMC
import Probability.Random
addSPRMoves rate tree = do
addMove rate $ sampleSPRAll tree
addMove (rate/2) $ sampleSPRFlat tree
addMove (rate/2) $ sampleSPRNodes tree
addTopologyMoves rate tree = do
addSPRMoves rate tree
addMove rate $ walkTreeSampleNNIandBranchLengths tree
addMove (2*rate) $ walkTreeSampleNNI tree -- if alignment is fixed this is really cheap -- increase weight?
addMove (0.5*rate) $ walkTreeSampleNNIandA tree
addLengthMoves rate tree = do
addMove rate $ walkTreeSampleBranchLengths tree
addTreeMoves rate tree = do
addTopologyMoves (rate*2) tree
addLengthMoves (rate*2) tree
|
