基于三个相关系数计算变量x总解释方差的方法咨询

1. 能否计算y和z对x的总解释方差？ #

可以。总解释方差即复相关系数的平方 ( R^2_{x.yz} )，完全可以通过给定的两两相关系数推导，不需要原始数据。

计算公式 #

当x作为因变量，y和z作为自变量时，总解释方差的推导公式为：
[
R^2_{x.yz} = \frac{r_{xy}^2 + r_{xz}^2 - 2r_{xy}r_{xz}r_{yz}}{1 - r_{yz}^2}
]
其中：

• ( r_{xy} ) = x与y的相关系数
• ( r_{xz} ) = x与z的相关系数
• ( r_{yz} ) = y与z的相关系数

代入题目中的数值（( r_{xy}=0.71 )，( r_{xz}=0.25 )，( r_{yz}=0.33 )）：

• 分子：( 0.71^2 + 0.25^2 - 2 \times 0.71 \times 0.25 \times 0.33 = 0.5041 + 0.0625 - 0.11715 = 0.44945 )
• 分母：( 1 - 0.33^2 = 1 - 0.1089 = 0.8911 )
• 最终结果：( R^2_{x.yz} \approx 0.44945 / 0.8911 \approx 0.504 )，也就是说y和z能共同解释x约50.4%的方差。

2. 用半偏相关方法得到大于1的结果的情况及公式 #

正常统计场景下，半偏相关系数的绝对值≤1，总解释方差 ( R^2 ) 也必然≤1。但如果给定的两两相关系数组合不满足相关矩阵的正定性（这种组合在现实中不存在对应的变量），代入公式会得到大于1的结果。

核心逻辑 #

三个变量的相关矩阵必须是正定的，即其行列式必须大于0：
[
\det(R) = 1 - r_{xy}^2 - r_{xz}^2 - r_{yz}^2 + 2r_{xy}r_{xz}r_{yz} > 0
]
若 ( \det(R) < 0 )，说明这个相关系数组合在统计上不成立，此时用半偏相关推导总解释方差的公式会得出大于1的结果。

对应计算公式 #

通过半偏相关推导总解释方差的公式和直接计算复相关系数平方的公式一致：
[
R^2_{x.yz} = r_{xy}^2 + r_{xz.y}^2
]
其中 ( r_{xz.y} ) 是z对x的半偏相关系数（控制y的影响后，z与x的相关），其公式为：
[
r_{xz.y} = \frac{r_{xz} - r_{xy}r_{yz}}{\sqrt{1 - r_{yz}^2}}
]
将半偏相关系数代入总解释方差公式，最终会等价于：
[
R^2_{x.yz} = \frac{r_{xy}^2 + r_{xz}^2 - 2r_{xy}r_{xz}r_{yz}}{1 - r_{yz}^2}
]

示例（得到大于1的结果） #

假设取 ( r_{xy}=0.9 )，( r_{xz}=0.9 )，( r_{yz}=0.1 )：

• 先计算行列式：( \det(R) = 1 - 0.9^2 - 0.9^2 - 0.1^2 + 2 \times 0.9 \times 0.9 \times 0.1 = 1 - 0.81 - 0.81 - 0.01 + 0.162 = -0.468 < 0 )，相关矩阵不定
• 代入总解释方差公式：( R^2_{x.yz} = \frac{0.9^2 + 0.9^2 - 2 \times 0.9 \times 0.9 \times 0.1}{1 - 0.1^2} = \frac{0.81 + 0.81 - 0.162}{0.99} = \frac{1.458}{0.99} = 1.47 )，结果大于1

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