期刊/特刊下一篇文章
非加性熵S及其在物理和其他方面的应用:一些评论
期刊上的上一篇文章
基于f-发散的地震数据盲反褶积
特刊上一篇文章
文件分类的Tsallis互信息
 
 
订购文章重印
字体类型:
宋体 佐治亚州 宋体,Verdana
字体大小:
澳大利亚 澳大利亚 澳大利亚
行距:
列宽:
背景:
第条

投影幂熵和最大Tsallis熵分布

通过
平口真道(Shinto Eguchi)
*,
小森Osamu Komori
加藤昭吾
日本东京塔川统计数学研究所,190-8562
*
信件应寄给的作者。
2011,13(10), 1746-1764;https://doi.org/10.3390/e13101746
收到的意见:2011年7月26日/修订日期:2011年9月20日/接受日期:2011年9月20日/发布日期:2011年9月26日
(本文属于特刊Tsallis熵)

摘要

:
我们讨论了具有幂指数的两个分布之间的广义交叉熵的单参数族,称为投影幂熵。如果两个分布相等,则交叉熵本质上被简化为Tsallis熵。与射影功率熵相关的统计和概率特性被广泛研究,包括一个表征问题,其中的条件唯一地确定射影功率熵值直至功率指数。熵和勒贝格空间的密切关系 L(左) 第页 和双重 L(左) 探索了护送分发与一个有趣的属性相关联的问题。当我们考虑均值向量和方差矩阵约束下的最大Tsallis熵分布时,模型变为多变量-具有椭圆轮廓的高斯模型,包括高斯和t分布模型。我们通过最小化与投影功率熵相关的经验损失来讨论统计估计。结果表明,最大熵模型下均值向量和方差矩阵的最小损失估计量是样本均值向量和样本方差矩阵。最大熵分布的伴随分布对推导起着关键作用。

1.简介

在经典统计物理和信息论中,与玻尔兹曼-香农熵的密切关系已被很好地建立起来,以提供基本而清晰的理解。Kullback-Leibler散度与极大似然直接相关,极大似然是统计学中最基本的概念之一。Tsallis为功率熵开辟了新的视角,以阐明统计物理中的非平衡态,这些都对非紧张和混沌现象的研究产生了强大的影响,囊性纤维变性。[1,2]. 提出了熵和散度的几种广义形式,囊性纤维变性。[3,4,5,6,7]. 我们考虑定义在有限质量密度函数空间上的广义熵和散度,
F类 = (f) : (f) ( x个 ) d日 x个 < , (f) ( x个 ) 0   几乎无处不在   x个
在Amari提出的信息几何框架中,囊性纤维变性。[8,9].
A功能性 D类 : F类 × F类 [ 0 , ) 称为散度,如果 D类 ( , (f) ) 0 等式当且仅当 = (f) 。如所示[10,11]在建模的流形中,任何散度都与黎曼度量和一对共轭连接相关联 F类 在温和的条件下。
我们从幂交叉熵的原始形式开始[12]使用索引β属于 R(右) 由定义
C类 β ( o个 ) ( , (f) ) = - 1 β ( x个 ) { (f) ( x个 ) β - 1 } d日 x个 + 1 1 + β (f) ( x个 ) 1 + β d日 x个
为所有人(f)在里面 F类 幂(对角线)熵
H(H) β ( o个 ) ( (f) ) = C类 β ( o个 ) ( (f) , (f) ) = - 1 β ( β + 1 ) (f) ( x个 ) 1 + β d日 x个 + 1 β (f) ( x个 ) d日 x个
请参见[13,14]用于独立成分分析和模式识别的信息几何和统计应用。注意,这是在概率密度函数的连续情况下定义的,但可以简化为离散情况,参见Tsallis[2]对于统计物理的广泛讨论。事实上,Tsallis熵
S公司 ( (f) ) = 1 - 1 1 - (f) ( x个 ) d日 x个
对于概率密度函数 (f) ( x个 ) 与功率熵和常数成正比 H(H) β ( o个 ) ( ) - 1 ,其中 = 1 + β 。功率发散由下式给出
D类 β ( o个 ) ( , (f) ) = C类 β ( o个 ) ( , (f) ) - H(H) β ( o个 ) ( )
通常,as由交叉熵和对角熵的差值定义。
在本文中,我们重点讨论了由
C类 γ ( , (f) ) = - 1 γ ( 1 + γ ) ( x个 ) (f) ( x个 ) γ d日 x个 (f) ( x个 ) 1 + γ d日 x个 γ 1 + γ
所以投影幂熵是
H(H) γ ( (f) ) = - 1 γ ( 1 + γ ) (f) ( x个 ) 1 + γ d日 x个 1 1 + γ
的日志表达式 C类 γ ( , (f) ) 由定义
C类 γ 日志 ( , (f) ) = - 1 γ 日志 { - γ ( 1 + γ ) C类 γ ( , (f) ) }
请参见[15,16]用于推导 C类 γ 日志 ,并详细讨论了 C类 β ( o个 ) ( , (f) ) C类 γ ( , (f) ) .投影幂交叉熵 C类 γ ( , (f) ) 满足以下方面的线性以及投影不变性,即 C类 γ ( , λ (f) ) = C类 γ ( , (f) ) 对于任何常数 λ > 0 。请注意 H(H) γ ( (f) ) 与……一一对应 S公司 ( (f) ) 由提供
H(H) γ ( (f) ) = - 1 ( - 1 ) { 1 - ( - 1 ) S公司 ( (f) ) } 1
哪里 = 1 + γ .投影功率散度为
D类 γ ( , (f) ) = C类 γ ( , (f) ) - H(H) γ ( )
这将在与H的密切关系上进行讨论 o个 ¨ 奥尔德不等式。散度定义为 C类 γ ( , (f) ) 满足
D类 γ 日志 ( , (f) ) = C类 γ 日志 ( , (f) ) - C类 γ 日志 ( , ) 0
为所有人γ属于 R(右) 如果存在积分 D类 γ 日志 ( , (f) ) .非负性导致
D类 γ ( , (f) ) 0
我们注意到幂指数的存在范围γ对于 C类 γ ( , (f) ) H(H) γ ( (f) ) 取决于样本空间(f)定义。如果样本空间紧凑 C类 γ ( , (f) ) H(H) γ 对所有人都有明确的定义 γ R(右) .如果样本空间不紧凑, C类 γ ( , (f) ) 为定义 γ 0 H(H) γ ( (f) ) 是用于 γ > - 1 。更准确地说,我们将探讨样本空间是 R(右) d日 在随后的讨论中以及力矩条件。通常我们观察到
极限 γ 0 D类 γ ( , (f) ) = D类 0 ( , (f) )
哪里 D类 0 ( , (f) ) 表示Kullback-Leibler散度,
D类 0 ( , (f) ) = ( x个 ) 日志 ( x个 ) (f) ( x个 ) d日 x个
请参见附录1用于推导(5).
{ x个 1 , , x个 n个 } 是概率密度函数分布的随机样本 ( x个 ) .统计模型 { (f) ( x个 , θ ) : θ Θ } 带参数θ假设充分接近潜在密度函数 ( x个 ) ,其中θ是参数空间。然后是与投影功率熵相关的损失函数 C类 γ ( , (f) ( · , θ ) ) 根据以下示例
L(左) γ ( θ ) = - 1 γ ( 1 + γ ) 1 n个 = 1 n个 k个 γ ( θ ) (f) ( x个 , θ ) γ
我们在其中呼叫
θ ^ γ 选择 argmin(最小值) θ Θ L(左) γ ( θ )
这个γ-估计器,其中
k个 γ ( θ ) = (f) ( x个 , θ ) 1 + γ d日 x个 - γ 1 + γ
我们注意到
E类 { L(左) γ ( θ ) } = C类 γ ( , (f) ( · , θ ) )
哪里 E类 表示关于以下方面的统计期望可以观察到,由于损失,0-估计器只不过是最大似然估计器(MLE) L(左) γ ( θ ) 收敛到负对数似然函数,
L(左) 0 ( θ ) = - 1 n个 = 1 n个 日志 (f) ( x个 , θ )
在这个意义上
极限 γ 0 L(左) γ ( θ ) + 1 γ ( 1 + γ ) = L(左) 0 ( θ )
如果潜在密度函数 ( x个 ) 属于具有均值的高斯模型μ和方差 σ 2 ,然后是MLEμ σ 2 是样本均值和样本方差。相反的语句如所示[17,18]. 我们将把这个理论推广到γ-估计值低于γ-模型。
第2节我们讨论了投影幂熵的特征。第3节具有Tsallis熵的最大熵分布 S公司 具有 = 1 + γ 在平均向量约束下μ并且考虑了方差矩阵∑。我们讨论了最大熵分布模型,称为γ-模型,其中0模型和2模型分别等于高斯模型和维格纳模型。然后我们显示γ-估计值μ和∑γ-模型是样本均值和样本方差。第4节给出结论和进一步的评论。

2.射影不变性

让我们看看 F类 用勒贝格的空间
L(左) 第页 = (f) ( x个 ) : | (f) ( x个 ) | 第页 d日 x个 <
哪里 第页 1 L(左) 第页 -规范 第页 由定义
(f) 第页 = | (f) ( x个 ) | 第页 d日 x个 1 第页
是的共轭指数第页令人满意的 1 / 第页 + 1 / = 1 ,其中第页可以表示为参数的函数 γ > 0 这样的话 第页 = 1 + γ - 1 = 1 + γ 。我们注意到等于索引Tsallis熵 S公司 在关系中 = 1 + γ .对于任何概率密度函数 (f) ( x个 ) 我们用概率密度函数定义护送分布,
e(电子) ( (f) ( x个 ) ) = (f) ( x个 ) (f) ( ) d日
囊性纤维变性。[2]供广泛讨论。我们讨论了投影交叉熵的一个有趣的关系(1)护送分发。根据护送分配的定义,
C类 γ ( , (f) ) = - 1 γ ( 1 + γ ) { e(电子) ( (f) ( x个 ) ) } 1 第页 ( x个 ) d日 x个
我们注意到 e(电子) ( (f) ) 1 第页 在的单位范围内 L(左) 第页 在代表中。投影幂对角熵(2)与 L(左) -也就是说,
H(H) γ ( (f) ) = - 1 γ ( 1 + γ ) (f)
其中H o个 ¨ 奥尔德不等式
( x个 ) (f) ( x个 ) γ d日 x个 (f) γ 第页
声称 C类 γ ( , (f) ) H(H) γ ( ) ,或同等
D类 γ ( , (f) ) 0
为所有人(f)在里面 F类 也是由 C类 γ ( o个 ) ( , (f) ) H(H) γ ( o个 ) ( ) .中的平等(10)仅当且仅当 (f) ( x个 ) = λ ( x个 ) 几乎无处不在x个,其中λ是一个正常数。权力转换表明了空间之间的相互作用 L(左) 第页 L(左) 根据关系,
(f) γ 第页 = (f) γ
取的极限γH中为0 o个 ¨ 奥尔德不等式(9)产生这样的结果
( x个 ) 日志 (f) ( x个 ) d日 x个 ( x个 ) 日志 ( x个 ) d日 x个
自从
极限 γ 0 ( x个 ) (f) ( x个 ) γ - 1 γ d日 x个 = ( x个 ) 日志 (f) ( x个 ) d日 x个
极限 γ 0 (f) γ 第页 - 1 γ = ( x个 ) 日志 ( x个 ) d日 x个
此限制涉及第页与另一个空间而不是 L(左) 比如说,空间就是所有具有有限Boltzmann-Shannon熵的密度函数的空间 L(左) 日志 .功率指数γ重新参数化勒贝格空间 L(左) 第页 和双重空间 L(左) 与关系 第页 = 1 + γ - 1 然而,要进行权力转换 (f) ( x个 ) γ 与勒贝格空间的普通讨论完全不同,因此对偶性收敛于 ( L(左) 日志 , L(左) 1 ) 如中所观察到的(11). 在信息几何中,这对 ( L(左) 日志 , L(左) 1 ) 对应于混合连接和指数连接,囊性纤维变性。[9]. 另请参见另一个参数化 L(左) 第页 空间[19].
我们现在讨论一个关于 C类 γ ( , (f) ) 如下定理所示。关于特性的一般讨论见[16]然而,推导过程相当复杂。这里我们假设一个关键条件,即交叉熵 Γ ( , (f) ) 是线性的给出一个基本的证明。Riesz表示定理表明
Γ ( , (f) ) = c(c) ( (f) ) ( x个 ) ψ ( (f) ( x个 ) ) d日 x个
哪里 c(c) ( (f) ) 是一个常数,取决于(f)因此,当我们为 c(c) ( (f) ) 以保证缩放不变性。
定理1。.定义函数 Γ : F类 × F类 R(右) 通过
Γ ( , (f) ) = φ ρ ( (f) ( x个 ) ) d日 x个 ( x个 ) ψ ( (f) ( x个 ) ) d日 x个
哪里φ,ρψ是可微的单调函数。假设
(i) 。 Γ ( , ) = 最小值 (f) F类 Γ ( , (f) ) 为所有人 F类 ,
还有那个
(ii)。 Γ ( , (f) ) = Γ ( , λ (f) ) 为所有人 λ > 0 以及所有 , (f) F类 .
然后就有了γ这样的话 Γ ( , (f) ) = C类 γ ( , (f) ) 高达常数因子,其中 C类 γ ( , (f) ) 是投影幂交叉熵,定义为(1).
证明。要求(ii)意味着
λ φ ρ ( λ (f) ( x个 ) ) d日 x个 ψ ( λ (f) ( x个 ) ) ( x个 ) d日 x个 = 0
这意味着,如果(f)是绝对连续的狄拉克的测量值是 x个 0 ,然后
ψ ˙ ( λ (f) ( x个 0 ) ) ψ ( λ (f) ( x个 0 ) ) λ (f) ( x个 0 ) = c(c) ( λ )
哪里
c(c) ( λ ) = - λ φ ˙ ρ ( λ (f) ( x个 ) ) d日 x个 ρ ˙ ( λ (f) ( x个 ) ) (f) ( x个 ) d日 x个 φ ρ ( λ (f) ( x个 ) ) d日 x个
因为我们可以取任意值 (f) ( x个 0 ) 对于任何固定λ,
ψ ˙ ( t吨 ) ψ ( t吨 ) = c(c) ( λ ) t吨 - 1
其唯一解为 ψ ( t吨 ) = t吨 γ 哪里 γ = c(c) ( λ ) 。接下来,我们考虑一个有限离散空间的情况, { x个 : 1 } 然后,因为 ψ ( (f) ) = (f) γ ,我们可以写
Γ ( , (f) ) = φ = 1 ρ ( (f) ) = 1 (f) γ
哪里 (f) = (f) ( x个 ) = ( x个 ) 要求(i)导致 ( / (f) j个 ) Γ ( , (f) ) | (f) = = 0 为所有人 j个 , 1 j个 ,这意味着
ρ ˙ ( j个 ) = - γ c(c) ( 1 , , ) j个 γ
哪里
c(c) ( 1 , , ) = φ = 1 ρ ( ) = 1 1 + γ φ ˙ = 1 ρ ( )
它源自(13)那个 c(c) ( 1 , , ) 必须是中的常量 1 , , ,说吧C类,以便我们解决(13)作为 ρ ( j个 ) = - γ C类 j个 1 + γ / ( 1 + γ ) 因此,方程式(14)由编写
φ ˙ ( t吨 ) φ ( t吨 ) = - γ 1 + γ t吨 - 1
这将导致 φ ( t吨 ) = t吨 - γ 1 + γ 我们的结论是 Γ ( , (f) ) C类 γ ( , (f) ) ,这就完成了证明。
备注1。上述证明基本上适用于积分(11)仅由二进制分布的求和给出的情况。在这个意义上,定理1的陈述并不严密,但是,统计推断是以统一的方式进行讨论的,这样分布要么是连续的,要么是离散的。在随后的讨论中,我们将重点讨论在 R(右) d日 .
备注2。我们看到的乘法分解 C类 γ ( , (f) ) 实现统计独立性。事实上,如果(f)分解为 (f) = (f) 1 (f) 2 , = 1 2 在同一分区中,则
C类 γ ( , (f) ) = C类 γ ( 1 , (f) 1 ) C类 γ ( 2 , (f) 2 )
这个性质也是基本的,但我们不假设这个可分解性是定理1的要求。

3.最大熵分布模型

我们将阐明上的最大熵模型之间的二元结构 H(H) γ ,定义于(2)最小交叉熵估计 C类 γ ,定义于(1). 在讨论之前,我们概述了经典的情况,其中最大似然估计在玻尔兹曼-香农熵的最大熵模型下表现良好,也就是说,如果我们考虑均值和方差约束,则为高斯模型。我们将使用传统的符号X(X)表示具有值的随机变量x个.让 { x个 1 , , x个 n个 } 是具有密度函数的高斯分布的随机样本
(f) 0 ( x个 , μ , ) = det(探测) ( 2 π ) - 1 2 经验 { - 1 2 ( x个 - μ ) T型 - 1 ( x个 - μ ) }
高斯密度函数由标准形式表示
( 2 π ) - d日 2 经验 - 1 2 ( x个 - μ ) T型 Ξ ( x个 - μ ) + 1 2 日志 det(探测) ( Ξ )
其中Ξ被称为由 - 1 .分化(15)上的μ和Ξ产量
E类 (f) 0 ( · , μ , ) ( X(X) ) = μ V(V) (f) 0 ( · , μ , ) ( X(X) ) =
哪里 E类 (f) V(V) (f) 表示关于概率密度函数的期望向量和方差矩阵 (f) ( x个 ) 分别是。
最大似然估计量由下式给出
( μ ^ 0 , ^ 0 ) = ( x个 ¯ , S公司 )
哪里 x个 ¯ S公司是样本均值向量和样本方差矩阵,
x个 ¯ = 1 n个 = 1 n个 x个 , S公司 = 1 n个 = 1 n个 ( x个 - x个 ¯ ) ( x个 - x个 ¯ ) T型
这是因为负log-likelihood函数是
L(左) 0 ( μ , ) = - 1 n个 = 1 n个 日志 (f) 0 ( x个 , μ , )
由以下人员编写
1 2 追踪 ( S公司 ( μ ) Ξ ) - 1 2 日志 det(探测) ( Ξ )
除了常数,其中
S公司 ( μ ) = 1 n个 = 1 n个 ( x个 - μ ) ( x个 - μ ) T型
因此,估算系统为
μ L(左) 0 ( μ , ) Ξ L(左) 0 ( μ , ) = Ξ ( x个 ¯ - μ ) 1 2 { S公司 ( μ ) - } = 0 O(运行)
它结束了表达式(16)MLE自 S公司 ( μ ) = S公司 + ( x个 ¯ - μ ) ( x个 ¯ - μ ) T型 。或者,我们还有另一条路线要展示(16)如下所示。中定义的Kullback-Leibler散度(6)由提供
D类 0 ( (f) 0 ( · , μ , ) , (f) 0 ( · , μ 1 , 1 ) )
= 1 2 ( μ - μ 1 ) T型 1 - 1 ( μ - μ 1 ) + 1 2 追踪 { ( - 1 ) 1 - 1 } - 1 2 日志 det(探测) ( 1 - 1 )
因此,我们观察到
L(左) 0 ( μ , ) - L(左) 0 ( x个 ¯ , S公司 ) = D类 0 ( (f) 0 ( · , x个 ¯ , S公司 ) , (f) 0 ( · , μ , ) )
当且仅当 ( μ , ) = ( x个 ¯ , S公司 ) 。这意味着(16).
在温和的规则性条件下,相反的语句成立,即位置和散布模型的MLE满足(16)当且仅当模型为高斯时,囊性纤维变性。[17,18]. 然而,即使我们对底层分布不做任何假设 ( x个 ) 、统计数据 x个 ¯ S公司是渐近一致的
μ = E类 ( X(X) ) = V(V) ( X(X) )
这是强大数定律的直接结果,中心极限定理导致这两个统计量的渐近正态性。从这个意义上说, ( x个 ¯ , S公司 ) 也是一个非参数估计量 ( μ , ) .
我们探讨了统计模型和估计方法之间的密切关系。我们考虑具有γ-熵 H(H) γ 在…的空间上d日-具有共同均值和方差的维度分布,
F类 ( μ , ) = { (f) F类 : (f) ( x个 ) d日 x个 = 1 , E类 (f) ( x个 ) = μ , V(V) (f) ( x个 ) = }
然后我们定义一个概率密度函数的分布
(f) γ ( x个 , μ , ) = c(c) γ det(探测) ( 2 π ) 1 2 1 - γ 2 + d日 γ + 2 γ ( x个 - μ ) T型 - 1 ( x个 - μ ) + 1 γ
哪里 ( ) + 表示积极的部分 c(c) γ 是归一化因子,
c(c) γ = 2 γ 2 + d日 γ + 2 γ d日 2 Γ ( 1 + d日 2 + 1 γ ) Γ ( 1 + 1 γ ) 如果   γ > 0 - 2 γ 2 + d日 γ + 2 γ d日 2 Γ ( - 1 γ ) Γ ( - 1 γ - d日 2 ) 如果   - 2 d日 + 2 < γ < 0
参见推导 c(c) γ 在里面附录2.如果尺寸d日等于1,则 (f) γ ( x个 , μ , ) 是一个-高斯分布 = γ + 1 。我们注意到
极限 γ 0 c(c) γ = 极限 γ 0 c(c) γ = 1
在哪儿 (f) γ ( x个 , μ , ) 减少为d日-变化高斯密度,当 γ = 0 .支持 (f) γ ( · , μ , ) 变为椭球体,定义为
x个 R(右) d日 : ( x个 - μ ) T型 - 1 ( x个 - μ ) 2 + d日 γ + 2 γ γ
如果 γ > 0 另一方面,如果 - 2 d日 + 2 < γ < 0 ,密度函数(21)写为
(f) γ ( x个 , μ , ) = det(探测) ( π τ ) - 1 2 Γ ( - 1 γ ) Γ ( - 1 γ - d日 2 ) 1 + 1 τ ( x个 - μ ) T型 - 1 ( x个 - μ ) 1 γ
哪里
τ = - 2 + ( d日 + 2 ) γ γ
这个d日-变量t分布定义为
ν ( x个 , μ , P(P) ) = det(探测) ( π ν P(P) ) - 1 2 Γ ( ν + d日 2 ) Γ ( ν 2 ) 1 + 1 ν ( x个 - μ ) T型 P(P) - 1 ( x个 - μ ) - ν + d日 2
囊性纤维变性。[20]进行广泛讨论。假设
ν + d日 2 = - 1 γ ν P(P) = τ
然后我们从(23)和(24)那个
(f) γ ( x个 , μ , ) = ν ( x个 , μ , P(P) )
因此,密度函数 (f) γ ( x个 , μ , ) 具有 - 2 d日 + 2 < γ < 0 是t分布。分布在欧氏空间上具有椭圆轮廓 R(右) d日 对于任何 γ > - 2 d日 + 2 ,如所示图1对于以下典型情况γ.
图1。t分布 ( γ = - 0 . 4 ) ,高斯 ( γ = 0 ) 和维格纳 ( γ = 2 ) 分配。
图1。t分布 ( γ = - 0 . 4 ) ,高斯 ( γ = 0 ) 和维格纳 ( γ = 2 ) 分配。
熵13 01746 g001a熵13 01746 g001b
M(M) γ = (f) γ ( x个 , μ , ) : μ R(右) d日 , S公司 d日
我们称之为γ-模型,其中 S公司 d日 表示所有对称正定矩阵的空间d日.我们确认了γ-模型如下。
柠檬.在模型下 M(M) γ 定义于(25)使用索引 γ > - 2 d日 + 2 ,
E类 (f) γ ( · , μ , ) ( X(X) ) = μ V(V) (f) γ ( · , μ , ) ( X(X) ) =
证明。我们需要考虑一系列护送分发。在模型中 M(M) γ 我们可以将护送分配定义为
e(电子) ( (f) γ ( x个 , μ , ) ) = c(c) γ * det(探测) ( ) 1 2 1 - γ 2 + d日 γ + 2 γ ( x个 - μ ) T型 - 1 ( x个 - μ ) + 1 + γ γ
哪里 = 1 + γ c(c) γ * 是归一化因子。因此,
e(电子) ( (f) γ ( x个 , μ , ) ) = c(c) γ * det(探测) ( ) - 1 2 γ 1 + γ - γ 2 + d日 γ + 2 γ ( x个 - μ ) T型 { det(探测) ( ) - 1 2 γ 1 + γ - 1 } ( x个 - μ ) + 1 + γ γ
这里我们定义了可选参数 Ξ γ 转换为原始参数∑
Ξ γ = det(探测) ( ) - 1 2 γ 1 + γ - 1
因此,逆变换由下式给出
= det(探测) ( Ξ γ ) γ d日 γ + 2 γ + 2 Ξ γ - 1
注意到了 det(探测) ( Ξ γ ) = det(探测) ( ) - 1 2 d日 γ + 2 γ + 2 1 + γ 因此,我们得到了(26)作为
e(电子) ( (f) γ ( x个 , μ , ) ) = c(c) γ * det(探测) ( Ξ γ ) γ 2 + d日 γ + 2 γ - γ 2 + d日 γ + 2 γ ( x个 - μ ) T型 Ξ γ ( x个 - μ ) + 1 + γ γ
通过对指数族的讨论进行类比,我们得到了中支撑项的以下表达式(30)作为
- 2 γ 2 + d日 γ + 2 γ 1 2 追踪 ( x个 x个 T型 Ξ γ ) - μ T型 Ξ γ x个 + 1 2 μ T型 Ξ γ μ - 2 + d日 γ + 2 γ 2 γ det(探测) ( Ξ γ ) γ d日 γ + 2 γ + 2
护送分布的一个特性建议了该分布的力矩公式(25)如下:我们有一个身份
μ c(c) γ * det(探测) ( Ξ γ ) γ d日 γ + 2 γ + 2 - γ 2 + d日 γ + 2 γ ( x个 - μ ) T型 Ξ γ ( x个 - μ ) + 1 + γ γ d日 x个 = 0
这意味着
det(探测) ( Ξ γ ) γ d日 γ + 2 γ + 2 - γ 2 + d日 γ + 2 γ ( x个 - μ ) T型 Ξ γ ( x个 - μ ) + 1 γ Ξ γ ( x个 - μ ) d日 x个 = 0
由此得出结论
E类 (f) γ ( · , μ , ) ( X(X) ) = μ
同样,
Ξ γ c(c) γ * det(探测) ( Ξ γ ) γ d日 γ + 2 γ + 2 - γ 2 + d日 γ + 2 γ ( x个 - μ ) T型 Ξ γ ( x个 - μ ) + 1 + γ γ d日 x个 = 0
哪个是
c(c) γ * det(探测) ( Ξ γ ) γ d日 γ + 2 γ + 2 - γ 2 + d日 γ + 2 γ ( x个 - μ ) T型 Ξ γ ( x个 - μ ) + 1 γ
× γ d日 γ + 2 γ + 2 det(探测) ( Ξ γ ) γ d日 γ + 2 γ + 2 Ξ γ - 1 - γ 2 + d日 γ + 2 γ ( x个 - μ ) ( x个 - μ ) T型 d日 x个 = 0
由此得出结论
V(V) (f) γ ( · , μ , ) ( X(X) ) =
因为 Ξ γ 和∑,如(29). 证据是完整的。 ☐
备注3.规范形式(30)护送分发的(26)在引理的证明中起着重要作用。基本上我们可以写出(21)然而,像指数族的情况一样,它与分布性质没有任何联系。
备注4.在方程式中(31)函数
φ ( Ξ ) = 1 2 ω det(探测) ( Ξ ) ω
被视为Fenchel凸对偶中的一个势函数,其中
ω = γ 2 + d日 γ + 2 γ
囊性纤维变性。[21,22]协方差结构模型。
从引理我们观察到 (f) γ ( · , μ , ) F类 ( μ , ) 下一步,我们展示了密度分布 (f) γ ( · , μ , ) 最大化γ-熵 H(H) γ 在空间上方 F类 ( μ , ) ,其中 H(H) γ 定义于(2).
定理2.
(i) ●●●●。如果 - 2 d日 + 2 < γ 0 ,然后
(f) γ ( · , μ , ) = argmax(最大值) (f) F类 ( μ , ) H(H) γ ( (f) )
哪里 F类 ( μ , ) 在中定义(20).
(ii)。如果 γ > 0 ,然后
(f) γ ( · , μ , ) = argmax(最大值) (f) F类 ( μ , ) ( γ ) H(H) γ ( (f) )
哪里
F类 ( μ , ) ( γ ) = { (f) F类 ( μ , ) : (f) ( x个 ) = 0   几乎无处不在   x个 B类 ( μ , ) }
具有 B类 ( μ , ) 存在 { x个 R(右) d日 : ( x个 - μ ) T型 - 1 ( x个 - μ ) > 2 + d日 γ + 2 γ γ } .
证明根据的定义 F类 ( μ , ) ,我们从引理中看到 (f) γ ( · , μ , ) F类 ( μ , ) 对于任何 γ ( - 2 d日 + 2 , 0 ) 。这导致
E类 (f) γ ( · , μ , ) { (f) γ ( X(X) , μ , ) γ } = E类 (f) { (f) γ ( X(X) , μ , ) γ }
对于任何(f)在里面 F类 ( μ , ) ,这意味着
H(H) γ ( (f) γ ( · , μ , ) ) = C类 γ ( (f) , (f) γ ( · , μ , ) )
因此
H(H) γ ( (f) γ ( · , μ , ) ) - H(H) γ ( (f) ) = D类 γ ( (f) , (f) γ ( · , μ , ) )
如中所述为非负(4). 这就结束了(34). 同样,我们观察到(36)保留任何 γ > 0 以及任何(f)在里面 F类 ( μ , ) ( γ ) 自从支持(f)包括 (f) ( · , μ , ) 。此结论(35).
我们想通过数据集的最小投影交叉熵来阐明统计推断的类似结构 { x个 1 , , x个 n个 } 假设遵循模型 M(M) γ .我们回忆起(8)从射影交叉熵与伴随分布的关系
C类 γ ( , (f) ) = - 1 γ ( 1 + γ ) e(电子) ( (f) ( x个 ) ) γ 1 + γ ( x个 ) d日 x个
当我们有数据时 { x个 1 , , x个 n个 } 安装在模型上 M(M) γ ,损失函数为
L(左) γ ( μ , ) = - 1 γ ( 1 + γ ) 1 n个 = 1 n个 e(电子) ( (f) γ ( x个 , μ , ) ) γ 1 + γ
哪里 (f) γ ( x个 , μ , ) 定义于(21). 这个γ-估计器定义为
( μ ^ γ , ^ γ ) = argmin(最小值) ( μ , ) L(左) γ ( μ , )
参见一般定义(7). 它遵循中定义的规范形式(30)使用规范参数 Ξ γ 定义于(28)那个
L(左) γ ( μ , ) = - 1 γ ( 1 + γ ) ( c(c) γ * ) γ γ + 1 det(探测) ( Ξ γ ) ω - ω { 追踪 ( Ξ γ S公司 ) + ( μ - x个 ¯ ) T型 Ξ γ ( μ - x个 ¯ ) }
哪里 ( x个 ¯ , S公司 ) ω定义于(17)和(33)、和 c(c) γ * 是中定义的标准化因子(27). 这里我们注意到,如果 γ > 0 ,然后是参数 ( μ , ) 必须假定处于 Θ n个 ,其中
Θ n个 = { ( μ , ) R(右) d日 × S公司 d日 : ( x个 - μ ) T型 - 1 ( x个 - μ ) < ω - 1 ( = 1 , , n个 ) }
我们注意到 L(左) γ ( μ , ) = C类 γ ( (f) ( · , x个 ¯ , S公司 ) , (f) ( · , μ , ) ) L(左) γ ( x个 ¯ , S公司 ) = H(H) γ ( (f) ( · , x个 ¯ , S公司 ) ) 自从
E类 (f) ( · , x个 ¯ , S公司 ) ( X(X) ) = x个 ¯ , V(V) (f) ( · , x个 ¯ , S公司 ) ( X(X) ) = S公司
因此,我们观察到的论点类似于(19)对于MLE。投影发散 D类 γ 定义于(3)等于γ-损失函数为
L(左) γ ( μ , ) - L(左) γ ( x个 ¯ , S公司 ) = D类 γ ( (f) γ ( · , x个 ¯ , S公司 ) , (f) γ ( · , μ , ) ) ,
当且仅当 ( μ , ) = ( x个 ¯ , S公司 ) 见方程式后的讨论(10). 这样,我们可以将上述讨论总结如下:
定理3.让 { x个 1 , , x个 n个 } 是来自γ-中定义的模型(21). 然后γ-中定义的估计器(7)的 ( μ , ) ( x个 ¯ , S公司 ) ,其中 ( x个 ¯ , S公司 ) 定义于(17).
证明。让我们再提供一个证据。估算系统如下所示
μ L(左) γ ( μ , ) Ξ γ L(左) γ ( μ , ) = Ξ γ ( x个 ¯ - μ ) ω { det(探测) ( Ξ γ ) ω Ξ γ - 1 - S公司 ( μ ) } = 0 O(运行)
相当于
μ - x个 ¯ - S公司 ( μ ) = 0 O(运行)
因为 Ξ γ 转换为∑,如(29). 因此,我们也得出了结论 ( μ ^ γ , ^ γ ) = x个 ¯ , S公司 . 这样,我们就得到了由(40)通过参数 Ξ γ 利用押运分布与损失函数的关系(37).
备注5考虑位置模型 { (f) γ ( · , μ , ) } 使用位置参数μ,其中∑在定理3中已知。然后我们很容易看到γ-估计量μ x个 ¯ 相反的说法呢?我们发现,如果γ-估计量μ x个 ¯ 根据样本大小 n个 3 ,则模型为γ-模型, { (f) γ ( · , μ , ) 已知∑。该证明与[17]. 事实上,我们得出的结论是模型密度函数 (f) ( x个 ) 满足这一点
{ (f) ( x个 - μ ) } γ = + b条 ( x个 - μ ) T型 - 1 ( x个 - μ )
哪里b条是常数。
备注6。如果我们共同研究定理2和3,那么
最小值 ( μ , ) R(右) d日 × S公司 d日 L(左) γ ( μ , ) = 最大值 (f) F类 ( x个 ¯ , S公司 ) H(H) γ ( (f) )
自从 L(左) γ ( x个 ¯ , S公司 ) = H(H) γ ( (f) γ ( · , x个 ¯ , S公司 ) ) .两侧(41)与不平等联系在一起(39)和(36)上的γ-单独讨论中的分歧。
备注7.的推导γ-定理3中的估计器是由正则参数提供的 Ξ γ 中给出的护送分配(28). 这里,我们直接计算损耗相对于∑的梯度,如下所示:
L(左) γ ( μ , ) = - 1 2 ( 1 + γ ) 2 det(探测) ( ) - 1 2 γ 1 + γ ( 1 - ω 追踪 { S公司 ( μ ) - 1 } ) - 1 + γ ( 1 + γ ) ω det(探测) ( ) - 1 2 γ 1 + γ - 1 S公司 ( μ ) - 1 = - 1 2 ( 1 + γ ) 2 det(探测) ( ) - 1 2 γ 1 + γ × ( 1 - ω 追踪 { S公司 ( μ ) - 1 } ) - 1 - 1 + γ 1 + 1 2 d日 γ + γ - 1 S公司 ( μ ) - 1
因此,我们观察到,如果我们 μ = x个 ¯ = α S公司 ( x个 ¯ ) ,然后
L(左) γ ( x个 ¯ , α S公司 ( x个 ¯ ) ) = - 1 2 γ 1 + γ det(探测) ( α S公司 ( x个 ¯ ) ) - 1 2 γ 1 + γ ( α S公司 ( x个 ¯ ) ) - 1
× ( 1 - ω 追踪 { S公司 ( x个 ¯ ) ( α S公司 ( x个 ¯ ) ) - 1 } ) α S公司 ( x个 ¯ ) - 1 + γ 1 + 1 2 d日 γ + γ S公司 ( x个 ¯ ) ( α S公司 ( x个 ¯ ) ) - 1
括号中的术语(42)由提供
α ( 1 - ω 追踪 { S公司 ( x个 ¯ ) ( α S公司 ( x个 ¯ ) ) - 1 } ) - 1 + γ 1 + 1 2 d日 γ + γ S公司 ( x个 ¯ ) = α - d日 γ 2 + d日 γ + 2 γ - 1 + γ 1 + 1 2 d日 γ + γ S公司 ( x个 ¯ )
结论是如果 α = 1 ,然后 ( / ) L(左) γ ( x个 ¯ , α S公司 ( x个 ¯ ) ) = 0 这是定理3的直接证明,但它将伴随着对替换 ( μ , ) 进入之内 ( x个 ¯ , α S公司 ( x个 ¯ ) ) .

4.结束语

我们探索了优雅的酒店(39),毕达哥拉斯经验关系γ-模型和γ-估计器,在这个意义上(39)直接给出定理3,无需任何微分。另一个优雅的表达是大自然和决策者之间的minimax游戏,参见[23]. 考虑一下空间 F类 ( μ , ) 定义于(20). The intersection of theγ-模型(21)和 F类 ( μ , ) 是独生子女 { (f) γ ( · , μ , ) } ,这是的极大极小解
最大值 F类 ( μ , ) 最小值 (f) F类 C类 γ ( , (f) ) = C类 γ ( (f) γ ( · , μ , ) , (f) γ ( · , μ , ) )
考虑不同的指数γ γ * 其中指定了γ-模型和 γ * -估计值。基本上 γ * -估计量在γ-适用于任何选择的模型γ γ * 。如果我们特别修复 γ = 0 对于模型,即高斯模型,则 γ * -该估计器对于任何 γ * > 0 ,请参阅[16]. 鲁棒性程度与 γ * 以效率为代价。这个 γ * -估计量 ( μ , ) 高斯模型的解是
μ = = 1 n个 (f) 0 ( x个 , μ , ) γ * x个 = 1 n个 (f) 0 ( x个 , μ , ) γ *
= ( 1 + γ * ) = 1 n个 (f) 0 ( x个 , μ , ) γ * ( x个 - μ ) ( x个 - μ ) T型 = 1 n个 (f) 0 ( x个 , μ , ) γ *
重量函数 (f) 0 ( x个 , μ , ) γ * 对于-第th次观测 x个 x个 是一个离群值。或者,经典稳健方法采用 γ * = 0 ,即指定错误的模型的MLE γ < 0 t吨-分布模型,请参见[24,25]. 因此,不同的指数γ γ * 以二元方式进行稳健的统计。
该属性是指数模型和MLE之间关联属性的扩展,然而,它在以下意义上是脆弱的(19)如果γ-模型和 γ * -估计值略有不同。在实践中,我们发现在求解最大似然估计的数值任务中存在一些困难γ-模型,带有 γ > 0 因为密度的支撑取决于参数和指数γ我们讨论了模型的统计和概率特性以及与特定交叉熵相关的估计。讨论的一部分性质仍然适用于更广泛的类中的任何交叉熵,这是从[13,26].

致谢

我们感谢两位裁判的宝贵意见和建设性建议。

工具书类

  1. Tsallis,C.Boltzmann-Gibbs统计的可能推广。J.统计。物理学。 1988,52, 479–487. [谷歌学者][交叉参考]
  2. C·查利斯。非扩展统计力学导论:走向复杂世界; Springer-Verlag:美国纽约州纽约市,2009年。[谷歌学者]
  3. Cichocki,A。;克鲁斯,S。;Amari,S.α-β-和γ-发散的家族:灵活而稳健的相似性度量。 2010,12, 1532–1568. [谷歌学者][交叉参考]
  4. Cichocki,A。;克鲁斯,S。;Amari,S.广义α-β发散及其在稳健非负矩阵分解中的应用。 2011,13, 134–170. [谷歌学者][交叉参考]
  5. Csiszár,I.概率分布差异和间接观测的信息型度量。匈牙利科学研究所 1967,2, 229–318. [谷歌学者]
  6. 关于熵和信息的度量。程序。伯克利第四交响乐团。数学。统计师。探针。 1961,1, 547–561. [谷歌学者]
  7. Tøpsoe,F.信息差异的一些不平等现象和相关的歧视措施。IEEE传输。通知。西奥。 2000,46,1602年至1609年。[谷歌学者]
  8. Amari,S.统计学中的微分几何方法。统计学课堂讲稿; 施普林格出版社:美国纽约州纽约市,1985年;第28卷。[谷歌学者]
  9. 阿玛里,S。;长冈,H.信息几何方法。数学专著的翻译; 美国数学学会:普罗维登斯,RI,美国,2000年;第191卷。[谷歌学者]
  10. Eguchi,S.曲线指数族中最小对比度估计的二阶效率。安。统计师。 1983,11, 793–803. [谷歌学者][交叉参考]
  11. Eguchi,S.最小对比度几何。广岛数学。J。 1992,22, 631–647. [谷歌学者]
  12. 巴苏,A。;I.R.哈里斯。;霍尔特,N.L。;Jones,M.C.通过最小化密度功率发散实现稳健高效的估计。生物特征 1988,85, 549–559. [谷歌学者][交叉参考]
  13. 信息散度几何及其在统计机器学习中的应用。信息论与统计学习; Emmert-Streib,F.,Dehmer,M.,编辑。;施普林格:美国纽约州纽约市,2008年;第309-332页。[谷歌学者]
  14. 米纳米。;Eguchi,S.通过β-发散实现鲁棒盲源分离。神经计算。 2002,14, 1859–1886. [谷歌学者]
  15. Eguchi,S。;与幂函数和统计应用相关的熵和散度。 2010,12, 262–274. [谷歌学者][交叉参考]
  16. 藤泽,H。;Eguchi,S.对严重污染有小偏差的稳健参数估计。《多元分析杂志》。 2008,99, 2053–2081. [谷歌学者][交叉参考]
  17. 阿扎里尼,A。;Genton,M.G.关于高斯正态分布的特征。伯努利 2007,13, 169–174. [谷歌学者][交叉参考]
  18. Teicher,H.分布的最大似然特征。安。数学。统计人员。 1961,32,1214–1222。[谷歌学者][交叉参考]
  19. 阿玛里,S。;Ohara,A.概率分布q指数族的几何。 2011,13,1170-1185。[谷歌学者][交叉参考]
  20. 科茨,S。;南卡罗来纳达拉亚。多元T分布及其应用; 剑桥大学出版社:英国剑桥,2004年。[谷歌学者]
  21. 基于对比泛函的统计推断的微分几何方法。广岛数学。J。 1985,15, 341–391. [谷歌学者]
  22. Wakaki,H。;Eguchi,S。;Fujikoshi,Y.一般协方差结构的一类检验。《多元分析杂志》。 1990,32, 313–325. [谷歌学者][交叉参考]
  23. 格伦瓦尔德,P.D。;Dawid,A.P.博弈论,最大熵,最小差异,稳健贝叶斯决策理论。安。统计师。 2004,32, 1367–1433. [谷歌学者]
  24. Kent,J.T。;Tyler,D.E.重新定义多元位置和散布的M估计。安。统计师。 1991,19, 2102–2119. [谷歌学者][交叉参考]
  25. Marrona,R.A.多元位置和散布的稳健M估计。安。统计师。 1976,4, 51–67. [谷歌学者][交叉参考]
  26. Eguchi,S.信息几何与统计模式识别。Sugaku博览会 2006,19, 197–216. [谷歌学者]

附录1

我们展示了(5). 它来自l'H o个 ^ 皮塔尔的规则是
极限 γ 0 D类 γ ( , (f) ) = γ ( x个 ) 1 + γ d日 x个 1 1 + γ - ( x个 ) (f) ( x个 ) γ d日 x个 (f) ( x个 ) 1 + γ d日 x个 γ 1 + γ γ = 0
写为
1 1 + γ ( x个 ) 1 + γ d日 x个 - γ 1 + γ ( x个 ) 1 + γ 日志 ( x个 ) d日 x个 - ( x个 ) (f) ( x个 ) γ 日志 (f) ( x个 ) d日 x个 (f) ( x个 ) 1 + γ d日 x个 γ 1 + γ + γ 1 + γ ( x个 ) (f) ( x个 ) γ d日 x个 (f) ( x个 ) 1 + γ d日 x个 1 + 2 γ 1 + γ (f) ( x个 ) 1 + γ 日志 (f) ( x个 ) d日 x个 γ = 0
减少到
( x个 ) 日志 ( x个 ) d日 x个 - ( x个 ) 日志 (f) ( x个 ) d日 x个
这就完成了(5).

附录2

首先,我们给出以下公式 c(c) γ 英寸(22)何时 γ > 0 .让
= 1 det(探测) ( 2 π ) 1 2 1 - ω ( x个 - μ ) T型 - 1 ( x个 - μ ) + 1 γ d日 x个
哪里 ω = γ 2 + d日 γ + 2 γ . 积分被重写为
= ( 2 π ω ) - d日 2 ( 1 - T型 ) + 1 γ d日
哪里 = ( ω ) 1 2 - 1 / 2 ( x个 - μ ) 。用极坐标表示为
= ( 2 π ω ) - d日 2 S公司 d日 - 1 0 1 ( 1 - 第页 2 ) 1 γ 第页 d日 - 1 d日 第页
哪里 S公司 d日 - 1 是单位球体的表面积 d日 - 1 维度,即,
S公司 d日 - 1 = 2 π d日 2 Γ ( d日 2 )
由于中的积分(43)用贝塔函数表示,我们有
c(c) γ = - 1 = ( 2 ω ) d日 2 Γ ( 1 + d日 2 + 1 γ ) Γ ( 1 + 1 γ )
其次,当 - 2 d日 + 2 < γ < 0 . 与上面类似的论点
= ( - 2 π ω ) - d日 2 ( 1 + T型 ) 1 γ d日
哪里 = ( - 2 π ω ) 1 / 2 - 1 / 2 ( x个 - μ ) 。以极坐标表示为
= ( - 2 π ω ) - d日 2 S公司 d日 - 1 0 ( 1 + 第页 2 ) 1 γ 第页 d日 - 1 d日 第页
这导致了
c(c) γ = ( - 2 ω ) d日 2 Γ ( - 1 γ ) Γ ( - 1 γ - d日 2 )

分享和引用

MDPI和ACS样式

Eguchi,S。;O.小森。;加藤,S。投影功率熵和最大Tsallis熵分布。 2011,13, 1746-1764.https://doi.org/10.3390/e13101746

AMA风格

Eguchi S、Komori O、Kato S。投影功率熵和最大Tsallis熵分布。. 2011; 13(10):1746-1764.https://doi.org/10.3390/e13101746

芝加哥/图拉宾风格

Eguchi、Shinto、Osamu Komori和Shogo Kato。2011.“投影功率熵和最大Tsallis熵分布”13,编号10:1746-1764。https://doi.org/10.3390/e13101746

文章指标

返回页首顶部