統計とR言語の勉強をしています。
「Rによるやさしい統計学」の写経で勉強。

実行環境

Windows10 Pro 64bit
R：3.4.3
RStudio：1.1.383
RStudioのConsoleを使って実行しています

基本計算

計算機のように数字と+の符号を打つだけで計算してくれます。
sum関数を使っても同じ結果です。

> 10+13+8+15+8
[1] 54
> sum(10,13,8,15,8)
[1] 54

変数に入れ込むことで計算が(計算機ではなく)プログラムらしくなります。

> testa <- c(10,13,8,15,8)
> testa
[1] 10 13  8 15  8
> sum(testa)
[1] 54

平均を求めるのであれば、以下の３つの方法があります(変数testaは上の続きです)。length関数はデータの個数を出力します。当然、mean関数を使うのが最もシンプルです。

> sum(testa)/5
[1] 10.8
> sum(testa)/length(testa)
[1] 10.8
> mean(testa)
[1] 10.8

median関数は中央値をtable関数は度数分布を出力します。

> median(testa)
[1] 10
> table(testa)
testa
 8 10 13 15 
 2  1  1  1

分散と標準偏差

地味に計算していきます。平均、偏差の二乗和、分散、標準偏差の順に計算しています。

> testaave <- mean(testa)
> testaave
[1] 10.8

# 平均との差
> testadev <- testa - testaave
> testaave
[1] 10.8
> testadev
[1] -0.8  2.2 -2.8  4.2 -2.8

# 偏差の2乗
> testadevsqr <- testadev^2
> testadevsqr
[1]  0.64  4.84  7.84 17.64  7.84

# 偏差の2乗和
> testadevsqrsum <- sum(testadevsqr)
> testadevsqrsum
[1] 38.8

# 分散(不偏分散でない)
> testavar <- testadevsqrsum/length(testa)
> testavar
[1] 7.76

# 標準偏差
> sqrt(testavar)
[1] 2.785678

varという関数があるのですが、不偏分散なので注意が必要です。上で算出した7.76と結果が異なります。

> var(testa)
[1] 9.7

標準偏差はsdです。こちらも不偏分散の平方根なので上で算出した2.785678と結果が異なります。

> sd(testa)
[1] 3.114482

その他

maxとminという関数もあります。最大値と最小値です。

> max(testa)
[1] 15
> min(testa)
[1] 8