Why Check an Accumulation Model at Endpoints?
An accumulation model predicts how a stock changes over time from its starting value and its rate of change. For example, a model might describe the amount of water in a reservoir, inventory in a warehouse, or material stored in a tank. Comparing its predictions with measured stock values at selected times is a direct way to check whether the model is consistent with observations.
A mismatch is useful evidence, but it does not automatically reveal the cause. A prediction may be off because the starting stock was recorded incorrectly, because the modeled rate is inaccurate, or because the measurement itself contains error. Comparing several endpoints—and comparing the changes between them—helps separate these possibilities.
Model Predictions and Endpoint Discrepancies
Suppose a stock \(S(t)\) has initial value \(S_0\) at time \(t_0\), and its rate of change is modeled by \(r(t)\). The accumulation model predicts
Here \(M(t)\) is the model-predicted stock. The integral accumulates the modeled rate from the initial time to the endpoint. This formula applies when the rate is integrable over the interval and the stock changes according to that rate.
If a measurement at time \(t_i\) is \(y_i\), define its endpoint discrepancy as observed stock minus predicted stock:
A positive discrepancy means the observed stock is greater than the model prediction; a negative discrepancy means it is smaller. A discrepancy of zero means the observation agrees with the prediction at that endpoint. The sign convention matters: reversing the subtraction reverses the interpretation.
Worked Example: A Persistent Baseline Offset
A warehouse model starts with \(S_0=100\) items and has a constant modeled net rate of \(5\) items per day. Compare its predictions with measurements of \(103\) items at day \(0\), \(113\) items at day \(2\), and \(128\) items at day \(5\).
The model is \(M(t)=100+\int_0^t 5\,du=100+5t\). At day \(0\), \(M(0)=100+5(0)=100\), so \(d_0=103-100=3\) items.
At day \(2\), \(M(2)=100+5(2)=110\), so \(d_2=113-110=3\) items. At day \(5\), \(M(5)=100+5(5)=125\), so \(d_5=128-125=3\) items.
The measured increase from day \(0\) to day \(2\) is \(113-103=10\) items; the model increase is \(110-100=10\) items. From day \(2\) to day \(5\), the measured increase is \(128-113=15\) items; the model increase is \(125-110=15\) items. Thus the interval changes agree, while every absolute measurement is \(3\) items above the prediction.
As a check, the shifted-baseline model \(103+5t\) gives \(103\), \(113\), and \(128\) items at days \(0\), \(2\), and \(5\), respectively. This pattern is consistent with a baseline offset, although measurement error could also contribute.
Why Comparing Endpoints Works
The distinction between a baseline discrepancy and an accumulated rate discrepancy follows from the accumulation equation itself. Let \(X(t)\) be the actual stock and suppose its actual rate is \(q(t)\), so that \(X(t)=X_0+\int_{t_0}^{t}q(u)\,du\). Assume a measurement at \(t_i\) is \(y_i=X(t_i)+\varepsilon_i\), where \(\varepsilon_i\) is its measurement error. The model prediction is \(M(t_i)=S_0+\int_{t_0}^{t_i}r(u)\,du\). Subtracting the prediction from the measurement gives
This identity separates three contributions: the initial-stock difference \(X_0-S_0\), the accumulated difference between actual and modeled rates, and measurement error. In particular, for two endpoints \(t_i<t_j\), subtracting their discrepancies cancels the baseline difference:
So, if measurement errors are small, a change in discrepancy between endpoints indicates a net rate mismatch over the intervening interval. It does not specify exactly when that mismatch occurred or what caused it. The endpoint comparison measures the accumulated difference over the interval, not the full rate history.
Compare Changes as Well as Stock Values
For consecutive observation times \(t_i\) and \(t_j\), compare the observed stock change \(y_j-y_i\) with the model-predicted change:
The initial value does not appear in this comparison. That is why interval changes help assess the rate component separately from the baseline. If the observed and predicted changes agree, but both endpoint discrepancies are similar, the mismatch may lie in the baseline. If the changes disagree, the model's accumulated rate over that interval does not match the observed change, subject to measurement uncertainty.
Worked Example: Discrepancies That Grow Across Endpoints
A stock model has initial value \(50\) units and constant modeled rate \(6\) units per day, so \(M(t)=50+6t\). Measurements are \(54\) units at day \(0\), \(70\) units at day \(2\), and \(94\) units at day \(5\).
At day \(0\), the prediction is \(M(0)=50\), giving \(d_0=54-50=4\) units. At day \(2\), \(M(2)=50+6(2)=62\), giving \(d_2=70-62=8\) units. At day \(5\), \(M(5)=50+6(5)=80\), giving \(d_5=94-80=14\) units.
From day \(0\) to day \(2\), the observed increase is \(70-54=16\) units, while the modeled increase is \(62-50=12\) units. The observed increase is \(4\) units greater. From day \(2\) to day \(5\), the observed increase is \(94-70=24\) units, while the modeled increase is \(80-62=18\) units. The observed increase is \(6\) units greater.
The discrepancies increase from \(4\) to \(8\) to \(14\) units. The measured changes are consistent with a rate \(2\) units per day greater than the modeled rate, together with a \(4\)-unit baseline difference: \(54+8(2)=70\) and \(54+8(5)=94\). This is a consistency check, not proof that these are the actual causes; measurement errors or a different varying rate could produce the same endpoints.
Check an Interval Before Naming the Cause
A discrepancy at a later endpoint may arise from a mismatch anywhere in the interval leading to that endpoint. Comparing successive endpoints can narrow down where the accumulated disagreement occurs. It cannot, by itself, reveal the detailed behavior of the rate between observations.
Worked Example: A Mismatch Concentrated in a Later Interval
A stock starts at \(200\) units, and its modeled rate is \(r(t)=10-2t\) units per day. Measurements are \(200\) units at day \(0\), \(216\) units at day \(2\), and \(220\) units at day \(4\).
An antiderivative of \(10-2t\) is \(10t-t^2\). Therefore, the prediction at day \(2\) is \(M(2)=200+[(10(2)-2^2)-(10(0)-0^2)]=200+(20-4)=216\) units.
At day \(4\), the prediction is \(M(4)=200+[(10(4)-4^2)-(10(0)-0^2)]=200+(40-16)=224\) units. The discrepancies are \(d_0=200-200=0\), \(d_2=216-216=0\), and \(d_4=220-224=-4\) units.
From day \(0\) to day \(2\), the observed increase is \(216-200=16\) units. The modeled increase is \((20-4)=16\) units, so they agree. From day \(2\) to day \(4\), the observed increase is \(220-216=4\) units. The modeled increase is \((40-16)-(20-4)=24-16=8\) units, a difference of \(-4\) units.
Thus the observations agree with the model through day \(2\), then fall \(4\) units short of its predicted increase by day \(4\). As a check, the model's predicted endpoint is \(224\) units, and the measured endpoint is \(220\) units, giving the same \(-4\)-unit discrepancy. The data identify an interval of disagreement, not the precise time or mechanism of a rate error.
A Common Mistake: Calling Every Difference a Rate Error
Suppose the first example is interpreted by saying, “The day-zero measurement is \(3\) items above prediction, so the modeled rate must be \(3\) items per day too low.” This confuses stock, measured in items, with rate, measured in items per day. It also ignores the later endpoint comparisons.
If the rate really were \(3\) items per day too low while the baseline remained \(100\), this interpretation would predict \(100+(5+3)(2)=116\) items at day \(2\), not the observed \(113\). The correction is to compare interval changes: the model and measurements both increase by \(10\) items over the first two days and \(15\) items over the next three days. The repeated \(3\)-item stock difference is consistent with a baseline offset, not evidence by itself of a rate error.
Limits of Endpoint Evidence
Endpoint measurements are valuable checks, but conclusions must remain proportionate to the evidence. A few endpoints do not determine the rate at every moment. A rate could be too high during part of an interval and too low during another part, with the two errors canceling in the total accumulation. Consequently, matching endpoints do not prove that the model's rate was accurate throughout the interval.
Measurement uncertainty also matters. If each measurement has a known tolerance, discrepancies smaller than that tolerance may not indicate a meaningful mismatch. When measurement errors are unknown, it is safest to describe the observed pattern and the explanations it is consistent with, rather than assert a unique cause.
Use the initial stock and the integral of the modeled rate to calculate stock at each observed time.
Calculate each discrepancy as observed stock minus predicted stock, keeping the units and sign.
Compare observed stock changes with the modeled rate accumulated over the same intervals.
Describe whether the evidence is consistent with baseline or rate mismatch, while allowing for measurement error and non-unique explanations.
Check Your Understanding
Use the definitions and comparisons in this tutorial to answer the following questions.
- If a measurement is \(7\) units below its model prediction, what is the endpoint discrepancy using the convention in this tutorial?
- Why does comparing stock changes between two endpoints remove the effect of a constant baseline offset?
- A model and measurements have the same increase over each observed interval, but every measured stock is \(5\) units higher. What kind of error is this pattern consistent with?
- What can a changing endpoint discrepancy indicate, and why does it not identify the exact rate history?
- Can matching predictions at two endpoints prove that the modeled rate was correct throughout the interval? Explain briefly.