使用dplyr基于多列双条件筛选数据框行的报错问题咨询

嗨，我看到你在尝试用dplyr筛选满足调整后p值<0.05且log2倍数变化>1.3的上调基因时遇到了报错，咱们来一步步解决这个问题~

先分析你的报错原因 #

你写的这段代码：

df %>% filter_at (vars (starts_with ("adj")), any_vars (. < 0.05) && vars (starts_with ("log2")), any_vars (. > 1.30))

问题出在语法逻辑上：filter_at只能针对一组变量应用筛选条件（比如只处理所有adj开头的列），你试图用&&把两组变量的条件硬凑在一起，而any_vars()返回的是dplyr专用的特殊对象，不能直接用逻辑运算符和vars()组合，这就导致了那个类型错误。

而且从dplyr 1.0.0版本开始，官方已经推荐用across()代替filter_at这类“作用域动词”，across()的写法更灵活清晰，完全能满足你的需求。

根据你的需求给出两种解决方案 #

首先明确你的核心需求：筛选上调基因，下面针对两种常见的筛选逻辑给出代码：

场景1：筛选存在至少一个对比组同时满足条件的基因（同一组的adj和log2都达标）

这是更符合“上调差异基因”定义的逻辑——某个时间点的对比中，基因既显著差异（adj_pval<0.05）又上调（log2_ratio>1.3）：

library(dplyr)
library(purrr)

# 先提取所有对比组的后缀（比如"hr6_vs_0"、"day1_vs_0"）
comparisons <- c("hr6_vs_0", "day1_vs_0", "day3_vs_0", "day7_vs_0")

# 遍历每个对比组，检查对应列的条件，只要有一组满足就保留该行
df_filtered <- df %>%
  filter(
    map_lgl(comparisons, ~ 
      (!!sym(paste0("adj_pval_", .x)) < 0.05) & 
      (!!sym(paste0("log2_ratio_", .x)) > 1.3)
    ) %>% any()
  )

场景2：筛选只要有任意一个adj列<0.05，且任意一个log2列>1.3的基因（不要求同一组）

如果你的需求是宽松一点的筛选（只要有一个显著p值，同时有一个上调的倍数，不管是不是同一个时间点），可以用更简洁的写法：

library(dplyr)

df_filtered <- df %>%
  filter(
    # 存在至少一个调整p值<0.05
    any(across(starts_with("adj_pval"), ~ . < 0.05)),
    # 存在至少一个log2倍数变化>1.3
    any(across(starts_with("log2_ratio"), ~ . > 1.3))
  )

验证你的示例数据 #

看你提供的head(df)数据，所有log2_ratio的最大值是0.651（A2M的day1对比），都小于1.3，所以用上面的代码筛选后结果会是空，这符合数据实际情况，但代码逻辑是完全正确的。

备注：内容来源于stack exchange，提问作者ginn