use moment instead of sum abs2

Co-authored-by: David Widmann <[email protected]>
JuliaStats · Oct 23, 2023 · 135e43a · 135e43a
1 parent 475498d
commit 135e43a
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/src/shapiro_wilk.jl b/src/shapiro_wilk.jl
@@ -72,8 +72,7 @@ end
 
 function unsafe_swstat(X::AbstractVector{<:Real}, A::AbstractVector{<:Real})
     AX = @inbounds dot(view(A, 1:length(X)), X)
-    m = mean(X)
-    S² = sum(x -> abs2(x - m), X)
+    S² = moment(X, 2)
     W = AX^2 / S²
     return min(W, one(W)) # to guard against numeric errors
 end