Maximum likelihood estimation of bivariate circular hidden Markov models from incomplete data