我正在尝试手动验证sd()函数,但是我的数学又回来了。 StdDev是方差的平方根。方差是减去均值并平方后的所有数字的均值。
y <- c(2,3,4,5)
m <- mean(y) # 3.5
a <- c()
a[1] <- (y[1] - m)^2
a[2] <- (y[2] - m)^2
a[3] <- (y[3] - m)^2
a[4] <- (y[4] - m)^2
v <- mean(a) # variance here
sqrt(v) # outputs 1.118034
sd(y) # outputs 1.290994