R循环报错(undefined columns selected)：前半段正常后半段异常

需求说明 #

需要根据原始分数Araw、性别（1=男，2=女）和年龄，计算标准化t分数：

• 男性的标准化分数存储在数据框male_a中，女性的存储在female_a中
• 两个数据框的第1列对应30岁以下的分数，第2列对应30岁及以上的分数

报错与异常现象 #

运行循环后出现报错：

Error in [.data.frame(female_a, ifelse(data$Araw[[i]] == 0, : undefined columns selected

进一步检查发现：

• tscore_A列表的前91行结果正确
• 从第91行开始，存储的是行号，结果异常
• 手动代入第91行数据到ifelse函数中能得到正确结果，但循环无法正常运行

循环代码 #

tscore_A <- as.list(seq(1, nrow(data))) # 创建空列表存储循环结果

for (i in 1:nrow(data)){
  
  tscore_A[[i]] <- ifelse(data$gender[[i]] == 1, 
          male_a[ifelse(data$Araw[[i]] == 0, 0, data$Araw[[i]] + 1), ifelse(data$age[[i]] < 30, 1, 2)], 
         female_a[ifelse(data$Araw[[i]] == 0, 0, data$Araw[[i]] + 1), ifelse(data$age[[i]] < 30, 1, 2)])

}

数据结构 #

主数据框data #

dput(data)
structure(list(gender = c(1, 2, 2, 1, 1, 2, 2, 1, 1, 1, 2, 
1, 2, 2, 2, 2, 2, 1, 1, 2, 1, 2, 2, 1, 1, 1, 2, NA, 2, 2, 2, 
2, 2, 1, 1, 2, 1, 1, 2, NA, 2, 1, 2, 2, 2, 1, 1, 1, 2, 1, 2, 
1, 2, NA, 2, NA, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 1, 1, 2, 2, 2, 
2, 1, 2, 2, 2, 2, 1, 2, 2, 1, 2, 2, 2, 1, 2, 2, 1, 2, 1, 2, 1, 
2, 2, 2, 1, 2, 2, 1, 2, 2, 2, 1, 2, 1, 1, 2, 1, 1, 1, 2, NA, 
1, 1, 1, 2, 2, 2, 1, 1, 1, 1, 2, 2, 1, 2, 1, 1, 2, 1, 2, 1, 2, 
1, 1, NA, NA, 1, 1, 1, 2, 2, 2, 1, 2, 1, 2, 2, 1, 4, 1, 2, 2, 
2, 1, 2, 1, 1, 2, 2, 2, 1, 2, 2, 2, 2, 1, 1, 1, 1, 1, 1, 2, 2, 
2, 1, 2, NA, 1, 2, 1, 2, 2, 1, 1, 1, 2, 2, NA, 2, 2, 1, 2, 1, 
1, 1, 1, 2, 2, NA, 1, 2, 1, 2, 2, 2, 2, 1, 2, 1, NA, NA, 2, 2, 
NA, 2, 2, NA, 2, 2, 1, 1, 2, 1, 2, 2, NA, NA, 2, 2, 1, NA, NA, 
NA, NA, 1, 1, 1, NA, 1, NA, NA, NA, NA, NA, 1), age = c(29, 30, 
34, 26, 38, 33, 29, 32, 23, 31, 30, 28, 40, 29, 31, 26, 32, 28, 
33, 27, 26, 28, 25, 23, 32, 27, 30, NA, 25, 32, 27, 26, 31, 29, 
30, 30, 34, 29, 28, 23, 27, 30, 21, 29, 29, 31, 36, 32, 33, 30, 
37, 29, 23, 20, 29, 24, 31, 33, 29, 31, 23, 35, 26, 25, 23, 28, 
31, 31, 30, 28, 25, 25, 31, 31, 29, 24, 35, 26, 27, 27, 25, 25, 
33, 34, 30, 25, 31, 31, 25, 31, NA, 29, 27, 30, 29, 32, 26, 24, 
25, 27, 30, 30, 33, 29, 27, 25, 29, 30, 27, 28, 29, 27, 29, 29, 
26, 32, 36, 30, 31, 29, 26, 30, 26, 27, 28, 33, 30, 27, 28, 27, 
28, 30, 26, 25, 27, 25, 22, 35, 28, 25, 23, 37, 26, 32, 35, 30, 
27, 34, 28, 23, 32, 26, 34, 27, 34, 30, 28, 29, 25, 29, 28, 35, 
27, 25, 26, 20, 36, 27, 25, 27, 29, NA, 28, 25, 23, 33, 22, 24, 
29, 28, 25, 28, 25, 26, 33, 28, 29, 29, 23, 43, 23, 31, 31, 38, 
28, 27, 26, 28, NA, NA, 28, 23, 26, 23, 20, 31, NA, 30, 24, 30, 
26, 34, 26, 27, 28, 35, 29, 28, 27, 26, 35, 30, 35, 25, 30, 36, 
36, 25, 27, 26, 26, 31, 27, 34, 31, 32, 31, 32, 23, 32, 33, 32, 
NA, NA, 29, 33), Araw = c(6, 0, 9, 7, 7, 6, 4, 12, 13, 
11, 1, 10, 4, 6, 5, 3, 9, 4, 6, 4, 10, 3, 4, 5, 5, 5, 0, 9, 6, 
8, 10, 7, 7, 4, 10, 5, 12, 6, 8, NA, 6, 5, 3, 5, 5, 5, 5, 7, 
0, 8, 4, 9, 3, 8, 5, 10, 2, 10, 3, 5, 10, 7, 8, 5, 5, 8, 6, 8, 
8, 5, 4, 5, 6, 0, 8, 6, 4, 5, 5, 5, 4, 5, 6, 13, 12, 4, 3, 0, 
7, 7, 5, 5, 1, 3, 3, 6, 1, 4, 4, 6, 6, 5, 6, 2, 9, 1, 5, 13, 
4, 7, 8, 12, 7, 6, 6, 5, 6, 4, 3, 2, 6, 11, 4, 6, 5, 3, 7, 2, 
6, 10, 4, 8, 8, 6, 9, NA, NA, 9, 6, 8, 8, 7, 4, 4, 4, 8, 11, 
9, 6, 9, 6, 0, 5, 5, 5, 5, 7, 0, 1, 7, 10, 7, 5, 4, 4, 11, 5, 
12, 2, 10, 13, 8, 6, 7, 5, 8, 8, NA, 3, 9, 1, 5, 9, 11, 11, 4, 
5, 7, NA, 12, 5, 6, 2, 7, 12, 4, 8, 1, 11, NA, 5, 8, 6, 7, 10, 
8, 9, 6, 5, 9, NA, NA, 4, 7, NA, 8, 6, NA, 11, 5, 7, 3, 1, 6, 
7, 6, 10, NA, 5, 6, 9, NA, NA, NA, NA, 6, 11, 7, 7, 6, NA, NA, 
NA, 4, NA, 7)), class = "data.frame", row.names = c(NA, -246L
))

男性标准化分数数据框male_a #

dput(male_a)
structure(list(X1 = c(35L, 37L, 39L, 42L, 45L, 48L, 51L, 54L, 
57L, 60L, 63L, 66L, 69L, 72L, 75L, 78L), X2 = c(36L, 40L, 42L, 
45L, 49L, 52L, 55L, 59L, 62L, 66L, 69L, 73L, 76L, 79L, 83L, 86L
)), class = "data.frame", row.names = 0:15)

女性标准化分数数据框female_a #

dput(female_a)
structure(list(X1 = c(36L, 39L, 41L, 44L, 47L, 51L, 54L, 57L, 
61L, 64L, 67L, 70L, 74L, 77L, 80L, 83L), X2 = c(36L, 39L, 42L, 
45L, 49L, 52L, 56L, 59L, 63L, 66L, 70L, 73L, 77L, 80L, 84L, 87L
)), class = "data.frame", row.names = 0:15)

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