diff --git a/DESCRIPTION b/DESCRIPTION
index 4d75b29..870160a 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: piar
 Title: Price Index Aggregation
-Version: 0.8.1.9003
+Version: 0.8.1.9004
 Authors@R: c(
     person("Steve", "Martin", role = c("aut", "cre", "cph"),
            email = "marberts@protonmail.com",
diff --git a/R/aggregate.piar_index.R b/R/aggregate.piar_index.R
index fde6db8..3e6f43a 100644
--- a/R/aggregate.piar_index.R
+++ b/R/aggregate.piar_index.R
@@ -38,15 +38,15 @@ aggregate_index <- function(x, pias, pias2, na.rm, contrib, r, include_ea, chain
     }
     res2 <- aggregate_(x, pias2, na.rm, has_contrib, -r, include_ea, chainable)
     if (has_contrib) {
-      res[[2L]] <- Map(
+      res$contrib <- Map(
         super_aggregate_contrib(),
-        res[[2L]],
-        res2[[2L]],
-        res[[1L]],
-        res2[[1L]]
+        res$contrib,
+        res2$contrib,
+        res$index,
+        res2$index
       )
     }
-    res[[1L]] <- Map(\(x, y) (x * y)^0.5, res[[1L]], res2[[1L]])
+    res$index <- Map(\(x, y) (x * y)^0.5, res$index, res2$index)
   }
   
   if (include_ea) {
@@ -55,7 +55,7 @@ aggregate_index <- function(x, pias, pias2, na.rm, contrib, r, include_ea, chain
     lev <- unlist(pias$levels[-length(pias$levels)], use.names = FALSE)
   }
   
-  piar_index(res[[1L]], res[[2L]], lev, x$time, chainable)
+  piar_index(res$index, res$contrib, lev, x$time, chainable)
 }
 
 aggregate_ <- function(x, pias, na.rm, has_contrib, r, include_ea, chainable) {
@@ -122,7 +122,7 @@ aggregate_ <- function(x, pias, na.rm, has_contrib, r, include_ea, chainable) {
     contrib[[t]] <- unlist(rev(con), recursive = FALSE, use.names = FALSE)
   }
   
-  list(index, contrib)
+  list(index = index, contrib = contrib)
 }
 
 #' Aggregate elemental price indexes
@@ -180,8 +180,10 @@ aggregate_ <- function(x, pias, na.rm, has_contrib, r, include_ea, chainable) {
 #' `contrib = FALSE`).
 #' 
 #' If two aggregation structures are given then the steps above are done for
-#' each aggregation structure and the resulting indexes are combined with a
-#' geometric mean to make a superlative quadratic mean of order `r+1` index.
+#' each aggregation structure, with the aggregation for `pias` done with a
+#' generalized mean of order `r` the aggregation for `pias2` done with a
+#' generalized mean of order `-r`. The resulting indexes are combined with a
+#' geometric mean to make a superlative quadratic mean of order `2*r` index.
 #' Percent-change contributions are combined using a generalized van IJzeren
 #' decomposition; see [`gpindex::nested_transmute()`] for details.
 #'
@@ -191,8 +193,8 @@ aggregate_ <- function(x, pias, na.rm, has_contrib, r, include_ea, chainable) {
 #' @param x A price index, usually made by [elemental_index()].
 #' @param pias A price index aggregation structure or something that can be
 #' coerced into one. This can be made with [aggregation_structure()].
-#' @param pias2 An optional secondary aggregation structure with current-period
-#' weights to make a superlative index.
+#' @param pias2 An optional secondary aggregation structure, usually with
+#' current-period weights, to make a superlative index. See details.
 #' @param na.rm Should missing values be removed? By default, missing values
 #' are not removed. Setting `na.rm = TRUE` is equivalent to overall mean
 #' imputation.
@@ -202,7 +204,7 @@ aggregate_ <- function(x, pias, na.rm, has_contrib, r, include_ea, chainable) {
 #' averaging indexes over subperiods), or -1 for a harmonic index (usually for
 #' a Paasche index). Other values are possible; see
 #' [gpindex::generalized_mean()] for details. If `pias2` is given then the
-#' index is aggregated with a quadratic mean of order `r+1`.
+#' index is aggregated with a quadratic mean of order `2*r`.
 #' @param contrib Aggregate percent-change contributions in `x` (if any)?
 #' @param include_ea Should indexes for the elemental aggregates be included
 #' along with the aggregated indexes? By default, all index values are returned.
diff --git a/man/aggregate.piar_index.Rd b/man/aggregate.piar_index.Rd
index c830654..e5ed704 100644
--- a/man/aggregate.piar_index.Rd
+++ b/man/aggregate.piar_index.Rd
@@ -36,8 +36,8 @@ coerced into one. This can be made with \code{\link[=aggregation_structure]{aggr
 
 \item{...}{Not currently used.}
 
-\item{pias2}{An optional secondary aggregation structure with current-period
-weights to make a superlative index.}
+\item{pias2}{An optional secondary aggregation structure, usually with
+current-period weights, to make a superlative index. See details.}
 
 \item{na.rm}{Should missing values be removed? By default, missing values
 are not removed. Setting \code{na.rm = TRUE} is equivalent to overall mean
@@ -51,7 +51,7 @@ arithmetic index (the default for aggregating elemental indexes and
 averaging indexes over subperiods), or -1 for a harmonic index (usually for
 a Paasche index). Other values are possible; see
 \code{\link[gpindex:generalized_mean]{gpindex::generalized_mean()}} for details. If \code{pias2} is given then the
-index is aggregated with a quadratic mean of order \code{r+1}.}
+index is aggregated with a quadratic mean of order \code{2*r}.}
 
 \item{include_ea}{Should indexes for the elemental aggregates be included
 along with the aggregated indexes? By default, all index values are returned.}
@@ -116,8 +116,10 @@ index values have been imputed (i.e., when \code{na.rm = TRUE} and
 \code{contrib = FALSE}).
 
 If two aggregation structures are given then the steps above are done for
-each aggregation structure and the resulting indexes are combined with a
-geometric mean to make a superlative quadratic mean of order \code{r+1} index.
+each aggregation structure, with the aggregation for \code{pias} done with a
+generalized mean of order \code{r} the aggregation for \code{pias2} done with a
+generalized mean of order \code{-r}. The resulting indexes are combined with a
+geometric mean to make a superlative quadratic mean of order \code{2*r} index.
 Percent-change contributions are combined using a generalized van IJzeren
 decomposition; see \code{\link[gpindex:transmute_weights]{gpindex::nested_transmute()}} for details.
 }