PooledStPD | R Documentation |
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>.
PooledStPD(
data,
Response,
Location,
Replication,
RowFactor,
ColumnFactor,
alpha,
Mult_Comp_Test
)
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:
|
A list containing the following components:
Individual_ANOVA: Summary of ANOVA results for each location or year.
Location_wise: Multiple comparisons of interaction of RowFactor and ColumnFactor within each location or year.
Bartlett_Test: Results of Bartlett's test for homogeneity of variances.
Pooled_ANOVA: Combined (pooled) ANOVA table across all locations or years.
Interaction_Comparison: Summary of pooled interaction of RowFactor and ColumnFactor comparisons using the selected multiple comparison test..
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>.
# 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.