Error in model_parameters.aov plot() method?

DESCRIPTION

@@ -1,7 +1,7 @@
 Type: Package
 Package: see
 Title: Model Visualisation Toolbox for 'easystats' and 'ggplot2'
-Version: 0.9.0.12
+Version: 0.9.0.13
 Authors@R:
     c(person(given = "Daniel",
              family = "Lüdecke",

NEWS.md

@@ -11,9 +11,12 @@
 - `plot()` for `model_parameters()` gets a `show_direction` argument, to turn
   off the direction of the effect in the plot.
 
+- `plot()` for `model_parameters()` now gives an informative error message when
+  trying to plot Anova tables (which does not work).
+
 - `plot()` for `simulate_parameters()` now better copes with models that have
   multiple response levels (e.g. multinomial models).
-  
+
 - Gains `{patchwork}` as a hard dependency given its importance for the package.
 
 ## Bug fixes

R/plot.parameters_model.R

@@ -47,7 +47,7 @@
 #' result
 #' plot(result)
 #' @export
 plot.see_parameters_model <- function(x,
                                       show_intercept = FALSE,
                                       size_point = 0.8,
                                       size_text = NA,
@@ -65,6 +65,12 @@
   # retrieve settings ----------------
   model_attributes <- attributes(x)[!names(attributes(x)) %in% c("names", "row.names", "class")]
 
+  # no plot for anova
+  model_class <- attributes(x)$model_class
+  if (!is.null(model_class) && any(model_class %in% c("aov", "anova"))) {
+    insight::format_error("Plots cannot be created for Anova tables. Please fit a (linear) regression model and try again.") # nolint
+  }
+
   # for GAMs, remove smooth terms
   if (!is.null(x$Component) && any(x$Component == "smooth_terms")) {
     x <- x[x$Component != "smooth_terms", ]
@@ -91,7 +97,7 @@
 
   # Clean up column names
   if (!any(grepl("Coefficient", colnames(x), fixed = TRUE))) {
     colnames(x)[which.min(match(colnames(x), c("Median", "Mean", "Map", "MAP", model_attributes$coefficient_name)))] <- "Coefficient"
   }
 
   if (!any(grepl("Parameter", colnames(x), fixed = TRUE))) {
@@ -161,7 +167,7 @@
     if (!is.null(cleaned_parameters)) {
       x <- merge(x, cleaned_parameters, sort = FALSE)
       x$Parameter <- x$Cleaned_Parameter
       if (all(x$Group == "")) {
         x$Group <- NULL
       } else {
         x <- x[!startsWith(x$Group, "SD/Cor"), , drop = FALSE]
@@ -180,7 +186,7 @@
   x <- .fix_facet_names(x)
 
   if (is_bayesian && "Group" %in% colnames(x)) {
     x$Effects[x$Group != ""] <- paste0(x$Effects[x$Group != ""], " (", x$Group[x$Group != ""], ")")
   }
 
   # remember components
@@ -196,7 +202,7 @@
   is_meta_bma <- !is.null(mc) && any(mc %in% c("meta_random", "meta_fixed", "meta_bma"))
 
   # minor fixes for Bayesian models
   if (!is.null(mc) && !is.null(cp) && any(mc %in% c("stanreg", "stanmvreg", "brmsfit")) && length(cp) == length(x$Parameter)) {
     x$Parameter <- cp
   }
 
@@ -208,9 +214,9 @@
 
     if (is_bootstrap || isTRUE(is_bayesian)) {
       if (is_bootstrap) {
         data <- model_attributes$boot_samples
       } else {
         data <- as.data.frame(.retrieve_data(x))
       }
 
       # MCMC or bootstrapped models
@@ -218,7 +224,7 @@
         insight::format_error("Could not retrieve parameter simulations.")
       }
 
       data <- datawizard::reshape_longer(
         data,
         names_to = "Parameter",
         rows_to = "Iteration",

tests/testthat/test-plot.parameters_model.R

@@ -14,6 +14,16 @@ test_that("`plot.see_parameters_model()` works", {
   )
 })
 
+
+test_that("`plot.see_parameters_model()` errors for anova tables", {
+  skip_if_not_installed("parameters")
+  data(iris)
+  m <- aov(Sepal.Length ~ Species, data = iris)
+  mp <- parameters::model_parameters(m)
+  expect_error(plot(mp), regex = "Plots cannot be created")
+})
+
+
 test_that("`plot.see_parameters_model()` random parameters works", {
   skip_if_not_installed("lme4")
   skip_if_not_installed("parameters")