理学療法士がまとめたノート

統計学備忘録 since2016

Rを使って統計学を勉強するブログです

ヒストグラムのビン調整

head(iris)

  Sepal.Length Sepal.Width Petal.Length Petal.Width Species
1          5.1         3.5          1.4         0.2  setosa
2          4.9         3.0          1.4         0.2  setosa
3          4.7         3.2          1.3         0.2  setosa
4          4.6         3.1          1.5         0.2  setosa
5          5.0         3.6          1.4         0.2  setosa
6          5.4         3.9          1.7         0.4  setosa

変数Speciesがsetosaのデータを抽出
列Speciesがsetosatである行を全て抜き出す

iris_setosa <- iris[iris$Species=="setosa",]

パッケージMASSをインストール

install.packages("MASS")
library("MASS")

ヒストグラムのビン調整(0.1, 0.2, 0.3, 0.4)

par(mfrow = c(2, 2))
truehist(iris_setosa$Sepal.Length, h=0.1, prob=F) 
truehist(iris_setosa$Sepal.Length, h=0.2, prob=F) 
truehist(iris_setosa$Sepal.Length, h=0.3, prob=F) 
truehist(iris_setosa$Sepal.Length, h=0.4, prob=F) 
par(mfrow = c(1, 1))   

f:id:yoshida931:20180716223044p:plain:w600