问询SPSS中ANCOVA分析的计算公式及F统计量推导方法

一、ANCOVA核心逻辑 #

ANCOVA是方差分析的延伸，核心是先扣除**协变量（连续型控制变量）**对因变量的线性影响，再分析分组因素的真实效应。最终通过比较「扣除协变量后的组间变异」与「组内变异」的比值，得到F统计量来判断分组效应是否显著。

二、手动计算F统计量的代码（无预构建ANCOVA模块） #

以下用纯NumPy实现，每一步对应ANCOVA的推导逻辑：

import numpy as np

# 模拟示例数据：3组，每组20个样本，1个协变量
np.random.seed(42)
n_groups = 3
n_per_group = 20
total_n = n_groups * n_per_group

# 协变量x
x = np.random.normal(50, 10, total_n)
# 分组变量（编码为0、1、2）
group = np.repeat([0,1,2], n_per_group)
# 因变量y：受分组、协变量影响，添加随机噪声
y = 2*group + 0.5*x + np.random.normal(0, 3, total_n)

# 1. 计算总平方和与总自由度
y_mean = np.mean(y)
SS_total = np.sum((y - y_mean)**2)
df_total = total_n - 1

# 2. 拟合仅含协变量的回归模型，计算协变量解释的平方和
X_cov = np.column_stack((np.ones(total_n), x))  # 加截距项
beta_cov = np.linalg.inv(X_cov.T @ X_cov) @ X_cov.T @ y
y_hat_cov = X_cov @ beta_cov
SS_cov = np.sum((y_hat_cov - y_mean)**2)

# 3. 拟合全模型（协变量+分组），计算总解释平方和
# 构造分组哑变量（避免多重共线性，取n_groups-1个）
group_dummies = np.zeros((total_n, n_groups-1))
for i in range(n_groups-1):
    group_dummies[:, i] = (group == i).astype(int)
X_full = np.column_stack((np.ones(total_n), x, group_dummies))
beta_full = np.linalg.inv(X_full.T @ X_full) @ X_full.T @ y
y_hat_full = X_full @ beta_full
SS_full = np.sum((y_hat_full - y_mean)**2)

# 4. 计算分组因素的平方和与均方
SS_group = SS_full - SS_cov  # 扣除协变量后，分组解释的变异
df_group = n_groups - 1
MS_group = SS_group / df_group

# 5. 计算误差项的平方和与均方
SS_error = SS_total - SS_full  # 未被模型解释的变异
df_error = total_n - n_groups - 1  # 总样本 - 组数 - 协变量个数
MS_error = SS_error / df_error

# 6. 计算F统计量
F_stat = MS_group / MS_error

print(f"分组因素的F统计量: {F_stat:.4f}")
print(f"自由度（组间/误差）: {df_group}/{df_error}")

三、SPSS中ANCOVA的计算公式 #

SPSS默认采用Type III调整平方和（适用于平衡/不平衡数据），核心公式与手动推导一致，具体如下：

1. 分组因素的调整平方和 #

SS(分组 | 协变量) = SS(全模型) - SS(仅协变量模型)
即扣除协变量的影响后，分组对因变量的解释变异，和手动计算的SS_group完全等价

2. 误差平方和 #

SS(误差) = SS(总) - SS(全模型)

3. 均方与F统计量 #

MS(分组) = SS(分组 | 协变量) / (k-1)，其中k为分组数
MS(误差) = SS(误差) / (N - k - p)，其中N为总样本量，p为协变量个数
F = MS(分组) / MS(误差)

SPSS的Type III平方和会确保在控制其他变量（协变量）的前提下，单独计算每个因素的效应，这是ANCOVA分析中的标准选择，能避免样本量不平衡带来的偏差。

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