R语言vtable包sumtable()函数不可复现错误排查求助

问题背景 #

使用vtable包的sumtable()函数时，添加factor.numeric=TRUE参数仅在部分数据集上报错，同事运行相同代码和数据无问题。已尝试卸载重装R、RStudio及所有包，切换数据格式（csv转dta）仍无法解决。

可正常运行的示例（mtcars数据集） #

# Attempt 1
names <- c('cyl', 'vs', 'am', 'gear', 'carb')
mtcars[,names] <- lapply(mtcars[,names] , factor) 
mtcars2<-mtcars
mtcars2$cyl<-factor(mtcars2$cyl)
sumtable(mtcars2, c('mpg','cyl','disp'), factor.numeric = TRUE)

报错示例（世界价值观调查数据集） #

# Attempt 2
WVS <- read_delim("[my path]/WV5_Data_csv_v20180912.csv", delim = ";", escape_double = FALSE, trim_ws = TRUE)
WVS <- replace(WVS, WVS < 0, NA)
WVS2<-WVS %>% select('V147', 'V91','V95','V123','V235','V255')
WVS2$V255<-factor(WVS2$V255)
sumtable(WVS2 %>% select(V255, V95), factor.numeric = TRUE)

报错信息：

Error in model.frame.default(object, data, xlev = xlev) : 
  invalid type (list) for variable 'va'

移除factor.numeric=TRUE参数后代码可正常运行：

# Attempt 3
sumtable(WVS2 %>% select(V255, V95))

排查步骤与解决方案 #

1. 核对包版本差异 #

vtable包的不同版本可能存在bug修复，先确认你和同事的vtable版本是否一致：

packageVersion("vtable")

若版本不同，升级到最新版（或和同事一致的版本）：

install.packages("vtable")

2. 检查变量类型与隐藏属性 #

报错提示变量为list类型，即使你已将V255转为factor，仍可能存在隐藏的list结构或属性：

• 查看变量结构：

  str(WVS2$V255)
  str(WVS2$V95)
  

  确认输出中无List of ...的标识，若存在需先清理list元素（比如用unlist()处理）。
• 移除tibble的附加属性（read_delim默认返回tibble，可能携带spec属性干扰函数运行）：

  WVS_simple <- as.data.frame(WVS2 %>% select(V255, V95))
  sumtable(WVS_simple, factor.numeric = TRUE)
  

3. 检查factor水平的异常值 #

若V255的factor水平包含NA、空值或非标准格式，可能导致函数转换失败：

• 查看并清理factor水平：

  levels(WVS2$V255)
  # 移除空水平或NA水平
  WVS2$V255 <- droplevels(WVS2$V255)
  # 重新尝试
  sumtable(WVS2 %>% select(V255, V95), factor.numeric = TRUE)
  

4. 核对R版本差异 #

不同R版本的底层函数行为可能存在差异，确认你和同事的R版本是否一致：

R.version.string

若版本差异较大，尝试升级/降级到和同事一致的R版本。

5. 临时替代方案 #

若上述步骤仍无法解决，可手动将factor转为数值后再生成汇总表：

WVS2$V255_num <- as.numeric(WVS2$V255)
sumtable(WVS2 %>% select(V255_num, V95))

内容的提问来源于stack exchange，提问作者J. Holiday