## [1] "DualEndpointEmax"
The relationships between dose and toxicity and between dose and PD biomarker will be modelled simultaneously.
A probit log normal model will describe the relationship between dose and toxicity: \[ \Phi^{-1}(Tox | d) = f(X = 1 | \theta, d) = \alpha + \beta \cdot log(d/d^*) \]where d* denotes a reference dose.
The prior for θ is given by\[ \boldsymbol\theta = \begin{bmatrix}\alpha \\ \beta\end{bmatrix}\sim N \left(\begin{bmatrix} 0.00 \\ 1.00\end{bmatrix} , \begin{bmatrix} 1.00 & 0.00 \\ 0.00 & 1.00\end{bmatrix} \right) \]
The reference dose will be 1.00.
The PD biomarker response w at dose d is
modelled as \[ w(d) \sim N(f(d), \sigma_w^2)
\]
where f(d) is a parametric Emax density function such that
\[ f(d) = E_0 + \frac{(E_{max} - E_0) \times \frac{d}{d^*}}{\text{ED}_{50} + \frac{d}{d^*}} \]
where d* is the reference dose, E0 and Emax are, respectively, the minimum and maximum levels of the biomarker and ED50 is the dose achieving half the maximum effect, 0.5 × Emax.
In this case,
\[ E_0 \sim U(0, 100) \]
\[ E_{max} \sim U(0, 500) \]
\[ \text{ED}_{50} \sim U(10, 200) \]
and
\[ d^* = 1000 \]