sa_contour <- function(df,
L = seq(-40, 40, by = 10),
title = "Spray Angle",
NCOL = 2){
if(is.data.frame(df) == TRUE) {
df <- list(df)
names(df) <- "Group"
}
N_df <- length(df)
if(is.list(df) == TRUE){
if(length(names(df)) == 0){
names(df) <- paste("Group", 1:N_df)
}
}
df_p <- NULL
for(j in 1:N_df){
df[[j]] %>%
setup_inplay() %>%
sa_gam_fit() %>%
grid_predict() %>%
mutate(Group = names(df)[j]) -> df_c
df_p <- rbind(df_p, df_c)
}
if(N_df == 1){
contour_graph(df_p, L, title)
} else {
contour_graph(df_p, L, title) +
facet_wrap(~ Group, ncol = NCOL)
}
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.