Skip to content

Commit 5934f31

Browse files
committed
Update figures and figure widths
1 parent 5ef94e9 commit 5934f31

6 files changed

+10
-8
lines changed

README.md

+10-8
Original file line numberDiff line numberDiff line change
@@ -13,6 +13,8 @@ subtitle: Prior selection and clock calibration using Influenza A data.
1313

1414

1515

16+
17+
1618
# Background
1719

1820
In the Bayesian analysis of sequence data, priors play an important role. When priors are not specified correctly, it may cause runs to take a long time to converge, not converge at all or cause a bias in the inferred trees and model parameters. Specifying proper priors and starting values is crucial and can be a difficult exercise in the beginning. It is not always easy to pick a proper model of tree generation (tree prior), substitution model, molecular clock model or the prior distribution for an unknown parameter.
@@ -125,7 +127,7 @@ Since we will be using the birth-death skyline model (**BDSKY**) {% cite Stadler
125127
126128
<figure>
127129
<a id="packageManage2"></a>
128-
<img style="width:75.0%;" src="figures/beast2_package_manager.png" alt="">
130+
<img src="figures/beast2_package_manager.png" alt="">
129131
<figcaption>Figure 2: The BEAST2 Package Manager.</figcaption>
130132
</figure>
131133
<br>
@@ -151,7 +153,7 @@ We will first analyse the alignment of sequences sampled through time (heterochr
151153
152154
<figure>
153155
<a id="importAlignment"></a>
154-
<img style="width:75.0%;" src="figures/beast2_import_alignment.png" alt="">
156+
<img src="figures/beast2_import_alignment.png" alt="">
155157
<figcaption>Figure 3: Importing the alignment into BEAUti.</figcaption>
156158
</figure>
157159
<br>
@@ -176,7 +178,7 @@ The sequence labels (headers in the FASTA file) contain sampling times specified
176178
177179
<figure>
178180
<a id="timeUnitsAndFlow"></a>
179-
<img style="width:100.0%;" src="figures/beast2_time_specification.png" alt="">
181+
<img src="figures/beast2_time_specification.png" alt="">
180182
<figcaption>Figure 4: Specifying time units and direction of time flow.</figcaption>
181183
</figure>
182184
<br>
@@ -192,13 +194,13 @@ You could specify the tip dates by hand, by clicking for each row (i.e. for each
192194
193195
<figure>
194196
<a id="guessDates"></a>
195-
<img style="width:50.0%;" src="figures/beast2_guess_dates.png" alt="">
197+
<img style="width:75.0%;" src="figures/beast2_guess_dates.png" alt="">
196198
<figcaption>Figure 5: Specifying tip dates.</figcaption>
197199
</figure>
198200
<br>
199201

200202

201-
You should now see that the tip ages have been filled in for all of the taxa with the **Date (raw value)** columns showing the date strings extracted from the taxon names, and the **Age/Height** column showing numbers on the order of 0.1 (the age in years of each tip relative to the most recent sample).
203+
You should now see that the tip ages have been filled in for all of the taxa with the **Date (raw value)** columns showing the date strings extracted from the taxon names, and the **Age/Height** column showing numbers on the order of 0.1 (the age in years of each tip relative to the most recent sample). If you did everything correctly, the sequence with the _most recent_ sampling date (2009/06/08) should have the samllest height (0.0 in this case).
202204

203205
Now we are done with the data specification and we are about to start specifying models and priors for the model parameters.
204206

@@ -234,7 +236,7 @@ Let us therefore choose the HKY model with 4 gamma rate categories for the subst
234236
235237
<figure>
236238
<a id="substitutionModel"></a>
237-
<img style="width:75.0%;" src="figures/beast2_substitution_model.png" alt="">
239+
<img src="figures/beast2_substitution_model.png" alt="">
238240
<figcaption>Figure 6: Specifying the substitution model.</figcaption>
239241
</figure>
240242
<br>
@@ -262,7 +264,7 @@ Since we are observing the sequence data from a single epidemic of H3N2 virus in
262264

263265
<figure>
264266
<a id="clockModel"></a>
265-
<img style="width:75.0%;" src="figures/beast2_clock_model.png" alt="">
267+
<img src="figures/beast2_clock_model.png" alt="">
266268
<figcaption>Figure 7: Specifying the clock model.</figcaption>
267269
</figure>
268270
<br>
@@ -453,7 +455,7 @@ It is rarely necessary to specify a strong prior for equilibrium frequencies. Eq
453455
<figure>
454456
<a id="freqPrior"></a>
455457
<img src="figures/beast2_prior_frequencies.png" alt="">
456-
<figcaption>Figure 16: Specifying the equilibrium nucleotide frequencies prior.</figcaption>
458+
<figcaption>Figure 16: Specifying the equilibrium nucleotide frequencies prior. As it is a multivariate distribution, BEAUti cannot plot it.</figcaption>
457459
</figure>
458460
<br>
459461

figures/beast2_homochronous_tMRCA.png

11.3 KB
Loading

figures/beast2_package_manager.png

18.6 KB
Loading

figures/beast2_prior_clockRate.png

26.3 KB
Loading

figures/beast2_prior_rho.png

27.1 KB
Loading
16 KB
Loading

0 commit comments

Comments
 (0)