Caps extreme outliers at specified percentiles. Crucially, this function allows for stratified winsorization (e.g., by region) to avoid over-trimming poor/rich areas, and it creates new `_w` suffixed columns to preserve raw data provenance.
Usage
ihs_winsorize(data, vars, by = NULL, probs = c(0.01, 0.99))