PooledStPD {CANE}R Documentation

Pooled Strip Plot Design Analysis

Description

This function conducts a pooled analysis of variance (ANOVA) using the strip plot design (StPD) for data collected across multiple locations or years. In this design, the interaction between factors (RowFactor and ColumnFactor) is estimated with higher precision. For more details see Dean et al. (2017)<doi:10.1007/978-3-319-52250-0> and Ruíz et al. (2024)<doi:10.1007/978-3-031-65575-3>.

Usage

PooledStPD(
  data,
  Response,
  Location,
  Replication,
  RowFactor,
  ColumnFactor,
  alpha,
  Mult_Comp_Test
)

Arguments

data

A data frame containing the experimental data.

Response

A numeric variable representing the dependent variable (response).

Location

A factor indicating different locations or years.

Replication

A factor indicating replications.

RowFactor

A factor used for horizontal strips.

ColumnFactor

A factor used for vertical strips.

alpha

A numeric value specifying the significance level for Bartlett’s test.

Mult_Comp_Test

An integer specifying the type of multiple comparison test:

  • 1 = Tukey's honestly significant difference (Tukey's HSD) test

  • 2 = Duncan's multiple range test (DMRT)

  • 3 = least significant difference (LSD) test

Value

A list containing the following components:

References

Dean A, Voss D, Draguljic D (2017)<doi:10.1007/978-3-319-52250-0>.

Ruíz JS, López OAM, Crossa J (2024)<doi:10.1007/978-3-031-65575-3>.

Examples

# Creating a sample dataset for Pooled Strip Plot Design (StPD)
df <- data.frame(
  Location = factor(rep(c("Londan", "Agumbe"), each = 12)),  # Locations
  Replication = factor(rep(c(1, 2), each = 6, times = 2)),  # Replications
  RowFactor = factor(rep(c(1, 2), each = 3, times = 4)),  # Row factor
  ColumnFactor = factor(rep(1:3, times = 8)),  # Column factor
  Yield = c(4940, 4810, 5150, 4900, 4920, 5070, 
            4830, 5110, 4920, 5020, 5110, 5230,
            4964, 4997, 5011, 5102, 4858, 4888, 
            5100, 5165, 4965, 5113, 5086, 5176)  # Yield values
)

# Running PooledStPD function on the dataset
out <- PooledStPD(df, "Yield", "Location", "Replication", "RowFactor", "ColumnFactor", 0.05, 1)

# Print results
print(out)

[Package CANE version 0.1.1 Index]