diff --git a/src/var_equality.jl b/src/var_equality.jl
index c80a8d69..7268fc9b 100644
--- 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(scores))/sum(Nᵢ)
     SStᵢ = Nᵢ .* (Z̄ᵢ .- Z̄).^2
     SSeᵢ = sum.( (z .- z̄).^2 for (z, z̄) in zip(scores, Z̄ᵢ) )
     (Nᵢ, SStᵢ, SSeᵢ)
diff --git a/test/var_equality.jl b/test/var_equality.jl
index e58383ce..98af052d 100644
--- a/test/var_equality.jl
+++ b/test/var_equality.jl
@@ -31,7 +31,8 @@ using DelimitedFiles
     t = OneWayANOVATest(groups...)
     @test nobs(t) == [7, 8, 8, 6]
     @test dof(t) == (3,25)
-    @test pvalue(t) ≈ 0.072 atol=1e-3
+    @test pvalue(t) ≈ 0.07276 atol=1e-6
+    @test HypothesisTests.teststatistic(t) ≈ 2.62311 atol=1e-6
     @test occursin("reject h_0", sprint(show, t))
 
     # http://www.real-statistics.com/one-way-analysis-of-variance-anova/homogeneity-variances/levenes-test/