在R Studio中使用read_q2metadata读取QIIME2元数据遇列数不匹配错误

问题背景 #

使用qiime2R包的read_q2metadata函数读取QIIME2格式的元数据文件时，触发列数不匹配错误，但肉眼观察元数据的列名和数据行列数一致。

错误信息 #

metadata_genevstemp <- read_q2metadata("sample_metadata_genevstemp.txt")
Error in read.table(file, header = F, col.names = coltitles, skip = 2,  : 
  more columns than column names

使用代码 #

install.packages(qiime2R)
library(qiime2R)
metadata_genevstemp <- read_q2metadata("sample_metadata_genevstemp.txt")

元数据文件示例 #

SampleID    Temperature Genotype    GxT                                         
#q2:types   Categorical Categorical Categorical                                         
CSkab18C1   18  III-M   III-M18                                         
CSkab18C2   18  III-M   III-M18                                         
# 其余行省略

问题原因 #

错误本质是read_q2metadata内部调用的read.table函数，将文件中行尾多余的空白字符（空格/制表符）解析成了额外列，导致列数和列名不匹配。QIIME2元数据通常为制表符分隔，但如果文件保存时自动添加了行尾空白，就会触发这个问题。

解决方案 #

方案1：传递分隔符和清理参数给read_q2metadata #

直接在read_q2metadata中指定制表符分隔，并开启空白清理：

library(qiime2R)
metadata_genevstemp <- read_q2metadata(
  "sample_metadata_genevstemp.txt",
  sep = "\t",       # 指定制表符为分隔符
  trim = TRUE       # 清理单元格前后空白
)

方案2：手动读取后转换为QIIME2元数据格式 #

如果方案1无效，可先手动读取数据，再转换为QIIME2元数据结构：

library(qiime2R)
# 读取元数据，跳过类型行，指定制表符分隔并清理空白
raw_metadata <- read.delim(
  "sample_metadata_genevstemp.txt",
  header = TRUE,    # 使用第一行作为列名
  skip = 1,         # 跳过第二行的#q2:types类型定义
  sep = "\t",
  trim = TRUE,
  stringsAsFactors = FALSE
)
# 转换为QIIME2元数据格式
metadata_genevstemp <- as_q2metadata(raw_metadata)

方案3：手动清理文件 #

直接打开元数据文件，删除所有行尾的空白字符，以及文件末尾的空行，保存后再用原代码读取。

内容的提问来源于stack exchange，提问作者Sohana Al Sanjee