The Central Idea

Suppose a forward model predicts

\[ \mathbf{y}_{\text{true}} = F(\boldsymbol{\theta}). \]

The actual observation is often better represented as

\[ \boxed{ \mathbf{y}_{\text{obs}} = \mathbf{y}_{\text{true}} + \boldsymbol{\varepsilon} } \]

where

\[ \boldsymbol{\varepsilon} \]

represents measurement error or observational noise.

Measurement Error Is Part of the Model

A complete inverse model has at least two components:

  1. a forward model describing the underlying system;
  2. an observation model describing how measurements differ from the underlying system.

Conceptually,

\[ \boldsymbol{\theta} \longrightarrow F(\boldsymbol{\theta}) \longrightarrow \mathbf{y}_{\text{obs}}. \]

Ignoring the second stage can produce misleading parameter estimates and uncertainty statements.

Random Error

Random error varies unpredictably from one measurement to another.

Repeated measurements of the same quantity may produce values scattered around a central value.

A common model is

\[ \varepsilon \sim \text{distribution with mean }0. \]

Random error contributes uncertainty but does not necessarily introduce systematic bias.

Systematic Error

Systematic error shifts measurements consistently in one direction or according to a repeatable pattern.

Examples include:

Bias

The bias of a measurement error is

\[ \boxed{ \operatorname{Bias} = \mathbb{E}[\varepsilon] }. \]

If

\[ \mathbb{E}[\varepsilon]=0, \]

the measurement error is unbiased in expectation.

If

\[ \mathbb{E}[\varepsilon]\ne0, \]

the observation process contains systematic bias.

Variance

The variance measures the spread of random measurement error:

\[ \boxed{ \operatorname{Var}(\varepsilon) = \mathbb{E} \left[ (\varepsilon-\mathbb{E}[\varepsilon])^2 \right] }. \]

Large variance means repeated measurements are widely scattered.

Small variance means the measurements are relatively precise.

Standard Deviation

The standard deviation is

\[ \sigma = \sqrt{ \operatorname{Var}(\varepsilon) }. \]

Because \(\sigma\) has the same units as the measured quantity, it is often easier to interpret than variance.

Accuracy Versus Precision

Accuracy concerns closeness to the true value.

Precision concerns repeatability or spread.

Measurements can be:

A Precise but Biased Sensor

Suppose a true value is

\[ 10. \]

A sensor repeatedly reports values near

\[ 10.8. \]

The readings may have very small variance while still being systematically wrong.

Precision alone does not guarantee accuracy.

Gaussian Noise

A common model for additive measurement error is

\[ \varepsilon \sim N(0,\sigma^2). \]

Then

\[ y_{\text{obs}} = y_{\text{true}} + \varepsilon. \]

This model is appropriate when many small independent effects combine to produce approximately symmetric noise.

Gaussian Likelihood

If the forward prediction is

\[ \mu_i(\boldsymbol{\theta}), \]

and

\[ y_i \sim N( \mu_i(\boldsymbol{\theta}), \sigma_i^2 ), \]

then the negative log-likelihood contains

\[ \frac{ \left[ y_i-\mu_i(\boldsymbol{\theta}) \right]^2 }{ 2\sigma_i^2 }. \]

Measurement variance therefore determines the natural weighting of residuals.

Homoscedastic Noise

Noise is homoscedastic when all observations share the same variance:

\[ \operatorname{Var}(\varepsilon_i) = \sigma^2. \]

Under independent Gaussian errors, ordinary least squares is then a natural estimator.

Heteroscedastic Noise

Noise is heteroscedastic when the measurement variance differs across observations:

\[ \operatorname{Var}(\varepsilon_i) = \sigma_i^2. \]

Variance may depend on:

Inverse-Variance Weighting

If observation \(i\) has known variance

\[ \sigma_i^2, \]

a common residual weight is

\[ w_i = \frac{1}{\sigma_i^2}. \]

The weighted least-squares objective becomes

\[ \boxed{ \Phi(\boldsymbol{\theta}) = \sum_i \frac{ \left[ y_i-F_i(\boldsymbol{\theta}) \right]^2 }{ \sigma_i^2 } }. \]

Why More Precise Measurements Receive More Weight

If one measurement has small uncertainty, a given residual is more surprising than the same residual from a noisy measurement.

Inverse-variance weighting formalizes this idea.

A residual of \(1\) unit may be severe if

\[ \sigma=0.1, \]

but modest if

\[ \sigma=10. \]

Covariance

Errors from different observations may vary together.

The covariance between errors \(i\) and \(j\) is

\[ \boxed{ \operatorname{Cov}( \varepsilon_i,\varepsilon_j ) = \mathbb{E} \left[ ( \varepsilon_i-\mu_i ) ( \varepsilon_j-\mu_j ) \right] }. \]

Nonzero covariance means the errors are statistically related.

The Covariance Matrix

For a noise vector

\[ \boldsymbol{\varepsilon} = (\varepsilon_1,\ldots,\varepsilon_n)^T, \]

the covariance matrix is

\[ \boxed{ \Sigma = \mathbb{E} \left[ ( \boldsymbol{\varepsilon} - \mathbb{E}[\boldsymbol{\varepsilon}] ) ( \boldsymbol{\varepsilon} - \mathbb{E}[\boldsymbol{\varepsilon}] )^T \right] }. \]

Diagonal and Off-Diagonal Terms

The diagonal entries

\[ \Sigma_{ii} \]

are error variances.

The off-diagonal entries

\[ \Sigma_{ij} \]

are error covariances.

A diagonal covariance matrix represents uncorrelated errors.

Correlated Measurement Error

Correlation can arise when several measurements share:

Treating correlated measurements as independent can overstate the amount of information in the data.

Generalized Least Squares

If

\[ \operatorname{Cov}( \boldsymbol{\varepsilon} ) = \Sigma, \]

then a natural quadratic residual measure is

\[ \boxed{ \Phi = \mathbf{r}^T \Sigma^{-1} \mathbf{r} }. \]

This accounts for both unequal variances and correlations.

Whitening

If

\[ \Sigma=LL^T, \]

define transformed residuals

\[ \tilde{\mathbf{r}} = L^{-1}\mathbf{r}. \]

Then

\[ \mathbf{r}^T \Sigma^{-1} \mathbf{r} = \tilde{\mathbf{r}}^T \tilde{\mathbf{r}}. \]

The transformation converts correlated, unequal-variance residuals into standardized residual coordinates.

Correlation

Covariance depends on measurement scale.

The correlation coefficient standardizes covariance:

\[ \rho_{ij} = \frac{ \operatorname{Cov}( \varepsilon_i,\varepsilon_j ) }{ \sigma_i\sigma_j }. \]

It lies between

\[ -1 \]

and

\[ 1. \]

Temporal Correlation

Time-series measurements often have residuals that are correlated across nearby times.

For example,

\[ \operatorname{Cov}( \varepsilon_t, \varepsilon_{t+\Delta} ) \]

may decrease gradually as the lag \(\Delta\) increases.

Ignoring this structure can lead to overly optimistic uncertainty estimates.

Spatial Correlation

Measurements taken at nearby locations may share environmental or physical disturbances.

A covariance model may therefore depend on distance:

\[ \operatorname{Cov}( \varepsilon(\mathbf{x}), \varepsilon(\mathbf{x}') ) = C( \|\mathbf{x}-\mathbf{x}'\| ). \]

Spatially correlated error is common in imaging, geophysics, climate data, and remote sensing.

Calibration Error

A sensor may follow

\[ y_{\text{obs}} = a + b\,y_{\text{true}} + \varepsilon. \]

Here:

Calibration parameters may need to be estimated alongside the physical parameters.

Additive Bias

An additive sensor bias can be modeled as

\[ y_{\text{obs}} = y_{\text{true}} + b + \varepsilon. \]

If \(b\) is ignored, the inverse method may incorrectly force physical model parameters to compensate for the sensor offset.

Multiplicative Error

Some instruments exhibit proportional error:

\[ y_{\text{obs}} = y_{\text{true}} (1+\varepsilon). \]

The absolute error then grows naturally with signal magnitude.

Relative-error or logarithmic observation models may be more appropriate than constant-variance additive noise.

Lognormal Error

For positive quantities with multiplicative variability, one may model

\[ \log y_{\text{obs}} = \log y_{\text{true}} + \varepsilon, \]

where

\[ \varepsilon \sim N(0,\sigma^2). \]

Then the original measurement follows a lognormal-type distribution.

Count Noise

Count measurements may naturally follow a Poisson model:

\[ Y \sim \operatorname{Poisson}(\lambda). \]

Then

\[ \mathbb{E}[Y]=\lambda \]

and

\[ \operatorname{Var}(Y)=\lambda. \]

The variance is therefore tied to the expected signal.

Quantization Error

Digital instruments record values on a finite grid.

If the measurement resolution is

\[ \Delta, \]

then continuous values are rounded or truncated to nearby representable values.

This introduces quantization error even when the underlying sensor is otherwise noise-free.

Finite Resolution

Some instruments do not measure a point value but average over a spatial, temporal, or spectral window.

The observation may therefore be better modeled as

\[ y_i = \int w_i(s) f(s) \,ds + \varepsilon_i. \]

Instrument resolution is then part of the observation operator, not merely an additive noise term.

Timing Error

Suppose a measurement is intended at time

\[ t_i \]

but actually occurs at

\[ t_i+\delta t_i. \]

Then the observation may be

\[ y_i = f( t_i+\delta t_i ) + \varepsilon_i. \]

In rapidly changing systems, timing uncertainty can produce substantial apparent measurement error.

Error in the Independent Variable

Standard least squares often assumes predictors are known exactly.

If both \(x\) and \(y\) contain measurement error, the problem changes.

For example,

\[ x_{\text{obs}} = x_{\text{true}} + \delta, \] \[ y_{\text{obs}} = f( x_{\text{true}}, \boldsymbol{\theta} ) + \varepsilon. \]

Ignoring predictor error can bias fitted parameters.

Errors-in-Variables Models

When explanatory variables are noisy, one may model their latent true values explicitly.

This leads to errors-in-variables formulations or total least-squares methods in suitable linear settings.

Outliers

An outlier is an observation unusually far from the pattern predicted by the model.

Possible causes include:

Outliers should be investigated rather than automatically deleted.

Heavy-Tailed Noise

Gaussian models assign rapidly decreasing probability to extreme residuals.

If large deviations occur more often than Gaussian theory predicts, a heavy-tailed distribution may be more appropriate.

Heavy-tailed likelihoods reduce the influence of extreme observations.

Robust Estimation

Robust methods reduce sensitivity to a small number of extreme residuals.

Instead of penalizing every residual quadratically, a robust loss may grow more slowly when

\[ |r| \]

becomes large.

This can improve stability when contamination or outliers are plausible.

Detection Limits

Instruments may not resolve values below a threshold.

A reported measurement may therefore mean only

\[ y rather than giving an exact numerical value.

Treating such a censored observation as if it were exactly zero can bias parameter estimates.

Censoring

If the instrument reports only that

\[ Y then the likelihood contribution is

\[ P( Y not a point-density evaluation.

Likelihood-based methods can incorporate this information directly.

Missing Measurements

Missing data should not automatically be treated as zero or replaced by an arbitrary value.

The correct treatment depends on why the data are missing and whether the missingness mechanism is related to the underlying system.

Preprocessing Error

Data often undergo processing before estimation:

These operations can alter both the mean structure and covariance of the measurement errors.

Smoothing Creates Correlation

Suppose raw measurement errors are independent.

If neighboring measurements are averaged, the processed values share some of the same raw observations.

The resulting errors become correlated even if the original noise was not.

Replication

Repeating a measurement under nominally identical conditions helps estimate random variability.

If measurements are

\[ y_1,\ldots,y_m, \]

the sample mean is

\[ \bar y = \frac{1}{m} \sum_{j=1}^{m}y_j. \]

Replicates provide direct information about repeatability.

Variance of an Average

If independent replicated measurements have common variance

\[ \sigma^2, \]

then

\[ \operatorname{Var}(\bar y) = \frac{\sigma^2}{m}. \]

Therefore, averaging independent replicates reduces random uncertainty.

Replication Does Not Remove Bias

If every measurement contains the same additive bias

\[ b, \]

then averaging many measurements does not eliminate it.

Random noise decreases through averaging, but systematic error remains.

Error Propagation Through a Function

Suppose

\[ z=g(x) \]

and \(x\) has small uncertainty.

Linearization gives

\[ \Delta z \approx g'(x)\Delta x. \]

Therefore,

\[ \boxed{ \operatorname{Var}(z) \approx [g'(x)]^2 \operatorname{Var}(x) }. \]

Multivariable Error Propagation

Let

\[ \mathbf{z} = g(\mathbf{x}). \]

If \(\mathbf{x}\) has covariance matrix

\[ \Sigma_x, \]

then local linearization gives

\[ \boxed{ \Sigma_z \approx J_g \Sigma_x J_g^T } \]

where \(J_g\) is the Jacobian of \(g\).

Inverse Error Propagation

In inverse problems, the mapping from data to parameters may be highly sensitive.

If

\[ \boldsymbol{\theta} = G(\mathbf{y}), \]

then approximately

\[ \Sigma_{\theta} \approx J_G \Sigma_y J_G^T. \]

Large entries in \(J_G\) indicate strong amplification of measurement uncertainty.

Linear Inverse Error Propagation

If

\[ \hat{\boldsymbol{\theta}} = A\mathbf{y}, \]

then

\[ \boxed{ \operatorname{Cov} ( \hat{\boldsymbol{\theta}} ) = A \Sigma_y A^T }. \]

This formula makes the effect of measurement covariance explicit.

Noise Amplification and Ill-Posedness

If the inverse operator \(A\) has very large gain in some directions, then modest measurement covariance can produce enormous parameter covariance.

This is the uncertainty form of ill-posedness.

Standardized Residuals

If a measurement has standard deviation

\[ \sigma_i, \]

a standardized residual is

\[ \boxed{ r_i^{\text{std}} = \frac{ y_i-F_i(\boldsymbol{\theta}) }{ \sigma_i } }. \]

This measures model-data mismatch in units of expected measurement uncertainty.

Interpreting Standardized Residuals

A raw residual of \(2\) units cannot be judged without knowing the noise scale.

If

\[ \sigma=0.1, \]

then \(2\) units is extremely large relative to the expected noise.

If

\[ \sigma=10, \]

the same raw residual is modest.

Residual Diagnostics

After fitting, residuals should be inspected for:

Residual structure often reveals an inadequate observation model.

Residual Mean

If residuals have a persistent nonzero average,

\[ \bar r \ne 0, \]

this may indicate:

Residual Variance Patterns

If residual spread increases with predicted magnitude, constant-variance Gaussian noise may be inappropriate.

A heteroscedastic or multiplicative error model may better represent the data-generating process.

Autocorrelation Diagnostics

For time-indexed residuals, examine whether

\[ r_t \]

is correlated with

\[ r_{t+k}. \]

Persistent autocorrelation can indicate unmodeled dynamics or correlated measurement noise.

Noise Model Misspecification

If the assumed noise model is wrong, several consequences can follow:

Estimating Noise Parameters

Noise parameters such as

\[ \sigma, \]

correlation length, or dispersion may be unknown.

They can sometimes be estimated jointly with the physical parameters.

This creates a larger inverse problem in which both signal and uncertainty structure are inferred.

Separating Process Noise and Measurement Noise

In dynamic systems, uncertainty may enter both the system evolution and the measurement process.

A state model may include

\[ \mathbf{x}_{k+1} = f( \mathbf{x}_k ) + \mathbf{w}_k, \]

while observations satisfy

\[ \mathbf{y}_k = h( \mathbf{x}_k ) + \mathbf{v}_k. \]

Here \(\mathbf{w}_k\) is process noise and \(\mathbf{v}_k\) is measurement noise.

Why the Distinction Matters

Measurement noise describes imperfect observation of a system state.

Process noise describes genuine uncertainty in how the system evolves.

Treating one as the other can distort state and parameter estimates.

Instrument Specifications

Real instruments may report uncertainty through:

These specifications should be translated carefully into the mathematical observation model.

Error Budgets

An error budget separates different uncertainty sources.

For example,

\[ \sigma_{\text{total}}^2 = \sigma_{\text{sensor}}^2 + \sigma_{\text{calibration}}^2 + \sigma_{\text{sampling}}^2 \]

when independent components can reasonably be added in variance.

Error budgets help identify which source dominates total uncertainty.

Repeated Calibration

Calibration uncertainty may change over time because instruments drift.

Repeated calibration allows analysts to distinguish:

Data Quality Versus Data Quantity

More measurements are not always equivalent to more information.

A large number of highly correlated, biased, or low-resolution measurements may be less useful than a smaller number of precise independent measurements.

Experimental Design for Noise Reduction

Measurement uncertainty can often be reduced through:

Noise and Regularization

Regularization strength should depend partly on measurement uncertainty.

If noise is large, fitting every detail of the observations is undesirable.

If noise is very small, stronger data fidelity may be justified.

Discrepancy Principle

If the expected noise magnitude is known, one may choose regularization so that

\[ \|\mathbf{y}-F(\hat{\boldsymbol{\theta}})\| \]

is comparable to the expected observational error.

Residuals far below the noise level may indicate overfitting.

Noise and Bayesian Estimation

In Bayesian inference, the noise model enters the likelihood:

\[ p( \mathbf{y} \mid \boldsymbol{\theta} ). \]

Larger assumed noise generally produces a broader likelihood and weaker data-driven posterior concentration.

Smaller assumed noise gives the data greater influence.

Underestimating Noise

If the analysis assumes

\[ \sigma_{\text{assumed}} < \sigma_{\text{true}}, \]

then residuals appear more informative than they really are.

Parameter uncertainty may be reported too narrowly.

Overestimating Noise

If

\[ \sigma_{\text{assumed}} > \sigma_{\text{true}}, \]

the analysis may underuse the information in the observations.

Parameter uncertainty can become unnecessarily broad.

Common Error: Assuming All Error Is Random

Repeated averaging can reduce random noise but cannot remove a persistent calibration bias.

Systematic and random errors should be modeled separately when possible.

Common Error: Assuming Independence

If observations share a sensor, preprocessing step, or environmental disturbance, their errors may be correlated.

Treating them as independent can overstate effective sample size.

Common Error: Using Equal Weights Automatically

Equal weighting is appropriate only when observations have comparable error scales and the objective is defined accordingly.

Large differences in precision should generally influence residual weighting.

Common Error: Deleting Outliers Without Investigation

An extreme observation may be erroneous, but it may also reflect genuine system behavior that the model fails to explain.

Outliers should be investigated before removal.

Common Error: Ignoring Error in Predictors

Standard regression assumptions may fail when explanatory variables are also measured with uncertainty.

Predictor error can bias slope and parameter estimates.

Common Error: Treating Instrument Resolution as Pure Noise

Finite resolution may represent averaging or blurring by the instrument.

That effect belongs in the forward observation operator, not only in an additive error term.

Common Error: Ignoring Preprocessing Effects

Smoothing, interpolation, and normalization can change error variance and induce correlation.

The statistical model should reflect the data actually used in estimation, not merely the raw instrument output.

A Reliable Measurement-Error Workflow

  1. Identify the physical quantity actually observed.
  2. Separate forward-model output from observation-model behavior.
  3. List random and systematic error sources.
  4. Estimate or justify measurement variance.
  5. Check whether variance changes with condition or signal level.
  6. Assess temporal, spatial, or instrument-induced correlation.
  7. Model calibration and predictor uncertainty when relevant.
  8. Translate uncertainty into residual weights or a likelihood.
  9. Inspect standardized residuals after fitting.
  10. Test for bias, heteroscedasticity, and correlation.
  11. Propagate measurement uncertainty into parameter uncertainty.
  12. Use replication or calibration experiments where possible.
  13. Document preprocessing and uncertainty assumptions.

Practice 1: Bias

If a measurement error satisfies

\[ \mathbb{E}[\varepsilon]=0.4, \]

is the measurement process unbiased?

Practice 2: Variance and Standard Deviation

If

\[ \operatorname{Var}(\varepsilon)=9, \]

what is the standard deviation?

Practice 3: Inverse-Variance Weighting

Two observations have standard deviations

\[ \sigma_1=1 \]

and

\[ \sigma_2=2. \]

What are their inverse-variance weights?

Practice 4: Replication

Independent repeated measurements each have variance

\[ 16. \]

If \(4\) such measurements are averaged, what is the variance of the average?

Practice 5: Systematic Error

A sensor always reads \(0.7\) units too high.

What type of measurement error is this?

Practice 6: Correlation

Why can smoothing independent raw measurements create correlated processed errors?

Practice 7: Standardized Residual

Suppose

\[ y=15, \qquad \hat y=13, \qquad \sigma=4. \]

Compute the standardized residual.

Practice 8: Error Propagation

Suppose

\[ z=3x \]

and

\[ \operatorname{Var}(x)=2. \]

Approximate \(\operatorname{Var}(z)\).

Practice 9: Multiplicative Noise

If the absolute measurement spread grows in proportion to the signal level, which may be more appropriate: constant-variance additive noise or a multiplicative error model?

Practice 10: Bias and Averaging

Will averaging many repeated measurements eliminate a fixed calibration offset?

Practice Answers

Practice 1

No. The measurement process is biased because

\[ \boxed{ \mathbb{E}[\varepsilon]=0.4\ne0 }. \]

Practice 2

\[ \sigma = \sqrt{9} = \boxed{3}. \]

Practice 3

The inverse-variance weights are

\[ w_1 = \frac{1}{1^2} = \boxed{1}, \]

and

\[ w_2 = \frac{1}{2^2} = \boxed{\frac14}. \]

Practice 4

For \(m=4\) independent measurements,

\[ \operatorname{Var}(\bar y) = \frac{16}{4} = \boxed{4}. \]

Practice 5

This is a systematic additive bias.

Practice 6

Neighboring processed values can share some of the same raw observations. Because the same random errors contribute to several processed points, their resulting errors become statistically correlated.

Practice 7

\[ r^{\text{std}} = \frac{15-13}{4} = \boxed{0.5}. \]

Practice 8

Since

\[ \frac{dz}{dx}=3, \]

local variance propagation gives

\[ \operatorname{Var}(z) \approx 3^2 \operatorname{Var}(x) = 9(2) = \boxed{18}. \]

Practice 9

A multiplicative error model is generally more appropriate because the noise scale grows with the signal.

Practice 10

No. Averaging reduces random variance but does not remove a fixed systematic calibration offset.

What to Remember

Measurement error determines how strongly observations should influence an inverse estimate and how much uncertainty remains afterward.

The central principle is:

\[ \boxed{ \text{an inverse estimate can be no more trustworthy than the uncertainty model attached to its data} }. \]

Where We Go Next

Once measurement uncertainty is represented explicitly, the next question is how to decide among competing model structures that may explain the same observations.

The next lesson, Model Selection, develops methods for comparing candidate models while balancing fit, complexity, predictive performance, and uncertainty.