Tutorials › Pharmacometric Workflows and Reproducible Modeling › Verifying Fixed Parameters and Estimation Controls

Quality Control and Independent Review · Tutorial 336 of 399

Verifying Fixed Parameters and Estimation Controls

Learn to distinguish fixed values from starting values, verify bounds and units, and compare estimation and numerical settings against a traceable analysis record.

Intermediate 12 min read

What You'll Learn

  • Distinguish a parameter fixed at a specified value from an estimated parameter with an initial value.
  • Compare submitted parameter bounds with the analysis record, including units and parameterization.
  • Verify that the estimation method and its required options match the stated analysis.
  • Review numerical controls using explicit, reproducible comparison rules.
  • Document discrepancies without inferring model performance from configuration settings.

Why Configuration Verification Matters

A submitted run can use the intended data and model equations yet still implement a different analysis if a parameter was estimated instead of fixed, a bound was transcribed incorrectly, or an optimizer option changed. Configuration review asks a narrow, important question: do the settings that govern estimation match the analysis record and the modeler’s stated intent? It is not a review of whether the resulting estimates are plausible or whether the outputs are internally consistent.

The earlier material on reproducing a submitted run and comparing model code with its written specification provides a foundation for this check. Here, the focus is on turning the analysis record into explicit, testable requirements, then comparing those requirements with the configuration actually used. A run control file alone is not necessarily authoritative: compare it with the approved analysis record, relevant model code, and any run-specific documentation, and resolve conflicts rather than silently choosing one source.

Definition: A fixed parameter is assigned a specified value and is not varied by the estimation procedure. An estimated parameter is varied by the procedure; its initial value is the starting point, not the final value. A bound specifies an allowed range for an estimated parameter. Estimation controls specify the method and its options, as well as numerical settings that govern how the procedure searches for a solution.

Translate Intent Into Checkable Requirements

Before inspecting the submitted settings, create a compact configuration record from the approved analysis documentation. For every parameter, record its name, role (fixed or estimated), value or initial value, lower and upper bounds where applicable, units, and parameterization. For the estimation procedure, record the method and any options that change its definition, such as whether a specified interaction option is enabled. Record numerical controls separately: examples include convergence tolerances, iteration or evaluation limits, and finite-difference step settings.

This separation prevents common category errors. A starting value is not a fixed value; an optimizer limit is not a parameter bound; and a convergence tolerance is not a criterion for statistical significance. A parameter may be estimated even when its starting value equals a bound, and that does not make it fixed. Conversely, a fixed parameter should not be treated as estimated merely because a bounds field happens to contain the same number at both ends. Check the software’s documented semantics and the actual parameter status.

Use the same parameter scale and units on both sides of the comparison. If one record uses a different but equivalent parameterization, document the conversion before comparing values. A value for clearance in L/h, for example, cannot be compared directly with a value in mL/min. Also distinguish a value represented on the natural scale from one represented on a transformed scale: the numerical entries can differ even when they describe the same quantity.

Worked Example: Fixed Value Versus Starting Value

Suppose the approved record for a hypothetical run specifies clearance \(CL\) fixed at \(5.0\) L/h, volume \(V\) estimated from an initial value of \(20\) L with bounds \(1\) to \(100\) L, and absorption rate \(k_a\) estimated from an initial value of \(1.2\) h\(^{-1}\) with bounds \(0.05\) to \(5.0\) h\(^{-1}\). A submitted configuration lists \(CL=5.0\) L/h as estimated, while its other entries match the record.

The numerical value of clearance matches, but its status does not. In the submitted run, the estimation procedure may vary \(CL\); in the approved analysis it must not. The discrepancy is therefore substantive even if an eventual output happens to report \(CL=5.0\) L/h.

For \(V\) and \(k_a\), compare the starting values and both bounds separately, along with their units. Their initial values are not estimates of the final parameter values. The review finding should identify the exact mismatch—“\(CL\) was configured as estimated rather than fixed at 5.0 L/h”—and should not claim anything about the resulting estimate or model fit.

Compare Values and Bounds Transparently

Values copied between systems can differ because of rounding or text representation. Define a comparison rule before deciding whether a numeric difference is material. For a recorded value \(r\) and submitted value \(s\), one common rule is to compute the absolute difference and compare it with a tolerance that combines absolute and relative scales:

$$ \Delta = |s-r|,\qquad \tau = \tau_{\mathrm{abs}}+\tau_{\mathrm{rel}}|r|. $$

Treat the values as matching under this rule when \(\Delta \leq \tau\). The tolerance must be chosen and documented for the particular audit; it is not a universal statistical threshold. For values expected to match exactly, such as a discrete method code or a fixed status flag, use exact comparison rather than a numeric tolerance. When \(r=0\), the absolute component is important because the relative component alone is zero.

For an estimated parameter, verify the initial value and each bound as distinct fields. Check that the initial value is permitted by the recorded bound convention. For inclusive bounds, the condition is \(L \leq x_0 \leq U\), where \(L\) is the lower bound, \(U\) is the upper bound, and \(x_0\) is the initial value. If the software imposes different endpoint rules, follow those documented rules. Do not assume that a bound is implemented on the natural scale if the software applies it on a transformed scale.

Worked Example: Converting and Comparing Bounds

In a hypothetical analysis record, clearance is estimated with initial value \(6.0\) L/h and bounds \(0.5\) to \(20.0\) L/h. A submitted configuration expresses clearance in mL/min. Use \(1\) L/h \(= 1000/60\) mL/min to translate all three quantities before comparing.

The converted initial value is \(6.0 \times 1000/60 = 100\) mL/min. The lower bound is \(0.5 \times 1000/60 = 8.333\) mL/min, and the upper bound is \(20.0 \times 1000/60 = 333.3\) mL/min, rounded to four significant figures.

Suppose the submitted entries are an initial value of \(100\) mL/min and bounds of \(8.33\) to \(300\) mL/min. The initial value and lower bound agree up to rounding, but the upper bound does not: the expected value is about \(333.3\) mL/min, not \(300\) mL/min. The discrepancy cannot be explained by the unit conversion or the displayed rounding. The review should flag the upper bound and verify whether the submitted value is a transcription error or reflects a documented change in intent.

Verify Estimation Method and Numerical Controls

Method names and options should be checked as a coordinated specification. First verify the named estimation method exactly as documented. Then verify options that alter how that method is applied, rather than assuming that a similar method label implies equivalent settings. Record whether options are enabled, disabled, or not applicable. If a setting is absent, determine whether the software supplies a documented default; an undocumented assumption that “missing means off” is not a reproducible check.

Numerical controls also require field-by-field comparison. A convergence tolerance determines a numerical stopping condition according to the implementation; an evaluation limit caps computational work; and a finite-difference step controls a numerical approximation when that approximation is used. These settings can affect whether and how an algorithm proceeds, but a configuration comparison by itself does not establish that a fit converged or that an estimate is reliable. Those questions require review of the estimation outputs and diagnostics.

Worked Example: Method Option and Evaluation Limit

A hypothetical analysis record specifies the FOCEI estimation method, the interaction option enabled, a convergence tolerance of \(1 \times 10^{-6}\), and a maximum of \(5000\) evaluations. The submitted run specifies FOCE without the interaction option, the same tolerance, and a maximum of \(500\) evaluations.

Compare the method first: FOCEI in the record does not match FOCE in the submission. Compare the method option next: interaction is required by the record but is not enabled in the submitted settings. The tolerance matches exactly, but the evaluation limit does not: \(500\) is one tenth of the specified \(5000\), a difference of \(4500\) evaluations.

The review finding should list the method, interaction setting, and evaluation limit as separate discrepancies. It should not state that the lower limit caused a failed fit or that the chosen method produced a particular estimate; those conclusions would require evidence from the run outputs. The settings alone establish only that the submitted configuration does not match the record.

A Reproducible Review Procedure

The following sequence makes the audit traceable and reduces the chance that a correct numerical value is accepted under the wrong status, scale, or setting.

1
Establish the comparison source.
Identify the approved analysis record and the submitted run configuration. If they conflict, document the conflict and obtain a resolved source of intent before declaring one configuration correct.
2
Build a parameter inventory.
For each parameter, record fixed or estimated status, specified or initial value, bounds, units, and parameterization. Include parameters that are fixed, not only those with bounds.
3
Normalize before comparison.
Convert units or parameter scales using documented relationships. Apply a pre-specified numeric tolerance only where approximate numeric comparison is appropriate; compare statuses and discrete method options exactly.
4
Check method and controls.
Compare the estimation method, its relevant options, and each numerical control as separate fields. Resolve omitted settings by checking the applicable software defaults rather than guessing.
5
Record evidence and disposition.
For each field, retain the expected setting, observed setting, comparison result, source or file location, and any resolution. Separate confirmed mismatches from unresolved documentation gaps.

Why the Comparison Rule Works—and Where It Can Mislead

The tolerance rule is a direct way to allow for harmless representation differences without hiding meaningful discrepancies. If a recorded value \(r\) is nonzero, \(\tau_{\mathrm{rel}}|r|\) scales with the size of the value, while \(\tau_{\mathrm{abs}}\) protects comparisons near zero. Thus \(\Delta \leq \tau\) means the observed difference is no larger than the declared allowance at that scale. This conclusion is only as defensible as the tolerance choice: a generous tolerance can incorrectly label a changed setting as a match.

For example, let \(r=2.00\), \(s=2.002\), \(\tau_{\mathrm{abs}}=0.0001\), and \(\tau_{\mathrm{rel}}=0.001\). Then \(\Delta=|2.002-2.00|=0.002\), while \(\tau=0.0001+0.001(2.00)=0.0021\). The comparison passes because \(0.002 \leq 0.0021\). If the relative tolerance had instead been \(0.0005\), then \(\tau=0.0011\), and the same values would fail. Neither setting is inherently right for every audit; the chosen tolerance must reflect the record’s precision and the consequences of changing that parameter.

A direct field-by-field comparison is usually preferable for independent review because it makes discrepancies easy to localize and report. Automated configuration diffs can help with scale, but should be checked for aliases, defaults, omitted fields, units, and formatting differences. A full rerun can establish what a submitted configuration does in practice, but it does not replace checking whether that configuration matches the approved intent. These approaches complement rather than substitute for a clear analysis record.

Diagnostics, Limitations, and Common Pitfalls

A useful diagnostic is a completed comparison table in which every required field has a recorded expected value, observed value, and disposition. Blank cells, unexplained defaults, inconsistent units, or status labels that do not map clearly to the software’s parameter behavior are signs that the review is incomplete. When a setting cannot be interpreted from the available documentation, classify it as unresolved; do not convert uncertainty into a pass.

  • Confusing a starting value with a fixed value. Matching starting values do not show that the parameter was held constant. Check the parameter’s estimation status and implementation.
  • Comparing bounds without units or scale. Equal-looking numbers can represent different physical values; different-looking numbers can be equivalent after a documented conversion.
  • Assuming defaults. An absent option is not automatically equivalent to the intended default. Verify the relevant software behavior and version.
  • Treating tolerances as universal. A tolerance is an audit rule, not a property of the parameter or a statistical significance threshold. State it and justify its use.
  • Inferring success from settings. Matching controls show configuration agreement, not that the estimation completed successfully, reached a suitable solution, or yielded precise estimates.

This review is limited to configuration fidelity. It does not establish that the model is scientifically appropriate, that the analysis record is itself correct, or that the estimation outputs are internally consistent. Those require their own evidence and checks. The purpose here is to detect and document whether the submitted run was set up as intended.

Key takeaway: Verify parameter status, values, bounds, units, estimation method, method options, and numerical controls as distinct, traceable fields. A matching number alone is not proof of a matching analysis configuration.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. A parameter is specified as estimated with initial value \(3.0\), and the submitted run lists it as fixed at \(3.0\). Do the settings match? Why?
  2. A recorded bound is \(0.2\) L/h and the submitted bound is \(3.33\) mL/min. Using \(1\) L/h \(=1000/60\) mL/min, are these approximately equivalent?
  3. The record gives a maximum of \(2000\) evaluations and the run configuration gives \(1000\). Can a reviewer conclude from these settings alone that the run failed?
  4. For \(r=10.0\), \(s=10.01\), \(\tau_{\mathrm{abs}}=0.001\), and \(\tau_{\mathrm{rel}}=0.001\), does the comparison rule \(\Delta \leq \tau_{\mathrm{abs}}+\tau_{\mathrm{rel}}|r|\) pass?

Solutions. (1) No. The value matches, but the status differs: the record requires estimation while the submission fixes the parameter. (2) Yes, approximately: \(0.2 \times 1000/60 = 3.333\) mL/min, which is close to \(3.33\) mL/min at the displayed precision; apply the audit’s stated tolerance. (3) No. The evaluation limit is a configuration mismatch, but settings alone do not establish whether the run failed. (4) Yes. \(\Delta=|10.01-10.0|=0.01\), and \(\tau=0.001+0.001(10.0)=0.011\), so \(0.01 \leq 0.011\).

Clinical Biostats