diff --git a/R/simStructure.R b/R/simStructure.R
index b7f7121..cc768a3 100644
--- a/R/simStructure.R
+++ b/R/simStructure.R
@@ -94,6 +94,11 @@ simStructure <- function(dataS, method=c("direct", "multinom", "distribution"),
     strata <- factor(rep(1, nrow(dataS@data)))
   }
 
+  # Check levels for multinom method
+  if ( method == "multinom" & length(levels(strata)) < 2 ) {
+    stop("Your strata region must have 2 or more levels!\n")
+  }
+
   ##### set up household structure
 
   # generate variables on household level (indicated by H)