correctly pass unbalanced param to HSD.test #62

computational-metabolomics · Jul 26, 2021 · 4c1e8a8 · 4c1e8a8
1 parent 448e85d
commit 4c1e8a8
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/HSD_class.R b/R/HSD_class.R
@@ -138,7 +138,7 @@ setMethod(f="model_apply",
 
             # for each combination of factors...
             out2=lapply(FF,function(x) {
-                A=agricolae::HSD.test(LM,x,group = FALSE)$comparison
+                A=agricolae::HSD.test(LM,x,group = FALSE,unbalanced=M$unbalanced)$comparison
                 if (ALIAS) {
                     A[!is.na(A)]=NA # replace with NA if alias are present
                 }