使用索引的嵌套for循环未覆盖数据帧，求助列值出现次数计算

嘿，我来帮你搞定这个嵌套for循环的问题！先让我重现你的数据框，确保我们在同一个语境下：

ID<-seq(1:5) #patient ID
snp1<-c("A","T","A","A","T")
snp2<-c("C","C","0","C","C")
snp3<-c("A","G","A","A","G")
snp4<-c("T","0","C","G","T")
snp5<-c("G","G","G","G","A")
dat<-data.frame(ID,snp1,snp2,snp3,snp4,snp5)
print(dat)

输出和你给出的完全一致：

ID snp1 snp2 snp3 snp4 snp5
1  1    A    C    A    T    G
2  2    T    C    G    0    G
3  3    A    0    A    C    G
4  4    A    C    A    G    G
5  5    T    C    G    T    A

你大概率踩了这些嵌套循环的坑 #

用嵌套for循环处理这类问题时，新手常碰到的问题不外乎这几个：

• 没选对列范围：比如不小心把ID列也包含进去了，或者列索引起始值搞错了
• 没提前准备存储结果的容器，导致每次循环都覆盖之前的计数
• 遍历值的时候逻辑出错，比如没拿到当前列的所有唯一值，或者计数时误把"0"当成了NA处理

修复后的嵌套for循环实现 #

假设你想计算每列中每个不同值的出现次数（包括那个代表缺失的"0"），下面是一个能正常跑通的嵌套循环写法：

# 先建个空列表存每列的计数结果
count_results <- list()

# 遍历从snp1到snp5的列（索引从2到6，因为第1列是ID）
for (col_idx in 2:ncol(dat)) {
  # 先拿到当前列的名字和数据
  col_name <- colnames(dat)[col_idx]
  col_data <- dat[[col_idx]]
  
  # 获取当前列的所有唯一值
  unique_vals <- unique(col_data)
  
  # 给当前列的计数建个空向量
  col_counts <- c()
  
  # 遍历每个唯一值，计算出现次数
  for (val in unique_vals) {
    count <- sum(col_data == val, na.rm = FALSE) # 这里不用na.rm=TRUE，因为你的"0"是字符，不是缺失值NA
    col_counts[val] <- count
  }
  
  # 把当前列的计数放进结果列表
  count_results[[col_name]] <- col_counts
}

# 看看最终结果
print(count_results)

运行后你会得到每个SNP列的各值计数：

$`snp1`
A T 
3 2 

$`snp2`
C 0 
4 1 

$`snp3`
A G 
3 2 

$`snp4`
T 0 C G 
2 1 1 1 

$`snp5`
G A 
4 1 

更简洁的替代方案（不用嵌套循环） #

其实在R里，完全没必要用嵌套for循环，用lapply结合table函数就能一行搞定，代码更简洁还不容易出错：

# 对ID以外的所有列应用table函数
count_results_simple <- lapply(dat[, -1], table)
print(count_results_simple)

这个结果和上面的嵌套循环完全一样，而且更符合R的编程风格。

如果只是想统计某个特定值的出现次数 #

比如你只想知道每列中"A"出现了多少次，那嵌套循环可以简化成这样：

target_val <- "A"
specific_counts <- c()

for (col_idx in 2:ncol(dat)) {
  col_name <- colnames(dat)[col_idx]
  count <- sum(dat[[col_idx]] == target_val, na.rm = FALSE)
  specific_counts[col_name] <- count
}

print(specific_counts)

输出结果：

snp1 snp2 snp3 snp4 snp5 
   3    0    3    0    1 

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