R/familyScorePctID.R

In serratus-bio/tantalus: R Data Exploration for Serratus

#' Plot a scatterplot/histogram of percentage identity vs family score for a given family
#'
#'
#' @param  family_df family data frame
#' @param  family_name name of the family
#'
#' @return ggplot
#'
#' @examples
#' 
#' familyScorPctID(ZOO, "Coronaviridae")
#' 
#' @export
#'
familyScorPctID <- function(family_df, family_name){
  family_df <- family_df %>% filter(family_name %in% c(family_name))
  
  p <- ggplot(data=family_df, aes(x = score, y = percent_identity)) +
    geom_jitter(color="blue", alpha = 0.4) + 
    theme_bw() +
    labs(x = "Family score", y = "Percentage identity") +
    labs(title = paste(family_name, "family", sep = " "))
  s <- ggExtra::ggMarginal(p, type = "histogram")
  return(s)
}