From a Precision Goal to a Data Plan
A precision goal is useful only if it can guide a design decision. For example, a team may require a 95% confidence interval for a pathway rate constant to have a half-width no greater than a specified value. The design question is then concrete: how many independent observations are needed, how many replicates should be collected, which sampling times are most informative, or how precise must the assay be?
Fisher information provides a local, model-based translation from a design to expected parameter uncertainty. Its approximation is most useful when the model is locally identifiable, the likelihood is sufficiently regular, and the estimate is not strongly affected by a boundary or severe nonlinearity. Profile likelihood provides a complementary check by evaluating likelihood support over a range of parameter values rather than relying only on local curvature.
Translate an Interval Target Into Required Information
For a scalar parameter \(\theta\), let \(I(\theta)\) denote the Fisher information for the planned observations, evaluated at a plausible parameter value. Under appropriate regularity conditions, an asymptotically efficient estimator (such as the maximum likelihood estimator) is approximately normal with variance \(1/I(\theta)\). A two-sided approximate 95% confidence interval therefore has half-width \(1.96/\sqrt{I(\theta)}\). If the desired half-width is at most \(h\), the local information requirement is:
This calculation is for a marginal interval for one parameter under the stated normal approximation. A one-sided target, a joint confidence region, or a simultaneous guarantee across several parameters requires an appropriate critical value or joint criterion; the number 1.96 should not be used automatically for every precision statement. For a relative-precision target, it is often clearer to work with a transformed parameter such as \(\log(\theta)\), provided \(\theta>0\), and define the desired uncertainty on that scale.
Information depends on the model, the assumed parameter values, the observation error, and the design. Thus a calculation does not establish that a design will meet its target in every realized data set. It predicts precision under the assumptions used to compute the information.
Why Replicates and Independent Observations Can Be Counted
For independent observations with the same design and error model, log-likelihood contributions add, so their Fisher information adds. If one observation contributes information \(I_1\), then \(n\) such observations contribute \(nI_1\). In the scalar case, approximate standard error consequently decreases as \(1/\sqrt{n}\), not as \(1/n\).
With multiple estimated parameters, the information is a matrix. The variance for a parameter of interest depends not only on its own diagonal information but also on its correlation with nuisance parameters. For a two-parameter information matrix with target \(\theta\) and nuisance parameter \(\phi\), the local information for \(\theta\) after accounting for uncertainty in \(\phi\) is the Schur complement \(I_{\theta\theta}-I_{\theta\phi}^2/I_{\phi\phi}\). Ignoring that adjustment can overstate precision.
Worked Example: Required Number of Independent Observations
Suppose a hypothetical model has a scalar rate parameter \(k\), measured in \(\mathrm{h}^{-1}\). Under a proposed fixed sampling design and assumed observation-error model, each independent observation is predicted to contribute 12 units of Fisher information, expressed in the inverse squared units of \(k\). The target is an approximate two-sided 95% interval with half-width no greater than \(0.0500\ \mathrm{h}^{-1}\).
The required total information is:
With \(n\) independent observations, the predicted information is \(12n\ \mathrm{h}^2\). Therefore:
The required integer sample size is \(n=129\). At that size, the predicted half-width is \(1.96/\sqrt{12(129)}=1.96/\sqrt{1548}\approx0.04982\ \mathrm{h}^{-1}\), just within the target. This result depends on the assumed information per observation being valid across the planned observations and on the observations being independent.
Replicate Counts, Assay Precision, and Sampling Times
When each replicate has the same independent contribution to information, the replicate count can be found by dividing required information by information per replicate and rounding up. This is not the same as simply multiplying the number of measurements if replicates share substantial within-subject or batch correlation. In that case, the joint likelihood and its dependence structure determine the information; treating correlated measurements as independent exaggerates the gain.
Worked Example: Replicates With a Nuisance Parameter
Consider a hypothetical two-parameter model with target \(\theta\) and nuisance parameter \(\phi\). A single independent replicate is expected to contribute the information matrix \(J=\begin{pmatrix}4&1.2\\1.2&1\end{pmatrix}\), with entries in the corresponding inverse parameter units. The target is a marginal 95% interval for \(\theta\) with half-width at most \(0.250\) units.
The information matrix is positive definite because its diagonal entries are positive and its determinant is \(4(1)-1.2^2=2.56>0\). Accounting for \(\phi\), the effective information about \(\theta\) from one replicate is:
The target requires \(I_{\mathrm{required}}=(1.96/0.250)^2=61.4656\). With \(m\) independent replicates, the effective information is \(2.56m\), so \(m\ge61.4656/2.56=24.01\). Rounding up gives \(m=25\). The predicted half-width is \(1.96/\sqrt{25(2.56)}=1.96/8=0.245\) units. Twenty-four replicates would give information \(61.44\), slightly below the requirement, so rounding down would not meet the stated target.
For independent normally distributed assay errors with known standard deviation \(\sigma\), a single observation with model mean \(f(\theta)\) contributes scalar information \(\{ \partial f(\theta)/\partial\theta\}^2/\sigma^2\). Thus, with a fixed number of observations and fixed sensitivities, information is proportional to \(1/\sigma^2\). This relationship makes it possible to translate a parameter target into an assay precision requirement.
Worked Example: Required Assay Precision
In a hypothetical one-parameter model, 40 independent measurements are planned. At the design point, the predicted sensitivity of the measurement mean to parameter \(\theta\) is \(0.800\) assay units per unit of \(\theta\). Assume additive, independent, normally distributed errors with a common known standard deviation \(\sigma\). The desired two-sided 95% interval half-width for \(\theta\) is \(0.100\) units.
The required information is \((1.96/0.100)^2=384.16\). The planned information is \(40(0.800)^2/\sigma^2=25.6/\sigma^2\). Meeting the target therefore requires:
Under these assumptions, the assay standard deviation must be no greater than approximately \(0.2581\) assay units. If measurement error is estimated rather than known, has a different distribution, or is correlated across observations, this calculation alone does not capture all uncertainty.
Sampling times change the sensitivities of model outputs to parameters, and therefore change the information contributed at each time. For a single parameter, a simple design comparison can sum the expected information at the selected times. With several parameters, add the information matrices first and then assess the target parameter’s uncertainty; nuisance-parameter correlations can make the best schedule different from the one suggested by one sensitivity alone.
Worked Example: Comparing Candidate Sampling Schedules
Suppose a hypothetical one-parameter model predicts effective information contributions of \(0.200\), \(1.10\), and \(0.500\) inverse squared parameter units per independent measurement at times 1, 4, and 8 hours, respectively. Four measurements are available, and the target is a 95% interval half-width no greater than 1.00 parameter unit. The required information is \((1.96/1.00)^2=3.8416\).
A schedule with two measurements at 4 hours and one each at 1 and 8 hours provides \(2(1.10)+0.200+0.500=2.90\) units of information. Its approximate half-width is \(1.96/\sqrt{2.90}\approx1.151\) units, so it misses the target. Four independent measurements at 4 hours provide \(4(1.10)=4.40\) units; the predicted half-width is \(1.96/\sqrt{4.40}\approx0.9344\) units, meeting the target. This simplified comparison assumes the one-parameter information contributions are additive and that all four measurements can be taken at the selected time. Real designs may impose timing, burden, or dependence constraints.
Check the Local Prediction With Profile Likelihood
A Fisher-information calculation uses local curvature near a parameter estimate or assumed design point. If the likelihood is asymmetric, flat in one direction, multimodal, or constrained by a parameter boundary, its normal-based interval may give a poor description of uncertainty. The profile likelihood addresses this by fixing the target parameter at candidate values and maximizing the likelihood over the remaining parameters. A profile-based interval consists of values whose likelihood-ratio statistic does not exceed the chosen cutoff.
For prospective planning, a practical approach is to specify a candidate design, generate or otherwise evaluate its expected data under plausible parameter values, fit the model, and calculate the target parameter’s profile interval. The design meets a maximum half-width target only if the profile interval satisfies the prespecified width criterion. Because a future data set is not known, repeated simulations can estimate how often the criterion is met. The likelihood-ratio cutoff and any small-sample calibration should match the intended inferential setting.
The two approaches answer related but different questions. Fisher information is efficient for screening many designs and showing how changes in sample size, replicate count, or error scale affect local precision. Profile likelihood checks the shape of likelihood support beyond the local approximation, but requires repeated fitting and profiling and can be computationally demanding. If a target must be met with a stated probability across future studies, simulation of the full design and analysis is more informative than either a single expected-information calculation or a single profile.
State the parameter or parameter combination, its scale and units, the interval type, and the maximum allowed uncertainty.
Choose sampling times, replicate structure, observation-error assumptions, and parameter values at which expected information will be evaluated.
Convert the target to an information threshold, calculate information for each candidate design, and account for nuisance parameters.
Inspect identifiability and numerical diagnostics, then use profile likelihood or simulation when local approximation or expected performance is not adequate.
Diagnostics, Limitations, and a Common Pitfall
A calculated data requirement is suspect if the information matrix is singular or nearly singular, if estimated parameters approach bounds, if fitted profiles are highly asymmetric, or if small changes in the assumed parameter values substantially change the required design. Numerical instability in optimization or profiling is another warning: an apparently narrow interval is not credible if the fit has not converged reliably. Sensitivity to the assumed observation-error variance also deserves attention because information scales inversely with variance in the simple Gaussian setting.
A common mistake is to calculate the required information for each parameter separately and conclude that every target is met once its diagonal information is large enough. This ignores covariance with other parameters. The correct marginal variance comes from the inverse of the full information matrix, not from the reciprocal of a diagonal entry alone. Another mistake is to treat a nominal expected-information target as a guarantee: the realized data can be less informative than expected, and the local approximation itself can fail. Evaluate the full matrix, inspect profiles where needed, and make clear whether the goal is expected precision or a probability of meeting a target.
These calculations also rely on the assumed model and design being meaningful. More observations at an uninformative time may add little information; improved assay precision cannot resolve structural confounding by itself. A precision target may be infeasible for the proposed measurement types and parameterization. In that case, the result is a reason to reconsider the design, the estimand, or the target—not to report a rounded sample size as if it solved the underlying identifiability problem.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- A scalar parameter has a 95% interval half-width target of 0.20 units. What information is required under the local normal approximation?
- Each independent observation contributes 5 units of scalar information. How many observations are needed to meet a half-width target of 0.40 units?
- Why can a large diagonal entry in a multi-parameter information matrix fail to imply precise estimation of that parameter?
- For fixed sensitivities and independent normal errors, what happens to Fisher information if the assay standard deviation is doubled?
- Name two reasons to check a Fisher-based precision prediction with profile likelihood or simulation.
Solutions. (1) \(I_{\mathrm{required}}=(1.96/0.20)^2=96.04\) inverse squared parameter units. (2) The required information is \((1.96/0.40)^2=24.01\), so \(5n\ge24.01\) and the minimum integer sample size is \(n=5\). (3) Correlation with nuisance parameters can reduce the effective information; the marginal variance is obtained from the inverse of the full information matrix. (4) Information is proportional to \(1/\sigma^2\), so doubling \(\sigma\) reduces information to one quarter. (5) For example, the profile may be asymmetric or flat, estimates may lie near boundaries, or the design’s realized precision may vary enough that an expected-information calculation does not answer the probability-of-success question.