fix per JuliaStats#242 (comment)

tbenst · Sep 10, 2021 · 9690adb · 9690adb
1 parent dd35ae0
commit 9690adb
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/var_equality.jl b/src/var_equality.jl
@@ -60,7 +60,7 @@ end
 function anova(scores::AbstractVector{<:Real}...)
     Nᵢ = [length(g) for g in scores]
     Z̄ᵢ = mean.(scores)
-    Z̄ = mean(Z̄ᵢ)
+    Z̄ = sum(Iterators.flatten(groups))/sum(Nᵢ)
     SStᵢ = Nᵢ .* (Z̄ᵢ .- Z̄).^2
     SSeᵢ = sum.( (z .- z̄).^2 for (z, z̄) in zip(scores, Z̄ᵢ) )
     (Nᵢ, SStᵢ, SSeᵢ)