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


## ----sim----------------------------------------------------------------------
s <-
  DiallelSim(
      Parents = 6
    , Method  = 1
    , Rep     = 4
    , Env     = 1
    , Mean    = 100
    , gca     = 5
    , sca     = 3
    , rec     = 1
    , gcaEnv  = 0
    , scaEnv  = 0
    , EnvSD   = 0
    , Block   = 2
    , Error   = 4
    , seed    = 123
    )
fs <-
  Griffing(
      y          = y
    , Rep        = Rep
    , Cross1     = Cross1
    , Cross2     = Cross2
    , data       = s
    , Method     = 1
    , Model      = 1
    , Reciprocal = c("combined", "split")[1]
    )
round(x = cbind(true = attr(x = s, which = "truth")$gca, estimated = fs$GCA), digits = 2)


## ----power--------------------------------------------------------------------
pw <-
  DiallelPower(
      Parents = 6
    , Method  = 2
    , Rep     = 2:5
    , Model   = 1
    , gca     = 1
    , sca     = 1.5
    , rec     = 0
    , Error   = 4
    , nsim    = 200
    , alpha   = 0.05
    , seed    = 1
    )
pw
plot(x = pw)