使用WRS2包pbad2way函数报错：'list'无法转为integer的求助

问题确认 #

你的猜想完全正确——报错确实是由因子组合存在空单元格导致的。pbad2way函数在处理不平衡设计（部分因子组合无观测值）时，内部的序列生成逻辑会因缺失分组映射抛出类型转换错误。

报错重现信息 #

你运行的代码：

pbad2way(avg_rad ~ antifungal + dose + antifungal:dose, data = y2, est = "mom", nboot = 5000) 

报错内容：

Error in sequence.default(piece, ...) : 'list' object cannot be coerced to type 'integer'

你的数据因子组合分布：

解决办法 #

方法1：移除无效因子组合 #

直接清理数据中无观测值的因子水平，让函数只处理有数据的分组：

# 自动删除无数据的因子水平
y2$antifungal <- droplevels(y2$antifungal)
y2$dose <- droplevels(y2$dose)

# 重新运行函数
pbad2way(avg_rad ~ antifungal + dose + antifungal:dose, data = y2, est = "mom", nboot = 5000)

方法2：替换为兼容不平衡设计的函数 #

如果需要保留完整因子水平，可以改用WRS2包中对不平衡设计更友好的bwtrim函数，它支持稳健的双向分析：

bwtrim(avg_rad ~ antifungal * dose, data = y2, est = "mom", nboot = 5000)

方法3：补全空单元格（应急方案，不推荐） #

若必须使用pbad2way，可给空组合添加虚拟数据（会引入偏差，谨慎使用）：

library(dplyr)
# 生成空组合的虚拟行
empty_groups <- expand.grid(antifungal = levels(y2$antifungal),
                           dose = levels(y2$dose)) %>%
  anti_join(y2, by = c("antifungal", "dose")) %>%
  mutate(avg_rad = min(y2$avg_rad) * 0.1) # 用极小值填充

# 合并数据并重新运行
y2_full <- bind_rows(y2, empty_groups)
pbad2way(avg_rad ~ antifungal + dose + antifungal:dose, data = y2_full, est = "mom", nboot = 5000)

内容的提问来源于stack exchange，提问作者mr.raccoon