From 8d891ce61725c2ff5799b9a53af39d1e0131b590 Mon Sep 17 00:00:00 2001
From: Michael Mumford <111483220+michaelhm-daf@users.noreply.github.com>
Date: Thu, 12 Feb 2026 11:46:48 +1000
Subject: [PATCH 1/5] Include line plot as an option in autoplot.mct

---
 .gitignore   | 3 +++
 R/autoplot.R | 7 +++++--
 2 files changed, 8 insertions(+), 2 deletions(-)

diff --git a/.gitignore b/.gitignore
index 6332a611..2d8302fd 100644
--- a/.gitignore
+++ b/.gitignore
@@ -11,3 +11,6 @@ biometryassist.Rproj
 *.Rdata
 *.Rds
 inst/doc
+.httr-oauth
+.DS_Store
+.quarto
diff --git a/R/autoplot.R b/R/autoplot.R
index 2b46139e..133cc58c 100644
--- a/R/autoplot.R
+++ b/R/autoplot.R
@@ -85,11 +85,14 @@ autoplot.mct <- function(object, size = 4, label_height = 0.1,
         ggplot2::theme(axis.text.x = ggplot2::element_text(angle = axis_rotation, vjust = 0.5, hjust = hjust_value, ...)) +
         ggplot2::labs(x = "", y = paste0("Predicted ", ylab))
 
-    if(tolower(type) == "point") {
+    if(tolower(type) %in% c("point","line")) {
         plot <- plot + ggplot2::geom_point(ggplot2::aes(y = {{ yval }}), colour = "black", shape = 16, size = 2) +
             ggplot2::geom_errorbar(aes(ymin = .data[["low"]], ymax = .data[["up"]]), width = 0.2)
+        if(tolower(type) == "line"){
+          plot <- plot + ggplot2::geom_line(ggplot2::aes(y = {{ yval }}), colour="black", linewidth=1)
+        }
     }
-    else if(tolower(type) %in% c("col", "column")) {
+    else if(tolower(type) %in% c(,"bar", "col", "column")) {
         plot <- plot + ggplot2::geom_col(ggplot2::aes(y = {{ yval }}), colour = "black", fill = "cornflowerblue", alpha = 0.75) +
             ggplot2::geom_errorbar(aes(ymin = .data[["low"]], ymax = .data[["up"]]), width = 0.2)
     }

From aea92e2205242fc3a9e581cac3dab8f4d768fd42 Mon Sep 17 00:00:00 2001
From: Michael Mumford <111483220+michaelhm-daf@users.noreply.github.com>
Date: Thu, 12 Feb 2026 12:57:19 +1000
Subject: [PATCH 2/5] Add options to exclude errorbars or group lettering from
 autoplot.mct

---
 R/autoplot.R    | 20 +++++++++++++-------
 man/autoplot.Rd |  6 ++++++
 2 files changed, 19 insertions(+), 7 deletions(-)

diff --git a/R/autoplot.R b/R/autoplot.R
index 133cc58c..24aba001 100644
--- a/R/autoplot.R
+++ b/R/autoplot.R
@@ -7,6 +7,8 @@
 #' @param axis_rotation Enables rotation of the x axis independently of the group labels within the plot.
 #' @param label_rotation Enables rotation of the treatment group labels independently of the x axis labels within the plot.
 #' @param type A string specifying the type of plot to display. The default of 'point' will display a point estimate with error bars. The alternative, 'column' (or 'col'), will display a column graph with error bars.
+#' @param include_errorbar A logical indicating whether to include errorbars when plotting the predicted values from a multiple comparisons test
+#' @param include_lettering A logical indicating whether to include group lettering when plotting the predicted values from a multiple comparisons test
 #' @param margin Logical (default `FALSE`). A value of `FALSE` will expand the plot to the edges of the plotting area i.e. remove white space between plot and axes.
 #' @param palette A string specifying the colour scheme to use for plotting or a vector of custom colours to use as the palette. Default is equivalent to "Spectral". Colour blind friendly palettes can also be provided via options `"colour blind"` (or `"colour blind"`, both equivalent to `"viridis"`), `"magma"`, `"inferno"`, `"plasma"`, `"cividis"`, `"rocket"`, `"mako"` or `"turbo"`. Other palettes from [scales::brewer_pal()] are also possible.
 #' @param row A variable to plot a column from `object` as rows.
@@ -40,7 +42,8 @@ ggplot2::autoplot
 #' autoplot(output, label_height = 0.5)
 autoplot.mct <- function(object, size = 4, label_height = 0.1,
                          rotation = 0, axis_rotation = rotation,
-                         label_rotation = rotation, type = "point", ...) {
+                         label_rotation = rotation, type = "point",
+                         include_errorbar=TRUE, include_lettering=TRUE, ...) {
     stopifnot(inherits(object, "mct"))
 
     rlang::check_dots_used()
@@ -86,18 +89,21 @@ autoplot.mct <- function(object, size = 4, label_height = 0.1,
         ggplot2::labs(x = "", y = paste0("Predicted ", ylab))
 
     if(tolower(type) %in% c("point","line")) {
-        plot <- plot + ggplot2::geom_point(ggplot2::aes(y = {{ yval }}), colour = "black", shape = 16, size = 2) +
-            ggplot2::geom_errorbar(aes(ymin = .data[["low"]], ymax = .data[["up"]]), width = 0.2)
+        plot <- plot + ggplot2::geom_point(ggplot2::aes(y = {{ yval }}), colour = "black", shape = 16, size = 2) #+
+            #ggplot2::geom_errorbar(aes(ymin = .data[["low"]], ymax = .data[["up"]]), width = 0.2)
         if(tolower(type) == "line"){
-          plot <- plot + ggplot2::geom_line(ggplot2::aes(y = {{ yval }}), colour="black", linewidth=1)
+          plot <- plot + ggplot2::geom_line(ggplot2::aes(y = {{ yval }}, group=1), colour="black", linewidth=0.4)
         }
     }
     else if(tolower(type) %in% c(,"bar", "col", "column")) {
-        plot <- plot + ggplot2::geom_col(ggplot2::aes(y = {{ yval }}), colour = "black", fill = "cornflowerblue", alpha = 0.75) +
-            ggplot2::geom_errorbar(aes(ymin = .data[["low"]], ymax = .data[["up"]]), width = 0.2)
+        plot <- plot + ggplot2::geom_col(ggplot2::aes(y = {{ yval }}), colour = "black", fill = "cornflowerblue", alpha = 0.75)
+    }
+    
+    if( ("low" %in% colnames(pred_df)) && (include_errorbar==TRUE) ){
+      plot <- plot + ggplot2::geom_errorbar(aes(ymin = .data[["low"]], ymax = .data[["up"]]), width = 0.2)
     }
 
-    if("groups" %in% colnames(pred_df)) {
+    if( ("groups" %in% colnames(pred_df)) && (include_lettering==TRUE) ) {
         # Calculate outside of aes()
         y_pos <- ifelse(pred_df$up > pred_df$low, pred_df$up, pred_df$low)
         nudge_val <- ifelse(abs(label_height) <= 1,
diff --git a/man/autoplot.Rd b/man/autoplot.Rd
index 8026d982..086218fc 100644
--- a/man/autoplot.Rd
+++ b/man/autoplot.Rd
@@ -16,6 +16,8 @@ autoplot(object, ...)
   axis_rotation = rotation,
   label_rotation = rotation,
   type = "point",
+  include_errorbar = TRUE,
+  include_lettering = TRUE,
   ...
 )
 
@@ -50,6 +52,10 @@ autoplot(object, ...)
 
 \item{type}{A string specifying the type of plot to display. The default of 'point' will display a point estimate with error bars. The alternative, 'column' (or 'col'), will display a column graph with error bars.}
 
+\item{include_errorbar}{A logical indicating whether to include errorbars when plotting the predicted values from a multiple comparisons test}
+
+\item{include_lettering}{A logical indicating whether to include group lettering when plotting the predicted values from a multiple comparisons test}
+
 \item{margin}{Logical (default \code{FALSE}). A value of \code{FALSE} will expand the plot to the edges of the plotting area i.e. remove white space between plot and axes.}
 
 \item{palette}{A string specifying the colour scheme to use for plotting or a vector of custom colours to use as the palette. Default is equivalent to "Spectral". Colour blind friendly palettes can also be provided via options \code{"colour blind"} (or \code{"colour blind"}, both equivalent to \code{"viridis"}), \code{"magma"}, \code{"inferno"}, \code{"plasma"}, \code{"cividis"}, \code{"rocket"}, \code{"mako"} or \code{"turbo"}. Other palettes from \code{\link[scales:pal_brewer]{scales::brewer_pal()}} are also possible.}

From 068faeb8cb5f27ff1b700c08597ad521b4b2ad62 Mon Sep 17 00:00:00 2001
From: Michael Mumford <111483220+michaelhm-daf@users.noreply.github.com>
Date: Thu, 12 Feb 2026 15:01:26 +1000
Subject: [PATCH 3/5] Add option to display average HSD errorbar instead of
 confidence interval in autoplot.mct

---
 R/autoplot.R    | 21 ++++++++++++++++-----
 man/autoplot.Rd |  3 +++
 2 files changed, 19 insertions(+), 5 deletions(-)

diff --git a/R/autoplot.R b/R/autoplot.R
index 24aba001..d733dddd 100644
--- a/R/autoplot.R
+++ b/R/autoplot.R
@@ -9,6 +9,7 @@
 #' @param type A string specifying the type of plot to display. The default of 'point' will display a point estimate with error bars. The alternative, 'column' (or 'col'), will display a column graph with error bars.
 #' @param include_errorbar A logical indicating whether to include errorbars when plotting the predicted values from a multiple comparisons test
 #' @param include_lettering A logical indicating whether to include group lettering when plotting the predicted values from a multiple comparisons test
+#' @param errorbar_type A character (default is "ci") that indicates what the errorbars in the plot represent. Current options are 95% confidence interval ("ci") or Tukeys (average) HSD value ("hsd")
 #' @param margin Logical (default `FALSE`). A value of `FALSE` will expand the plot to the edges of the plotting area i.e. remove white space between plot and axes.
 #' @param palette A string specifying the colour scheme to use for plotting or a vector of custom colours to use as the palette. Default is equivalent to "Spectral". Colour blind friendly palettes can also be provided via options `"colour blind"` (or `"colour blind"`, both equivalent to `"viridis"`), `"magma"`, `"inferno"`, `"plasma"`, `"cividis"`, `"rocket"`, `"mako"` or `"turbo"`. Other palettes from [scales::brewer_pal()] are also possible.
 #' @param row A variable to plot a column from `object` as rows.
@@ -43,6 +44,7 @@ ggplot2::autoplot
 autoplot.mct <- function(object, size = 4, label_height = 0.1,
                          rotation = 0, axis_rotation = rotation,
                          label_rotation = rotation, type = "point",
+                         errorbar_type="ci",
                          include_errorbar=TRUE, include_lettering=TRUE, ...) {
     stopifnot(inherits(object, "mct"))
 
@@ -89,20 +91,29 @@ autoplot.mct <- function(object, size = 4, label_height = 0.1,
         ggplot2::labs(x = "", y = paste0("Predicted ", ylab))
 
     if(tolower(type) %in% c("point","line")) {
-        plot <- plot + ggplot2::geom_point(ggplot2::aes(y = {{ yval }}), colour = "black", shape = 16, size = 2) #+
-            #ggplot2::geom_errorbar(aes(ymin = .data[["low"]], ymax = .data[["up"]]), width = 0.2)
+        plot <- plot + ggplot2::geom_point(ggplot2::aes(y = {{ yval }}), colour = "black", shape = 16, size = 2)
         if(tolower(type) == "line"){
           plot <- plot + ggplot2::geom_line(ggplot2::aes(y = {{ yval }}, group=1), colour="black", linewidth=0.4)
         }
     }
-    else if(tolower(type) %in% c(,"bar", "col", "column")) {
+    else if(tolower(type) %in% c("bar", "col", "column")) {
         plot <- plot + ggplot2::geom_col(ggplot2::aes(y = {{ yval }}), colour = "black", fill = "cornflowerblue", alpha = 0.75)
     }
     
-    if( ("low" %in% colnames(pred_df)) && (include_errorbar==TRUE) ){
+    if( (tolower(errorbar_type)=="ci") && (include_errorbar==TRUE) ){
       plot <- plot + ggplot2::geom_errorbar(aes(ymin = .data[["low"]], ymax = .data[["up"]]), width = 0.2)
     }
-
+    else if( (tolower(errorbar_type)=="hsd") && (include_errorbar==TRUE) ){
+      subset_df <- pred_df[1,]
+      plot <- plot + ggplot2::geom_errorbar(data=subset_df, 
+                        aes(x = {{ classify }},
+                            ymin = .data[["predicted.value"]] - 0.5 * object$hsd , 
+                            ymax = .data[["predicted.value"]] + 0.5 * object$hsd ), 
+                        width = 0.2#, 
+                        #position = ggplot2::position_dodge(width = 0.5) # does not work for some reason
+                        )
+    }
+    
     if( ("groups" %in% colnames(pred_df)) && (include_lettering==TRUE) ) {
         # Calculate outside of aes()
         y_pos <- ifelse(pred_df$up > pred_df$low, pred_df$up, pred_df$low)
diff --git a/man/autoplot.Rd b/man/autoplot.Rd
index 086218fc..8975f9a3 100644
--- a/man/autoplot.Rd
+++ b/man/autoplot.Rd
@@ -16,6 +16,7 @@ autoplot(object, ...)
   axis_rotation = rotation,
   label_rotation = rotation,
   type = "point",
+  errorbar_type = "ci",
   include_errorbar = TRUE,
   include_lettering = TRUE,
   ...
@@ -52,6 +53,8 @@ autoplot(object, ...)
 
 \item{type}{A string specifying the type of plot to display. The default of 'point' will display a point estimate with error bars. The alternative, 'column' (or 'col'), will display a column graph with error bars.}
 
+\item{errorbar_type}{A character (default is "ci") that indicates what the errorbars in the plot represent. Current options are 95\% confidence interval ("ci") or Tukeys (average) HSD value ("hsd")}
+
 \item{include_errorbar}{A logical indicating whether to include errorbars when plotting the predicted values from a multiple comparisons test}
 
 \item{include_lettering}{A logical indicating whether to include group lettering when plotting the predicted values from a multiple comparisons test}

From d71ef6fc3c9a81a660a6cb236566f23ed73c239d Mon Sep 17 00:00:00 2001
From: Michael Mumford <111483220+michaelhm-daf@users.noreply.github.com>
Date: Thu, 12 Feb 2026 16:25:31 +1000
Subject: [PATCH 4/5] Add option to display predicted values on the transformed
 scale in autoplot.mct

---
 R/autoplot.R    | 53 +++++++++++++++++++++++++++++++++++++++++++------
 man/autoplot.Rd |  7 +++++--
 2 files changed, 52 insertions(+), 8 deletions(-)

diff --git a/R/autoplot.R b/R/autoplot.R
index d733dddd..0a647c3a 100644
--- a/R/autoplot.R
+++ b/R/autoplot.R
@@ -7,9 +7,10 @@
 #' @param axis_rotation Enables rotation of the x axis independently of the group labels within the plot.
 #' @param label_rotation Enables rotation of the treatment group labels independently of the x axis labels within the plot.
 #' @param type A string specifying the type of plot to display. The default of 'point' will display a point estimate with error bars. The alternative, 'column' (or 'col'), will display a column graph with error bars.
-#' @param include_errorbar A logical indicating whether to include errorbars when plotting the predicted values from a multiple comparisons test
-#' @param include_lettering A logical indicating whether to include group lettering when plotting the predicted values from a multiple comparisons test
+#' @param include_errorbar Logical (default 'TRUE') indicating whether to include errorbars when plotting the predicted values from a multiple comparisons test
+#' @param include_lettering Logical (default 'TRUE') indicating whether to include group lettering when plotting the predicted values from a multiple comparisons test
 #' @param errorbar_type A character (default is "ci") that indicates what the errorbars in the plot represent. Current options are 95% confidence interval ("ci") or Tukeys (average) HSD value ("hsd")
+#' @param trans_scale Logical (default 'FALSE') that indicates whether the predicted values should be displayed on the transformed scale.
 #' @param margin Logical (default `FALSE`). A value of `FALSE` will expand the plot to the edges of the plotting area i.e. remove white space between plot and axes.
 #' @param palette A string specifying the colour scheme to use for plotting or a vector of custom colours to use as the palette. Default is equivalent to "Spectral". Colour blind friendly palettes can also be provided via options `"colour blind"` (or `"colour blind"`, both equivalent to `"viridis"`), `"magma"`, `"inferno"`, `"plasma"`, `"cividis"`, `"rocket"`, `"mako"` or `"turbo"`. Other palettes from [scales::brewer_pal()] are also possible.
 #' @param row A variable to plot a column from `object` as rows.
@@ -45,10 +46,19 @@ autoplot.mct <- function(object, size = 4, label_height = 0.1,
                          rotation = 0, axis_rotation = rotation,
                          label_rotation = rotation, type = "point",
                          errorbar_type="ci",
-                         include_errorbar=TRUE, include_lettering=TRUE, ...) {
+                         include_errorbar=TRUE, include_lettering=TRUE,
+                         trans_scale=FALSE, ...) {
     stopifnot(inherits(object, "mct"))
 
     rlang::check_dots_used()
+    
+    # Force trans_scale = TRUE if errorbar_type is "hsd"
+    if (tolower(errorbar_type) == "hsd" && include_errorbar == TRUE) {
+      if (!trans_scale) {
+        trans_scale <- TRUE  # Force trans_scale to TRUE
+        warning("The error bar is an average HSD value and not a confidence interval.")
+      }
+    }
 
     # Extract the predictions data frame from the mct object
     # For new structure: object is a list with $predictions
@@ -73,7 +83,7 @@ autoplot.mct <- function(object, size = 4, label_height = 0.1,
     # Get ylab as attribute (works for both old and new structure)
     ylab <- attributes(object)$ylab
 
-    yval <- ifelse("PredictedValue" %in% colnames(pred_df), "PredictedValue", "predicted.value")
+    yval <- ifelse( ("PredictedValue" %in% colnames(pred_df)) && (trans_scale==FALSE) , "PredictedValue", "predicted.value")
     yval <- rlang::ensym(yval)
 
     # Calculate hjust based on axis rotation
@@ -101,7 +111,13 @@ autoplot.mct <- function(object, size = 4, label_height = 0.1,
     }
     
     if( (tolower(errorbar_type)=="ci") && (include_errorbar==TRUE) ){
-      plot <- plot + ggplot2::geom_errorbar(aes(ymin = .data[["low"]], ymax = .data[["up"]]), width = 0.2)
+      if(trans_scale==TRUE){
+        plot <- plot + ggplot2::geom_errorbar(aes(ymin = .data[["predicted.value"]] -2*.data[["std.error"]], 
+                                                  ymax = .data[["predicted.value"]] +2*.data[["std.error"]]), 
+                                              width = 0.2)
+      } else {
+        plot <- plot + ggplot2::geom_errorbar(aes(ymin = .data[["low"]], ymax = .data[["up"]]), width = 0.2)
+      }
     }
     else if( (tolower(errorbar_type)=="hsd") && (include_errorbar==TRUE) ){
       subset_df <- pred_df[1,]
@@ -133,7 +149,32 @@ autoplot.mct <- function(object, size = 4, label_height = 0.1,
     else if(exists("classify2")) {
         plot <- plot + ggplot2::facet_wrap(stats::as.formula(paste("~", classify2)))
     }
-    return(plot)
+    
+    # Add secondary y-axis using the PredictedValue column
+    if (trans_scale) {
+      # Ensure the PredictedValue column exists in the data
+      if (!"PredictedValue" %in% colnames(pred_df)) {
+        stop("The column 'PredictedValue' does not exist in the data.")
+      }
+      
+      # Calculate the range of the secondary axis
+      primary_range <- range(pred_df[["predicted.value"]], na.rm = TRUE)
+      secondary_range <- range(pred_df[["PredictedValue"]], na.rm = TRUE)
+      
+      # Define a linear transformation between the primary and secondary axes
+      scale_factor <- diff(secondary_range) / diff(primary_range)
+      offset <- secondary_range[1] - primary_range[1] * scale_factor
+      
+      # Add the secondary y-axis
+      plot <- plot + ggplot2::scale_y_continuous(
+        sec.axis = ggplot2::sec_axis(
+          trans = ~ . * scale_factor + offset,  # Linear transformation
+          name = "Back-transformed Scale"  # Label for the secondary y-axis
+        )
+      )
+    }
+    
+  return(plot)
 }
 
 
diff --git a/man/autoplot.Rd b/man/autoplot.Rd
index 8975f9a3..708a0cb8 100644
--- a/man/autoplot.Rd
+++ b/man/autoplot.Rd
@@ -19,6 +19,7 @@ autoplot(object, ...)
   errorbar_type = "ci",
   include_errorbar = TRUE,
   include_lettering = TRUE,
+  trans_scale = FALSE,
   ...
 )
 
@@ -55,9 +56,11 @@ autoplot(object, ...)
 
 \item{errorbar_type}{A character (default is "ci") that indicates what the errorbars in the plot represent. Current options are 95\% confidence interval ("ci") or Tukeys (average) HSD value ("hsd")}
 
-\item{include_errorbar}{A logical indicating whether to include errorbars when plotting the predicted values from a multiple comparisons test}
+\item{include_errorbar}{Logical (default 'TRUE') indicating whether to include errorbars when plotting the predicted values from a multiple comparisons test}
 
-\item{include_lettering}{A logical indicating whether to include group lettering when plotting the predicted values from a multiple comparisons test}
+\item{include_lettering}{Logical (default 'TRUE') indicating whether to include group lettering when plotting the predicted values from a multiple comparisons test}
+
+\item{trans_scale}{Logical (default 'FALSE') that indicates whether the predicted values should be displayed on the transformed scale.}
 
 \item{margin}{Logical (default \code{FALSE}). A value of \code{FALSE} will expand the plot to the edges of the plotting area i.e. remove white space between plot and axes.}
 

From 8086c149427ee56e3b90ed7f2975ea36d5252415 Mon Sep 17 00:00:00 2001
From: Michael Mumford <111483220+michaelhm-daf@users.noreply.github.com>
Date: Fri, 13 Feb 2026 13:36:46 +1000
Subject: [PATCH 5/5] Add function to format skeleton ANOVA for strip-plot
 designs.

---
 R/satab.R | 68 +++++++++++++++++++++++++++++++++++++++++++++++--------
 1 file changed, 58 insertions(+), 10 deletions(-)

diff --git a/R/satab.R b/R/satab.R
index ca07ce01..fd9abd9d 100644
--- a/R/satab.R
+++ b/R/satab.R
@@ -309,6 +309,9 @@ format_satab <- function(anova_structure, design_type) {
   if (design_type == "split") {
     return(format_satab_split(anova_structure))
   }
+  else if (design_type == "strip") {
+    return(format_satab_strip(anova_structure))
+  }
 
   # Standard formatting
   output <- c(
@@ -344,11 +347,13 @@ format_satab_split <- function(anova_structure) {
   names <- anova_structure$names
 
   # Determine width based on df magnitude
-  width1 <- ifelse(df[1] > 10, 44, 45)
-  width2 <- ifelse(df[2] > 10, 35, 36)
-  width3 <- ifelse(df[4] > 10, 35, 36)
-  width4 <- ifelse(df[5] > 10, 35, 36)
-  width5 <- ifelse(df[7] > 10, 44, 45)
+  width1 <- ifelse(df[1] > 9, 44, 45)
+  width2 <- ifelse(df[2] > 9, 35, 36)
+  width3 <- ifelse(df[3] > 9, 44, 45)
+  width4 <- ifelse(df[4] > 9, 35, 36)
+  width5 <- ifelse(df[5] > 9, 35, 36)
+  width6 <- ifelse(df[6] > 9, 44, 45)
+  width7 <- ifelse(df[7] > 9, 44, 45)
 
   output <- c(
     paste0(format("Source of Variation", width = 45), "df", "\n"),
@@ -357,20 +362,63 @@ format_satab_split <- function(anova_structure) {
     "--------------------------------------------------\n",
     "Whole plot stratum\n",
     paste0(format(" ", width = 9), format(sources[2], width = width2), df[2], "\n"),
-    paste0(format(sources[3], width = 45), df[3], "\n"),
+    paste0(format(sources[3], width = width3), df[3], "\n"),
     "==================================================\n",
     "Subplot stratum\n",
-    paste0(format(" ", width = 9), format(sources[4], width = width3), df[4], "\n"),
-    paste0(format(" ", width = 9), format(sources[5], width = width4), df[5], "\n"),
-    paste0(format(" ", width = 9), format(sources[6], width = 35), df[6], "\n"),
+    paste0(format(" ", width = 9), format(sources[4], width = width4), df[4], "\n"),
+    paste0(format(" ", width = 9), format(sources[5], width = width5), df[5], "\n"),
+    paste0(format(sources[6], width = width6), df[6], "\n"),
     "==================================================\n",
-    paste0(format("Total", width = width5), df[7], "\n")
+    paste0(format("Total", width = width7), df[7], "\n")
   )
 
   class(output) <- c("satab", class(output))
   return(output)
 }
 
+
+#' Format SATAB for Split Plot (special case)
+#' @noRd
+format_satab_strip <- function(anova_structure) {
+  sources <- anova_structure$sources
+  df <- anova_structure$df
+  names <- anova_structure$names
+  
+  # Determine width based on df magnitude
+  width1 <- ifelse(df[1] > 9, 44, 45)
+  width2 <- ifelse(df[2] > 9, 35, 36)
+  width3 <- ifelse(df[3] > 9, 44, 45)
+  width4 <- ifelse(df[4] > 9, 35, 36)
+  width5 <- ifelse(df[5] > 9, 44, 45)
+  width6 <- ifelse(df[6] > 9, 35, 36)
+  width7 <- ifelse(df[7] > 9, 44, 45)
+  width8 <- ifelse(df[8] > 9, 44, 45)
+  
+  output <- c(
+    paste0(format("Source of Variation", width = 45), "df", "\n"),
+    "==================================================\n",
+    paste0(format(sources[1], width = width1), df[1], "\n"),
+    "--------------------------------------------------\n",
+    "Row strip stratum\n",
+    paste0(format(" ", width = 9), format(sources[2], width = width2), df[2], "\n"),
+    paste0(format(sources[3], width = width3), df[3], "\n"),
+    "==================================================\n",
+    "Column strip stratum\n",
+    paste0(format(" ", width = 9), format(sources[4], width = width4), df[4], "\n"),
+    paste0(format(sources[5], width = width5), df[5], "\n"),
+    "==================================================\n",
+    "Observational unit stratum\n",
+    paste0(format(" ", width = 9), format(sources[6], width = width6), df[6], "\n"),
+    paste0(format(sources[7], width = width7), df[7], "\n"),
+    "==================================================\n",
+    paste0(format("Total", width = width8), df[8], "\n")
+  )
+  
+  class(output) <- c("satab", class(output))
+  return(output)
+}
+
+
 #' @noRd
 #' @method print satab
 #' @export