［R］女子大生200人のバストの度数分布表とそのヒストグラム（「統計解析のはなし」（東京図書）、p.21）

> dtf <- read.csv("table1_1_1.csv", header = TRUE)
> bust <- dtf$bust
> w <- 3
> sumn <- length(which(bust > 71 & bust <= 98))
> cumn <- cumrf <- 0.0
> cat("女子大生199人のバストの度数分布表\n")
女子大生199人のバストの度数分布表
> for (i in seq(71, 98 - w, w)) {
+ n <- length(which(bust > i & bust <= i + w))
+ rf <- n / sumn
+ cumn <- cumn + n
+ cumrf <- cumrf + rf
+ cat(sprintf("%d～%d %2d %4.1f%% %3d %5.1f%%\n", i, i + w, n, rf * 100, cumn, cumrf * 100))
+ }
71～74  0  0.0%   0   0.0%
74～77  8  4.0%   8   4.0%
77～80 88 44.2%  96  48.2%
80～83 53 26.6% 149  74.9%
83～86 35 17.6% 184  92.5%
86～89  6  3.0% 190  95.5%
89～92  4  2.0% 194  97.5%
92～95  3  1.5% 197  99.0%
95～98  2  1.0% 199 100.0%
> cat(sprintf("計 %d\n", sumn))
計 199
> hist(dtf$bust, breaks = seq(71, 98, length = 10))