powered by
Perform threshold analysis
perform_threshold_analysis(water_sf, variable, thresholds, verbose)