From 6699dad277e07816d15ed4e3dff33e3228293e2e Mon Sep 17 00:00:00 2001 From: Jannes Reddig Date: Thu, 4 Jul 2019 10:03:42 +0200 Subject: [PATCH] Update Fisher Skewness Previously incorrect formula. x should be replaced by (x - mean(x)) --- R/skewness.R | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/R/skewness.R b/R/skewness.R index cd9915b5..fa28ab62 100644 --- a/R/skewness.R +++ b/R/skewness.R @@ -94,7 +94,7 @@ skewness <- if (n < 3) skewness = NA else - skewness = ((sqrt(n*(n-1))/(n-2))*(sum(x^3)/n))/((sum(x^2)/n)^(3/2)) + skewness = ((sqrt(n*(n-1))/(n-2))*(sum((x-mean(x))^3)/n)) / ((sum((x-mean(x))^2)/n)^(3/2)) } if (method == "sample") { skewness = sum((x-mean(x))^3/sqrt(var(x)*(n-1)/n)^3)*n/((n-1)*(n-2))