diff --git a/NAMESPACE b/NAMESPACE
index 8cf346e8..12244cea 100644
--- a/NAMESPACE
+++ b/NAMESPACE
@@ -10,6 +10,7 @@ export(plot_miss)
 export(plot_pattern)
 export(plot_pred)
 export(plot_trace)
+export(plot_variance)
 export(stripplot)
 export(xyplot)
 importFrom(magrittr,"%>%")
diff --git a/R/plot_miss.R b/R/plot_miss.R
index 0ac1f898..53b935e0 100644
--- a/R/plot_miss.R
+++ b/R/plot_miss.R
@@ -4,6 +4,7 @@
 #' @param vrb String, vector, or unquoted expression with variable name(s), default is "all".
 #' @param border Logical indicating whether borders should be present between tiles.
 #' @param row.breaks Optional numeric input specifying the number of breaks to be visualized on the y axis.
+#' @param ordered Logical indicating whether rows should be ordered according to their pattern.
 #'
 #' @return An object of class [ggplot2::ggplot].
 #'
@@ -15,7 +16,8 @@ plot_miss <-
   function(data,
            vrb = "all",
            border = FALSE,
-           row.breaks = nrow(data)) {
+           row.breaks = nrow(data),
+           ordered = FALSE) {
     # input processing
     if (is.matrix(data) && ncol(data) > 1) {
       data <- as.data.frame(data)
@@ -35,9 +37,23 @@ plot_miss <-
         )
       )
     }
-    # Create missingness indicator matrix
-    na.mat <- purrr::map_df(data[,vrb], function(y) as.numeric(is.na(y)))
+    if(ordered){
+      # extract md.pattern matrix
+      mdpat <- mice::md.pattern(data, plot = FALSE) %>%
+        head(., -1)
+      # save frequency of patterns
+      freq.pat <- rownames(mdpat) %>%
+        as.numeric()
 
+      na.mat <- mdpat %>%
+        as.data.frame() %>%
+        dplyr::select(-ncol(.)) %>%
+        dplyr::mutate(nmis = freq.pat) %>%
+        tidyr::uncount(nmis)
+    } else{
+    # Create missingness indicator matrix
+    na.mat <- purrr::map_df(data[,vrb], function(y) as.numeric(!is.na(y)))
+    }
     # extract pattern info
     vrb <- colnames(na.mat)
     rws <- nrow(na.mat)
@@ -87,11 +103,18 @@ plot_miss <-
         fill = "",
         alpha = ""
       ) +
+      ggplot2::coord_cartesian(expand = FALSE) +
       theme_minimice()
     if(border){
       gg <- gg + ggplot2::geom_tile(color = "black")
     } else{
       gg <- gg + ggplot2::geom_tile()
     }
+    if(ordered){
+      gg <- gg +
+        ggplot2::theme(axis.text.y = ggplot2::element_blank(),
+                       axis.ticks.y = ggplot2::element_blank()
+      )
+    }
     return(gg)
   }
diff --git a/man/plot_miss.Rd b/man/plot_miss.Rd
index b4a2e122..246d1d4b 100644
--- a/man/plot_miss.Rd
+++ b/man/plot_miss.Rd
@@ -4,7 +4,13 @@
 \alias{plot_miss}
 \title{Plot missingness in a dataset}
 \usage{
-plot_miss(data, vrb = "all", border = FALSE, row.breaks = nrow(data))
+plot_miss(
+  data,
+  vrb = "all",
+  border = FALSE,
+  row.breaks = nrow(data),
+  ordered = FALSE
+)
 }
 \arguments{
 \item{data}{An incomplete dataset of class \code{data.frame} or \code{matrix}.}
@@ -14,6 +20,8 @@ plot_miss(data, vrb = "all", border = FALSE, row.breaks = nrow(data))
 \item{border}{Logical indicating whether borders should be present between tiles.}
 
 \item{row.breaks}{Optional numeric input specifying the number of breaks to be visualized on the y axis.}
+
+\item{ordered}{Logical indicating whether rows should be ordered according to their pattern.}
 }
 \value{
 An object of class \link[ggplot2:ggplot]{ggplot2::ggplot}.
diff --git a/tests/testthat/test-plot_miss.R.R b/tests/testthat/test-plot_miss.R.R
new file mode 100644
index 00000000..2fd60495
--- /dev/null
+++ b/tests/testthat/test-plot_miss.R.R
@@ -0,0 +1,26 @@
+# create test objects
+dat <- mice::nhanes
+
+# tests
+test_that("plot_miss produces plot", {
+  expect_s3_class(plot_miss(dat), "ggplot")
+  expect_s3_class(plot_miss(dat), "ggplot")
+  expect_s3_class(plot_miss(cbind(dat, "testvar" = NA)), "ggplot")
+})
+
+test_that("plot_miss works with different inputs", {
+  expect_s3_class(plot_miss(dat, c("age", "bmi")), "ggplot")
+  expect_s3_class(plot_miss(dat, c(age, bmi)), "ggplot")
+  expect_s3_class(plot_miss(data.frame(age = dat$age, testvar = NA)), "ggplot")
+  expect_s3_class(plot_miss(cbind(dat, "with space" = NA)), "ggplot")
+})
+
+
+test_that("plot_miss with incorrect argument(s)", {
+  expect_output(plot_miss(na.omit(dat)))
+  expect_error(plot_miss("test"))
+  expect_error(plot_miss(dat, vrb = "test"))
+  expect_error(plot_miss(dat, cluster = "test"))
+  expect_error(plot_miss(cbind(dat, .x = NA)))
+  expect_error(plot_miss(dat, npat = "test"))
+})