diff --git a/R/modMAarea.R b/R/modMAarea.R
index cd6588f7..77b0495d 100644
--- a/R/modMAarea.R
+++ b/R/modMAarea.R
@@ -52,6 +52,8 @@
 #' only one domain, rowvar = domain variable. If more than one domain, include
 #' colvar. If no domain, rowvar = NULL.
 #' @param colvar String. Name of the column domain variable in cond or tree.
+#' @param sumunits Logical. If TRUE, estimation units are summed and returned
+#' in one table.
 #' @param bootstrap Logical. If TRUE, returns bootstrap variance estimates,
 #' otherwise uses Horvitz-Thompson estimator under simple random sampling
 #' without replacement.
@@ -250,6 +252,7 @@ modMAarea <- function(MApopdat,
                       pcfilter = NULL, 
                       rowvar = NULL, 
                       colvar = NULL, 
+                      sumunits = FALSE,
                       bootstrap = FALSE,
                       returntitle = FALSE, 
                       savedata = FALSE, 
@@ -283,7 +286,6 @@ modMAarea <- function(MApopdat,
   esttype="AREA"
   parameters <- FALSE
   returnlst <- list()
-  sumunits <- FALSE
   rawdata <- TRUE
   
   ## Set global variables
diff --git a/R/modMAratio.R b/R/modMAratio.R
index f575c52d..d1acbd14 100644
--- a/R/modMAratio.R
+++ b/R/modMAratio.R
@@ -230,7 +230,7 @@ modMAratio <- function(MApopdat,
                        FIA = TRUE,
                        rowvar = NULL, 
                        colvar = NULL, 
-                       sumunits = TRUE, 
+                       sumunits = FALSE, 
                        returntitle = FALSE, 
                        savedata = FALSE, 
                        table_opts = NULL, 
diff --git a/R/modMAtree.R b/R/modMAtree.R
index c4bdd752..aae81c3b 100644
--- a/R/modMAtree.R
+++ b/R/modMAtree.R
@@ -67,6 +67,8 @@
 #' @param colvar String. Optional. If rowvar != NULL, name of domain variable
 #' to group estvar by for columns in table output. Colvar must be included in
 #' an input data frame (i.e., plt, cond, tree).
+#' @param sumunits Logical. If TRUE, estimation units are summed and returned
+#' in one table.
 #' @param prednames String vector. Name(s) of predictor variables to include in
 #' model.
 #' @param modelselect Logical. If TRUE, an elastic net regression model is fit 
@@ -284,6 +286,7 @@ modMAtree <- function(MApopdat,
                       pcfilter = NULL, 
                       rowvar = NULL, 
                       colvar = NULL, 
+                      sumunits = FALSE,
                       prednames = NULL, 
                       modelselect = FALSE, 
                       FIA = TRUE,
@@ -319,7 +322,6 @@ modMAtree <- function(MApopdat,
   esttype="TREE"
   parameters <- FALSE
   returnlst <- list()
-  sumunits <- FALSE
   
   ## Set global variables
   ONEUNIT=n.total=n.strata=strwt=TOTAL <- NULL
diff --git a/man/modMAarea.Rd b/man/modMAarea.Rd
index 74b3c27b..82c1f820 100644
--- a/man/modMAarea.Rd
+++ b/man/modMAarea.Rd
@@ -14,6 +14,7 @@ modMAarea(
   pcfilter = NULL,
   rowvar = NULL,
   colvar = NULL,
+  sumunits = FALSE,
   bootstrap = FALSE,
   returntitle = FALSE,
   savedata = FALSE,
@@ -54,6 +55,9 @@ colvar. If no domain, rowvar = NULL.}
 
 \item{colvar}{String. Name of the column domain variable in cond or tree.}
 
+\item{sumunits}{Logical. If TRUE, estimation units are summed and returned
+in one table.}
+
 \item{bootstrap}{Logical. If TRUE, returns bootstrap variance estimates,
 otherwise uses Horvitz-Thompson estimator under simple random sampling
 without replacement.}
diff --git a/man/modMAratio.Rd b/man/modMAratio.Rd
index c8f156c6..91fe1f62 100644
--- a/man/modMAratio.Rd
+++ b/man/modMAratio.Rd
@@ -19,7 +19,7 @@ modMAratio(
   FIA = TRUE,
   rowvar = NULL,
   colvar = NULL,
-  sumunits = TRUE,
+  sumunits = FALSE,
   returntitle = FALSE,
   savedata = FALSE,
   table_opts = NULL,
diff --git a/man/modMAtree.Rd b/man/modMAtree.Rd
index 3843d474..0ed62c5c 100644
--- a/man/modMAtree.Rd
+++ b/man/modMAtree.Rd
@@ -15,6 +15,7 @@ modMAtree(
   pcfilter = NULL,
   rowvar = NULL,
   colvar = NULL,
+  sumunits = FALSE,
   prednames = NULL,
   modelselect = FALSE,
   FIA = TRUE,
@@ -65,6 +66,9 @@ for the total estimation unit. Include colvar for grouping by 2 variables.}
 to group estvar by for columns in table output. Colvar must be included in
 an input data frame (i.e., plt, cond, tree).}
 
+\item{sumunits}{Logical. If TRUE, estimation units are summed and returned
+in one table.}
+
 \item{prednames}{String vector. Name(s) of predictor variables to include in
 model.}