The last columns with cal_ variable names were calibrated with fsqca and show a marked difference in the box plot.

This plot shows why GNIPC is the major predictor of Productivity since the distributions are closely related to each other.

Using clustering (findTh) on each column

Used hierarchical cluster function findTh() in R to set thresholds in R function calibrate