repo.f: Automatic report for a factorial experiment

In reyzaguirre/pepa: Package for the Execution of Pre Cooked Analysis

Automatic report for a factorial experiment

Description

Produces an automatic report for selected variables in a factorial experiment with a CRD or RCBD.

Usage

repo.f(
  dfr,
  vars,
  factors,
  rep = NULL,
  maxp = 0.1,
  pe = FALSE,
  se = FALSE,
  title = "Automatic report for a factorial experiment",
  subtitle = NULL,
  author = "International Potato Center",
  format = c("html", "word", "pdf")
)

Arguments

dfr	The name of the data frame containing the data.
vars	The variables to analize.
factors	The factors.
rep	The replications or blocks, default is NULL for a CRD.
maxp	Maximum allowed proportion of missing values to estimate, default is 10%.
pe	Logical. If "pe = TRUE" multiple comparison tests for principal effects are included even if interaction is significat, only in the case of 2 factors. Default to "pe = FALSE".
se	Logical. If "se = TRUE" multiple comparison tests for simple effects are included even if interaction is not significat, only in the case of 2 factors. Default to "se = 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 factorial experiment with a CRD or RCBD for the selected variables. If data is unbalanced, missing values are estimated up to an specified maximum proportion, 10% by default. All factors are considered as fixed.

Value

It returns an automatic report for the factorial experiment with a CRD or RCBD fitted model.

Author(s)

Raul Eyzaguirre.

Examples

## Example 1: Two factors factorial
repo.f(asc, c("asc.dw", "asc.fw"), c("geno", "treat"))

## Example 2: Three factors factorial
# Create design with 3 factors
fnames <- c('A', 'B', 'C')
flevels <- list(c('a1', 'a2'), c('b1', 'b2'), c('c1', 'c2', 'c3'))
tmp <- cr.f(fnames, flevels, 'rcbd', 2, 10)$book
# Simulate random data
tmp$y <- rnorm(24)
# Run report
repo.f(tmp, 'y', c('A', 'B', 'C'), 'block')

reyzaguirre/pepa documentation built on March 29, 2025, 9:56 p.m.