如何在R中从ordiR2step结果（列表的列表）提取调整后R2？

在使用vegan包的ordiR2step做模型选择时，控制台会输出调整后R²（R2.adjusted），但这个统计量并没有直接存放在test对象的列表结构里。要提取并整理成指定格式的数据框，可以通过捕获并解析test的打印输出来实现，具体步骤如下：

1. 捕获打印输出 #

用capture.output()函数把print(test)的控制台输出内容保存下来：

output <- capture.output(print(test))

2. 解析输出内容，提取变量和对应R² #

通过正则表达式匹配和提取目标内容，整理成数据框：

# 筛选包含R2.adjusted的行
r2_lines <- grep("R2.adjusted", output, value = TRUE)
# 筛选包含变量信息的行（Start/Step/All variables）
var_lines <- grep("^(Start|Step|<All variables>)", output, value = TRUE)

# 处理变量名称
variables <- sapply(var_lines, function(line) {
  if (grepl("Start", line)) {
    "<none>"
  } else if (grepl("Step", line)) {
    # 提取Step行里的变量名，转为"+ 变量名"格式
    sub("Step \\d+:.*~ 1 \\+ (.*)", "+ \\1", line)
  } else if (grepl("<All variables>", line)) {
    "<All variables>"
  }
})

# 提取R2.adjusted的数值
r2_adjusted <- as.numeric(sub(".*R2.adjusted= ([0-9.]+)", "\\1", r2_lines))

# 合并成目标数据框
result_df <- data.frame(
  variables = variables,
  R2.adjusted = r2_adjusted,
  stringsAsFactors = FALSE
)

# 调整行顺序，匹配示例格式
result_df <- result_df[c(2, 3, 4, 1), ]

运行这段代码后，result_df就是你需要的结构：

variables R2.adjusted
2          + env1  0.08636035
3 <All variables>  0.08220187
4          + env2  0.00842866
1           <none>  0.00000000

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