配对t检验报错：Error in complete.cases(x,y):参数长度不一致求助

报错原因分析 #

你遇到的问题核心在于配对t检验的前提条件不满足：配对t检验要求用于对比的两组数据必须是一一对应的配对观测，也就意味着两组的样本量必须完全相等。

当你用公式形式 t.test(Diarrhea ~ Max_by_90_H2SPos_12, data = noch4, paired=T) 时，R会自动将Diarrhea拆分成两组：一组是Max_by_90_H2SPos_12=0的180个观测，另一组是Max_by_90_H2SPos_12=1的57个观测。这两个组的长度（180 vs 57）不相等，直接触发了complete.cases(x, y)的报错——因为配对检验要求每一个观测都有对应的配对项。

针对性解决方法 #

根据你的研究设计，分两种情况处理：

情况1：你确实是配对设计（同一研究对象的两次测量）

如果Max_by_90_H2SPos_12=0和1是同一研究对象的两种状态测量（比如干预前/后），那你的数据结构大概率是错的：

• 正确的结构应该是宽格式：每个研究对象占一行，包含两个变量（比如Diarrhea_0对应状态0的测量，Diarrhea_1对应状态1的测量），此时直接对两个变量做配对t检验：

  t.test(noch4$Diarrhea_0, noch4$Diarrhea_1, paired = TRUE)
  

• 如果你的数据是长格式（每个研究对象两行数据，分别对应两种状态），那说明存在部分对象只有单状态的测量数据，需要先筛选出配对完整的样本：

  # 假设你的数据有ID列标记每个研究对象
  library(dplyr)
  # 只保留同时有两种状态测量的对象
  matched_data <- noch4 %>%
    group_by(ID) %>%
    filter(n() == 2) %>%
    ungroup()
  # 再运行配对t检验
  t.test(Diarrhea ~ Max_by_90_H2SPos_12, data = matched_data, paired = TRUE)
  

情况2：你混淆了配对t检验和独立样本t检验

如果Max_by_90_H2SPos_12=0和1是两组独立的研究对象（不是同一对象的两次测量），那你应该使用独立样本t检验，去掉paired=T参数即可：

t.test(Diarrhea ~ Max_by_90_H2SPos_12, data = noch4)

验证步骤 #

你可以先手动确认两组数据的长度，验证报错的直接原因：

# 提取两组数据
group0 <- noch4$Diarrhea[noch4$Max_by_90_H2SPos_12 == 0]
group1 <- noch4$Diarrhea[noch4$Max_by_90_H2SPos_12 == 1]
# 查看长度
length(group0) # 输出180
length(group1) # 输出57

这两个长度不等，就是导致报错的直接原因。

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