Pelo que sei kurtosis() do pacote e1071 que é em relação a 0, é calcular e somar 3. Mas do pacote momentes kurtosis() é em relação a 3 direto. Vide o código. Fiz os cálculos à mão devido ter este monte de estimadores. Os dois cálculos abaixo é usando definição de momentos.
============================== Código R ======================================================
require(e1071)
x=c(72, 70, 77, 60, 67, 69, 68, 66, 65, 71, 69)
kurtosis(x,type=1)
+3skewness(x,type=1)
=============================Código R========================================================
require(moments)
x=c(72, 70, 77, 60, 67, 69, 68, 66, 65, 71, 69)
skewness(x)
kurtosis(x)
agostino.test(x, alternative = "two.sided")
anscombe.test(x, alternative="two.sided")
#two sided = dois lados
#less = abaixo
#greater = acima
===========================================================================================
André Oliveira Souza.
Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES
Em Quarta-feira, 17 de Dezembro de 2014 22:55, Luiz Roberto Martins Pinto <luizroberto.uesc@gmail.com> escreveu:
Olá André,
Agradeço a sua sugestão de leitura.
De fato, a kurtose média na distribuição normal é 3. Ocorre que nesta função kurtosis(x), da biblioteca agricolae, a saída é um desvio desta curtose média, por isto kurtosis(x)=0.
Se eu estiver enganado, por favor me informe.
Luiz Roberto.