［R］共分散（「まずはこの一冊から 意味が分かる統計解析」（ペレ出版）、p.214）

以下のp.152の表の数値をファイルtable_p214.csvに保存して、カレントディレクトリに置いておく。

 i,     h,    w
 1, 147.9, 41.7
 2, 163.5, 60.2
 3, 159.8, 47.0
 4, 155.1, 53.2
 5, 163.3, 48.3
 6, 158.7, 55.2
 7, 172.0, 58.5
 8, 161.2, 49.0
 9, 153.9, 46.7
10, 161.6, 52.5

計算する。

> dtf <- read.csv("table_p214.csv")
> n <- nrow(dtf)
> h <- dtf$h
> w <- dtf$w
> mh <- mean(h)
> mw <- mean(w)
> sxy <- sum((h - mh) * (w - mw)) / n
> spxy <- sum((h - mh) * (w - mw)) / (n - 1)
> cat(sprintf("身長の平均 = %f\n", mh))
身長の平均 = 159.700000
> cat(sprintf("体重の平均 = %f\n", mw))
体重の平均 = 51.230000
> cat(sprintf("sxy = %f\n", sxy))
sxy = 23.730000
> cat(sprintf("cov関数の結果 = %f\n", cov(h, w)))
cov関数の結果 = 26.366667
> cat(sprintf("s'xy = %f\n", spxy))
s'xy = 26.366667