# Boxplot data
#
# Function to subset a data file containing only boxplot data
# As data is lognormal, we use logaritmic plotting. Therefore all data which value = 0 is omitted.
#' Function to obtain boxplot data
#' @import data.table
#'
#' @export
bxdata <- function(file) {
bx <- file[P>0]
bx<- within(bx,
STN <- factor(STN))
bx[, ':='(Qm = mean(P)), by=STN]
bx[, ':='(Q50 = quantile(P, probs=c(0.50))), by=STN]
bx[, ':='(Q99 = quantile(P, probs=c(0.99)), Q99.9 = quantile(P, probs=c(0.999))), by=STN]
bx <- subset(bx, select=c("STN", "P", "Qm", "Q99", "Q99.9"))
return(bx)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.