Bug fix

yufongpeng · Feb 28, 2024 · 93b3ee5 · 93b3ee5
1 parent dcc2dc5
commit 93b3ee5
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/anova.jl b/src/anova.jl
@@ -86,7 +86,7 @@ function anova(::Type{FTest}, aovm::FullModel{<: TRM_LM})
     elseif aovm.type == 2
         fstat = ntuple(last(fullasgn) - offset) do fix
             s1 = sort!(collect(select_super_interaction(fullpred, fix + offset)))
-            s2 = setdiff(select1, fix + offset)
+            s2 = setdiff(s1, fix + offset)
             select1 = findall(in(s1), fullasgn)
             select2 = findall(in(s2), fullasgn)
             (β[select1]' * (varβ[select1, select1] \ β[select1]) - β[select2]' * (varβ[select2, select2] \ β[select2])) / df[fix]