D is matrix of lagged Y (covariates) Y is observations
update_MLE_asymmetric(D, U, V, W, Z, DDT, DTD, DYT, DTY)