R语言mutate函数报错：f需为因子或字符向量求助

报错里的f是什么？ #

报错里的f是fct_relevel()函数的第一个参数的官方名称，在这里就是你传入的type变量。函数报错时会用参数名指代你传入的具体变量，所以它说的f就是你刚生成的type字段。

为什么会触发这个错误？ #

核心问题出在type的生成与后续处理的衔接上：

• 你用嵌套ifelse()生成的type是字符向量，虽然fct_relevel()理论上支持字符向量，但如果当前数据集中type的取值没有覆盖你指定的所有重排类别（比如你的数据里没有bad-hit [evalue > 0.01]这个情况），就会触发这个看似“类型不匹配”的误报（本质是指定的重排水平不存在于当前向量中）。
• 另外，嵌套ifelse()的逻辑如果出现意外的NA值，也可能导致type的类型混乱，不过结合你之前数据集能正常运行的情况，更大概率是当前数据集缺少了某个type类别。

修复代码的两种实用方式 #

方式1：先显式转为因子，再重排 #

先把type转为包含所有可能水平的因子，再进行重排（其实指定levels后顺序已经符合要求，重排步骤可省略）：

B1000_SEQ <- BrachCDS %>% 
  filter(gene %in% B1000$V1) %>% 
  mutate(length=nchar(seq)) %>% 
  left_join(blasts, by=c('gene'='qseqid')) %>% 
  mutate(
    type=ifelse(is.na(evalue),'non-hit', ifelse(evalue <= 0.01,'hit','bad-hit [evalue > 0.01]')),
    # 显式指定所有可能的水平，确保覆盖重排需要的类别
    type=factor(type, levels=c('hit','bad-hit [evalue > 0.01]','non-hit')),
    # 若已经指定levels，这一步重排可以省略
    type=fct_relevel(type,c('hit','bad-hit [evalue > 0.01]','non-hit'))
  )

方式2：用case_when()替代嵌套ifelse()，逻辑更清晰 #

case_when()的可读性更强，也能更稳定地生成符合要求的字符向量：

B1000_SEQ <- BrachCDS %>% 
  filter(gene %in% B1000$V1) %>% 
  mutate(length=nchar(seq)) %>% 
  left_join(blasts, by=c('gene'='qseqid')) %>% 
  mutate(
    type=case_when(
      evalue <= 0.01 ~ 'hit',
      evalue > 0.01 ~ 'bad-hit [evalue > 0.01]',
      is.na(evalue) ~ 'non-hit'
    ),
    type=fct_relevel(type,c('hit','bad-hit [evalue > 0.01]','non-hit'))
  )

为什么之前的数据集能正常运行？ #

因为之前的数据集里，type的三个类别（hit/bad-hit [evalue > 0.01]/non-hit）都存在，ifelse()生成的字符向量包含了所有你要重排的水平，fct_relevel()可以正常识别并处理。而当前数据集缺少了其中某个类别，导致字符向量的水平不全，触发了报错。

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