give.n <- function(x){
return(c(y = -0.05, label = length(x)))
}
#' Courtship boxplots
#' A function to plot courtship conditioning results as boxplots with pairwise Wilcoxon Rank Sum tests
#'
#' @param dat A data frame generated from loadData() followed by findOutlier()
#' @return A boxplot with pairiwise Wilcoxon Rank Sum test p-values
#' @export
ccBoxplots <- function(dat){
stat <- dat %>% dplyr::group_by(genotype) %>%
rstatix::wilcox_test(CI ~ condition, paired=F) %>%
rstatix::add_xy_position(x="condition")
ggplot2::ggplot(dat,ggplot2::aes(y=CI,x=condition)) + ggplot2::geom_boxplot() + ggplot2::facet_grid(cols=ggplot2::vars(genotype)) +
ggpubr::stat_pvalue_manual(stat, tip.length = 0) + ggplot2::stat_summary(fun.data = give.n,geom="text") +
ggplot2::scale_y_continuous(expand = c(0,0.1)) + ggplot2::labs(x=NULL)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.