Rと入門統計学

栗原伸一「入門統計学ー検定から多変量解析・実験計画法まで」の章末問題をRで解いてみる。

第1章

問2

myDataFrame <- data.frame(fid = 1:20, sales=c(400,15,480,993,600,150,115,50,0,130,3000,500,200,55,2200,1,900,1000,450,400), area=c(60,30,365,190,136,15,37,100,170,70,783,560,50,35,595,200,300,356,155,250))

hist(myDataFrame$area)

解答にある対数のヒストグラム

hist(log(myDataFrame$area,10))

算術平均
> mean(myDataFrame$sales)
[1] 581.95

分散
> var(myDataFrame$sales)
[1] 596000.5

標準偏差
> sd(myDataFrame$sales)
[1] 772.0107

変動係数
sd(myDataFrame$sales)/mean(myDataFrame$sales)
[1] 1.326593

問3
算術平均

> mean(c(2,4,5,7))
[1] 4.5
> mean(c(2,4,5,70))
[1] 20.25

幾何平均

> exp(mean(log(abs(c(2,4,5,7)))))
[1] 4.090623