## ----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)


## ----data---------------------------------------------------------------------
data(ge_data)
str(ge_data)
with(data = ge_data, expr = table(Env, Rep))[1:3, ]


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


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


## ----reg----------------------------------------------------------------------
reg <-
  stab_reg(
      .data = ge_data
    , .y    = Yield
    , .rep  = Rep
    , .gen  = Gen
    , .env  = Env
    )
head(x = reg$StabIndvReg[order(-Mean)], n = 8)
plot(x = reg, which = "slope")


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


## ----gge, fig.height = 6------------------------------------------------------
p <-
  gge_biplot(
      .data  = ge_data
    , .y     = Yield
    , .rep   = Rep
    , .gen   = Gen
    , .env   = Env
    , .which = "which-won-where"
    , .svp   = "symmetric"
    )
p


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