Skip to content

Commit 6f46aeb

Browse files
Reformatted commandlines for mcmctree doc.
1 parent 8964655 commit 6f46aeb

1 file changed

Lines changed: 45 additions & 22 deletions

File tree

doc/Dating.md

Lines changed: 45 additions & 22 deletions
Original file line numberDiff line numberDiff line change
@@ -8,8 +8,20 @@ icon: info-circle
88
doctype: tutorial
99
tags:
1010
- tutorial
11-
description: "Building time tree with node dates on phylogenetic trees."
11+
description: "Inferring time trees"
1212
sections:
13+
- name: Bayesian dating with MCMCtree
14+
url: bayesian-dating-with-mcmctree
15+
- name: IQ2MC workflow for time tree inference
16+
url: iq2mc-workflow-for-time-tree-inference
17+
- name: Estimating the gradients and the Hessian for MCMCtree dating
18+
url: estimating-the-gradients-and-the-hessian-for-mcmctree-dating
19+
- name: Using partitions and Mixture models for approximate likelihood dating
20+
url: partitions-and-mixture-models-for-approximate-likelihood-dating
21+
- name: How to run MCMCtree
22+
url: how-to-run-mcmctree
23+
- name: Least Square Dating (LSD2)
24+
url: least-square-dating
1325
- name: Inferring time tree with tip dates
1426
url: inferring-time-tree-with-tip-dates
1527
- name: Calibrating tree using ancestral dates
@@ -77,7 +89,10 @@ approximate likelihood dating.
7789
If the alignment file is called `example.phy` and the rooted tree file is called
7890
`example_tree.nwk`,
7991

80-
iqtree -s example.phy -m GTR+G4 -te example_tree.nwk --dating mcmctree --prefix example
92+
```
93+
iqtree -s example.phy -m GTR+G4 -te example_tree.nwk --dating mcmctree --prefix example
94+
```
95+
8196

8297
Note that, Here we generate the Hessian file for a fixed rooted tree. You can
8398
directly input the rooted tree which already contains fossil/tip calibration
@@ -102,7 +117,9 @@ input to MCMCtree to save compute.
102117
You can specify more parameters in the workflow to generate the control file
103118
accurately for the analysis with IQ-TREE.
104119

105-
iqtree -s example.phy -m GTR+G4 -te example_tree.nwk --dating mcmctree --mcmc-iter 20000,200,50000 --mcmc-bds 1,1,0.5 --mcmc-clock IND
120+
```
121+
iqtree -s example.phy -m GTR+G4 -te example_tree.nwk --dating mcmctree --mcmc-iter 20000,200,50000 --mcmc-bds 1,1,0.5 --mcmc-clock IND
122+
```
106123

107124
* `--mcmc-iter burnin,samplefreq,nsample` : use to set number of burin samples,
108125
sample frequency and number of MCMC samples in the control file. In the above
@@ -117,7 +134,7 @@ Currently supported clocks models are EQUAL: global clock with equal rates, IND:
117134
independent rates model with independent rates across lineages and CORR:
118135
correlated clock model with auto-correlated rates across the lineages.
119136

120-
Using partitions and Mixture models for approximate likelihood dating.
137+
Using partitions and Mixture models for approximate likelihood dating
121138
---------------------------------------------------------------------
122139

123140
IQ-TREE supports three partition models for approximate likelihood dating. Under
@@ -129,25 +146,28 @@ Hessian as branches are shared across partitions.
129146
Since IQ-TREE supports RAxML and NEXUS style partitions input file, you can use
130147
partitions defined in the following format.
131148

132-
DNA, part1 = 1-100
133-
DNA, part2 = 101-450
134-
149+
```
150+
DNA, part1 = 1-100
151+
DNA, part2 = 101-450
152+
```
135153
If your partition file is called `example.nex`,
136154

137-
iqtree -s example.phy -Q example.nex -m GTR+G4 -te example_tree.nwk --dating mcmctree
138-
155+
```
156+
iqtree -s example.phy -Q example.nex -m GTR+G4 -te example_tree.nwk --dating mcmctree
157+
```
139158

140159
Here, IQ-TREE generates the Hessian file using the `GTR+G4` model for all
141160
partitions. If you need to use different models for each partition, you need to
142161
create a more flexible NEXUS file like the following.
143162

144-
#nexus
145-
begin sets;
146-
charset part1 = 1-100;
147-
charset part2 = 101-450;
148-
charpartition mine = GTR+G4:part1, HKY:part2;
149-
end;
150-
163+
```
164+
#nexus
165+
begin sets;
166+
charset part1 = 1-100;
167+
charset part2 = 101-450;
168+
charpartition mine = GTR+G4:part1, HKY:part2;
169+
end;
170+
```
151171
Here, IQ-TREE uses `GTR+G4` model for partition 1, and `HKY` model for partition
152172
2 respectively. Using `-q` and `-p` options, you can generate the Hessian file
153173
which considers `edge-linked equal branch partition models` and `edge-linked
@@ -156,12 +176,14 @@ proportional branch length models` respectively.
156176
IQ-TREE also supports mixture models for the Hessian file generation. You can
157177
simply specify DNA or Amino Acid Mixture model as following,
158178

159-
iqtree -s example.phy -m "MIX{GTR,HKY}+G4" -te example_tree.nwk –-dating mcmctree
160-
179+
```
180+
iqtree -s example.phy -m "MIX{GTR,HKY}+G4" -te example_tree.nwk –-dating mcmctree
181+
```
161182
If you need to use an Amino Acid profile mixture model such as C60 model,
162183

163-
iqtree -s example.phy -m LG+G4+C60 -te example_tree.nwk –-dating mcmctree
164-
184+
```
185+
iqtree -s example.phy -m LG+G4+C60 -te example_tree.nwk –-dating mcmctree
186+
```
165187
If you are using ModelFinder or MixtureFinder, you need to follow a two-step
166188
approach. First, you can estimate the best-fit model for the data using
167189
ModelFinder or MixtureFinder. Then, the Hessian file can be generated using
@@ -173,8 +195,9 @@ How to run MCMCtree
173195
You can directly run MCMCtree from the control file generated by IQ-TREE in step
174196
2. The command to run MCMCtree with the control file is,
175197

176-
mcmctree example.mcmctree.ctl
177-
198+
```
199+
mcmctree example.mcmctree.ctl
200+
```
178201

179202

180203
The control file generated by IQ-TREE has the following format. You can simply

0 commit comments

Comments
 (0)