Rで標本分散、標準偏差、CVを計算する
Rでvarやsdで計算されるのは標本分散、標本標準偏差ではなくそれぞれ不偏分散と不偏標準偏差になります。 をサンプル数、を平均値とします。
標本分散
標本標準偏差
不偏分散
不偏標準偏差
標本分散と標準偏差、ついでにCV(coefficient of variation, 標準偏差/平均)を求める関数は下のようになります。
# 標本分散 sample_var <- function(x,...){ return(var(x,...) * (length(x)-1) / length(x)) } # 標本標準偏差 sample_sd <-function(x,...){ return(sd(x,...) * sqrt((length(x)-1)/length(x))) } # 変動係数 cv <- function(x,...){ return(sample_sd(x,...)/mean(x)) }