This paper considers the maximum likelihood estimation (MLE) of a class of stationary and invertible vector autoregressive fractionally integrated moving-average (VARFIMA) processes considered in (26) of Luceño [1] or Model A of Lobato [2] where each component yi,t is a fractionally integrated process of order di, i = 1, . . . , r. Under the conditions outlined in Assumption 1 of this paper, th...