## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(collapse = FALSE, comment = "", fig.width = 7, fig.height = 4.6,
                      dev = if (knitr::is_latex_output()) "cairo_pdf" else
                        if (requireNamespace("svglite", quietly = TRUE)) "svglite" else "png")
if (requireNamespace("stability", quietly = TRUE)) library(stability) else
  pkgload::load_all(export_all = FALSE, helpers = FALSE, quiet = TRUE)


## ----indiv--------------------------------------------------------------------
ia <-
  indiv_anova(
      .data = ge_data
    , .y    = Yield
    , .rep  = Rep
    , .gen  = Gen
    , .env  = Env
    )
ia$Summary
ia$Homogeneity
plot(x = ia, which = "error")


## ----indiv-int----------------------------------------------------------------
stab_interpret(
    x             = ia
  , higher_better = TRUE
  , top           = 5
  )


## ----add----------------------------------------------------------------------
aa <-
  add_anova(
      .data = ge_data
    , .y    = Yield
    , .rep  = Rep
    , .gen  = Gen
    , .env  = Env
    )
aa$anova
aa$Summary


## ----er-----------------------------------------------------------------------
er <-
  er_anova(
      .data = ge_data
    , .y    = Yield
    , .rep  = Rep
    , .gen  = Gen
    , .env  = Env
    )
er$er_anova[c(1:6, nrow(er$er_anova)), ]
head(x = er$Parameters[order(-Mean)], n = 8)


## ----er-int-------------------------------------------------------------------
stab_interpret(
    x             = er
  , higher_better = TRUE
  , top           = 5
  )