歪度と同じく面倒
# 尖度
kurt <- function(x, method = c("SPSS", "ordinary"))
{
method <- match.arg(method)
n <- length(x)
if (method == "SPSS") {
n*(n+1)*sum(scale(x)^4)/(n-1)/(n-2)/(n-3)-3*(n-1)^2/(n-2)/(n-3)
}
else {
sum(((x-mean(x))/sqrt((n-1)*var(x)/n))^4)/n-3
}
}
# 通常の（?）定義（これ以外にもある）
> kurt(c(1,3,2,5,7,9), method="ordinary")
[1] -1.289086
# SPSS や Excel による結果
# C. I. Bliss, Statistics in Biology, Vol. 1, 1967, McGraw-Hill. による
> kurt(c(1,3,2,5,7,9), method="SPSS")
[1] -1.259834
# こちらをデフォールトにしてあります
> kurt(c(1,3,2,5,7,9))
[1] -1.259834