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


## ----random-------------------------------------------------------------------
r1 <-
  Griffing(
      y          = Yield
    , Rep        = Rep
    , Cross1     = Cross1
    , Cross2     = Cross2
    , data       = GriffingData1
    , Method     = 1
    , Model      = 2
    , Reciprocal = c("combined", "split")[1]
    )
knitr::kable(x = r1$ANOVA, digits = 3)
knitr::kable(x = r1$Genetic.Components, digits = 3)


## ----all----------------------------------------------------------------------
comp <-
  sapply(
      X   = 1:4
    , FUN = function(m) {
        f <-
          Griffing(
              y          = Yield
            , Rep        = Rep
            , Cross1     = Cross1
            , Cross2     = Cross2
            , data       = get(paste0("GriffingData", m))
            , Method     = m
            , Model      = 2
            , Reciprocal = c("combined", "split")[1]
            )
        f$Genetic.Components[c("gca", "sca", "Error"), "Components"]
      }
    )
dimnames(comp) <- list(c("gca", "sca", "Error"), paste("Method", 1:4))
round(x = comp, digits = 2)


## ----ad-----------------------------------------------------------------------
r1$Genetic.Components[c("Additive", "Dominant", "gca/sca"), , drop = FALSE]