为何在R中残差为0时Ljung-Box检验的Q统计量和p值返回NA

核心原因：残差全0导致自相关系数计算失效 #

Ljung-Box检验的Q统计量完全依赖残差序列的样本自相关系数，其计算公式为：
Q = n(n+2) * Σ(ρ_k²/(n-k))
其中ρ_k是残差的k阶样本自相关系数，求和范围为指定的滞后阶数。

而样本自相关系数的计算公式是：
ρ_k = Σ[(e_t - ē)(e_{t-k} - ē)] / Σ[(e_t - ē)²]
当残差全部为0时，残差均值ē=0，此时分子为Σ(0*0)=0，分母为Σ0²=0——分母为0的除法无意义，直接导致所有ρ_k的计算结果为NA。

Q统计量与p值的连锁失效 #

由于Q统计量完全基于ρ_k计算，当所有ρ_k都是NA时，Q统计量自然无法生成有效数值，返回NA。而p值是通过将Q统计量与对应自由度的卡方分布对比得到的，没有有效的Q值，就无法完成分布映射，最终p值也只能返回NA。

为何不会出现p=0.99999这类结果？ #

p值的计算前提是Q统计量能正常取值并符合卡方分布假设。但残差全0的场景下，统计量的计算基础（非零的残差方差）不存在，根本无法进入卡方分布的映射环节，因此不可能得到接近1的p值。

内容的提问来源于stack exchange，提问作者Mattymat