Tutorials › Optimal Design and Model-Informed Dosing › Quantization and Rounding in Titration Algorithms

Dose Regimens and Titration Algorithms · Tutorial 268 of 497

Quantization and Rounding in Titration Algorithms

Translate continuous dose adjustments into feasible dose strengths and determine how rounding changes target attainment and convergence.

Advanced 13 min read

What You'll Learn

  • Define a dose quantization rule on a specified set of available doses
  • Calculate the dose and response discrepancy caused by rounding
  • Derive the residual-error band in which a rounded update can stop changing the dose
  • Distinguish target rounding from rounding each step of an iterative rule
  • Diagnose dose cycles and other effects that rounding can introduce

From a Continuous Recommendation to a Usable Dose

A titration rule may calculate a dose such as 52 mg/day, while the available tablet strengths or administration constraints permit only doses on a coarser grid. The prescribed dose must then be selected from the feasible options. This is not merely a matter of presentation: rounding can change the achieved exposure or response, stop an iterative titration before its continuous target is reached, or contribute to oscillation between dose levels.

The key is to distinguish two operations: rounding a final, continuously calculated target dose once, and rounding after every update in an iterative algorithm. They need not produce the same trajectory. The analysis below uses a deterministic dose-response model and dose updates; practical decisions also need to account for measurement error, patient safety, and the actual available formulations.

Definition: Let \(G\) be the set of feasible doses, including any minimum and maximum dose constraints. A quantization rule \(Q_G(x)\) maps a continuous dose recommendation \(x\) to a dose in \(G\). Under nearest-dose quantization, \(Q_G(x)\) is a feasible dose minimizing \(|x-d|\) over \(d\in G\). If two doses are equally near, a specified tie rule is required.

For a uniform dose grid, write \(G=\{d_0+jh:j\text{ is an allowed integer}\}\), where \(d_0\) is the grid origin and \(h\) is the spacing. The spacing alone does not identify the available doses: both \(d_0\) and \(h\), as well as the bounds, matter. If there are irregular combinations of tablets or other feasibility restrictions, the set \(G\) may not be a uniform grid; the nearest feasible dose must then be found from the actual set.

Assumptions and the Effect of One Rounding Step

Suppose a dose-response relationship \(R(D)\) is known and sufficiently accurate for the decision, and a target response \(R^*\) corresponds to a continuous target dose \(D^*\). For the simple calculations below, assume \(R(D)\) is increasing and locally linear, \(R(D)=a+bD\), with \(b>0\). Thus \(D^*=(R^*-a)/b\). The model is treated as deterministic: the examples do not include response measurement noise, parameter uncertainty, or a delay between changing the dose and measuring the response.

If the continuous target is rounded once to the nearest dose on an unbounded uniform grid of spacing \(h\), the dose error is at most half a grid step. This follows because every real number lies within \(h/2\) of a grid point. With the linear response assumption, the corresponding response error is bounded by the slope times that dose error.

$$ |Q_G(D^*)-D^*|\leq \frac{h}{2}, \qquad |R(Q_G(D^*))-R^*| =b|Q_G(D^*)-D^*| \leq \frac{bh}{2}. $$

These bounds assume the nearest grid dose is feasible. At a minimum or maximum dose, a bound can still be calculated by comparing the target with the feasible set, but it need not be \(h/2\): the target may lie outside the allowable range. With a nonuniform set of doses, the relevant error is the actual distance to the selected feasible dose, not half of a presumed spacing.

Worked Example: Rounding to Explicitly Specified Strengths

Suppose a linear exposure-response calculation gives a target dose of 73 mg/day. Available doses are multiples of 20 mg/day, so the feasible doses include 60 and 80 mg/day. Let the hypothetical response model be \(R(D)=1.2D\), with response measured in a consistent response unit. What dose is nearest, and what response difference does this rounding imply?

The two candidates surrounding 73 mg/day are 60 and 80 mg/day. Their distances from the target are \( |73-60|=13\) mg/day and \( |80-73|=7\) mg/day, respectively. Therefore 80 mg/day is the unique nearest dose. Specifying that the doses are multiples of 20 mg/day matters: the spacing alone would not establish that 60 and 80 mg/day are available.

At the continuous target, the model predicts \(R(73)=1.2(73)=87.6\). At the selected dose it predicts \(R(80)=1.2(80)=96\). The dose rounding error is \(80-73=7\) mg/day, so the response difference is \(1.2(7)=8.4\) response units. This agrees with \(96-87.6=8.4\). The response is higher than the target-model response because the selected dose is above the continuous target.

Rounding at Every Iteration

An iterative algorithm can be affected differently because it does not necessarily calculate \(D^*\) and round that value once. Consider the proportional dose update \(D_{k+1}^{\mathrm{cont}}=D_k+\lambda(D^*-D_k)\), where \(D_k\) is the current dose, \(D_{k+1}^{\mathrm{cont}}\) is the continuous next-dose recommendation, and \(\lambda\) is a positive adjustment gain. If the dose actually administered must be feasible, the implemented rule is \(D_{k+1}=Q_G(D_k+\lambda(D^*-D_k))\).

Before quantization, the dose error \(e_k=D^*-D_k\) obeys a simple recurrence. Substituting the continuous update gives \(e_{k+1}=D^*-D_k-\lambda(D^*-D_k)=(1-\lambda)e_k\). In the absence of rounding, the error therefore shrinks in magnitude when \(0<\lambda<2\), remains unchanged in magnitude at the boundary \(\lambda=2\), and grows for \(\lambda>2\). This continuous result alone does not establish convergence of the rounded rule.

For a uniform grid, write the quantization discrepancy at update \(k+1\) as \(q_{k+1}=D_{k+1}-(D_k+\lambda e_k)\). If the recommendation is within the grid bounds and nearest-dose quantization is used, \(|q_{k+1}|\leq h/2\). The rounded error consequently satisfies \(e_{k+1}=(1-\lambda)e_k-q_{k+1}\), and hence \(|e_{k+1}|\leq |1-\lambda||e_k|+h/2\). This inequality shows that rounding can leave a residual error even when the continuous rule would continue making progress. It is a bound on one update, not a proof that the rounded sequence converges.

A particularly important case occurs when an update rounds back to the current grid dose. If \(D_k\) is on a uniform grid, the continuous recommendation moves by \(\lambda(D^*-D_k)\). When the magnitude of this move is less than \(h/2\), nearest-dose rounding returns \(D_k\), assuming the current dose is feasible and the recommendation is not affected by a bound. Thus an iterative rule can stop changing while the continuous target remains unreached. At exact midpoints, the result depends on the tie rule.

$$ Q_G(D_k+\lambda(D^*-D_k))=D_k \quad\text{if}\quad |\lambda(D^*-D_k)|<\frac{h}{2}. $$

For \(\lambda>0\), this no-change region corresponds to a dose error smaller than \(h/(2\lambda)\). Under the linear response model, its response discrepancy from the target is smaller than \(bh/(2\lambda)\). This bound describes a possible stopping band for the iterative rule, not the error from rounding the final target once. In particular, a smaller gain can make the stopping band wider.

Worked Example: A Rounded Update Stops Short of the Continuous Target

Consider feasible doses in 5 mg/day increments, with a target dose \(D^*=52\) mg/day. Start at \(D_0=30\) mg/day and use the proportional rule with \(\lambda=0.5\), rounding each recommendation to the nearest feasible dose. Assume no dose bounds interfere.

From 30 mg/day, the continuous recommendation is \(30+0.5(52-30)=41\) mg/day. The nearest feasible dose is 40 mg/day, so \(D_1=40\). From there the recommendation is \(40+0.5(52-40)=46\) mg/day, which rounds to \(D_2=45\). The next recommendation is \(45+0.5(52-45)=48.5\) mg/day; 50 mg/day is nearer than 45 mg/day, so \(D_3=50\).

At 50 mg/day, the next recommendation is \(50+0.5(52-50)=51\) mg/day. This rounds back to 50 mg/day, because it is 1 mg/day away while 55 mg/day is 4 mg/day away. The algorithm therefore stays at 50 mg/day. The target dose is not feasible; among 50 and 55 mg/day, 50 mg/day is the nearer available dose, with a residual dose difference of 2 mg/day.

The no-change threshold is \(h/(2\lambda)=5/(2(0.5))=5\) mg/day. The 2 mg/day residual lies inside that band, as expected. If the response were \(R(D)=0.8D\), the target response would be \(0.8(52)=41.6\), whereas the response at the stopping dose would be \(0.8(50)=40\), a difference of \(-1.6\) response units. This example illustrates why the stopping tolerance should be judged in the clinically relevant response scale as well as in dose units.

When Rounding Produces a Dose Cycle

Rounding does not always make an iterative update settle at one dose. If the continuous update gain is large enough to cross the target repeatedly, quantization can turn the resulting movements into a repeating cycle. A sequence that alternates between two feasible doses does not converge to a single dose, even if the unrounded error recurrence would shrink for that gain.

Worked Example: A Two-Dose Cycle on a 5 mg/day Grid

Use feasible doses in 5 mg/day increments, target dose \(D^*=52\) mg/day, and nearest-dose rounding with \(\lambda=1.4\). Start at \(D_0=40\) mg/day. Assume the dose range includes all doses calculated here and that no measurements or other rules modify the update.

At 40 mg/day, the recommendation is \(40+1.4(52-40)=56.8\) mg/day, which rounds to 55 mg/day: it is 1.8 mg/day from 56.8, compared with 3.2 mg/day for 60. At 55 mg/day, the recommendation is \(55+1.4(52-55)=50.8\) mg/day, which rounds to 50 mg/day. At 50 mg/day, it is \(50+1.4(52-50)=52.8\) mg/day, which rounds to 55 mg/day. The sequence then repeats: 55 mg/day again produces a recommendation of 50.8 mg/day and returns to 50 mg/day.

Thus the rounded trajectory is \(40,55,50,55,50,\ldots\) mg/day. The target is between 50 and 55 mg/day, and 50 mg/day is the nearer of those two doses, but the algorithm repeatedly visits 55 mg/day as well. The unrounded error factor is \(1-\lambda=-0.4\), whose magnitude is less than 1; without quantization, the error would alternate in sign while shrinking. That fact is not enough to guarantee convergence after rounding.

Diagnostics and Practical Choices

When implementing or reviewing a titration algorithm, inspect the actual sequence of feasible doses, not only the continuous recommendations. For each update, record the current dose, the continuous recommendation, the selected feasible dose, and the predicted target metric. The following checks help reveal when quantization is materially changing the rule:

  • Verify the feasible set. List dose strengths, permitted combinations, dose bounds, and any limits on splitting or administration. A grid spacing without an origin or a complete feasible set is insufficient to determine a rounded dose.
  • Check the rounding rule. State whether doses are rounded to nearest, always down, always up, or chosen by another clinically justified rule. Specify how exact ties are handled.
  • Look for repeated doses and cycles. A dose that remains unchanged despite a nonzero target error may indicate that updates fall inside a rounding deadband. Alternation or repeated dose patterns may indicate quantization-induced cycling.
  • Evaluate the consequence in response units. Convert the difference between selected and target doses through the response or exposure model. A small dose difference is not necessarily clinically small if the response is sensitive to dose.
  • Reassess near bounds and with uncertain inputs. Clipping at a minimum or maximum dose changes the available choices. Measurement noise or uncertainty in the model can also make recommendations jump between neighboring strengths, so deterministic calculations should not be treated as guarantees about individual outcomes.

The main alternatives answer different design questions. Rounding the continuous target once by nearest dose gives a simple bound of \(h/2\) in dose units on a uniform grid, but does not ensure that every interim dose follows a prescribed stepwise titration. Rounding each iterative recommendation preserves the update structure but can create a stopping band or cycle. Directional rounding, such as always rounding toward the target, can avoid some small steps being rounded away, but may systematically bias the dose and does not by itself prevent overshoot. A finer feasible grid can reduce rounding error, but may not be available or appropriate. These choices should be evaluated against dose safety constraints and the target metric, rather than selected solely for mathematical convenience.

Quantization analysis has important limits. The simple half-step and deadband results require nearest-dose rounding on a uniform grid and, for response-error calculations, a locally linear response relationship. Nonlinear response requires evaluating the model at the actual selected doses; one cannot infer a response difference from dose spacing alone. If the response is noisy, delayed, or uncertain, observed changes may not follow the deterministic trajectory. A discrete update rule should therefore be checked in its full operating range, including boundaries and plausible measurement variation.

Common pitfall: It is incorrect to infer that a rounded iterative algorithm converges just because the continuous update converges, or to assume that a grid spacing identifies the nearest available dose. Specify the actual feasible-dose set and tie rule, then examine the implemented sequence and its predicted response.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. A target dose is 46 mg/day. Available doses are multiples of 10 mg/day from 0 mg/day upward. Which dose is selected by nearest-dose rounding, and what is the dose rounding error?
  2. On a 4 mg/day grid, an iterative rule has \(\lambda=0.5\). What is the no-change dose-error threshold under nearest rounding, away from bounds and ties?
  3. For a linear response with slope \(b=0.6\) response units per mg/day and grid spacing \(h=10\) mg/day, what is the maximum response error from rounding a feasible continuous target once to its nearest grid dose?
  4. Why does convergence of an unrounded update with \(0<\lambda<2\) not prove convergence of the rounded update?
  5. A dose recommendation lies exactly halfway between two feasible doses. What additional information is needed to determine the selected dose?

Solutions. (1) The neighboring doses are 40 and 50 mg/day, each 6 and 4 mg/day from the target, respectively. Select 50 mg/day; the error is \(+4\) mg/day. (2) The threshold is \(h/(2\lambda)=4/(2(0.5))=4\) mg/day, with exact midpoint behavior depending on the tie rule. (3) The bound is \(bh/2=0.6(10)/2=3\) response units. (4) Rounding adds a dose discrepancy at each update; it can leave a nonzero error unchanged or produce a repeating discrete cycle. (5) The tie-breaking rule, such as choosing the lower dose, must be specified.

Clinical Biostats