## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(collapse = TRUE, comment = "#>", fig.width = 6.8,
                      fig.height = 4.4, dpi = 120)

## ----setup--------------------------------------------------------------------
library(trialSizing)

## ----basic--------------------------------------------------------------------
reps <- calc_replicates(
  treatments  = c(3, 10, 50),
  cv_percent  = 9.25,
  lsd_percent = c(10, 20, 30),
  design      = "CRD"
)
reps

## ----summary------------------------------------------------------------------
summary(reps)

## ----access-------------------------------------------------------------------
head(reps$data[, c("Treatments", "LSD_percent", "r_continuous",
                   "r_optimal", "df_error", "q_tukey")])

## ----design-------------------------------------------------------------------
rbind(
  CRD  = calc_replicates(c(3, 50), 9.25, 10, design = "CRD")$data$r_continuous,
  RCBD = calc_replicates(c(3, 50), 9.25, 10, design = "RCBD")$data$r_continuous
)

## ----alpha--------------------------------------------------------------------
rbind(
  `alpha = 0.05` = calc_replicates(10, 9.25, 20, alpha = 0.05)$data$r_continuous,
  `alpha = 0.01` = calc_replicates(10, 9.25, 20, alpha = 0.01)$data$r_continuous
)

## ----workflow, message = FALSE------------------------------------------------
grid1 <- as.matrix(uniformity_trial[uniformity_trial$trial == "T1",
                                    grep("^col", names(uniformity_trial))])
cv_tab <- calc_cv_shapes(grid1)

lrp  <- fit_lrp(cv_tab, x = "x", cv = "cv", step = 0.01)
cvxo <- unname(lrp$parameters["Breakpoint_Response"])
cvxo

calc_replicates(treatments = c(5, 10, 20), cv_percent = cvxo,
                lsd_percent = c(10, 20), design = "RCBD")$data[
                  , c("Treatments", "LSD_percent", "r_continuous", "r_optimal")]

## ----plot---------------------------------------------------------------------
reps <- calc_replicates(treatments = 3:50, cv_percent = 9.25,
                        lsd_percent = c(10, 20, 30), design = "CRD")

plot(reps, title = "Replications needed (black oat)")

## ----plot-continuous----------------------------------------------------------
plot(reps, y_var = "r_continuous", title = "Replications (continuous)")

## ----save, eval = FALSE-------------------------------------------------------
# plot(reps, title = "Replications needed",
#      save = TRUE, file = "replications.tiff", format = "tiff", dpi = 300)