R语言中summary()函数为何无法正确显示数值？

问题场景 #

我有一个含10000个数值的向量，分布偏斜且存在极值。调用base R的summary()函数时，部分统计量显示为0，但实际存储的数值并非0。

实际输出 #

> summary(vec)
     Min.   1st Qu.    Median      Mean   3rd Qu.      Max. 
0.000e+00 0.000e+00 0.000e+00 1.244e+16 0.000e+00 1.225e+20

> str(summary(vec))
'summaryDefault' Named num [1:6] 3.69e-207 1.73e-01 2.84e-01 1.24e+16 4.52e-01 ...
 - attr(*, "names")= chr [1:6] "Min." "1st Qu." "Median" "Mean" ...

期望输出 #

Min.   1st Qu.    Median      Mean   3rd Qu.      Max. 
3.69e-207 1.73e-01   2.84e-01 1.244e+16  4.52e-01 1.225e+20

可复现示例（仅最小值显示为0） #

> x = seq(from=3.69e-207, to=1.23e+20, length.out=10000)

> summary(x)
     Min.   1st Qu.    Median      Mean   3rd Qu.      Max. 
0.000e+00 3.075e+19 6.150e+19 6.150e+19 9.225e+19 1.230e+20 

> str(summary(x))
 'summaryDefault' Named num [1:6] 3.69e-207 3.08e+19 6.15e+19 6.15e+19 9.22e+19 ...
 - attr(*, "names")= chr [1:6] "Min." "1st Qu." "Median" "Mean" ...

> summary(x)["Min."]
     Min. 
3.69e-207 

问题原因 #

summary()默认打印时会使用统一的科学计数法精度，当部分统计量的数值远小于其他极值时，会被格式化为0，但实际数值仍正确存储在结果对象中。

解决方法 #

1. 自定义打印函数 #

手动提取summary结果，指定灵活的格式输出：

print_summary <- function(x) {
  s <- summary(x)
  # 打印表头
  cat(sprintf("%-8s %-10s %-10s %-10s %-10s %-10s\n", names(s)))
  # 打印数值，按需求调整精度
  cat(sprintf("%-8.3e %-10.3e %-10.3e %-10.4e %-10.3e %-10.3e\n", s))
}

# 调用示例
print_summary(vec)

2. 直接提取目标统计量 #

如果只需查看特定统计量，直接索引提取即可：

summary(vec)[c("Min.", "1st Qu.", "Median", "3rd Qu.")]

3. 用format()自定义输出 #

通过format()设置足够的精度，避免小数值被截断：

format(summary(vec), scientific = TRUE, digits = 3)

内容的提问来源于Stack Exchange，提问作者datadon