id:bob3にコメントをいただけました。ありがとうございます。
で、修正しました。
Rのデフォルトにある関数はこんな感じ。
統計量 | 関数名 |
平均 | mean |
分散 | var |
標準偏差 | sd |
中央値 | median |
最小値 | min |
最大値 | max |
範囲 | range |
MAD(median absolute deviation) | mad |
> x <- sample(10,1000,replace=TRUE) > mean(x) # 平均 [1] 5.618 > var(x) # 分散 [1] 7.689766 > sd(x) # 標準偏差 [1] 2.773043 > median(x) # 中央値 [1] 6 > min(x) # 最小値 [1] 1 > max(x) # 最大値 [1] 10 > range(x) # 最大値と最小値を返します。 [1] 1 10 > mad(x) [1] 2.9652
で、歪度と尖度はe1071ライブラリ、MAD(median absolute deviation) はgdataパッケージ、CI(confidence intervals) はgmodelsパッケージ にある。
統計量 | 関数名 |
歪度 | kurtosis |
尖度 | skewness |
MAD | mad |
CI | ci |
> library(e1071) 要求されたパッケージ class をロード中です > kurtosis(x) # 歪度 [1] -1.148149 > skewness(x) # 尖度 [1] -0.01719176 > #> library(gdata) # madは基本パッケージにありました。 #> mad(x) #[1] 2.9652 > > library(gmodels) > ci(x) Estimate CI lower CI upper Std. Error 5.61800000 5.44591971 5.79008029 0.08769131
で、最頻値はどのパッケージにあるかしらないので自作する。
こんな感じかな?
> as.numeric(names(table(x))[table(x)==max(table(x))]) [1] 4 6