We consider the lossless compression bound of any single data sequence. If we fit the data by a parametric model, the entropy quantity $nH({\hat \theta}_n)$ obtained by plugging in the maximum likelihood estimate is an underestimate of the bound, where $n$ is the number of words. Shtarkov showed that the normalized maximum likelihood (NML) distribution or code length is optimal in a minimax sense for any parametric family. We show by the local asymptotic normality that the NML code length for the exponential families is $nH(\hat \theta_n) +\frac{d}{2}\log \, \frac{n}{2\pi} +\log \int_{\Theta} |I(\theta)|^{1/2}\, d\theta+o(1)$, where $d$ is the model dimension or dictionary size, and $|I(\theta)|$ is the determinant of the Fisher information matrix. We also demonstrate that sequentially predicting the optimal code length for the next word via a Bayesian mechanism leads to the mixture code, whose pathwise length is given by $nH({\hat \theta}_n) +\frac{d}{2}\log \, \frac{n}{2\pi} +\log \frac{|\, I({\hat \theta}_n)|^{1/2}}{w({\hat \theta}_n)}+o(1) $, where $w(\theta)$ is a prior. The asymptotics apply to not only discrete symbols but also continuous data if the code length for the former is replaced by the description length of the latter. The analytical result is exemplified by calculating compression bounds of protein-encoding DNA sequences under different parsing models. Typically, the highest compression is achieved when the parsing is in phase of the amino acid codons. On the other hand, the compression rates of pseudo-random sequences are larger than 1 regardless parsing models. These model-based results are in consistency with that random sequences are incompressible as asserted by the Kolmogorov complexity theory. The empirical lossless compression bound is particularly more accurate when dictionary size is relatively large.
翻译:我们考虑任意单一数据序列的无损压缩界限。若采用参数模型拟合数据,则通过代入最大似然估计得到的熵量 $nH({\hat \theta}_n)$ 会低估该界限,其中 $n$ 为数据词条数目。Shtarkov 证明,归一化最大似然分布或码长在极小化极大意义上对于任意参数族是最优的。本文借助局部渐近正态性证明,指数族的 NML 码长为 $nH(\hat \theta_n) +\frac{d}{2}\log \, \frac{n}{2\pi} +\log \int_{\Theta} |I(\theta)|^{1/2}\, d\theta+o(1)$,其中 $d$ 为模型维度或词典大小,$|I(\theta)|$ 为 Fisher 信息矩阵的行列式。我们还论证了通过贝叶斯机制对下一个词条的最优码长进行序贯预测会得到混合码,其路径码长为 $nH({\hat \theta}_n) +\frac{d}{2}\log \, \frac{n}{2\pi} +\log \frac{|\, I({\hat \theta}_n)|^{1/2}}{w({\hat \theta}_n)}+o(1)$,其中 $w(\theta)$ 为先验分布。该渐近性不仅适用于离散符号数据——当将离散符号的码长替换为连续数据的描述长度时,同样适用于连续型数据。我们通过计算不同解析模型下蛋白质编码 DNA 序列的压缩界限对分析结果进行例证。典型情况下,当解析处于氨基酸密码子相位时压缩率最高;而伪随机序列的压缩率无论采用何种解析模型均大于1。这些基于模型的结果与柯尔莫哥洛夫复杂度理论中"随机序列不可压缩"的论断一致。该经验性无损压缩界限在词典尺寸较大时尤为准确。