SSPD3: Simulated data for a split-split-plot experimental design.
In modeldiagramR: Generate Model Diagrams for Linear Mixed Effect Models

SSPD3

R Documentation

Simulated data for a split-split-plot experimental design.

Description

The data is simulated using simulate_SSPD3_data()

Usage

SSPD3

Format

a data.frame with 120 observations on 13 variables:

ID: Observation/split-split-plot identification number, numeric
LOCATION: One of three locations "A", "B", and "C", a blocking variable, categorical
PLOT: Plot identification variable, one for each location and whole plot combination, numeric
REP: Replication number, all 1 as there was no replication, integer
WHOLE_PLOT: Treatment, whether or not the plot was irrigated "IRR_NO" for not irrigated and "IRR_YES" for irrigated, categorical
SPLIT_PLOT: Treatment, fungicide applied to split-plot, 4 fungicides and one control, "Fung1", ..., "Fung4", and "NFung", categorical
SPLIT_SPLIT_PLOT: Treatment, variety of bean seeded in split-split-plot, 4 varieties, "Beans1", ..., "Beans4", categorical
TRT_COMB: Treatments combined over WHOLE_PLOT, SPLIT_PLOT, and SPLIT_SPLIT_PLOT
RESP: Simulated response, not based on a meaningful input values, numeric
LocationF: Factor version of LOCATION
Irrigation: Factor version of WHOLE_PLOT treatment
Fungicide: Factor version of SPLIT_PLOT treatment
Variety: Factor version of SPLIT_SPLIT_PLOT treatment

References

Murillo D, Gezan S (2024). FielDHub: A Shiny App for Design of Experiments in Life Sciences. R package version 1.4.2, https://CRAN.R-project.org/package=FielDHub.

Examples

data("SSPD3")

# Simple model for testing
library(nlme)
lme_simple <- lme(RESP ~ Fungicide + Variety, random=~1|LOCATION/WHOLE_PLOT,
                     data=SSPD3)
model_diagram(lme_simple)


# Intended full model
aov_results <- aov(RESP ~ Irrigation*Fungicide*Variety + Error(LOCATION/WHOLE_PLOT/SPLIT_PLOT),
                   data=SSPD3)
summary(aov_results)
library(lme4)
library(modeldiagramR)
lmer_results <- lmer(RESP ~ Irrigation*Fungicide*Variety + (1|LOCATION/WHOLE_PLOT/SPLIT_PLOT),
                     data=SSPD3)
summary(lmer_results)

anova(lmer_results)

model_diagram(lmer_results)

model_diagram(lmer_results, width = 800, height=400, orientation="horizontal",
              shiftFixed = 2, shiftRandom = 5, scaleFontSize = 2)

modeldiagramR documentation built on April 15, 2026, 5:07 p.m.