Tutorials › Optimal Design and Model-Informed Dosing › Recomputing Information After Interim Data

Robust and Adaptive Design · Tutorial 135 of 497

Recomputing Information After Interim Data

Learn to combine accrued information with posterior-averaged information from candidate future samples, and check that the resulting comparison is valid.

Advanced 12 min read

What You'll Learn

  • Define interim-conditional expected Fisher information for a candidate future design.
  • Add accrued and future information on the same likelihood-information scale.
  • Average parameter-dependent information over the interim parameter distribution.
  • Compare candidate schedules without treating a plug-in estimate as an uncertainty average.
  • Diagnose double counting, invalid additivity, and model or design mismatch.

Why Recompute Information at an Interim Analysis?

At an interim analysis, some observations have already been collected, and the current parameter distribution may differ from the one used to plan the original design. A useful comparison of candidate future samples therefore needs two pieces: the information contributed by the data already collected, and the expected information from each proposed future design under the updated parameter uncertainty. Comparing future samples alone can give the wrong impression of how much total information a design will provide.

The calculation is conditional on the interim data and on a specified model. It is not a new parameter update: the interim parameter distribution is an input. The task is to evaluate the information in each candidate future design across that distribution, then account for the accrued information once.

Definition: Let \(D\) denote the interim data and \(\theta\) the model parameter vector. For a candidate future design \(d\), the interim-conditional expected total Fisher information is \(E_{\theta\mid D}\{I_A(\theta)+I_F(\theta;d)\}\), where \(I_A(\theta)\) is the Fisher information from the accrued design and \(I_F(\theta;d)\) is the Fisher information expected from the future design, both evaluated at the same \(\theta\). For a vector parameter, these are matrices; for a scalar parameter, they are numbers.

Information to Add, and What the Average Means

For a likelihood with score vector \(U(\theta)=\partial \ell(\theta)/\partial\theta\), the Fisher information is the expected outer product of the score, \(I(\theta)=E_\theta\{U(\theta)U(\theta)^\mathsf{T}\}\), under regularity conditions. Equivalently, when differentiation and expectation can be interchanged, it is the negative expected Hessian of the log likelihood. The information is a function of the parameter, even when the observations themselves have already been collected.

For fixed, conditionally independent observations, log likelihoods add, so their scores add. Each score has expectation zero under the model. Consequently, the expected cross-products between scores from distinct observations are zero, and Fisher information adds. If \(A\) denotes the accrued sampling plan and \(d\) the proposed future plan, then

$$ I_{\mathrm{total}}(\theta;d)=I_A(\theta)+I_F(\theta;d). $$

The interim distribution is then used to average this total. Linearity of expectation gives an equivalent and often convenient calculation: average the accrued and future contributions separately, then add them. The accrued contribution is common to all candidates if every candidate continues from the same interim design. It still belongs in the reported total, although it does not change a ranking based only on a scalar criterion that increases with total information.

1
Fix the interim information.
Specify the accrued design and the interim parameter distribution. Keep the latter conditional on the observed interim data.
2
Evaluate each contribution at a parameter value.
Compute the accrued-design information and the candidate future-design information using the same model and parameter value.
3
Average and add.
Average the total information over the interim parameter distribution, or average the two components separately and add them.
4
Compare on a stated scale.
For scalar information, compare the expected totals. For matrices, state the matrix summary or design criterion used; a matrix comparison need not give a unique ranking.

Worked Examples

Worked Example: Comparing Two Future Sampling Times

Consider a hypothetical scalar model with parameter \(k>0\). Each independent measurement at time \(t\) is normally distributed as \(Y(t)\sim N(e^{-kt},1)\), with known variance one. The derivative of the mean with respect to \(k\) is \(-t e^{-kt}\), so the Fisher information from one measurement is \(I(t;k)=t^2e^{-2kt}\). Suppose 25 measurements have already been taken at \(t=1\), and the interim distribution assigns probability 0.75 to \(k=0.2\) and probability 0.25 to \(k=0.8\). Compare 25 additional measurements at \(t=2\) with 25 at \(t=4\), using expected total scalar information.

The accrued information is \(25e^{-2k}\). At \(k=0.2\), it is \(25e^{-0.4}=16.758\); at \(k=0.8\), it is \(25e^{-1.6}=5.047\). These values are the same for both future candidates.

At \(t=2\), the future information is \(25(2^2)e^{-4k}=100e^{-4k}\). At \(k=0.2\), this is \(100e^{-0.8}=44.9329\), giving total information \(16.7580+44.9329=61.6909\). At \(k=0.8\), it is \(100e^{-3.2}=4.07622\), giving total \(5.04741+4.07622=9.12363\). Thus the posterior-weighted mean total is

$$ 0.75(61.6909)+0.25(9.12363)=48.5491. $$

At \(t=4\), the future information is \(25(4^2)e^{-8k}=400e^{-8k}\). At \(k=0.2\), this is \(400e^{-1.6}=80.7586\), giving total \(16.7580+80.7586=97.5166\). At \(k=0.8\), it is \(400e^{-6.4}=0.664623\), giving total \(5.04741+0.664623=5.71204\). The posterior-weighted mean total is

$$ 0.75(97.5166)+0.25(5.71204)=74.5655. $$

Under this particular criterion and these hypothetical assumptions, \(t=4\) has the larger posterior-mean total scalar information: 74.57 rather than 48.55. This is a comparison of expected Fisher information, not a guarantee about the information realized in a particular future sample.

Worked Example: Adding Accrued and Future Information Matrices

Suppose a hypothetical linear model has independent observations \(Y(t)\sim N(\beta_0+\beta_1t,1)\), where the variance is known. For one observation, the information matrix for \((\beta_0,\beta_1)\) is \(x(t)x(t)^\mathsf{T}\), with \(x(t)=(1,t)^\mathsf{T}\). Six observations have already been collected at \(t=0\) and four at \(t=1\). Consider adding three observations at \(t=2\).

The accrued matrix is the sum of the contributions at the two observed design points:

$$ I_A=6 \begin{pmatrix}1&0\\0&0\end{pmatrix} +4 \begin{pmatrix}1&1\\1&1\end{pmatrix} = \begin{pmatrix}10&4\\4&4\end{pmatrix}. $$

For the three future observations, \(x(2)=(1,2)^\mathsf{T}\), so

$$ I_F=3 \begin{pmatrix}1&2\\2&4\end{pmatrix} = \begin{pmatrix}3&6\\6&12\end{pmatrix}, \qquad I_{\mathrm{total}}= \begin{pmatrix}13&10\\10&16\end{pmatrix}. $$

In this model the information does not depend on \((\beta_0,\beta_1)\). Averaging either matrix over any interim parameter distribution leaves it unchanged. As a basic check, the determinant of the total matrix is \(13(16)-10^2=108>0\), so the matrix is positive definite. The example also shows why accrued and future information must be expressed in compatible parameter order and units before being added.

Worked Example: Averaging Across Three Interim Parameter Values

Return to the scalar model \(I(t;k)=t^2e^{-2kt}\) for one measurement. Suppose two measurements have been accrued at \(t=1\), and the interim distribution places probabilities 0.20, 0.50, and 0.30 on \(k=0.1,0.5,\) and \(0.9\), respectively. A candidate adds one measurement at \(t=2\). Compute the expected total information.

At \(k=0.1\), accrued information is \(2e^{-0.2}=1.63746\), and future information is \(4e^{-0.4}=2.68128\), for a total of 4.31874. At \(k=0.5\), the corresponding values are \(2e^{-1}=0.735759\) and \(4e^{-2}=0.541341\), for a total of 1.27710. At \(k=0.9\), they are \(2e^{-1.8}=0.330598\) and \(4e^{-3.6}=0.109268\), for a total of 0.439866.

Weighting each total by its interim probability gives

$$ 0.20(4.31874)+0.50(1.27710)+0.30(0.439866)=1.63426. $$

As a check, the expected accrued information is \(0.20(1.63746)+0.50(0.735759)+0.30(0.330598)=0.794551\), and the expected future information is \(0.20(2.68128)+0.50(0.541341)+0.30(0.109268)=0.839707\). Their sum is 1.63426, matching the weighted total.

Assumptions, Diagnostics, and Alternatives

The addition-and-average calculation requires a clearly specified model and design. In particular, the candidate future schedule is treated as fixed conditional on the interim data; its information is evaluated under the assumed observation model; and the likelihood contributions are conditionally independent, or have their dependence correctly represented. The posterior distribution used for averaging must correspond to the same parameterization as the information calculation. If the future sampling plan changes in response to future outcomes, its information must account for that outcome-dependent rule rather than treating the schedule as fixed.

Several checks can reveal errors before a design comparison is reported:

  • Check the accounting. Each accrued observation should appear once, and future observations should not be included in the accrued component. Candidate totals should share the same accrued component when they start from the same interim design.
  • Check dimensions and units. Information matrices for different candidates must use the same parameter order. Scalar information has units reciprocal to the squared parameter unit; matrix entries have units determined by the corresponding parameter pair.
  • Check numerical behavior. Information matrices should be symmetric and positive semidefinite up to numerical rounding. Recalculate exponentials, weights, and matrix sums at adequate precision; posterior weights should be nonnegative and sum to one.
  • Check the model and support. If the posterior places weight in parameter regions where the model is invalid or the information calculation is numerically unstable, the resulting average is not trustworthy. Confirm that every parameter value used is admissible.

A common mistake is to compute the future information at a single interim point estimate and then add the accrued information at that same point, calling the result an uncertainty-averaged total. That calculation is a plug-in evaluation, not an average over the interim parameter distribution. The correction is to evaluate \(I_A(\theta)+I_F(\theta;d)\) across the distribution and average those values, or equivalently average the two terms separately. Another mistake is to add prior precision as though it were accrued likelihood information. If a criterion uses posterior precision, prior curvature may be relevant, but it must be identified separately and not counted as data-derived Fisher information.

The main alternative summaries answer different questions. A nominal calculation evaluates information at one designated parameter value; it is simpler, but does not average over interim uncertainty. A worst-case calculation takes the least favorable information over a specified parameter set, emphasizing protection rather than posterior average performance. For matrix information, averaging the matrix first and then applying a criterion is generally not equivalent to applying that criterion at each parameter value and then averaging. The chosen operation should match the design objective rather than being treated as interchangeable.

Expected Fisher information is a model-based planning quantity, not realized information, a posterior covariance, or a direct measure of clinical benefit. Its reliability depends on the model, parameter distribution, and proposed schedule. A posterior average can also conceal a poor result in a low-probability but important region; when that matters, report additional summaries of performance across parameter values. For a longer adaptive sequence, future information may depend on interim outcomes and the policy that selects later samples, requiring an expectation over those outcomes as well as parameter uncertainty.

Key takeaway: For each candidate future design, calculate accrued plus future Fisher information at each relevant parameter value, then average that total using the interim parameter distribution. Keep the data-information calculation distinct from prior information, posterior covariance, and realized future outcomes.

Check Your Understanding

Work through these questions, then compare your reasoning with the solutions.

  1. A design has accrued scalar information 7 at every parameter value. A candidate contributes expected future information 3.5 under the interim distribution. What is the expected total information?
  2. The interim distribution assigns probabilities 0.25 and 0.75 to two parameter values. A candidate's future information at those values is 8 and 2. What is its expected future information?
  3. Why should a prior precision matrix not automatically be added to accrued Fisher information when reporting total data information?
  4. In the scalar model \(I(t;k)=t^2e^{-2kt}\), how much information is contributed by two independent measurements at \(t=1\) when \(k=0.5\)?

Solutions. (1) Add the common accrued component once: \(7+3.5=10.5\). (2) The weighted average is \(0.25(8)+0.75(2)=2+1.5=3.5\). (3) Prior precision comes from the prior distribution, not from the accrued likelihood; including it changes the quantity from data Fisher information to a posterior-precision-type summary, so it must be stated and handled separately. (4) Each measurement contributes \(1^2e^{-2(0.5)(1)}=e^{-1}=0.367879\); two independent measurements contribute \(2e^{-1}=0.735759\), or 0.7358 to four significant figures.

Clinical Biostats