Identifiability from linear observations (source code)

= Identifiability from linear observations

\b[Identifiability from linear observations] asks which features of an unknown vector $\theta\in\mathbb R^p$ are fixed by observing $y=A\theta$. Two parameter vectors are observationally equivalent exactly when their difference lies in the <kernel of a linear map> $A$:
$$
A\theta_1=A\theta_2\quad\Longleftrightarrow\quad
\theta_1-\theta_2\in\ker A.
$$
Thus the full vector is identifiable precisely when $\ker A=\{0\}$, equivalently when $A$ has full column <matrix rank>. Repeating the same rows can improve precision in a noisy measurement without increasing this rank.

A linear target $g(\theta)=c^{\mathsf T}\theta$ requires less information. It is identifiable if and only if $c$ is in the <row space> of $A$. If $c=A^{\mathsf T}d$, then $g(\theta)=d^{\mathsf T}y$, so the target is recovered directly. Conversely, an identifiable target must vanish on every observationally invisible difference: $c^{\mathsf T}h=0$ for every $h\in\ker A$. Orthogonality to the kernel is exactly row-space membership. This is a concrete instance of <quantity-specific identifiability>; the primary identification literature also treats identification of individual parameter components by rank conditions.<Identification in Parametric Models by Thomas J. Rothenberg>[§5, printed p.586, Theorem8]{ref}

For an original example, suppose
$$
y_1=\theta_1+\theta_2=5,\qquad
y_2=\theta_2+\theta_3=7.
$$
Both $(1,4,3)$ and $(2,3,4)$ fit perfectly. Their difference is a multiple of $(1,-1,1)$, which spans the kernel. The target $\theta_1-\theta_3=y_1-y_2=-2$ is identified, whereas $\theta_1+\theta_2+\theta_3$ takes values 8 and 9 and is not. A new independent measurement of $\theta_2$ restores full identification; another copy of $y_1$ does not.

The same distinction applies to a <statistical model> $Y\sim N(A\theta,\Sigma)$ with known positive-definite <covariance matrix> $\Sigma$: two parameters give the same observable law precisely when their means agree. A <Moore-Penrose pseudoinverse> selects one solution when the system is underdetermined, but that numerical convention does not show that the selected vector is the unique physical parameter. Additional restrictions must be justified as assumptions, not silently supplied by a solver.