Implements a framework for creating boxplots where the
whisker lengths are determined by formal multiple testing procedures,
making them adaptive to sample size and data characteristics. The function bh_boxplot()
generates boxplots that control the False Discovery
Rate (FDR) via the Benjamini-Hochberg procedure, and the function holm_boxplot() generates
boxplots that control the Family-Wise Error Rate (FWER) via the Holm procedure.
The methods are based on the framework in Gang, Lin, and Tong (2025)