repo.rcbd: Automatic report for a Randomized Complete Block Design...
In reyzaguirre/pepa: Package for the Execution of Pre Cooked Analysis

repo.rcbd

R Documentation

Automatic report for a Randomized Complete Block Design (RCBD)

Description

Produces an automatic report for selected traits in an experiment with a RCBD.

Usage

repo.rcbd(
  traits,
  trt,
  trt.lab = "treatment",
  rep,
  eu = NULL,
  dfr,
  maxp = 0.1,
  mc = FALSE,
  title = "Automatic report for a Randomized Complete Block Design (RCBD)",
  subtitle = NULL,
  author = "International Potato Center",
  format = c("html", "word", "pdf")
)

Arguments

`traits`	The traits to analize.
`trt`	The treatments.
`trt.lab`	The labels for treatments.
`rep`	The replications.
`eu`	The experimental unit. Must be defined in case of subsamples.
`dfr`	The name of the data frame.
`maxp`	Maximum allowed proportion of missing values to estimate, default is 10%.
`mc`	Logical. If `"mc = TRUE"` multiple comparison tests are included even if the factor effect is not significat. Default to `"mc = FALSE"`.
`title`	Report title.
`subtitle`	Report subtitle.
`author`	Report author.
`format`	The output file format for the report, `"html"` by default. Other options are `"word"` and `"pdf"`.

Details

It fits a linear model for a RCBD and explains the results.

Under the assumption of fixed effects an ANOVA table is computed with missing values estimated up to a specified percentage (10% by default). In case of subsampling the ANOVA table is computed over the means for each experimental unit, and a second ANOVA table is computed considering subsampling only if there are no missing values. If the ANOVA results in a significant value for treatments then the Tukey HSD method for pairwise differences is applied. Assumptions of the model are evaluated with residual plots.

Under the assumption of random effects the model is estimated using REML and the variance components are shown. Missing values are not estimated in this case.

Value

It returns an explanation about the RCBD fitted model.

Author(s)

Raul Eyzaguirre.

Examples

repo.rcbd(c("trw", "vw"), "geno", "genotype", "rep", dfr = pjpz09)
# With a small data set
temp <- pjpz09[1:18, ]
repo.rcbd(c("trw", "vw", "crw"), "geno", "genotype", "rep", dfr = temp)

reyzaguirre/pepa documentation built on Sept. 24, 2023, 7:08 p.m.