Why the Input Is Part of the Identifiability Question
A signaling model may contain parameters that are distinguishable in principle, yet a particular experiment may provide no way to tell their values apart. This can happen even when the output is observed continuously, without measurement error, and at arbitrarily fine time resolution. The reason is that identifiability is assessed for a specified model and a specified input and initial state. If the input does not engage a modeled process, or if the experiment only observes a steady state, some parameter effects may be absent from the output trajectory.
To make this precise, consider a simple signaling state \(x(t)\), such as an active signaling pool, driven by an input \(u(t)\). The state increases in response to input and is removed at a rate proportional to its current level. The examples below use this model to separate exact ambiguity from weaker warning signs based on sensitivities.
These assumptions are deliberately favorable. In practice, observations are finite and noisy, but here the question is whether the exact output trajectory uniquely determines a parameter or parameter combination at all. A failure under ideal observation cannot be repaired merely by collecting the same output more precisely under the same conditions.
What Zero Input Leaves Unobserved
With zero input throughout the observation interval, \(u(t)=0\), the model becomes \(\dot{x}=-\gamma x\). Its solution is \(x(t)=x_0e^{-\gamma t}\). The input-dependent parameter \(\beta\) has disappeared from the output equation. If \(x_0\) is known and positive, the decay trajectory can contain information about \(\gamma\), but no observation of this trajectory can distinguish values of \(\beta\).
Worked Example: A Zero Input Hides the Activation Gain
Suppose a signaling state begins at \(x_0=0.60\) signal units and is observed for 10 minutes with \(u(t)=0\). The decay parameter is \(\gamma=0.40\ \mathrm{min}^{-1}\), but the activation gain \(\beta\) is unknown.
Under zero input, the predicted trajectory is \( y(t)=x(t)=0.60e^{-0.40t} \) signal units. For example, at \(t=2\) minutes, \( y(2)=0.60e^{-0.80}=0.2696 \) signal units.
Changing \(\beta\) from \(0.20\) to \(2.00\) signal units per minute leaves this entire trajectory unchanged, since \(\beta u(t)=0\) at every time. Thus \(\beta\) is exactly non-identifiable for this experiment, even with ideal observations. The nonzero known initial state allows the decay rate to affect the trajectory; the absence of input prevents the activation gain from affecting it.
The initial state also matters. If \(x_0=0\) and \(u(t)=0\), then \(y(t)=0\) for every time, for every admissible \(\beta\) and \(\gamma\). In that case neither parameter affects the observed output. If the initial state is unknown, conclusions must account for that additional uncertainty; known and unknown initial states are distinct experimental assumptions, not interchangeable conveniences.
Constant Input: A Transient Can Still Be Informative
A constant, nonzero input does not automatically make parameters indistinguishable. With \(u(t)=c>0\), the model solution is \(x(t)=\frac{\beta c}{\gamma}+(x_0-\frac{\beta c}{\gamma})e^{-\gamma t}\). When the state changes over time, the trajectory contains a transient as well as an equilibrium level. Under the model assumptions, the transient rate identifies \(\gamma\); the input-dependent level then identifies \(\beta\), since \(c\) is known and nonzero.
Worked Example: Recovering Parameters From a Constant-Input Transient
Suppose a constant input \(c=2\) is applied, the known initial state is \(x_0=0.20\) signal units, and the model parameters are \(\beta=0.30\) signal units per minute and \(\gamma=0.50\ \mathrm{min}^{-1}\). The equilibrium level implied by the model is \( \frac{\beta c}{\gamma}=\frac{(0.30)(2)}{0.50}=1.20 \) signal units.
The resulting trajectory is \( x(t)=1.20+(0.20-1.20)e^{-0.50t}=1.20-e^{-0.50t} \) signal units. At 1 minute, \(x(1)=1.20-e^{-0.50}=0.5935\) signal units; at 2 minutes, \(x(2)=1.20-e^{-1.00}=0.8321\) signal units. The state is changing, so this is not a steady-state-only observation.
To see how the two parameters are distinguished, write the differential equation as \( \dot{x}=\beta c-\gamma x \). At two times \(t_1,t_2\) where \(x(t_1)\ne x(t_2)\), subtract the equations: \( \dot{x}(t_1)-\dot{x}(t_2)=-\gamma[x(t_1)-x(t_2)] \). Thus the ideal trajectory determines \( \gamma=-\frac{\dot{x}(t_1)-\dot{x}(t_2)}{x(t_1)-x(t_2)}=0.50\ \mathrm{min}^{-1} \). Then \( \beta=[\dot{x}(t)+\gamma x(t)]/c \), giving \(\beta=0.30\) signal units per minute. The distinct state values ensure the subtraction contains information about \(\gamma\), and \(c\ne0\) permits recovery of \(\beta\).
This argument depends on the transient being present and on the input magnitude being known. If \(c\) were unknown as well, the equation would contain the product \(\beta c\), so the stated recovery of \(\beta\) would not follow without further information. Similarly, if the output were observed only after the transient had ended, the data would not contain the changing state values used in this argument.
Steady-State Observation Can Hide a Parameter Direction
For a constant input \(c>0\), steady state satisfies \(0=\beta c-\gamma x_{\mathrm{ss}}\), hence \(x_{\mathrm{ss}}=\beta c/\gamma\). A steady-state output therefore determines the ratio \(\beta/\gamma\), given known \(c\), but does not in general determine \(\beta\) and \(\gamma\) separately. This is an exact ambiguity: every admissible pair with the same ratio produces the same constant output when the initial state is already at that equilibrium.
Worked Example: Distinct Rate Pairs Give the Same Steady Signal
Let the known constant input be \(c=1.5\), and let the known initial state be \(x_0=0.80\) signal units. Consider two parameter pairs: \((\beta,\gamma)=(0.40\ \mathrm{signal\ units/min},\,0.75\ \mathrm{min}^{-1})\) and \((\beta,\gamma)=(0.80\ \mathrm{signal\ units/min},\,1.50\ \mathrm{min}^{-1})\).
For the first pair, the equilibrium is \(x_{\mathrm{ss}}=(0.40)(1.5)/0.75=0.80\) signal units. For the second, it is \(x_{\mathrm{ss}}=(0.80)(1.5)/1.50=0.80\) signal units. Since \(x_0=x_{\mathrm{ss}}\) in both cases, the solution is \(x(t)=0.80\) signal units at every observation time for either parameter pair.
The complete ideal outputs are identical although both parameters differ. The data identify the combination \(\beta c/\gamma\), or equivalently \(\beta/\gamma\) because \(c\) is known, but not the two rates separately. Continuous, error-free observation does not help here: there is no transient in the output to reveal its rate.
Derivation: Exact Equivalence Versus a Sensitivity Warning
The steady-state ambiguity follows directly from the solution. If \(x_0=\beta c/\gamma\), then the coefficient of the exponential term is zero:
Now multiply both \(\beta\) and \(\gamma\) by the same positive constant \(q\). The ratio \(\beta c/\gamma\) remains unchanged, and the initial state remains equal to that equilibrium. Thus the output remains exactly constant at the same value for every \(q>0\). This constructs a continuum of distinct parameter pairs with identical outputs under the stated experiment, establishing non-identifiability there.
Sensitivity analysis is useful for finding possible trouble, but it is not a substitute for this exact check. For example, let a scalar output be \(y(\theta)=\theta^3\). At \(\theta=0\), its derivative is zero, yet the map is one-to-one: if \(\theta_1^3=\theta_2^3\), then \(\theta_1=\theta_2\). A zero first derivative in a parameter direction therefore establishes first-order insensitivity at that point; it may indicate local non-identifiability, but does not by itself establish that distinct nearby parameter values give identical outputs. To establish non-identifiability, check equality of complete outputs for distinct values, or apply an appropriate identifiability result under its assumptions.
For a dynamical model, a sensitivity direction is a direction in parameter space along which the first derivative of the predicted output with respect to parameter changes vanishes. Rank deficiencies in local sensitivity calculations can flag directions for closer examination. But the conclusion concerns exact identifiability only if exact equality or a valid structural argument establishes it. Conversely, a nonzero sensitivity may still be very small, making parameters difficult to estimate in finite noisy data; that is a practical-information issue, distinct from exact structural non-identifiability.
Diagnostics, Alternatives, and Limitations
When an input appears uninformative, analyze the model under the actual input rather than relying on the model's general form. Substitute the input into the equations, solve or simplify the resulting system, and ask whether distinct parameter values can yield the same observed output. In a signaling experiment, useful warning signs include a modeled input that is identically zero, an output that remains at equilibrium throughout observation, or a parameter that appears only multiplied by a zero or unknown input quantity.
Record the input profile, observed output, observation interval, and whether the initial state is known. Identifiability conclusions apply to this full specification.
Set an identically zero input to zero or substitute the actual constant value. Check which parameters remain in the output dynamics.
Construct distinct admissible parameter values and verify that they produce the same complete output, or show that the output equations force equality of the parameters.
A vanishing derivative or rank deficiency is a first-order diagnostic. Do not call it exact non-identifiability without checking the output map beyond that derivative.
An alternative to analyzing only equilibrium output is to retain observations during the approach to equilibrium. In the example above, the transient supplies a time scale through \(\gamma\), while the input-dependent level supplies information about \(\beta\). More generally, input variation can expose effects that disappear under a single restricted profile. The design and analysis of multiple input profiles are separate considerations; the essential point here is that conclusions from one profile should not be generalized to all possible experiments.
These arguments presume the stated model is adequate, the input is known, and the output corresponds directly to \(x(t)\). If the input is uncertain, the initial state is unknown, or the measured output is transformed or noisy, additional ambiguities or practical estimation difficulties may arise. Ideal structural analysis isolates whether the specified equations and experiment permit unique recovery; it does not guarantee precise estimation from real data.
Common Mistakes
- Assuming ideal observation guarantees identifiability. It removes noise and sampling limitations, but cannot reveal a parameter that has no effect on the output under the chosen input and initial condition.
- Claiming that every constant input hides parameters. A constant nonzero input can produce an informative transient. In the model here, a non-equilibrium initial state makes the changing trajectory informative about both rates.
- Treating a zero sensitivity as proof of ambiguity. It proves first-order insensitivity at the point examined, not exact equality of outputs at distinct parameter values. Check the output map itself.
- Confusing a parameter combination with separate parameters. A steady-state level may identify \(\beta c/\gamma\) without identifying \(\beta\) and \(\gamma\) individually.
- Ignoring the initial state. The same input can yield a transient from one initial state and a constant output from an equilibrium initial state, leading to different identifiability conclusions.
Check Your Understanding
Answer each question, then compare with the brief solutions below.
- Under \(u(t)=0\), \(x_0=0.5\) signal units, and \(\gamma=0.2\ \mathrm{min}^{-1}\), does changing \(\beta\) alter the ideal output? What is \(x(3)\)?
- For constant \(c=1\), \(\beta=0.6\) signal units per minute, and \(\gamma=0.3\ \mathrm{min}^{-1}\), calculate the steady-state signal.
- Explain why observing that steady-state signal alone does not distinguish \(\beta\) from \(\gamma\), if \(c\) is known.
- Does a zero derivative of an output with respect to a parameter at one value prove local non-identifiability? State why or why not.
Solutions. (1) No. Under zero input, \(x(t)=x_0e^{-\gamma t}\), so \(x(3)=0.5e^{-0.6}=0.2744\) signal units; \(\beta\) does not enter. (2) \(x_{\mathrm{ss}}=\beta c/\gamma=(0.6)(1)/(0.3)=2.00\) signal units. (3) The output supplies the ratio \(\beta/\gamma=2\) for this known input, not the two rates separately; for example, multiplying both rates by the same positive factor preserves the ratio. (4) No. A zero derivative establishes first-order insensitivity only. The output map must be checked for distinct nearby parameter values with exactly equal outputs.