diff --git a/articles/biokinetics.html b/articles/biokinetics.html
index 7992cc2..242e193 100644
--- a/articles/biokinetics.html
+++ b/articles/biokinetics.html
@@ -120,9 +120,9 @@ <h2 id="figure-2">Figure 2<a class="anchor" aria-label="anchor" href="#figure-2"
 population trajectories using the fitted population parameters.</p>
 <div class="sourceCode" id="cb2"><pre class="downlit sourceCode r">
 <code class="sourceCode R"><span><span class="va">res</span> <span class="op">&lt;-</span> <span class="va">mod</span><span class="op">$</span><span class="fu">simulate_population_trajectories</span><span class="op">(</span><span class="op">)</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:19] Summarising fits</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:22] Adjusting by regression coefficients</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:25] Summarising into quantiles</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:26] Summarising fits</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:29] Adjusting by regression coefficients</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:31] Summarising into quantiles</span></span>
 <span><span class="fu"><a href="https://rdrr.io/r/utils/head.html" class="external-link">head</a></span><span class="op">(</span><span class="va">res</span><span class="op">)</span></span>
 <span><span class="co">#&gt;    time_since_last_exp       me        lo       hi titre_type infection_history</span></span>
 <span><span class="co">#&gt;                  &lt;int&gt;    &lt;num&gt;     &lt;num&gt;    &lt;num&gt;     &lt;char&gt;            &lt;char&gt;</span></span>
@@ -191,11 +191,11 @@ <h2 id="figure-4">Figure 4<a class="anchor" aria-label="anchor" href="#figure-4"
 simulate estimates for these from the fitted model:</p>
 <div class="sourceCode" id="cb4"><pre class="downlit sourceCode r">
 <code class="sourceCode R"><span><span class="va">res</span> <span class="op">&lt;-</span> <span class="va">mod</span><span class="op">$</span><span class="fu">population_stationary_points</span><span class="op">(</span>n_draws <span class="op">=</span> <span class="fl">2000</span><span class="op">)</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:26] Extracting parameters</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:27] Adjusting by covariates</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:27] Calculating peak and switch titre values</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:27] Recovering covariate names</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:27] Calculating medians</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:32] Extracting parameters</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:33] Adjusting by covariates</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:33] Calculating peak and switch titre values</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:33] Recovering covariate names</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:33] Calculating medians</span></span>
 <span><span class="fu"><a href="https://rdrr.io/r/utils/head.html" class="external-link">head</a></span><span class="op">(</span><span class="va">res</span><span class="op">)</span></span>
 <span><span class="co">#&gt;    infection_history titre_type     mu_p     mu_s rel_drop_me  mu_p_me  mu_s_me</span></span>
 <span><span class="co">#&gt;               &lt;char&gt;     &lt;char&gt;    &lt;num&gt;    &lt;num&gt;       &lt;num&gt;    &lt;num&gt;    &lt;num&gt;</span></span>
@@ -261,15 +261,15 @@ <h2 id="figure-5">Figure 5<a class="anchor" aria-label="anchor" href="#figure-5"
 maximum number of draws to just 250 for speed.</p>
 <div class="sourceCode" id="cb6"><pre class="downlit sourceCode r">
 <code class="sourceCode R"><span><span class="va">res</span> <span class="op">&lt;-</span> <span class="va">mod</span><span class="op">$</span><span class="fu">simulate_individual_trajectories</span><span class="op">(</span>n_draws <span class="op">=</span> <span class="fl">250</span><span class="op">)</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:28] Extracting parameters</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:31] Simulating individual trajectories</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:36] Recovering covariate names</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:36] Calculating exposure days. Adjusting exposures by 0 days</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:39] Resampling</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:34] Extracting parameters</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:37] Simulating individual trajectories</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:41] Recovering covariate names</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:41] Calculating exposure days. Adjusting exposures by 0 days</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:45] Resampling</span></span>
 <span><span class="co">#&gt; Registered S3 method overwritten by 'mosaic':</span></span>
 <span><span class="co">#&gt;   method                           from   </span></span>
 <span><span class="co">#&gt;   fortify.SpatialPolygonsDataFrame ggplot2</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:49] Summarising into population quantiles</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:54] Summarising into population quantiles</span></span>
 <span><span class="fu"><a href="https://rdrr.io/r/utils/head.html" class="external-link">head</a></span><span class="op">(</span><span class="va">res</span><span class="op">)</span></span>
 <span><span class="co">#&gt;    calendar_day titre_type       me       lo       hi time_shift</span></span>
 <span><span class="co">#&gt;          &lt;IDat&gt;     &lt;char&gt;    &lt;num&gt;    &lt;num&gt;    &lt;num&gt;      &lt;num&gt;</span></span>
@@ -348,72 +348,72 @@ <h2 id="figure-5">Figure 5<a class="anchor" aria-label="anchor" href="#figure-5"
 <span>  <span class="va">shift</span> <span class="op">&lt;-</span> <span class="va">time_shift_values</span><span class="op">[</span><span class="va">index</span><span class="op">]</span></span>
 <span>  <span class="va">mod</span><span class="op">$</span><span class="fu">simulate_individual_trajectories</span><span class="op">(</span>n_draws <span class="op">=</span> <span class="fl">50</span>, time_shift <span class="op">=</span> <span class="va">shift</span><span class="op">)</span></span>
 <span><span class="op">}</span><span class="op">)</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:50] Extracting parameters</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:53] Simulating individual trajectories</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:54] Recovering covariate names</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:54] Calculating exposure days. Adjusting exposures by -75 days</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:55] Resampling</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:57] Summarising into population quantiles</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 16:59:57] Extracting parameters</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:01] Simulating individual trajectories</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:01] Recovering covariate names</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:01] Calculating exposure days. Adjusting exposures by -60 days</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:02] Resampling</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:04] Summarising into population quantiles</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:04] Extracting parameters</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:08] Simulating individual trajectories</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:08] Recovering covariate names</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:08] Calculating exposure days. Adjusting exposures by -45 days</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:09] Resampling</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:10] Summarising into population quantiles</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:11] Extracting parameters</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:14] Simulating individual trajectories</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:15] Recovering covariate names</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:15] Calculating exposure days. Adjusting exposures by -30 days</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:15] Resampling</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:17] Summarising into population quantiles</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:18] Extracting parameters</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:21] Simulating individual trajectories</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:22] Recovering covariate names</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:22] Calculating exposure days. Adjusting exposures by -15 days</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:22] Resampling</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:24] Summarising into population quantiles</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:25] Extracting parameters</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:28] Simulating individual trajectories</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:29] Recovering covariate names</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:29] Calculating exposure days. Adjusting exposures by 0 days</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:29] Resampling</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:31] Summarising into population quantiles</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:32] Extracting parameters</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:35] Simulating individual trajectories</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:36] Recovering covariate names</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:36] Calculating exposure days. Adjusting exposures by 15 days</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:37] Resampling</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:39] Summarising into population quantiles</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:39] Extracting parameters</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:43] Simulating individual trajectories</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:43] Recovering covariate names</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:43] Calculating exposure days. Adjusting exposures by 30 days</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:44] Resampling</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:46] Summarising into population quantiles</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:46] Extracting parameters</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:50] Simulating individual trajectories</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:50] Recovering covariate names</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:50] Calculating exposure days. Adjusting exposures by 45 days</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:51] Resampling</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:53] Summarising into population quantiles</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:53] Extracting parameters</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:57] Simulating individual trajectories</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:57] Recovering covariate names</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:57] Calculating exposure days. Adjusting exposures by 60 days</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:58] Resampling</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:00:59] Summarising into population quantiles</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:01:00] Extracting parameters</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:01:03] Simulating individual trajectories</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:01:04] Recovering covariate names</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:01:04] Calculating exposure days. Adjusting exposures by 75 days</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:01:04] Resampling</span></span>
-<span><span class="co">#&gt; INFO [2024-10-24 17:01:06] Summarising into population quantiles</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:56] Extracting parameters</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:38:59] Simulating individual trajectories</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:00] Recovering covariate names</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:00] Calculating exposure days. Adjusting exposures by -75 days</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:00] Resampling</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:02] Summarising into population quantiles</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:03] Extracting parameters</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:06] Simulating individual trajectories</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:06] Recovering covariate names</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:06] Calculating exposure days. Adjusting exposures by -60 days</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:07] Resampling</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:09] Summarising into population quantiles</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:09] Extracting parameters</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:12] Simulating individual trajectories</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:13] Recovering covariate names</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:13] Calculating exposure days. Adjusting exposures by -45 days</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:13] Resampling</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:15] Summarising into population quantiles</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:16] Extracting parameters</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:19] Simulating individual trajectories</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:20] Recovering covariate names</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:20] Calculating exposure days. Adjusting exposures by -30 days</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:20] Resampling</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:22] Summarising into population quantiles</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:23] Extracting parameters</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:26] Simulating individual trajectories</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:26] Recovering covariate names</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:26] Calculating exposure days. Adjusting exposures by -15 days</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:27] Resampling</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:29] Summarising into population quantiles</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:29] Extracting parameters</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:33] Simulating individual trajectories</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:33] Recovering covariate names</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:33] Calculating exposure days. Adjusting exposures by 0 days</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:34] Resampling</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:36] Summarising into population quantiles</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:36] Extracting parameters</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:39] Simulating individual trajectories</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:40] Recovering covariate names</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:40] Calculating exposure days. Adjusting exposures by 15 days</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:40] Resampling</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:42] Summarising into population quantiles</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:43] Extracting parameters</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:46] Simulating individual trajectories</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:46] Recovering covariate names</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:46] Calculating exposure days. Adjusting exposures by 30 days</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:47] Resampling</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:49] Summarising into population quantiles</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:49] Extracting parameters</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:52] Simulating individual trajectories</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:53] Recovering covariate names</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:53] Calculating exposure days. Adjusting exposures by 45 days</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:53] Resampling</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:55] Summarising into population quantiles</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:56] Extracting parameters</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:39:59] Simulating individual trajectories</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:40:00] Recovering covariate names</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:40:00] Calculating exposure days. Adjusting exposures by 60 days</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:40:00] Resampling</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:40:02] Summarising into population quantiles</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:40:02] Extracting parameters</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:40:06] Simulating individual trajectories</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:40:06] Recovering covariate names</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:40:06] Calculating exposure days. Adjusting exposures by 75 days</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:40:07] Resampling</span></span>
+<span><span class="co">#&gt; INFO [2024-11-04 17:40:09] Summarising into population quantiles</span></span>
 <span></span>
 <span><span class="va">combined_data</span> <span class="op">&lt;-</span> <span class="fu">data.table</span><span class="fu">::</span><span class="fu"><a href="https://rdatatable.gitlab.io/data.table/reference/data.table.html" class="external-link">data.table</a></span><span class="op">(</span><span class="fu">data.table</span><span class="fu">::</span><span class="fu"><a href="https://rdatatable.gitlab.io/data.table/reference/rbindlist.html" class="external-link">rbindlist</a></span><span class="op">(</span><span class="va">results_list</span><span class="op">)</span><span class="op">)</span></span></code></pre></div>
 <p>Plotting the median values:</p>
diff --git a/index.html b/index.html
index 54d4988..178afed 100644
--- a/index.html
+++ b/index.html
@@ -12,7 +12,7 @@
 <link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/v4-shims.min.css" integrity="sha256-wZjR52fzng1pJHwx4aV2AO3yyTOXrcDW7jBpJtTwVxw=" crossorigin="anonymous">
 <!-- clipboard.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/clipboard.js/2.0.6/clipboard.min.js" integrity="sha256-inc5kl9MA1hkeYUt+EC3BhlIgyp/2jDIyBLS6k3UxPI=" crossorigin="anonymous"></script><!-- headroom.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/headroom.min.js" integrity="sha256-AsUX4SJE1+yuDu5+mAVzJbuYNPHj/WroHuZ8Ir/CkE0=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/jQuery.headroom.min.js" integrity="sha256-ZX/yNShbjqsohH1k95liqY9Gd8uOiE1S4vZc+9KQ1K4=" crossorigin="anonymous"></script><!-- pkgdown --><link href="pkgdown.css" rel="stylesheet">
 <script src="pkgdown.js"></script><meta property="og:title" content="Biomarker Kinetics Modelling">
-<meta property="og:description" content="Fit kinetic curves to biomarker data, using a Bayesian hierarchical model">
+<meta property="og:description" content="Fit kinetic curves to biomarker data, using a Bayesian hierarchical model.">
 <!-- mathjax --><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/MathJax.js" integrity="sha256-nvJJv9wWKEm88qvoQl9ekL2J+k/RWIsaSScxxlsrv8k=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/config/TeX-AMS-MML_HTMLorMML.js" integrity="sha256-84DKXVJXs0/F8OTMzX4UR909+jtl4G7SPypPavF+GfA=" crossorigin="anonymous"></script><!--[if lt IE 9]>
 <script src="https://oss.maxcdn.com/html5shiv/3.7.3/html5shiv.min.js"></script>
 <script src="https://oss.maxcdn.com/respond/1.4.2/respond.min.js"></script>
diff --git a/pkgdown.yml b/pkgdown.yml
index 19465e5..c2b05f9 100644
--- a/pkgdown.yml
+++ b/pkgdown.yml
@@ -5,4 +5,4 @@ articles:
   biokinetics: biokinetics.html
   data: data.html
   model: model.html
-last_built: 2024-10-24T16:55Z
+last_built: 2024-11-04T17:34Z
diff --git a/reference/biokinetics.html b/reference/biokinetics.html
index 02c42f1..456b7c4 100644
--- a/reference/biokinetics.html
+++ b/reference/biokinetics.html
@@ -69,6 +69,8 @@ <h2>Methods</h2>
 <h3 id="public-methods">Public methods<a class="anchor" aria-label="anchor" href="#public-methods"></a></h3>
 
 <ul><li><p><a href="#method-biokinetics-new"><code>biokinetics$new()</code></a></p></li>
+<li><p><a href="#method-biokinetics-plot_prior_predictive"><code>biokinetics$plot_prior_predictive()</code></a></p></li>
+<li><p><a href="#method-biokinetics-plot_model_inputs"><code>biokinetics$plot_model_inputs()</code></a></p></li>
 <li><p><a href="#method-biokinetics-get_stan_data"><code>biokinetics$get_stan_data()</code></a></p></li>
 <li><p><a href="#method-biokinetics-get_covariate_lookup_table"><code>biokinetics$get_covariate_lookup_table()</code></a></p></li>
 <li><p><a href="#method-biokinetics-fit"><code>biokinetics$fit()</code></a></p></li>
@@ -127,39 +129,79 @@ <h4 id="returns">Returns<a class="anchor" aria-label="anchor" href="#returns"></
 <p>An epikinetics::biokinetics object.</p>
 </div>
 
+</div><p></p><hr><a id="method-biokinetics-plot_prior_predictive"></a><div class="section">
+<h3 id="method-plot-prior-predictive-">Method <code>plot_prior_predictive()</code><a class="anchor" aria-label="anchor" href="#method-plot-prior-predictive-"></a></h3>
+<p>Plot the kinetics trajectory predicted by the model priors.
+Note that this is on a log scale, regardless of whether the data was provided
+on a log or a natural scale.</p><div class="section">
+<h4 id="usage-1">Usage<a class="anchor" aria-label="anchor" href="#usage-1"></a></h4>
+<p></p><div class="r"><div class="sourceCode"><pre><code><span><span class="va">biokinetics</span><span class="op">$</span><span class="fu">plot_prior_predictive</span><span class="op">(</span>tmax <span class="op">=</span> <span class="fl">150</span>, n_draws <span class="op">=</span> <span class="fl">2000</span><span class="op">)</span></span></code></pre></div><p></p></div>
+</div>
+
+<div class="section">
+<h4 id="arguments-1">Arguments<a class="anchor" aria-label="anchor" href="#arguments-1"></a></h4>
+<p></p><div class="arguments"><dl><dt><code>tmax</code></dt>
+<dd><p>Integer. The number of time points in each simulated trajectory. Default 150.</p></dd>
+
+
+<dt><code>n_draws</code></dt>
+<dd><p>Integer. The number of trajectories to simulate. Default 2000.</p></dd>
+
+
+</dl><p></p></div>
+</div>
+<div class="section">
+<h4 id="returns-1">Returns<a class="anchor" aria-label="anchor" href="#returns-1"></a></h4>
+<p>A ggplot2 object.</p>
+</div>
+
+</div><p></p><hr><a id="method-biokinetics-plot_model_inputs"></a><div class="section">
+<h3 id="method-plot-model-inputs-">Method <code>plot_model_inputs()</code><a class="anchor" aria-label="anchor" href="#method-plot-model-inputs-"></a></h3>
+<p>Plot model input data with a smoothing function. Note that
+this plot is on a log scale, regardless of whether data was provided on a
+log or a natural scale.</p><div class="section">
+<h4 id="usage-2">Usage<a class="anchor" aria-label="anchor" href="#usage-2"></a></h4>
+<p></p><div class="r"><div class="sourceCode"><pre><code><span><span class="va">biokinetics</span><span class="op">$</span><span class="fu">plot_model_inputs</span><span class="op">(</span><span class="op">)</span></span></code></pre></div><p></p></div>
+</div>
+
+<div class="section">
+<h4 id="returns-2">Returns<a class="anchor" aria-label="anchor" href="#returns-2"></a></h4>
+<p>A ggplot2 object.</p>
+</div>
+
 </div><p></p><hr><a id="method-biokinetics-get_stan_data"></a><div class="section">
 <h3 id="method-get-stan-data-">Method <code>get_stan_data()</code><a class="anchor" aria-label="anchor" href="#method-get-stan-data-"></a></h3>
 <p>View the data that is passed to the stan model, for debugging purposes.</p><div class="section">
-<h4 id="usage-1">Usage<a class="anchor" aria-label="anchor" href="#usage-1"></a></h4>
+<h4 id="usage-3">Usage<a class="anchor" aria-label="anchor" href="#usage-3"></a></h4>
 <p></p><div class="r"><div class="sourceCode"><pre><code><span><span class="va">biokinetics</span><span class="op">$</span><span class="fu">get_stan_data</span><span class="op">(</span><span class="op">)</span></span></code></pre></div><p></p></div>
 </div>
 
 <div class="section">
-<h4 id="returns-1">Returns<a class="anchor" aria-label="anchor" href="#returns-1"></a></h4>
+<h4 id="returns-3">Returns<a class="anchor" aria-label="anchor" href="#returns-3"></a></h4>
 <p>A list of arguments that will be passed to the stan model.</p>
 </div>
 
 </div><p></p><hr><a id="method-biokinetics-get_covariate_lookup_table"></a><div class="section">
 <h3 id="method-get-covariate-lookup-table-">Method <code>get_covariate_lookup_table()</code><a class="anchor" aria-label="anchor" href="#method-get-covariate-lookup-table-"></a></h3>
 <p>View the mapping of human readable covariate names to the model variable p.</p><div class="section">
-<h4 id="usage-2">Usage<a class="anchor" aria-label="anchor" href="#usage-2"></a></h4>
+<h4 id="usage-4">Usage<a class="anchor" aria-label="anchor" href="#usage-4"></a></h4>
 <p></p><div class="r"><div class="sourceCode"><pre><code><span><span class="va">biokinetics</span><span class="op">$</span><span class="fu">get_covariate_lookup_table</span><span class="op">(</span><span class="op">)</span></span></code></pre></div><p></p></div>
 </div>
 
 <div class="section">
-<h4 id="returns-2">Returns<a class="anchor" aria-label="anchor" href="#returns-2"></a></h4>
+<h4 id="returns-4">Returns<a class="anchor" aria-label="anchor" href="#returns-4"></a></h4>
 <p>A data.table mapping the model variable p to human readable covariates.</p>
 </div>
 
 </div><p></p><hr><a id="method-biokinetics-fit"></a><div class="section">
 <h3 id="method-fit-">Method <code>fit()</code><a class="anchor" aria-label="anchor" href="#method-fit-"></a></h3>
 <p>Fit the model and return CmdStanMCMC fitted model object.</p><div class="section">
-<h4 id="usage-3">Usage<a class="anchor" aria-label="anchor" href="#usage-3"></a></h4>
+<h4 id="usage-5">Usage<a class="anchor" aria-label="anchor" href="#usage-5"></a></h4>
 <p></p><div class="r"><div class="sourceCode"><pre><code><span><span class="va">biokinetics</span><span class="op">$</span><span class="fu">fit</span><span class="op">(</span><span class="va">...</span><span class="op">)</span></span></code></pre></div><p></p></div>
 </div>
 
 <div class="section">
-<h4 id="arguments-1">Arguments<a class="anchor" aria-label="anchor" href="#arguments-1"></a></h4>
+<h4 id="arguments-2">Arguments<a class="anchor" aria-label="anchor" href="#arguments-2"></a></h4>
 <p></p><div class="arguments"><dl><dt><code>...</code></dt>
 <dd><p>Named arguments to the <code><a href="https://rdrr.io/r/base/sample.html" class="external-link">sample()</a></code> method of CmdStan model.
 objects: <a href="https://mc-stan.org/cmdstanr/reference/model-method-sample.html" class="external-link">https://mc-stan.org/cmdstanr/reference/model-method-sample.html</a></p></dd>
@@ -168,24 +210,24 @@ <h4 id="arguments-1">Arguments<a class="anchor" aria-label="anchor" href="#argum
 </dl><p></p></div>
 </div>
 <div class="section">
-<h4 id="returns-3">Returns<a class="anchor" aria-label="anchor" href="#returns-3"></a></h4>
+<h4 id="returns-5">Returns<a class="anchor" aria-label="anchor" href="#returns-5"></a></h4>
 <p>A CmdStanMCMC fitted model object: <a href="https://mc-stan.org/cmdstanr/reference/CmdStanMCMC.html" class="external-link">https://mc-stan.org/cmdstanr/reference/CmdStanMCMC.html</a></p>
 </div>
 
 </div><p></p><hr><a id="method-biokinetics-extract_population_parameters"></a><div class="section">
 <h3 id="method-extract-population-parameters-">Method <code>extract_population_parameters()</code><a class="anchor" aria-label="anchor" href="#method-extract-population-parameters-"></a></h3>
 <p>Extract fitted population parameters</p><div class="section">
-<h4 id="usage-4">Usage<a class="anchor" aria-label="anchor" href="#usage-4"></a></h4>
+<h4 id="usage-6">Usage<a class="anchor" aria-label="anchor" href="#usage-6"></a></h4>
 <p></p><div class="r"><div class="sourceCode"><pre><code><span><span class="va">biokinetics</span><span class="op">$</span><span class="fu">extract_population_parameters</span><span class="op">(</span></span>
-<span>  n_draws <span class="op">=</span> <span class="fl">2500</span>,</span>
+<span>  n_draws <span class="op">=</span> <span class="fl">2000</span>,</span>
 <span>  human_readable_covariates <span class="op">=</span> <span class="cn">TRUE</span></span>
 <span><span class="op">)</span></span></code></pre></div><p></p></div>
 </div>
 
 <div class="section">
-<h4 id="arguments-2">Arguments<a class="anchor" aria-label="anchor" href="#arguments-2"></a></h4>
+<h4 id="arguments-3">Arguments<a class="anchor" aria-label="anchor" href="#arguments-3"></a></h4>
 <p></p><div class="arguments"><dl><dt><code>n_draws</code></dt>
-<dd><p>Numeric</p></dd>
+<dd><p>Integer. Default 2000.</p></dd>
 
 
 <dt><code>human_readable_covariates</code></dt>
@@ -195,25 +237,25 @@ <h4 id="arguments-2">Arguments<a class="anchor" aria-label="anchor" href="#argum
 </dl><p></p></div>
 </div>
 <div class="section">
-<h4 id="returns-4">Returns<a class="anchor" aria-label="anchor" href="#returns-4"></a></h4>
+<h4 id="returns-6">Returns<a class="anchor" aria-label="anchor" href="#returns-6"></a></h4>
 <p>A data.table</p>
 </div>
 
 </div><p></p><hr><a id="method-biokinetics-extract_individual_parameters"></a><div class="section">
 <h3 id="method-extract-individual-parameters-">Method <code>extract_individual_parameters()</code><a class="anchor" aria-label="anchor" href="#method-extract-individual-parameters-"></a></h3>
 <p>Extract fitted individual parameters</p><div class="section">
-<h4 id="usage-5">Usage<a class="anchor" aria-label="anchor" href="#usage-5"></a></h4>
+<h4 id="usage-7">Usage<a class="anchor" aria-label="anchor" href="#usage-7"></a></h4>
 <p></p><div class="r"><div class="sourceCode"><pre><code><span><span class="va">biokinetics</span><span class="op">$</span><span class="fu">extract_individual_parameters</span><span class="op">(</span></span>
-<span>  n_draws <span class="op">=</span> <span class="fl">2500</span>,</span>
+<span>  n_draws <span class="op">=</span> <span class="fl">2000</span>,</span>
 <span>  include_variation_params <span class="op">=</span> <span class="cn">TRUE</span>,</span>
 <span>  human_readable_covariates <span class="op">=</span> <span class="cn">TRUE</span></span>
 <span><span class="op">)</span></span></code></pre></div><p></p></div>
 </div>
 
 <div class="section">
-<h4 id="arguments-3">Arguments<a class="anchor" aria-label="anchor" href="#arguments-3"></a></h4>
+<h4 id="arguments-4">Arguments<a class="anchor" aria-label="anchor" href="#arguments-4"></a></h4>
 <p></p><div class="arguments"><dl><dt><code>n_draws</code></dt>
-<dd><p>Numeric</p></dd>
+<dd><p>Integer. Default 2000.</p></dd>
 
 
 <dt><code>include_variation_params</code></dt>
@@ -227,23 +269,23 @@ <h4 id="arguments-3">Arguments<a class="anchor" aria-label="anchor" href="#argum
 </dl><p></p></div>
 </div>
 <div class="section">
-<h4 id="returns-5">Returns<a class="anchor" aria-label="anchor" href="#returns-5"></a></h4>
+<h4 id="returns-7">Returns<a class="anchor" aria-label="anchor" href="#returns-7"></a></h4>
 <p>A data.table</p>
 </div>
 
 </div><p></p><hr><a id="method-biokinetics-simulate_population_trajectories"></a><div class="section">
 <h3 id="method-simulate-population-trajectories-">Method <code>simulate_population_trajectories()</code><a class="anchor" aria-label="anchor" href="#method-simulate-population-trajectories-"></a></h3>
 <p>Process the model results into a data table of titre values over time.</p><div class="section">
-<h4 id="usage-6">Usage<a class="anchor" aria-label="anchor" href="#usage-6"></a></h4>
+<h4 id="usage-8">Usage<a class="anchor" aria-label="anchor" href="#usage-8"></a></h4>
 <p></p><div class="r"><div class="sourceCode"><pre><code><span><span class="va">biokinetics</span><span class="op">$</span><span class="fu">simulate_population_trajectories</span><span class="op">(</span></span>
 <span>  t_max <span class="op">=</span> <span class="fl">150</span>,</span>
 <span>  summarise <span class="op">=</span> <span class="cn">TRUE</span>,</span>
-<span>  n_draws <span class="op">=</span> <span class="fl">2500</span></span>
+<span>  n_draws <span class="op">=</span> <span class="fl">2000</span></span>
 <span><span class="op">)</span></span></code></pre></div><p></p></div>
 </div>
 
 <div class="section">
-<h4 id="arguments-4">Arguments<a class="anchor" aria-label="anchor" href="#arguments-4"></a></h4>
+<h4 id="arguments-5">Arguments<a class="anchor" aria-label="anchor" href="#arguments-5"></a></h4>
 <p></p><div class="arguments"><dl><dt><code>t_max</code></dt>
 <dd><p>Integer. Maximum number of time points to include.</p></dd>
 
@@ -255,13 +297,13 @@ <h4 id="arguments-4">Arguments<a class="anchor" aria-label="anchor" href="#argum
 
 
 <dt><code>n_draws</code></dt>
-<dd><p>Integer. Maximum number of samples to include. Default 2500.</p></dd>
+<dd><p>Integer. Maximum number of samples to include. Default 2000.</p></dd>
 
 
 </dl><p></p></div>
 </div>
 <div class="section">
-<h4 id="returns-6">Returns<a class="anchor" aria-label="anchor" href="#returns-6"></a></h4>
+<h4 id="returns-8">Returns<a class="anchor" aria-label="anchor" href="#returns-8"></a></h4>
 <p>A data.table containing titre values at time points. If summarise = TRUE, columns are time_since_last_exp,
 me, lo, hi, titre_type, and a column for each covariate in the hierarchical model. If summarise = FALSE, columns are
 time_since_last_exp, .draw, t0_pop, tp_pop, ts_pop, m1_pop, m2_pop, m3_pop, beta_t0, beta_tp, beta_ts, beta_m1, beta_m2, beta_m3, mu
@@ -272,20 +314,20 @@ <h4 id="returns-6">Returns<a class="anchor" aria-label="anchor" href="#returns-6
 </div><p></p><hr><a id="method-biokinetics-population_stationary_points"></a><div class="section">
 <h3 id="method-population-stationary-points-">Method <code>population_stationary_points()</code><a class="anchor" aria-label="anchor" href="#method-population-stationary-points-"></a></h3>
 <p>Process the stan model results into a data.table.</p><div class="section">
-<h4 id="usage-7">Usage<a class="anchor" aria-label="anchor" href="#usage-7"></a></h4>
-<p></p><div class="r"><div class="sourceCode"><pre><code><span><span class="va">biokinetics</span><span class="op">$</span><span class="fu">population_stationary_points</span><span class="op">(</span>n_draws <span class="op">=</span> <span class="fl">2500</span><span class="op">)</span></span></code></pre></div><p></p></div>
+<h4 id="usage-9">Usage<a class="anchor" aria-label="anchor" href="#usage-9"></a></h4>
+<p></p><div class="r"><div class="sourceCode"><pre><code><span><span class="va">biokinetics</span><span class="op">$</span><span class="fu">population_stationary_points</span><span class="op">(</span>n_draws <span class="op">=</span> <span class="fl">2000</span><span class="op">)</span></span></code></pre></div><p></p></div>
 </div>
 
 <div class="section">
-<h4 id="arguments-5">Arguments<a class="anchor" aria-label="anchor" href="#arguments-5"></a></h4>
+<h4 id="arguments-6">Arguments<a class="anchor" aria-label="anchor" href="#arguments-6"></a></h4>
 <p></p><div class="arguments"><dl><dt><code>n_draws</code></dt>
-<dd><p>Integer. Maximum number of samples to include. Default 2500.</p></dd>
+<dd><p>Integer. Maximum number of samples to include. Default 2000.</p></dd>
 
 
 </dl><p></p></div>
 </div>
 <div class="section">
-<h4 id="returns-7">Returns<a class="anchor" aria-label="anchor" href="#returns-7"></a></h4>
+<h4 id="returns-9">Returns<a class="anchor" aria-label="anchor" href="#returns-9"></a></h4>
 <p>A data.table of peak and set titre values. Columns are tire_type, mu_p, mu_s, rel_drop_me, mu_p_me,
 mu_s_me, and a column for each covariate. See the data vignette for details:
 <code><a href="../articles/data.html">vignette("data", package = "epikinetics")</a></code></p>
@@ -295,16 +337,16 @@ <h4 id="returns-7">Returns<a class="anchor" aria-label="anchor" href="#returns-7
 <h3 id="method-simulate-individual-trajectories-">Method <code>simulate_individual_trajectories()</code><a class="anchor" aria-label="anchor" href="#method-simulate-individual-trajectories-"></a></h3>
 <p>Simulate individual trajectories from the model. This is
 computationally expensive and may take a while to run if n_draws is large.</p><div class="section">
-<h4 id="usage-8">Usage<a class="anchor" aria-label="anchor" href="#usage-8"></a></h4>
+<h4 id="usage-10">Usage<a class="anchor" aria-label="anchor" href="#usage-10"></a></h4>
 <p></p><div class="r"><div class="sourceCode"><pre><code><span><span class="va">biokinetics</span><span class="op">$</span><span class="fu">simulate_individual_trajectories</span><span class="op">(</span></span>
 <span>  summarise <span class="op">=</span> <span class="cn">TRUE</span>,</span>
-<span>  n_draws <span class="op">=</span> <span class="fl">2500</span>,</span>
+<span>  n_draws <span class="op">=</span> <span class="fl">2000</span>,</span>
 <span>  time_shift <span class="op">=</span> <span class="fl">0</span></span>
 <span><span class="op">)</span></span></code></pre></div><p></p></div>
 </div>
 
 <div class="section">
-<h4 id="arguments-6">Arguments<a class="anchor" aria-label="anchor" href="#arguments-6"></a></h4>
+<h4 id="arguments-7">Arguments<a class="anchor" aria-label="anchor" href="#arguments-7"></a></h4>
 <p></p><div class="arguments"><dl><dt><code>summarise</code></dt>
 <dd><p>Boolean. If TRUE, average the individual trajectories to get lo, me and
 hi values for the population, disaggregated by titre type. If FALSE return the indidivudal trajectories.
@@ -312,7 +354,7 @@ <h4 id="arguments-6">Arguments<a class="anchor" aria-label="anchor" href="#argum
 
 
 <dt><code>n_draws</code></dt>
-<dd><p>Integer. Maximum number of samples to draw. Default 2500.</p></dd>
+<dd><p>Integer. Maximum number of samples to draw. Default 2000.</p></dd>
 
 
 <dt><code>time_shift</code></dt>
@@ -322,7 +364,7 @@ <h4 id="arguments-6">Arguments<a class="anchor" aria-label="anchor" href="#argum
 </dl><p></p></div>
 </div>
 <div class="section">
-<h4 id="returns-8">Returns<a class="anchor" aria-label="anchor" href="#returns-8"></a></h4>
+<h4 id="returns-10">Returns<a class="anchor" aria-label="anchor" href="#returns-10"></a></h4>
 <p>A data.table. If summarise = TRUE columns are calendar_date, titre_type, me, lo, hi, time_shift.
 If summarise = FALSE, columns are pid, draw, time_since_last_exp, mu, titre_type, exposure_day, calendar_day, time_shift
 and a column for each covariate in the regression model. See the data vignette for details:
diff --git a/reference/epikinetics-package.html b/reference/epikinetics-package.html
index 91f2415..ba85cce 100644
--- a/reference/epikinetics-package.html
+++ b/reference/epikinetics-package.html
@@ -1,5 +1,5 @@
 <!DOCTYPE html>
-<!-- Generated by pkgdown: do not edit by hand --><html lang="en"><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8"><meta charset="utf-8"><meta http-equiv="X-UA-Compatible" content="IE=edge"><meta name="viewport" content="width=device-width, initial-scale=1.0"><title>epikinetics: Biomarker Kinetics Modelling — epikinetics-package • epikinetics</title><!-- jquery --><script src="https://cdnjs.cloudflare.com/ajax/libs/jquery/3.7.1/jquery.min.js" integrity="sha512-v2CJ7UaYy4JwqLDIrZUI/4hqeoQieOmAZNXBeQyjo21dadnwR+8ZaIJVT8EE2iyI61OV8e6M8PP2/4hpQINQ/g==" crossorigin="anonymous" referrerpolicy="no-referrer"></script><!-- Bootstrap --><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.4.1/css/bootstrap.min.css" integrity="sha256-bZLfwXAP04zRMK2BjiO8iu9pf4FbLqX6zitd+tIvLhE=" crossorigin="anonymous"><script src="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.4.1/js/bootstrap.min.js" integrity="sha256-nuL8/2cJ5NDSSwnKD8VqreErSWHtnEP9E7AySL+1ev4=" crossorigin="anonymous"></script><!-- bootstrap-toc --><link rel="stylesheet" href="../bootstrap-toc.css"><script src="../bootstrap-toc.js"></script><!-- Font Awesome icons --><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/all.min.css" integrity="sha256-mmgLkCYLUQbXn0B1SRqzHar6dCnv9oZFPEC1g1cwlkk=" crossorigin="anonymous"><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/v4-shims.min.css" integrity="sha256-wZjR52fzng1pJHwx4aV2AO3yyTOXrcDW7jBpJtTwVxw=" crossorigin="anonymous"><!-- clipboard.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/clipboard.js/2.0.6/clipboard.min.js" integrity="sha256-inc5kl9MA1hkeYUt+EC3BhlIgyp/2jDIyBLS6k3UxPI=" crossorigin="anonymous"></script><!-- headroom.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/headroom.min.js" integrity="sha256-AsUX4SJE1+yuDu5+mAVzJbuYNPHj/WroHuZ8Ir/CkE0=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/jQuery.headroom.min.js" integrity="sha256-ZX/yNShbjqsohH1k95liqY9Gd8uOiE1S4vZc+9KQ1K4=" crossorigin="anonymous"></script><!-- pkgdown --><link href="../pkgdown.css" rel="stylesheet"><script src="../pkgdown.js"></script><meta property="og:title" content="epikinetics: Biomarker Kinetics Modelling — epikinetics-package"><meta property="og:description" content="Fit kinetic curves to biomarker data, using a Bayesian hierarchical model"><!-- mathjax --><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/MathJax.js" integrity="sha256-nvJJv9wWKEm88qvoQl9ekL2J+k/RWIsaSScxxlsrv8k=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/config/TeX-AMS-MML_HTMLorMML.js" integrity="sha256-84DKXVJXs0/F8OTMzX4UR909+jtl4G7SPypPavF+GfA=" crossorigin="anonymous"></script><!--[if lt IE 9]>
+<!-- Generated by pkgdown: do not edit by hand --><html lang="en"><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8"><meta charset="utf-8"><meta http-equiv="X-UA-Compatible" content="IE=edge"><meta name="viewport" content="width=device-width, initial-scale=1.0"><title>epikinetics: Biomarker Kinetics Modelling — epikinetics-package • epikinetics</title><!-- jquery --><script src="https://cdnjs.cloudflare.com/ajax/libs/jquery/3.7.1/jquery.min.js" integrity="sha512-v2CJ7UaYy4JwqLDIrZUI/4hqeoQieOmAZNXBeQyjo21dadnwR+8ZaIJVT8EE2iyI61OV8e6M8PP2/4hpQINQ/g==" crossorigin="anonymous" referrerpolicy="no-referrer"></script><!-- Bootstrap --><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.4.1/css/bootstrap.min.css" integrity="sha256-bZLfwXAP04zRMK2BjiO8iu9pf4FbLqX6zitd+tIvLhE=" crossorigin="anonymous"><script src="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.4.1/js/bootstrap.min.js" integrity="sha256-nuL8/2cJ5NDSSwnKD8VqreErSWHtnEP9E7AySL+1ev4=" crossorigin="anonymous"></script><!-- bootstrap-toc --><link rel="stylesheet" href="../bootstrap-toc.css"><script src="../bootstrap-toc.js"></script><!-- Font Awesome icons --><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/all.min.css" integrity="sha256-mmgLkCYLUQbXn0B1SRqzHar6dCnv9oZFPEC1g1cwlkk=" crossorigin="anonymous"><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/v4-shims.min.css" integrity="sha256-wZjR52fzng1pJHwx4aV2AO3yyTOXrcDW7jBpJtTwVxw=" crossorigin="anonymous"><!-- clipboard.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/clipboard.js/2.0.6/clipboard.min.js" integrity="sha256-inc5kl9MA1hkeYUt+EC3BhlIgyp/2jDIyBLS6k3UxPI=" crossorigin="anonymous"></script><!-- headroom.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/headroom.min.js" integrity="sha256-AsUX4SJE1+yuDu5+mAVzJbuYNPHj/WroHuZ8Ir/CkE0=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/jQuery.headroom.min.js" integrity="sha256-ZX/yNShbjqsohH1k95liqY9Gd8uOiE1S4vZc+9KQ1K4=" crossorigin="anonymous"></script><!-- pkgdown --><link href="../pkgdown.css" rel="stylesheet"><script src="../pkgdown.js"></script><meta property="og:title" content="epikinetics: Biomarker Kinetics Modelling — epikinetics-package"><meta property="og:description" content="Fit kinetic curves to biomarker data, using a Bayesian hierarchical model."><!-- mathjax --><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/MathJax.js" integrity="sha256-nvJJv9wWKEm88qvoQl9ekL2J+k/RWIsaSScxxlsrv8k=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/config/TeX-AMS-MML_HTMLorMML.js" integrity="sha256-84DKXVJXs0/F8OTMzX4UR909+jtl4G7SPypPavF+GfA=" crossorigin="anonymous"></script><!--[if lt IE 9]>
 <script src="https://oss.maxcdn.com/html5shiv/3.7.3/html5shiv.min.js"></script>
 <script src="https://oss.maxcdn.com/respond/1.4.2/respond.min.js"></script>
 <![endif]--></head><body data-spy="scroll" data-target="#toc">
@@ -56,7 +56,7 @@ <h1>epikinetics: Biomarker Kinetics Modelling</h1>
     </div>
 
     <div class="ref-description">
-    <p>Fit kinetic curves to biomarker data, using a Bayesian hierarchical model</p>
+    <p>Fit kinetic curves to biomarker data, using a Bayesian hierarchical model.</p>
     </div>
 
 
diff --git a/reference/index.html b/reference/index.html
index 21fc473..f1fa73b 100644
--- a/reference/index.html
+++ b/reference/index.html
@@ -73,6 +73,14 @@ <h2 id="all-functions">All functions <a href="#all-functions" class="anchor" ari
           <p><code><a href="convert_log2_scale_inverse.html">convert_log2_scale_inverse()</a></code> </p>
         </td>
         <td><p>Invert base 2 log scale conversion</p></td>
+      </tr><tr><td>
+          <p><code><a href="plot.biokinetics_population_trajectories.html">plot(<i>&lt;biokinetics_population_trajectories&gt;</i>)</a></code> </p>
+        </td>
+        <td><p>Plot method for "biokinetics_population_trajectories" class</p></td>
+      </tr><tr><td>
+          <p><code><a href="plot.biokinetics_priors.html">plot(<i>&lt;biokinetics_priors&gt;</i>)</a></code> </p>
+        </td>
+        <td><p>Simulate biomarker kinetics predicted by the given biokinetics priors and optionally compare to a dataset.</p></td>
       </tr></tbody></table></div>
 
   <div class="col-md-3 hidden-xs hidden-sm" id="pkgdown-sidebar">
diff --git a/reference/plot.biokinetics_population_trajectories.html b/reference/plot.biokinetics_population_trajectories.html
new file mode 100644
index 0000000..0eecbbb
--- /dev/null
+++ b/reference/plot.biokinetics_population_trajectories.html
@@ -0,0 +1,109 @@
+<!DOCTYPE html>
+<!-- Generated by pkgdown: do not edit by hand --><html lang="en"><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8"><meta charset="utf-8"><meta http-equiv="X-UA-Compatible" content="IE=edge"><meta name="viewport" content="width=device-width, initial-scale=1.0"><title>Plot method for "biokinetics_population_trajectories" class — plot.biokinetics_population_trajectories • epikinetics</title><!-- jquery --><script src="https://cdnjs.cloudflare.com/ajax/libs/jquery/3.7.1/jquery.min.js" integrity="sha512-v2CJ7UaYy4JwqLDIrZUI/4hqeoQieOmAZNXBeQyjo21dadnwR+8ZaIJVT8EE2iyI61OV8e6M8PP2/4hpQINQ/g==" crossorigin="anonymous" referrerpolicy="no-referrer"></script><!-- Bootstrap --><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.4.1/css/bootstrap.min.css" integrity="sha256-bZLfwXAP04zRMK2BjiO8iu9pf4FbLqX6zitd+tIvLhE=" crossorigin="anonymous"><script src="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.4.1/js/bootstrap.min.js" integrity="sha256-nuL8/2cJ5NDSSwnKD8VqreErSWHtnEP9E7AySL+1ev4=" crossorigin="anonymous"></script><!-- bootstrap-toc --><link rel="stylesheet" href="../bootstrap-toc.css"><script src="../bootstrap-toc.js"></script><!-- Font Awesome icons --><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/all.min.css" integrity="sha256-mmgLkCYLUQbXn0B1SRqzHar6dCnv9oZFPEC1g1cwlkk=" crossorigin="anonymous"><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/v4-shims.min.css" integrity="sha256-wZjR52fzng1pJHwx4aV2AO3yyTOXrcDW7jBpJtTwVxw=" crossorigin="anonymous"><!-- clipboard.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/clipboard.js/2.0.6/clipboard.min.js" integrity="sha256-inc5kl9MA1hkeYUt+EC3BhlIgyp/2jDIyBLS6k3UxPI=" crossorigin="anonymous"></script><!-- headroom.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/headroom.min.js" integrity="sha256-AsUX4SJE1+yuDu5+mAVzJbuYNPHj/WroHuZ8Ir/CkE0=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/jQuery.headroom.min.js" integrity="sha256-ZX/yNShbjqsohH1k95liqY9Gd8uOiE1S4vZc+9KQ1K4=" crossorigin="anonymous"></script><!-- pkgdown --><link href="../pkgdown.css" rel="stylesheet"><script src="../pkgdown.js"></script><meta property="og:title" content="Plot method for " biokinetics_population_trajectories class plot.biokinetics_population_trajectories><meta property="og:description" content='Plot method for "biokinetics_population_trajectories" class'><!-- mathjax --><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/MathJax.js" integrity="sha256-nvJJv9wWKEm88qvoQl9ekL2J+k/RWIsaSScxxlsrv8k=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/config/TeX-AMS-MML_HTMLorMML.js" integrity="sha256-84DKXVJXs0/F8OTMzX4UR909+jtl4G7SPypPavF+GfA=" crossorigin="anonymous"></script><!--[if lt IE 9]>
+<script src="https://oss.maxcdn.com/html5shiv/3.7.3/html5shiv.min.js"></script>
+<script src="https://oss.maxcdn.com/respond/1.4.2/respond.min.js"></script>
+<![endif]--></head><body data-spy="scroll" data-target="#toc">
+
+
+    <div class="container template-reference-topic">
+      <header><div class="navbar navbar-default navbar-fixed-top" role="navigation">
+  <div class="container">
+    <div class="navbar-header">
+      <button type="button" class="navbar-toggle collapsed" data-toggle="collapse" data-target="#navbar" aria-expanded="false">
+        <span class="sr-only">Toggle navigation</span>
+        <span class="icon-bar"></span>
+        <span class="icon-bar"></span>
+        <span class="icon-bar"></span>
+      </button>
+      <span class="navbar-brand">
+        <a class="navbar-link" href="../index.html">epikinetics</a>
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="">0.0.0.9000</span>
+      </span>
+    </div>
+
+    <div id="navbar" class="navbar-collapse collapse">
+      <ul class="nav navbar-nav"><li>
+  <a href="../reference/index.html">Reference</a>
+</li>
+<li class="dropdown">
+  <a href="#" class="dropdown-toggle" data-toggle="dropdown" role="button" data-bs-toggle="dropdown" aria-expanded="false">
+    Articles
+
+    <span class="caret"></span>
+  </a>
+  <ul class="dropdown-menu" role="menu"><li>
+      <a href="../articles/biokinetics.html">biokinetics case study</a>
+    </li>
+    <li>
+      <a href="../articles/data.html">A guide to biokinetics input and output data</a>
+    </li>
+    <li>
+      <a href="../articles/model.html">Summary of the Stan Model</a>
+    </li>
+  </ul></li>
+      </ul><ul class="nav navbar-nav navbar-right"></ul></div><!--/.nav-collapse -->
+  </div><!--/.container -->
+</div><!--/.navbar -->
+
+
+
+      </header><div class="row">
+  <div class="col-md-9 contents">
+    <div class="page-header">
+    <h1>Plot method for "biokinetics_population_trajectories" class</h1>
+
+    <div class="hidden name"><code>plot.biokinetics_population_trajectories.Rd</code></div>
+    </div>
+
+    <div class="ref-description">
+    <p>Plot method for "biokinetics_population_trajectories" class</p>
+    </div>
+
+    <div id="ref-usage">
+    <div class="sourceCode"><pre class="sourceCode r"><code><span><span class="co"># S3 method for class 'biokinetics_population_trajectories'</span></span>
+<span><span class="fu"><a href="https://rdrr.io/r/graphics/plot.default.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">x</span>, <span class="va">...</span>, data <span class="op">=</span> <span class="cn">NULL</span><span class="op">)</span></span></code></pre></div>
+    </div>
+
+    <div id="arguments">
+    <h2>Arguments</h2>
+
+
+<dl><dt id="arg-x">x<a class="anchor" aria-label="anchor" href="#arg-x"></a></dt>
+<dd><p>An object of class "biokinetics_population_trajectories". These are
+generated by running biokinetics$simulate_population_trajectories(). See
+<a href="../../epikinetics/html/biokinetics.html#method-biokinetics-simulate_population_trajectories"><code>biokinetics$simulate_population_trajectories()</code></a></p></dd>
+
+
+<dt id="arg--">...<a class="anchor" aria-label="anchor" href="#arg--"></a></dt>
+<dd><p>Further arguments passed to the method.</p></dd>
+
+
+<dt id="arg-data">data<a class="anchor" aria-label="anchor" href="#arg-data"></a></dt>
+<dd><p>Optional data.table containing raw data as provided to the biokinetics model.</p></dd>
+
+</dl></div>
+
+  </div>
+  <div class="col-md-3 hidden-xs hidden-sm" id="pkgdown-sidebar">
+    <nav id="toc" data-toggle="toc" class="sticky-top"><h2 data-toc-skip>Contents</h2>
+    </nav></div>
+</div>
+
+
+      <footer><div class="copyright">
+  <p></p><p>Developed by Timothy Russell, Alex Hill.</p>
+</div>
+
+<div class="pkgdown">
+  <p></p><p>Site built with <a href="https://pkgdown.r-lib.org/" class="external-link">pkgdown</a> 2.1.1.</p>
+</div>
+
+      </footer></div>
+
+
+
+
+
+
+  </body></html>
+
diff --git a/reference/plot.biokinetics_priors.html b/reference/plot.biokinetics_priors.html
new file mode 100644
index 0000000..ea93582
--- /dev/null
+++ b/reference/plot.biokinetics_priors.html
@@ -0,0 +1,121 @@
+<!DOCTYPE html>
+<!-- Generated by pkgdown: do not edit by hand --><html lang="en"><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8"><meta charset="utf-8"><meta http-equiv="X-UA-Compatible" content="IE=edge"><meta name="viewport" content="width=device-width, initial-scale=1.0"><title>Simulate biomarker kinetics predicted by the given biokinetics priors and optionally compare to a dataset. — plot.biokinetics_priors • epikinetics</title><!-- jquery --><script src="https://cdnjs.cloudflare.com/ajax/libs/jquery/3.7.1/jquery.min.js" integrity="sha512-v2CJ7UaYy4JwqLDIrZUI/4hqeoQieOmAZNXBeQyjo21dadnwR+8ZaIJVT8EE2iyI61OV8e6M8PP2/4hpQINQ/g==" crossorigin="anonymous" referrerpolicy="no-referrer"></script><!-- Bootstrap --><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.4.1/css/bootstrap.min.css" integrity="sha256-bZLfwXAP04zRMK2BjiO8iu9pf4FbLqX6zitd+tIvLhE=" crossorigin="anonymous"><script src="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.4.1/js/bootstrap.min.js" integrity="sha256-nuL8/2cJ5NDSSwnKD8VqreErSWHtnEP9E7AySL+1ev4=" crossorigin="anonymous"></script><!-- bootstrap-toc --><link rel="stylesheet" href="../bootstrap-toc.css"><script src="../bootstrap-toc.js"></script><!-- Font Awesome icons --><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/all.min.css" integrity="sha256-mmgLkCYLUQbXn0B1SRqzHar6dCnv9oZFPEC1g1cwlkk=" crossorigin="anonymous"><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/v4-shims.min.css" integrity="sha256-wZjR52fzng1pJHwx4aV2AO3yyTOXrcDW7jBpJtTwVxw=" crossorigin="anonymous"><!-- clipboard.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/clipboard.js/2.0.6/clipboard.min.js" integrity="sha256-inc5kl9MA1hkeYUt+EC3BhlIgyp/2jDIyBLS6k3UxPI=" crossorigin="anonymous"></script><!-- headroom.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/headroom.min.js" integrity="sha256-AsUX4SJE1+yuDu5+mAVzJbuYNPHj/WroHuZ8Ir/CkE0=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/jQuery.headroom.min.js" integrity="sha256-ZX/yNShbjqsohH1k95liqY9Gd8uOiE1S4vZc+9KQ1K4=" crossorigin="anonymous"></script><!-- pkgdown --><link href="../pkgdown.css" rel="stylesheet"><script src="../pkgdown.js"></script><meta property="og:title" content="Simulate biomarker kinetics predicted by the given biokinetics priors and optionally compare to a dataset. — plot.biokinetics_priors"><meta property="og:description" content="Simulate trajectories by drawing random samples from the given
+priors for each parameter in the biokinetics model."><!-- mathjax --><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/MathJax.js" integrity="sha256-nvJJv9wWKEm88qvoQl9ekL2J+k/RWIsaSScxxlsrv8k=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/config/TeX-AMS-MML_HTMLorMML.js" integrity="sha256-84DKXVJXs0/F8OTMzX4UR909+jtl4G7SPypPavF+GfA=" crossorigin="anonymous"></script><!--[if lt IE 9]>
+<script src="https://oss.maxcdn.com/html5shiv/3.7.3/html5shiv.min.js"></script>
+<script src="https://oss.maxcdn.com/respond/1.4.2/respond.min.js"></script>
+<![endif]--></head><body data-spy="scroll" data-target="#toc">
+
+
+    <div class="container template-reference-topic">
+      <header><div class="navbar navbar-default navbar-fixed-top" role="navigation">
+  <div class="container">
+    <div class="navbar-header">
+      <button type="button" class="navbar-toggle collapsed" data-toggle="collapse" data-target="#navbar" aria-expanded="false">
+        <span class="sr-only">Toggle navigation</span>
+        <span class="icon-bar"></span>
+        <span class="icon-bar"></span>
+        <span class="icon-bar"></span>
+      </button>
+      <span class="navbar-brand">
+        <a class="navbar-link" href="../index.html">epikinetics</a>
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="">0.0.0.9000</span>
+      </span>
+    </div>
+
+    <div id="navbar" class="navbar-collapse collapse">
+      <ul class="nav navbar-nav"><li>
+  <a href="../reference/index.html">Reference</a>
+</li>
+<li class="dropdown">
+  <a href="#" class="dropdown-toggle" data-toggle="dropdown" role="button" data-bs-toggle="dropdown" aria-expanded="false">
+    Articles
+
+    <span class="caret"></span>
+  </a>
+  <ul class="dropdown-menu" role="menu"><li>
+      <a href="../articles/biokinetics.html">biokinetics case study</a>
+    </li>
+    <li>
+      <a href="../articles/data.html">A guide to biokinetics input and output data</a>
+    </li>
+    <li>
+      <a href="../articles/model.html">Summary of the Stan Model</a>
+    </li>
+  </ul></li>
+      </ul><ul class="nav navbar-nav navbar-right"></ul></div><!--/.nav-collapse -->
+  </div><!--/.container -->
+</div><!--/.navbar -->
+
+
+
+      </header><div class="row">
+  <div class="col-md-9 contents">
+    <div class="page-header">
+    <h1>Simulate biomarker kinetics predicted by the given biokinetics priors and optionally compare to a dataset.</h1>
+
+    <div class="hidden name"><code>plot.biokinetics_priors.Rd</code></div>
+    </div>
+
+    <div class="ref-description">
+    <p>Simulate trajectories by drawing random samples from the given
+priors for each parameter in the biokinetics model.</p>
+    </div>
+
+    <div id="ref-usage">
+    <div class="sourceCode"><pre class="sourceCode r"><code><span><span class="co"># S3 method for class 'biokinetics_priors'</span></span>
+<span><span class="fu"><a href="https://rdrr.io/r/graphics/plot.default.html" class="external-link">plot</a></span><span class="op">(</span><span class="va">x</span>, <span class="va">...</span>, tmax <span class="op">=</span> <span class="fl">150</span>, n_draws <span class="op">=</span> <span class="fl">2000</span>, data <span class="op">=</span> <span class="cn">NULL</span><span class="op">)</span></span></code></pre></div>
+    </div>
+
+    <div id="arguments">
+    <h2>Arguments</h2>
+
+
+<dl><dt id="arg-x">x<a class="anchor" aria-label="anchor" href="#arg-x"></a></dt>
+<dd><p>A named list of type 'biokinetics_priors'.</p></dd>
+
+
+<dt id="arg--">...<a class="anchor" aria-label="anchor" href="#arg--"></a></dt>
+<dd><p>Further arguments passed to the method.</p></dd>
+
+
+<dt id="arg-tmax">tmax<a class="anchor" aria-label="anchor" href="#arg-tmax"></a></dt>
+<dd><p>Integer. The number of time points in each simulated trajectory. Default 150.</p></dd>
+
+
+<dt id="arg-n-draws">n_draws<a class="anchor" aria-label="anchor" href="#arg-n-draws"></a></dt>
+<dd><p>Integer. The number of trajectories to simulate. Default 2000.</p></dd>
+
+
+<dt id="arg-data">data<a class="anchor" aria-label="anchor" href="#arg-data"></a></dt>
+<dd><p>Optional data.frame with columns time_since_last_exp and value. The raw data to compare to.</p></dd>
+
+</dl></div>
+    <div id="value">
+    <h2>Value</h2>
+    <p>A ggplot2 object.</p>
+    </div>
+
+  </div>
+  <div class="col-md-3 hidden-xs hidden-sm" id="pkgdown-sidebar">
+    <nav id="toc" data-toggle="toc" class="sticky-top"><h2 data-toc-skip>Contents</h2>
+    </nav></div>
+</div>
+
+
+      <footer><div class="copyright">
+  <p></p><p>Developed by Timothy Russell, Alex Hill.</p>
+</div>
+
+<div class="pkgdown">
+  <p></p><p>Site built with <a href="https://pkgdown.r-lib.org/" class="external-link">pkgdown</a> 2.1.1.</p>
+</div>
+
+      </footer></div>
+
+
+
+
+
+
+  </body></html>
+
diff --git a/sitemap.xml b/sitemap.xml
index decf323..f2c48cf 100644
--- a/sitemap.xml
+++ b/sitemap.xml
@@ -13,5 +13,7 @@
 <url><loc>/reference/convert_log2_scale_inverse.html</loc></url>
 <url><loc>/reference/epikinetics-package.html</loc></url>
 <url><loc>/reference/index.html</loc></url>
+<url><loc>/reference/plot.biokinetics_population_trajectories.html</loc></url>
+<url><loc>/reference/plot.biokinetics_priors.html</loc></url>
 </urlset>