可辨识性

在 统计学中,可辨识是一个能够更为准确推断的模型必须满足的属性。 一个模型是可辨识的,如果它在理论上能通过无限的观测结果学习到的真正该模型背后参数的真实值。 在数学上,这相当于说基于这些观测结果的不同的参数值必须产生不同的概率分布。 通常情况下,模型只是在某些情况下是可识别的,这些情况的限定条件被称为识别条件

一个模型是不可识别的,如果:两个或两个以上的参数化是观察等价的。 在某些情况下,即使一个模型是不可识别的,它仍然可能学习到某些特定模型参数子集的真实值。 在这种情况下,我们称该模型是部分地可识别的。 在其他情况下,模型可能可以学习到参数空间中一定有限区域的真的参数值,在这种情况下,该模型是集合可识别的。

除了严格的理论探索模型的属性,当使用可识别性分析使用实验数据集检验模型时,可识别性可以在一个更宽泛的范围内被提及。

定义
令 \mathcal{P}=\{P_\theta:\theta\in\Theta\} 为一个 统计模型, 其中参数空间 \Theta 可以是有限或无限维。 我们说 \mathcal{P} 是可识别的,如果映射 \theta\mapsto P_\theta 是 一一映射:

这个定义意味着不同值的 θ 应当对应于不同的概率分布:如果 θ1≠θ2,那么也有 Pθ1≠Pθ2。 如果分布是以概率密度的函数(pdf)方式定义的,那么这两个概率密度函数只有在它们对于一个非零测度集合表现不同时被认为是不同的(例如两个函数ƒ1(x) = 10 ≤ x 和ƒ2(x) = 10 ≤ x ≤ 1 不同之处仅在一个单一点 x = 1—一个测度为零的集合--因此不能被视为不同的概率密度函数)。

模型的可辨识性在映射 \theta\mapsto P_\theta的可逆性的意义上等价于能够在模型无限长的观察后学习模型的真实的参数值。事实上,如果{Xt} ⊆ S 是模型的观测序列,那么根据大数定律,

:
\frac 1 T \sum_{t=1}^T \mathbf{1}_{\{X_t\in A\}} \ \xrightarrow{\text{a.s.}}\ \Pr[X_t\in A],

对于每个可测量的集合A  ⊆   S (此处1 {...}是指示函数 )。 因此,通过无限数量的观察,我们将能够在模型中找到真实概率分布P 0 ,并且由于上述可识别性条件需要映射\theta\mapsto P_\theta是可逆的,我们也能够找到产生给定分布P 0 的真实参数值。

例子
例1
令 \mathcal{P} 是正态位置尺度族:

:
\mathcal{P} = \Big\{\ f_\theta(x) = \tfrac{1}{\sqrt{2\pi}\sigma} e^{ -\frac{1}{2\sigma^2}(x-\mu)^2 }\ \Big|\ \theta=(\mu,\sigma): \mu\in\mathbb{R}, \,\sigma\!>0 \ \Big\}.

那么

:
\begin{align}
& f_{\theta_1}=f_{\theta_2} \\[6pt]
\Longleftrightarrow {} & \frac 1 {\sqrt{2\pi}\sigma_1} \exp\left( -\frac 1 {2\sigma_1^2} (x-\mu_1)^2 \right) = \frac 1 {\sqrt{2\pi}\sigma_2} \exp\left( -\frac 1 {2\sigma_2^2}(x-\mu_2)^2 \right) \\[6pt]
\Longleftrightarrow {} & \frac 1 {\sigma_1^2}(x-\mu_1)^2 + \ln \sigma_1 = \frac 1 {\sigma_2^2}(x-\mu_2)^2 + \ln \sigma_2 \\[6pt]
\Longleftrightarrow {} & x^2 \left(\frac 1 {\sigma_1^2}-\frac 1 {\sigma_2^2}\right) - 2x\left(\frac{\mu_1}{\sigma_1^2}-\frac{\mu_2}{\sigma_2^2} \right) + \left(\frac{\mu_1^2}{\sigma_1^2}-\frac{\mu_2^2}{\sigma_2^2}+\ln\sigma_1-\ln\sigma_2\right) = 0
\end{align}

对于几乎所有的 x 只有当其所有系数都等于零,该公式为零,唯一可能的情况是|σ1|=|σ2|且 μ1 = μ2。 由于在尺度参数 σ 是限制大于零的,我们得出结论,该模型是可辨识的:ƒθ1 = ƒθ2 ⇔ θ1 = θ2。

例2
令 \mathcal{P} 为标准 线性回归模型:

:
y = \beta'x + \varepsilon, \quad \mathrm{E}[\,\varepsilon\mid x\,]=0

(其中,'表示矩阵转置)。 参数 β 是可辨识的,当且仅当矩阵 \mathrm{E}[xx'] 是可逆的。 因此,这是该模型的可辨识条件

例3
假设 \mathcal{P} 是经典的变量误差线性模型:

: \begin{cases}
y = \beta x^* + \varepsilon, \\
x = x^* + \eta,
\end{cases}

其中,(ε,η,x) 是联合正态独立随机变量,其期望为零,方差未知,只有变量(x,y)是观察到的。 那么这个模型是不可识别的, 只有积βσ2∗ (其中σ²∗是差异的潜在回归量 x*)。这也是一个集合可识别的模式的例子:虽然确切的 β 值无法被学习到,我们可以保证,它一定在 (βy,1÷βx-y) 区间中的某处,其中, βyx 是y关于x 的普通最小二乘法 回归的系数,并且 βxy 也是 x 关于 y* 的普通最小二乘法回归的系数。

如果我们放弃正态假设并且要求 x 不是正态分布,仅保留独立的条件 ε ⊥ η ⊥ x**,那么该模型成为可以识别的。

软件
在可部分地观察的动力系统的参数估计情况下, 似然函数也可以被用于结构性和实际可识别性分析。 关于 [http://www.fdmold.uni-freiburg.de/~araue/pmwiki.php/Projects/ProfileLikelihoodApproach] 的一个实现可以在MATLAB工具箱 PottersWheel中获取。

参考

  • 可觀測性
  • 系統識別

参考文献
引文
来源
*
*
*
*
*

进一步阅读
*

评论 (0)

  • 还没有评论,来抢沙发吧。