adjusted_boxplot: Adjusted Boxplot

In ChristianGoueguel/specProc: Preprocessing Tools For Laser-Induced Breakdown Spectroscopy (LIBS) Spectra

Adjusted Boxplot

Description

This function generates the adjusted boxplot, which is a robust graphical method for visualizing skewed data distributions. It provides a more accurate representation of the data's spread and skewness compared to standard boxplot, especially in the presence of outliers.

Usage

adjusted_boxplot(
  x,
  plot = TRUE,
  xlabels.angle = 90,
  xlabels.vjust = 1,
  xlabels.hjust = 1,
  box.width = 0.5,
  notch = FALSE,
  notchwidth = 0.5,
  staplewidth = 0.5
)

Arguments

x	A numeric data frame or tibble.
plot	A logical value indicating whether to plot the adjusted boxplot (default is TRUE).
xlabels.angle	A numeric value specifying the angle (in degrees) for x-axis labels (default is 90).
xlabels.vjust	A numeric value specifying the vertical justification of x-axis labels (default is 1).
xlabels.hjust	A numeric value specifying the horizontal justification of x-axis labels (default is 1).
box.width	A numeric value specifying the width of the boxplot (default is 0.5).
notch	A logical value indicating whether to display a notched boxplot (default is FALSE).
notchwidth	A numeric value specifying the width of the notch relative to the body of the boxplot (default is 0.5).
staplewidth	A numeric value specifying the width of staples at the ends of the whiskers.

Details

The function is based on the medcouple (MC) measure computed on the data and which robustly measures skewness. This measure is bounded between −1 and 1. The medcouple is equal to zero when the observed distribution is symmetric, whereas a positive (resp. negative) value of MC corresponds to a right (resp. left) tailed distribution. It worth noting that this method is more appropriate for distributions that are not excessively skewed i.e., for |\text{MC}| \leq 0.6.

Value

• If plot = TRUE, returns a ggplot2 object containing the adjusted boxplot.

• If plot = FALSE, returns a list of tibbles with the adjusted boxplot statistics and potantial outliers.

Author(s)

Christian L. Goueguel

References

The adjusted boxplot is based on the methodology described in:

• Brys, G., Hubert, M., Struyf, A., (2004). A Robust Measure of Skewness. Journal of Computational and Graphical Statistics, 13(4):996-1017

• Hubert, M., Vandervieren, E., (2008). An adjusted boxplot for skewed distributions. Computational Statistics and Data Analysis, 52(12):5186-5201

Examples

set.seed(123)
data <- data.frame(
  normal = rnorm(100),
  skewed = rexp(100, rate = 0.5),
  heavy_tailed = rt(100, df = 3)
)

# Plot the adjusted boxplot
adjusted_boxplot(data)

# Retrieve the adjusted boxplot statistics
adjusted_boxplot(data, plot = FALSE)

ChristianGoueguel/specProc documentation built on Nov. 9, 2024, 3:23 p.m.