在使用R语言的过程中,会存在许多细节问题,再此mark一下。
1.var函数
var是R计算方差(variance)的函数,在概率论中,方差var(X)=E{[X-E(X)]^2},用R语言表示的计算方式为:
varpro <- function(x)
{
varp <- sum((x - mean(x))^2) / length(x)
print(varp)
}
在数理统计中,样本方差var(X)的计算方式为:
varstat <- function(x)
{
vars <- sum((x - mean(x))^2) / (length(x) - 1)
print(vars)
}