在R语言中合并APOE e4状态评分以生成交互项的方法咨询

针对你需要将APOE e4评分（0=无e4等位基因，1=1个e4等位基因，2=2个e4等位基因）重编码为二分类变量（0=无e4，1=存在e4）的需求，以下是几种简洁的实现方式：

方法1：Base R的ifelse()函数 #

这是最直观的条件赋值方法：

# 假设原变量为apoe_e4，创建新二分类变量apoe_e4_binary
apoe_e4_binary <- ifelse(apoe_e4 == 0, 0, 1)

逻辑：原变量等于0时保留0，所有非0值（1或2）统一转为1，同时自动保留原变量中的缺失值（NA）。

方法2：逻辑表达式转数值 #

利用R中TRUE对应1、FALSE对应0的特性，用逻辑判断直接生成二分类变量，代码更简洁：

apoe_e4_binary <- as.integer(apoe_e4 != 0)

逻辑：判断原变量是否不等于0，得到布尔结果后转为整数类型，同样会保留原数据中的NA。

方法3：tidyverse的case_when()（适合数据框操作） #

如果你用dplyr等tidyverse工具处理数据框，case_when()可以明确指定规则，还能处理异常值：

library(dplyr)

# 假设数据框为df，新增二分类变量
df <- df %>%
  mutate(apoe_e4_binary = case_when(
    apoe_e4 == 0 ~ 0,          # 无e4等位基因
    apoe_e4 %in% c(1, 2) ~ 1, # 存在1个或2个e4等位基因
    TRUE ~ NA_real_            # 处理其他异常值（如NA或无效数值）
  ))

测试示例 #

你可以用测试数据验证效果：

test_apoe <- c(0, 1, 2, 0, 2, 1, NA)
# 方法1结果
ifelse(test_apoe == 0, 0, 1)
# 输出：0 1 1 0 1 1 NA

# 方法2结果
as.integer(test_apoe != 0)
# 输出：0 1 1 0 1 1 NA

# 方法3结果
case_when(
  test_apoe == 0 ~ 0,
  test_apoe %in% c(1,2) ~1,
  TRUE ~ NA_real_
)
# 输出：0 1 1 0 1 1 NA

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