Skip to content

Commit d585de6

Browse files
seabbs-botseabbs
andcommitted
docs: fix notation consistency in marginal model section
Use \mathbb{P} instead of \mathrm{Pr} to match the rest of the document, and lowercase p_{L,i}, p_{R,i} in the likelihood denominator to match the latent model convention. Co-authored-by: Sam Abbott <contact@samabbott.co.uk>
1 parent 9806bfe commit d585de6

1 file changed

Lines changed: 2 additions & 2 deletions

File tree

vignettes/model.Rmd

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -184,14 +184,14 @@ This approach uses the primary event censored distribution implemented in the [`
184184

185185
Under the assumption that the forward distribution is shift-invariant (i.e. $f_x = f$ for all $x$), the double censoring probability from Section \@ref(interval-censoring) simplifies to
186186
$$
187-
\mathrm{Pr}(S_L < S < S_R \mid P_L < P < P_R) = \int_{P_L}^{P_R} g_P(x \mid P_L, P_R) \left[F(S_R - x) - F(S_L - x)\right] \text{d}x.
187+
\mathbb{P}(S_L < S < S_R \mid P_L < P < P_R) = \int_{P_L}^{P_R} g_P(x \mid P_L, P_R) \left[F(S_R - x) - F(S_L - x)\right] \text{d}x.
188188
$$
189189
For common delay and primary event distributions, such as gamma or lognormal delays with uniform primary events, `primarycensored` provides closed-form analytical solutions to this integral.
190190
For other combinations, numerical integration is used.
191191

192192
Right truncation at time $T$ is handled by normalising the likelihood as in the latent model:
193193
$$
194-
\mathcal{L}(\mathbf{Y} \mid \mathbf{\theta}) = \prod_i \frac{\mathrm{Pr}(S_{L,i} < S_i < S_{R,i} \mid P_{L,i} < P_i < P_{R,i})}{\int_{P_{L,i}}^{P_{R,i}} g_P(z \mid P_{L,i}, P_{R,i}) F(T - z) \, \text{d}z}.
194+
\mathcal{L}(\mathbf{Y} \mid \mathbf{\theta}) = \prod_i \frac{\mathbb{P}(S_{L,i} < S_i < S_{R,i} \mid P_{L,i} < P_i < P_{R,i})}{\int_{P_{L,i}}^{P_{R,i}} g_P(z \mid p_{L,i}, p_{R,i}) F(T - z) \, \text{d}z}.
195195
$$
196196

197197
Since no latent variables are required, identical observations can be aggregated and the likelihood computed once per unique combination of delay, censoring windows, and covariates.

0 commit comments

Comments
 (0)