mul2.errorbar(
    range(1, n_nonzero_coefs + 1),
    rmse_uniform3.mean(1),
    yerr=rmse_uniform3.std(1),