Evaluating exact VARMA likelihood and its gradient when data are incomplete