diff --git a/R/HSD_class.R b/R/HSD_class.R
index 6bc0410..c2fd0c0 100644
--- 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
                 }