Code now works, only use sigma in observation likelihood covariance