« [R]女子大生200人の身長の度数分布表とそのヒストグラム(「統計解析のはなし」(東京図書)、pp.17-19) | トップページ | [R]女子大生200人のバストの度数分布表とそのヒストグラム(「統計解析のはなし」(東京図書)、p.21) »

2022年11月11日 (金)

[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))

Fig1_1_5

« [R]女子大生200人の身長の度数分布表とそのヒストグラム(「統計解析のはなし」(東京図書)、pp.17-19) | トップページ | [R]女子大生200人のバストの度数分布表とそのヒストグラム(「統計解析のはなし」(東京図書)、p.21) »

R(本の計算を再現)」カテゴリの記事

コメント

コメントを書く

コメントは記事投稿者が公開するまで表示されません。

(ウェブ上には掲載しません)

« [R]女子大生200人の身長の度数分布表とそのヒストグラム(「統計解析のはなし」(東京図書)、pp.17-19) | トップページ | [R]女子大生200人のバストの度数分布表とそのヒストグラム(「統計解析のはなし」(東京図書)、p.21) »