uni.errorbar(
    range(1, n_nonzero_coefs + 1),
    rmse_gaussian1.mean(1),
    yerr=rmse_gaussian1.std(1),