R/plotPTC.R

In newrickey/mrmctools: Statistical Tools to Support Multi-Reader, Multi-Case (MRMC) Study Designs

#' Plot rating histogram.
#' 
#' This function plots a stacked histogram indicating the lesion localizations and non-lesion localizations by confidence rating. 
#' The function returns a ggplot object.
#' 
#' @param JAFROCfilename The is the excel file produced by the makeJAFROCfile function. This file is the source of most JAFROC analyses. 
#' @return \code{ggplot2} figure object.
#' @export


plotPTC <- function(JAFROCfilename){
datafilename <- JAFROCfilename

TPsaved <- readxl::read_excel(datafilename, "TP")

TPsaved$Source <- "Lesion Localization"

TPsaved <- dplyr::rename(TPsaved, Rating = TP_Rating)


FPsaved <- readxl::read_excel(datafilename,  "FP")
FPsaved$Source <- "Non-lesion Localization"
FPsaved <- dplyr::rename(FPsaved, Rating = FP_Rating)


datamarks <- dplyr::bind_rows(TPsaved, FPsaved)



ratinghist <- ggplot(aes(x = Rating ) , data = datamarks) + 
  geom_histogram(aes(fill = Source ), binwidth=5, colour="grey20", lwd=0.2) +
  #  stat_bin(binwidth=5, geom="text", colour="black", size=3.5, aes(label=..count.., fill=Source, y=.75*(..count..))) +
  scale_x_continuous(breaks=seq(0,100, 5)) + ylab("Frequency")+
  theme_bw() +
  theme( legend.position = c(0.40, 0.85), legend.text=element_text(size=12))  +
  guides(fill=guide_legend(title="Reader Mark Classification"))

return(ratinghist)
}