histval is broken in models with lag > 1
initialization doesn't take auxiliary variables into account initial matrix has wrong dimension