R语言报错求助：仅0可与负下标混合（rmean10函数问题）

报错核心原因 #

这个错误的本质是rmean10函数的逻辑存在边界漏洞：当某个被试的有效试次（即函数参数cx的长度）小于10时，length(cx)-10会变成负数，R不允许下标同时混合正数和负数，因此触发报错。

全量200试次数据能正常运行，是因为每个被试的试次远大于10；但拆分到最后45个试次后，可能存在个别被试的实际有效试次不足10个（比如数据缺失、试次记录不全），刚好触发了这个逻辑问题。

解决步骤 #

1. 先排查数据问题
   先统计每个被试的试次数量，确认是否存在试次不足10的情况：

   # 假设你的数据框是df，被试ID列名为subject_id
   table(df$subject_id)
   

   用dplyr可以更直观筛选出异常被试：

   library(dplyr)
   df %>% 
     group_by(subject_id) %>% 
     summarise(trial_count = n()) %>% 
     filter(trial_count < 10)
   

2. 修改rmean10函数，添加边界判断
   给函数增加前置检查，当试次不足10时直接返回提示，避免报错：

   rmean10 <- function(cx) {
     trial_num <- length(cx)
     # 试次不足10时返回NA并给出警告
     if(trial_num < 10) {
       warning(paste("被试试次仅", trial_num, "个，无法计算10次移动均值"))
       return(NA)
     }
     # 保留原计算逻辑
     means <- sapply(1:(trial_num - 10), function(i) mean(cx[i:(i+9)]))
     return(means)
   }
   

3. 预处理数据（可选）
   如果不需要保留试次不足10的被试，可以先过滤掉这类数据：

   filtered_df <- df %>%
     group_by(subject_id) %>%
     filter(n() >= 10) %>%
     ungroup()
   

   再用过滤后的数据运行代码即可。

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