R Tips：不偏分散と標本分散

Rで不偏分散ではなく、標本分散を計算しようとしたら、不偏分散の関数しかないようです。
それほど難しくないので、下記のような関数を作れば、標本分散を求めることが出来ます。

基本情報として、サンプルサイズをnとして
nで割るのが、標本分散
n-1で割るのが不偏分散

ということは、不偏分散の値を var(x)とすると

標本分散 = var(x) * (n-1) / n
となります。

# 不偏分散を求める関数 var
# 標本分散を求める関数 varp
varp <- function(x) { var(x) * (length(x)-1) / length(x) }