在R语言中将死亡率转换为概率时遇报错，求解决方案

核心问题 #

报错Error in -r : invalid argument to unary operator说白了就是你要取负的对象不是数值型数据，结合你提供的数据显示情况（向量里存在缺失值，或数据类型不符合要求），多半是以下两种情况：

• 1. v_r_HDage不是数值型向量：比如它是字符型、因子型，甚至是列表结构（而非纯原子向量）
• 2. 向量中混入非数值内容：比如存在乱码字符、格式错误的缺失值

排查与解决步骤 #

• 先检查数据类型：运行class(v_r_HDage)，如果结果是character或factor，需要转换为数值型：

  v_r_HDage <- as.numeric(as.character(v_r_HDage))
  

  （因子类型必须先转字符再转数值，避免因子编码被误当成数值）
• 检查缺失值与异常值：运行summary(v_r_HDage)或is.na(v_r_HDage)，查看是否存在NA或非数值内容。如果有NA，转概率时可选择保留NA（默认行为），或用na.rm=TRUE处理（若rate_to_prob支持该参数，比如epiR包中的版本）：

  v_p_HDage <- 1 - exp(-v_r_HDage)
  # 使用rate_to_prob的写法：
  v_p_HDage <- rate_to_prob(v_r_HDage, na.rm = TRUE)
  

• 确认向量结构：运行is.vector(v_r_HDage)，如果返回FALSE，说明它可能是数据框的列或列表，需转换为纯原子向量：

  # 示例：如果是数据框df中的col列
  v_r_HDage <- as.vector(df$col)
  

快速验证方法 #

先取向量前几个值测试，缩小问题范围：

test_vec <- v_r_HDage[1:5]
class(test_vec)
1 - exp(-test_vec)

如果这段代码能正常运行，说明问题出在向量后续的某个异常值上，逐一排查即可。

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