R语言如何在for循环中获取列标签并输出标准差统计结果表

问题场景 #

现有如下示例数据（实际业务场景中数据列数超过12列）：

dat<-read.table (text=" d1  d2  d3  d4  d5  d6  d7  d8  d9  d10 d11 d12
    100 100 60  50  100 100 20  60  100 100 100 20
    75  75  60  50  75  75  20  100 75  75  75  20
    75  40  60  50  41  41  20  100 75  75  75  75
    78  49  60  50  45  41  20  100 75  75  75  75
    75  50  60  50  50  48  20  60  75  75  75  75
    75  48  60  50  44  49  20  40  10  75  75  10
    20  10  60  50  100 100 20  100 100 100 100 20
    75  75  60  50  75  75  20  50  75  39  75  10
    75  75  60  50  75  75  20  90  75  75  75  40
    22  10  60  50  75  75  20  90  75  10  75  10
    ", header=TRUE)

需要输出包含列名、对应列标准差的结果表，目标格式如下：

D   sd
d1  25.44711
d2  29.01647
d3  0
d4  0
d5  22.11586
d6  22.26831
d7  0
d8  23.78141
d9  24.61368
d10 26.89672
d11 10.54093
d12 28.62109

原有循环代码仅能输出标准差数值，无法同步匹配列名：

for(i in seq_along(dat)) {
  k<- dat[[i]]
  if(class(k) %in%  c("numeric", "integer"))  
    print(sd(k))
  else print("F")
}

实现方案 #

方案1：修改原有循环逻辑 #

循环时通过列索引同步取对应列名，提前初始化结果表存储每一步的计算结果，最终直接输出结构化表格：

# 初始化空结果表
result <- data.frame(
  D = character(),
  sd = numeric(),
  stringsAsFactors = FALSE
)

for(i in seq_along(dat)) {
  # 按当前索引取列名
  col_name <- names(dat)[i]
  k <- dat[[i]]
  if(class(k) %in% c("numeric", "integer")){
    col_sd <- sd(k)
    result <- rbind(result, data.frame(D = col_name, sd = col_sd))
  } else {
    # 非数值列可按需标记，这里示例填NA
    result <- rbind(result, data.frame(D = col_name, sd = NA))
  }
}

# 输出结果，取消行名即可匹配目标格式
print(result, row.names = FALSE)

方案2：向量化计算（更简洁高效） #

R中优先用向量化函数替代显式循环，代码更短，列数较多时运行效率更高：

# 批量遍历所有列计算标准差
sd_calc <- sapply(dat, function(col) {
  ifelse(class(col) %in% c("numeric", "integer"), sd(col), NA)
})

# 组装为目标格式数据框
result <- data.frame(
  D = names(sd_calc),
  sd = as.numeric(sd_calc),
  row.names = NULL
)

两种方案运行后输出结果和目标格式完全一致，其中d3、d4、d7列所有取值相同，标准差计算结果为0，和示例预期匹配。

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