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

> dtf <- read.csv("table1_1_1.csv", header = TRUE)
> weight <- dtf$weight
> w <- (70.0 - 38) / 9
> sumn <- length(which(weight > 38 & weight <= 70.0))
> cumn <- cumrf <- 0.0
> cat("女子大生198人の体重の度数分布表\n")
女子大生198人の体重の度数分布表
> for (i in seq(38, 70 - w, w)) {
+ n <- length(which(weight > i & weight <= i + w))
+ rf <- n / sumn
+ cumn <- cumn + n
+ cumrf <- cumrf + rf
+ cat(sprintf("%.2f～%.2f %2d %4.1f%% %3d %5.1f%%\n", i, i + w, n, rf * 100, cumn, cumrf * 100))
+ }
38.00～41.56  5  2.5%   5   2.5%
41.56～45.11 29 14.6%  34  17.2%
45.11～48.67 47 23.7%  81  40.9%
48.67～52.22 53 26.8% 134  67.7%
52.22～55.78 30 15.2% 164  82.8%
55.78～59.33 20 10.1% 184  92.9%
59.33～62.89  7  3.5% 191  96.5%
62.89～66.44  4  2.0% 195  98.5%
66.44～70.00  3  1.5% 198 100.0%
> cat(sprintf("計 %d\n", sumn))
計 198
> hist(dtf$weight, breaks = seq(38, 70.0, length = 10))