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

本書では200人分の測定値をまとめているが、pp.10-11で階級に属するのは145～148であれば不等式で145＜x≦148ということと約束しており、測定値に143が1つあり、その約束どおりにすると143は143～146の階級には含まれない。本計算では、その約束に従い計算している。

> dtf <- read.csv("table1_1_1.csv", header = TRUE)
> height <- dtf$height
> w <- (170 - 143) / 9
> sumn <- length(which(height > 143 & height <= 170))
> cumn <- cumrf <- 0.0
> cat("女子大生199人の身長の度数分布表\n")
女子大生199人の身長の度数分布表
> for (i in seq(143, 170 - w, w)) {
+ n <- length(which(height > i & height <= 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))
+ }
143～146  2  1.0%   2   1.0%
146～149  5  2.5%   7   3.5%
149～152 17  8.5%  24  12.1%
152～155 41 20.6%  65  32.7%
155～158 47 23.6% 112  56.3%
158～161 37 18.6% 149  74.9%
161～164 30 15.1% 179  89.9%
164～167 12  6.0% 191  96.0%
167～170  8  4.0% 199 100.0%
> cat(sprintf("計 %d\n", sumn))
計 199
> hist(dtf$height, breaks = seq(143, 173, by = 3))