Implement filter_covariance option in Bayesian estimation

Currently only possible in ML and calibrated smoother