Run the Shiny App for performing one-way ANOVA on data from a Randomized block design
Source:R/twoway_rbd.R
twoway_rbd.RdApp allows user to import data through excel/.csv files or through clipboard and user can select the post-hoc test method and download the report which contains anova results and plots
Value
This function runs a local instance of the Shiny app in your default web browser. The app interface allows users to upload data, select analysis method, and download outputs.
Details
This Shiny app is part of the visvaR package and is designed for
analysis of variance randomized block design (two factor) and user can download the report in word format.To use custom fonts, please install the extrafont package and run extrafont::font_import() and extrafont::loadfonts().
References
Fisher, R. A. (1925). Statistical Methods for Research Workers. Oliver and Boyd, Edinburgh. Scheffe, H. (1959). The Analysis of Variance. John Wiley & Sons, New York. R Core Team (2024). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/
Examples
# Example 1: Basic usage
if(interactive()) {
twoway_rbd()
}
# Example 2: Sample workflow with factorial fertilizer experiment
if(interactive()) {
# Prepare sample data with two factors: Nitrogen and Phosphorus levels
fertilizer_data <- data.frame(
Nitrogen = rep(c("N0", "N30", "N60"), each = 12),
Phosphorus = rep(rep(c("P0", "P30", "P60", "P90"), each = 3), times = 3),
Replication = rep(1:3, times = 12),
Grain_Yield = c(
4.2, 4.0, 4.1, # N0-P0
4.8, 4.6, 4.7, # N0-P30
5.1, 5.0, 5.2, # N0-P60
5.0, 4.9, 5.1, # N0-P90
5.5, 5.3, 5.4, # N30-P0
6.2, 6.0, 6.1, # N30-P30
6.8, 6.6, 6.7, # N30-P60
6.7, 6.5, 6.6, # N30-P90
6.0, 5.8, 5.9, # N60-P0
6.8, 6.6, 6.7, # N60-P30
7.5, 7.3, 7.4, # N60-P60
7.4, 7.2, 7.3 # N60-P90
),
Protein_Content = c(
9.0, 8.8, 8.9, # N0-P0
9.5, 9.3, 9.4, # N0-P30
9.8, 9.6, 9.7, # N0-P60
9.7, 9.5, 9.6, # N0-P90
10.5, 10.3, 10.4, # N30-P0
11.2, 11.0, 11.1, # N30-P30
11.8, 11.6, 11.7, # N30-P60
11.7, 11.5, 11.6, # N30-P90
12.0, 11.8, 11.9, # N60-P0
12.8, 12.6, 12.7, # N60-P30
13.5, 13.3, 13.4, # N60-P60
13.4, 13.2, 13.3 # N60-P90
)
)
# Save as Excel file
write.xlsx(fertilizer_data, "fertilizer_data.xlsx")
# Launch the app
twoway_rbd()
# Instructions for users:
# 1. Click "Choose .xlsx or .csv" and select fertilizer_data.xlsx
# 2. Select post-hoc test method (e.g., "LSD" or "Tukey HSD")
# 3. Customize plot appearance:
# - Select font style
# - Adjust font size
# 4. Set output filename
# 5. Click "Analyze"
# 6. View results in different tabs
# 7. Download comprehensive Word report
# Clean up
unlink("fertilizer_data.xlsx")
}
# Example 3: Using clipboard data
if(interactive()) {
# Copy this to clipboard:
# Temperature,Light,Replication,Growth_Rate,Chlorophyll
# Low,Dark,1,2.1,15.2
# Low,Dark,2,2.0,15.0
# Low,Dark,3,2.2,15.1
# Low,Medium,1,2.8,16.5
# Low,Medium,2,2.7,16.3
# Low,Medium,3,2.9,16.4
# Low,High,1,3.2,17.8
# Low,High,2,3.1,17.6
# Low,High,3,3.3,17.7
# High,Dark,1,2.5,14.8
# High,Dark,2,2.4,14.6
# High,Dark,3,2.6,14.7
# High,Medium,1,3.5,16.2
# High,Medium,2,3.4,16.0
# High,Medium,3,3.6,16.1
# High,High,1,4.2,17.5
# High,High,2,4.1,17.3
# High,High,3,4.3,17.4
twoway_rbd()
# Click "Use Clipboard Data" after copying data
}
# Example 4: Multiple response variables with variety trial
if(interactive()) {
# Create data with multiple responses
variety_trial <- data.frame(
Irrigation = rep(c("Full", "Deficit"), each = 18),
Variety = rep(rep(c("V1", "V2", "V3"), each = 6), times = 2),
Replication = rep(1:6, times = 6),
Yield = rnorm(36, mean = rep(c(5.5, 5.0, 4.8, 4.2, 3.8, 3.5), each = 6), sd = 0.2),
WUE = rnorm(36, mean = rep(c(2.2, 2.4, 2.3, 2.8, 3.0, 2.9), each = 6), sd = 0.1),
Biomass = rnorm(36, mean = rep(c(12.5, 11.8, 11.2, 10.2, 9.8, 9.5), each = 6), sd = 0.5)
)
# Save as Excel file
write.xlsx(variety_trial, "variety_trial.xlsx")
# Launch the app
twoway_rbd()
# Clean up
unlink("variety_trial.xlsx")
}