SPSS中用WRS2做重复测量ANOVA时F值与p值为NA/NaN求助

针对你在SPSS 29中使用WRS2包做单因素重复测量稳健ANOVA时遇到的rmanova输出F值和p值为NA/NaN的问题，可按以下步骤排查解决：

1. 检查重复测量数据的完整性 #

稳健重复测量ANOVA要求每个被试拥有所有时间点的观测值，若存在缺失（部分被试缺少t0/t1/t2中某一个值），会直接导致计算失败。在SPSS的R代码块中加入以下代码检查：

Begin Program R.
# 交叉表查看每个被试的各时间点数据是否存在
missing_table <- table(df$id, df$variable)
print(missing_table)
# 筛选出数据不完整的被试ID
incomplete_subjects <- rownames(missing_table)[rowSums(missing_table == 0) > 0]
print("数据不完整的被试ID：")
print(incomplete_subjects)
End Program.

如果存在不完整被试，可删除这些案例后重新分析：

Begin Program R.
# 保留数据完整的被试
df <- df[!df$id %in% incomplete_subjects, ]
# 重新执行分析
rmanova(df$value, df$variable, df$id, tr = 0.2)
End Program.

2. 调整截尾比例（tr参数） #

你的样本量为27，设置tr=0.2意味着每个时间点会截掉双侧各20%的极端值（约5个样本/侧），剩余有效样本量过少，可能导致计算时的方差矩阵奇异。尝试降低截尾比例，比如tr=0.1（双侧各截10%，约2个样本/侧）：

rmanova(df$value, df$variable, df$id, tr = 0.1)

若必须使用tr=0.2，可尝试使用bootstrap版本的稳健重复测量ANOVA函数rmanova_pb，该函数通过重抽样方法计算p值，对小样本更友好：

rmanova_pb(df$value, df$variable, df$id, tr = 0.2)

3. 排查数据分布与极端值问题 #

当tr=0时输出NaN，可能是数据存在极端值、某时间点所有观测值相同（标准差为0）或方差极度不齐。加入以下代码查看数据分布：

Begin Program R.
# 箱线图查看各组分布与极端值
boxplot(value ~ variable, data = df, main = "PHQ9得分分布")
# 输出各组描述统计
library(dplyr)
summary_stats <- df %>% 
  group_by(variable) %>% 
  summarise(
    样本量 = n(),
    均值 = mean(value, na.rm = TRUE),
    标准差 = sd(value, na.rm = TRUE)
  )
print(summary_stats)
End Program.

如果某组标准差为0，需检查该组数据是否存在录入错误；若存在极端值，可考虑保留合适的截尾比例（而非tr=0）来降低影响。

4. 确保WRS2包为最新版本 #

旧版本WRS2可能存在小样本计算的bug，执行以下代码更新包：

Begin Program R.
install.packages("WRS2", repos = "https://cloud.r-project.org")
library(WRS2)
print("当前WRS2版本：")
print(packageVersion("WRS2"))
End Program.

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