add **kwargs to diagnostics.plot_recovery

[marvinschmitt]

allows modifying the key part of the plot, i.e., the scatter (with or without error bars).

Example:

_ = plot_recovery(post_samples, test_data['parameters'], param_names=parameter_names, fig_size=(16, 4), rasterized=True)