Join cluster-level statistics into proc dataframe
add_stat.RdMerges the output of `stat_anom()` into the anomaly dataframe based on `cluster` ID. The added columns include `t_cen`, `N_anom`, and `t_lag`, which are relocated to appear immediately after the `cluster` column.