I am creating a grouped boxplot in plotly using predefined quantiles. I want to color the lines around the boxes based on a separate variable. I can't seem to do this directly in my plotly call. There is a nice solution here which involves changing the line colours as a post-processing step using plotly_build. The example in that link works well but the structure of the data is different when using pre-defined quantiles, and I can't seem to access the data as in that example. Perhaps there is a way but I can't figure it out.
My attempted solution is the following, which involves adding a new trace, but with transparent fill as follows:
library(plotly)
library(dplyr)
# CREATE DUMMY DATA
set.seed(123) # Set seed for reproducibility
site_name <- rep(paste0("site_", 1:5), each = 40) # Create the site_name column with 5 different site names, each with 20 rows
site_type <- rep(c("A", "B"), each = 20, times = 5) # Create the site_type column with 10 'A's and 10 'B's for each site
value <- runif(100, min = 0, max = 200) # Create the value column with random numbers (between 0 and 100).
# Combine into a data frame
df <- data.frame(site_name, site_type, value)
# Generate site_status randomly for each combination of site_name and site_type
unique_combinations <- unique(df[c("site_name", "site_type")])
unique_combinations$site_status <- sample(c("Good", "Bad"), nrow(unique_combinations), replace = TRUE)
# Merge site_status back to the original df
df <- df %>%
left_join(unique_combinations, by = c("site_name", "site_type"))
# Display the first few rows of the dataset
head(df, 20)
# MAKE SUMMARY DATA
# Group by site_name and site_type, then calculate summary statistics
stats_df <- df %>%
group_by(site_name, site_type) %>%
summarise(
lower_fence = quantile(value, probs = c(0.05), type = 5, na.rm = TRUE),
q1 = quantile(value, probs = c(0.25), type = 5, na.rm = TRUE),
median = quantile(value, probs = c(0.5), type = 5, na.rm = TRUE),
mean = mean(value, na.rm = TRUE),
q3 = quantile(value, probs = c(0.75), type = 5, na.rm = TRUE),
upper_fence = quantile(value, probs = c(0.95), type = 5, na.rm = TRUE),
sd = sd(value, na.rm = TRUE),
site_status = unique(site_status),
.groups = 'drop'
)
# PLOTTING CODE
# make box plot
fig <- plot_ly(data = stats_df,
x = ~site_name,
color = ~site_type, # boxes
colors = c("blue","red"),
type = "box",
lowerfence = ~lower_fence,
q1 = ~q1,
median = ~median,
q3 = ~q3,
upperfence = ~upper_fence) %>%
layout(boxmode = "group", boxgap = 1/5)
# Filter out the boxes to be drawn with green boxes around them
bad_data <- stats_df %>% filter(site_status == "Bad")
# Add green boxes
fig <- fig %>% plotly::add_trace(
x = factor(bad_data$site_name),
color = factor(bad_data$site_type),
colors = c("blue","red"),
type = "box",
lowerfence = bad_data$lower_fence,
q1 = bad_data$q1,
median = bad_data$median,
q3 = bad_data$q3,
upperfence = bad_data$upper_fence,
line = list(color = "green"),
fillcolor = "rgba(255,0,0,0.0)", # Red with transparency
boxmean = FALSE, # Avoid adding box means
showlegend = TRUE,
inherit = FALSE
)
# Show the figure
fig
This works like a charm when the data aren't grouped, but unfortunately when using a grouped boxplot it creates the box outlines but as separate grouped items and changes the colors in the existing boxes as in this image:
I'm not sure if the grouping attribute can be forced somehow. TBH I'm not really sure if this approach is viable.
I would love to be able to change the line colors using the following kind of approach as per the example referenced above, but it doesn't seem to work with pre-defined quantiles:
built_fig <- plotly_build(built_fig)
lapply(1:length(stats_df$site_status),
function(i){
nm = stats_df$site_status[i]
cr = ifelse(nm == "Good",
"#66FF66", "black")
built_fig$x$data[[i]]$line$color <<- cr # change graph by age
}
)
Any suggestions greatly appreciated.
