
Combined Split-Split-Plot Analysis over Environments
Source:R/design_split_plot.R
anova_sspd_pooled.RdCombined (pooled) analysis of a split-split-plot experiment repeated in
several environments, with environments fixed and replications random
within environments. See anova_spd_pooled.
Usage
anova_sspd_pooled(
data,
response,
main_plot,
sub_plot,
sub_sub_plot,
environment,
replication,
verbose = TRUE,
posthoc = "lsd",
alpha = 0.05
)Arguments
- data
Data frame containing the combined data.
- response
Name of the response variable.
- main_plot
Name of the main-plot factor.
- sub_plot
Name of the sub-plot factor.
- sub_sub_plot
Name of the sub-sub-plot factor.
- environment
Name of the environment column.
- replication
Name of the replication column.
- verbose
Logical. If TRUE (default), prints the report.
- posthoc
Method for the letters (default
"lsd").- alpha
Significance level. Default 0.05.
Value
An object of class aridagri_sspd_pooled; see
anova_spd_pooled. Also contains subsub_means.
Examples
data <- expand.grid(env = c("E1", "E2"), rep = 1:3, A = c("A1", "A2"),
B = c("B1", "B2"), C = c("C1", "C2"))
set.seed(8)
data$yield <- rnorm(nrow(data), 1200, 90)
res <- anova_sspd_pooled(data, "yield", "A", "B", "C", "env", "rep",
verbose = FALSE)
res$anova_table
#> Df Sum Sq Mean Sq F value Pr(>F)
#> env 1 22934.2419 22934.2419 1.6537586 0.26784180
#> rep(env) 4 55471.7992 13867.9498 NA NA
#> A 1 1997.9245 1997.9245 0.1718856 0.69971131
#> env:A 1 6097.5511 6097.5511 0.5245849 0.50898738
#> Error(a) 4 46494.2975 11623.5744 NA NA
#> B 1 1952.9792 1952.9792 0.1809301 0.68178614
#> env:B 1 19087.5742 19087.5742 1.7683328 0.22025261
#> A:B 1 3435.7088 3435.7088 0.3182948 0.58809762
#> env:A:B 1 1359.1013 1359.1013 0.1259114 0.73188212
#> Error(b) 8 86352.8608 10794.1076 NA NA
#> C 1 2193.2727 2193.2727 0.3390368 0.56850077
#> env:C 1 1189.5573 1189.5573 0.1838822 0.67377558
#> A:C 1 1296.2490 1296.2490 0.2003746 0.66041981
#> B:C 1 17548.5227 17548.5227 2.7126564 0.11905144
#> env:A:C 1 1083.4479 1083.4479 0.1674797 0.68778748
#> env:B:C 1 12333.2535 12333.2535 1.9064784 0.18633913
#> A:B:C 1 24610.3775 24610.3775 3.8042802 0.06886361
#> env:A:B:C 1 223.5847 223.5847 0.0345618 0.85485336
#> Error(c) 16 103506.0563 6469.1285 NA NA