On a method of estimating parameters in non-negative multivariate autoregressive processes