diff --git a/NAMESPACE b/NAMESPACE
index 32dee53..2c0ff2f 100644
--- a/NAMESPACE
+++ b/NAMESPACE
@@ -78,6 +78,7 @@ export(position_nudge_keep)
 export(position_nudge_line)
 export(position_nudge_to)
 export(position_stack_keep)
+export(position_stack_minmax)
 export(position_stacknudge)
 export(scale_npcx_continuous)
 export(scale_npcy_continuous)
diff --git a/NEWS.md b/NEWS.md
index e4d522d..3353a46 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -8,8 +8,9 @@ editor_options:
 
 # ggpp 0.5.5
 
--   Add data set `birch.df` for birch seedlings.
--   Move transcriptomics example data from 'ggpmisc' to 'ggpp'.
+-   Add data sets `birch.df`and `birch_dw.df` with data for 350 birch seedlings.
+-   Move transcriptomics data sets `quadrant_example.df` and `volcano_example.df` from 'ggpmisc' to 'ggpp'.
+-   Add `position_stack_minmax()` to correctly position error bars in stacked plots.
 
 # ggpp 0.5.4
 
diff --git a/R/position-nudge-stack.R b/R/position-nudge-stack.R
index c9e08b1..a695d1b 100644
--- a/R/position-nudge-stack.R
+++ b/R/position-nudge-stack.R
@@ -5,16 +5,37 @@
 #' \code{\link[ggplot2]{position_nudge}}, \code{position_stacknudge()} returns
 #' in \code{data} both the original coordinates and the nudged coordinates.
 #'
-#' This position function is backwards compatible with
+#' \code{position_fillnudge()} is useful when labelling plots such as filled
+#' bars, filled columns, filled lines, etc. In contrast to
+#' \code{\link[ggplot2]{position_nudge}}, \code{position_fillnudge()} returns
+#' in \code{data} both the original coordinates and the nudged coordinates.
+#'
+#' The wrapper \code{position_nudge_keep()} has the same signature and
+#' behaviour as \code{\link[ggplot2]{position_nudge}} nad provides an easier to
+#' remember name when the need is only to have access to both the original and
+#' nudged coordinates.
+#'
+#' These position functions are backwards compatible with
 #' \code{\link[ggplot2]{position_nudge}} but extends it by adding support for
 #' stacking and for geometries that make use of the original position to draw
 #' connecting segments or arrows.
 #'
-#' The wrapper \code{position_nudge_keep()} with exactly the same signature and
-#' behaviour as \code{\link[ggplot2]{position_nudge}} provides an easier to
-#' remember name when the desire is only to have access to both the original and
+#' The wrapper \code{position_stack_keep()} has the same signature and
+#' behaviour as \code{\link[ggplot2]{position_stack}} and provides an easier to
+#' remember name when the need is only to have access to both the original and
+#' nudged coordinates.
+#'
+#' The wrapper \code{position_fill_keep()} has the same signature and
+#' behaviour as \code{\link[ggplot2]{position_fill}} and provides an easier to
+#' remember name when the need is only to have access to both the original and
 #' nudged coordinates.
 #'
+#' The wrapper \code{position_stack_minmax()} has the same signature and
+#' behaviour as \code{\link[ggplot2]{position_stack}} but stacks y, ymin and
+#' ymax in parallel, making it possible to stack summaries with error bars,
+#' works correctly with \code{geom_pointrange()}, \code{geom_linerange()} and
+#' \code{geom_errorbar()}.
+#'
 #' @family position adjustments
 #'
 #' @param vjust Vertical adjustment for geoms that have a position (like points
@@ -101,6 +122,17 @@
 #'     vjust = "bottom") +
 #'   theme(legend.position = "none")
 #'
+#' ggplot(birch_dw.df,
+#'        aes(y = dry.weight * 1e-3, x = Density, fill = Part)) +
+#'    stat_summary(geom = "col", fun = mean, na.rm = TRUE,
+#'                 position = "stack", alpha = 0.7, width = 0.67) +
+#'    stat_summary(geom = "linerange", fun.data = mean_cl_normal, na.rm = TRUE,
+#'                 position = position_stack_minmax()) +
+#'   labs(y = "Seedling dry mass (g)") +
+#'   scale_fill_grey(start = 0.7, end = 0.3) +
+#'   facet_wrap(facets = vars(Container))
+#'
+
 position_stacknudge <-
   function(vjust = 1,
            reverse = FALSE,
@@ -190,9 +222,27 @@ PositionStackAndNudge <-
                      )
                    },
 
+                   setup_data = function(self, data, params) {
+                     data <- flip_data(data, params$flipped_aes)
+                     if (is.null(params$var)) {
+                       return(data)
+                     }
+
+                     data$ymax <- switch(params$var,
+                                         y = data$y,
+                                         ymax = as.numeric(ifelse(data$ymax == 0, data$ymin, data$ymax))
+                     )
+
+                     data <- remove_missing(
+                       data,
+                       vars = c("x", "xmin", "xmax", "y"),
+                       name = "position_stack"
+                     )
+                     flip_data(data, params$flipped_aes)
+                   },
+
                    compute_layer = function(self, data, params, layout) {
                      x_orig <- data$x
-                     y_orig <- data$y
 
                      # operate on the stacked positions (updated in August 2020)
                      data = ggplot2::ggproto_parent(ggplot2::PositionStack, self)$compute_layer(data, params, layout)
@@ -216,7 +266,6 @@ PositionStackAndNudge <-
                                                            function(y) y + params$nudge_y * params$.fun_y(y))
                      }
 
-                     # add original position
                      if (params$kept.origin == "stacked") {
                        data$x_orig <- x_stacked
                        data$y_orig <- y_stacked
@@ -271,3 +320,93 @@ position_fill_keep <-
                        direction = "none",
                        kept.origin = "original")
   }
+
+#' @rdname ggpp-ggproto
+#' @format NULL
+#' @usage NULL
+#' @noRd
+PositionStackMinMax <-
+  ggplot2::ggproto("PositionStackMinMax", ggplot2::PositionStack,
+                   var = "y",
+
+                   setup_params = function(self, data) {
+                     c(
+                       list(kept.origin = self$kept.origin, var = self$var, fill = FALSE),
+                       ggplot2::ggproto_parent(ggplot2::PositionStack, self)$setup_params(data)
+                     )
+                   },
+
+                   setup_data = function(self, data, params) {
+                     data <- flip_data(data, params$flipped_aes)
+                     if (!(exists("ymax", data) && exists("ymin", data) && exists("y", data))) {
+                       stop("position_stack_minmax() requires y, ymin and ymax mappings in data")
+                     }
+                     flip_data(data, params$flipped_aes)
+                   },
+
+                   compute_layer = function(self, data, params, layout) {
+                     x_orig <- data$x
+                     if (exists("xmin", data)) {
+                       xmin_delta <-  data$xmin - data$x
+                     } else {
+                       xmin_delta <-  0
+                     }
+                     if (exists("xmax", data)) {
+                       xmax_delta <-  data$xmax - data$x
+                     } else {
+                       xmax_delta <-  0
+                     }
+
+                     y_orig <- data$y
+                     if (exists("ymin", data)) {
+                       ymin_delta <-  data$ymin - data$y
+                     } else {
+                       ymin_delta <-  0
+                     }
+                     if (exists("ymax", data)) {
+                       ymax_delta <-  data$ymax - data$y
+                     } else {
+                       ymax_delta <-  0
+                     }
+
+                     # operate on the stacked positions (updated in August 2020)
+                     data = ggplot2::ggproto_parent(ggplot2::PositionStack, self)$compute_layer(data, params, layout)
+                     x_stacked <- data$x
+                     y_stacked <- data$y
+                     if (exists("xmin", data)) data$xmin <- data$x + xmin_delta
+                     if (exists("xmax", data)) data$xmax <- data$x + xmax_delta
+                     if (exists("ymin", data)) data$ymin <- data$y + ymin_delta
+                     if (exists("ymax", data)) data$ymax <- data$y + ymax_delta
+
+                     # add original position
+                     if (params$kept.origin == "original") {
+                       data$x_orig <- x_orig
+                       data$y_orig <- y_orig
+                     }
+
+                     data
+                   },
+
+                   compute_panel = function(self, data, params, scales) {
+                     ggplot2::ggproto_parent(PositionStack, self)$compute_panel(data, params, scales)
+                   }
+  )
+
+#' @rdname position_stacknudge
+#'
+#' @export
+#'
+position_stack_minmax <-
+  function(vjust = 1,
+           reverse = FALSE,
+           kept.origin = c("none", "original")) {
+
+    kept.origin <- rlang::arg_match(kept.origin)
+
+    ggplot2::ggproto(NULL, PositionStackMinMax,
+                     kept.origin = kept.origin,
+                     vjust = vjust,
+                     reverse = reverse
+    )
+  }
+
diff --git a/README.Rmd b/README.Rmd
index 415493e..e730400 100644
--- a/README.Rmd
+++ b/README.Rmd
@@ -116,6 +116,7 @@ keep in the `data` object the original coordinates.
 | `position_nudge_keep()`   | nudge          | x, y (fixed distance)           | data labels                 |
 | `position_jitter_keep()`  | jitter         | x, y (random)                   | dot plots                   |
 | `position_stack_keep()`   | stack          | vertical (absolute)             | column and bar plots        |
+| `position_stack_minmax()` | stack          | vertical (absolute)             | error bars                  |
 | `position_fill_keep()`    | fill           | vertical (relative, fractional) | column plots                |
 | `position_dodge_keep()`   | dodge          | sideways (absolute)             | column and bar plots        |
 | `position_dosge2_keep()`  | dodge2         | sideways (absolute)             | box plots                   |
diff --git a/README.md b/README.md
index 880f0b5..baeffe8 100644
--- a/README.md
+++ b/README.md
@@ -105,6 +105,7 @@ keep in the `data` object the original coordinates.
 | `position_nudge_keep()`   | nudge          | x, y (fixed distance)           | data labels                 |
 | `position_jitter_keep()`  | jitter         | x, y (random)                   | dot plots                   |
 | `position_stack_keep()`   | stack          | vertical (absolute)             | column and bar plots        |
+| `position_stack_minmax()` | stack          | vertical (absolute)             | error bars                  |
 | `position_fill_keep()`    | fill           | vertical (relative, fractional) | column plots                |
 | `position_dodge_keep()`   | dodge          | sideways (absolute)             | column and bar plots        |
 | `position_dosge2_keep()`  | dodge2         | sideways (absolute)             | box plots                   |
diff --git a/man/position_stacknudge.Rd b/man/position_stacknudge.Rd
index 3d0730c..f354174 100644
--- a/man/position_stacknudge.Rd
+++ b/man/position_stacknudge.Rd
@@ -5,6 +5,7 @@
 \alias{position_fillnudge}
 \alias{position_stack_keep}
 \alias{position_fill_keep}
+\alias{position_stack_minmax}
 \title{Combined positions stack and nudge}
 \source{
 \url{https://github.com/slowkow/ggrepel/issues/161}.
@@ -31,6 +32,12 @@ position_fillnudge(
 position_stack_keep(vjust = 1, reverse = FALSE)
 
 position_fill_keep(vjust = 1, reverse = FALSE)
+
+position_stack_minmax(
+  vjust = 1,
+  reverse = FALSE,
+  kept.origin = c("original", "none")
+)
 }
 \arguments{
 \item{vjust}{Vertical adjustment for geoms that have a position (like points
@@ -63,15 +70,36 @@ bars, stacked columns, stacked lines, etc. In contrast to
 in \code{data} both the original coordinates and the nudged coordinates.
 }
 \details{
-This position function is backwards compatible with
+\code{position_fillnudge()} is useful when labelling plots such as filled
+bars, filled columns, filled lines, etc. In contrast to
+\code{\link[ggplot2]{position_nudge}}, \code{position_fillnudge()} returns
+in \code{data} both the original coordinates and the nudged coordinates.
+
+The wrapper \code{position_nudge_keep()} has the same signature and
+behaviour as \code{\link[ggplot2]{position_nudge}} nad provides an easier to
+remember name when the need is only to have access to both the original and
+nudged coordinates.
+
+These position functions are backwards compatible with
 \code{\link[ggplot2]{position_nudge}} but extends it by adding support for
 stacking and for geometries that make use of the original position to draw
 connecting segments or arrows.
 
-The wrapper \code{position_nudge_keep()} with exactly the same signature and
-behaviour as \code{\link[ggplot2]{position_nudge}} provides an easier to
-remember name when the desire is only to have access to both the original and
+The wrapper \code{position_stack_keep()} has the same signature and
+behaviour as \code{\link[ggplot2]{position_stack}} and provides an easier to
+remember name when the need is only to have access to both the original and
+nudged coordinates.
+
+The wrapper \code{position_fill_keep()} has the same signature and
+behaviour as \code{\link[ggplot2]{position_fill}} and provides an easier to
+remember name when the need is only to have access to both the original and
 nudged coordinates.
+
+The wrapper \code{position_stack_minmax()} has the same signature and
+behaviour as \code{\link[ggplot2]{position_stack}} but stacks y, ymin and
+ymax in parallel, making it possible to stack summaries with error bars,
+works correctly with \code{geom_pointrange()}, \code{geom_linerange()} and
+\code{geom_errorbar()}.
 }
 \examples{
 
@@ -129,6 +157,16 @@ ggplot(data = df, aes(x2, x1, group = grp)) +
     vjust = "bottom") +
   theme(legend.position = "none")
 
+ggplot(birch_dw.df,
+       aes(y = dry.weight * 1e-3, x = Density, fill = Part)) +
+   stat_summary(geom = "col", fun = mean, na.rm = TRUE,
+                position = "stack", alpha = 0.7, width = 0.67) +
+   stat_summary(geom = "linerange", fun.data = mean_cl_normal, na.rm = TRUE,
+                position = position_stack_minmax()) +
+  labs(y = "Seedling dry mass (g)") +
+  scale_fill_grey(start = 0.7, end = 0.3) +
+  facet_wrap(facets = vars(Container))
+
 }
 \seealso{
 \code{\link[ggplot2]{position_nudge}},
diff --git a/tests/testthat/Rplots.pdf b/tests/testthat/Rplots.pdf
index 10f8b44..e96ecdb 100644
Binary files a/tests/testthat/Rplots.pdf and b/tests/testthat/Rplots.pdf differ
diff --git a/vignettes/grammar-extensions.Rmd b/vignettes/grammar-extensions.Rmd
index 3b7776a..ee62a26 100644
--- a/vignettes/grammar-extensions.Rmd
+++ b/vignettes/grammar-extensions.Rmd
@@ -316,6 +316,7 @@ of the original coordinates kept in data.
 | `position_dodge_keep()`                     | `ggplot2::position_dodge()`       |
 | `position_dodge2_keep()`                    | `ggplot2::position_dodge2()`      |
 | `position_stack_keep()`                     | `ggplot2::position_stack()`       |
+| `position_stack_minmax()`                   |                                   |
 | `position_fill_keep()`                      | `ggplot2::position_fill()`        |
 | `position_jitterdodge_keep()` *(planned)*   | `ggplot2::position_jitterdodge()` |
 | `position_jitternudge()`                    |                                   |