Vector autoregressive moving-average (VARMA) processes are suitable models for producing linear forecasts of sets of time series variables. They provide parsimonious representations of linear data generation processes (DGPs). The setup for these processes in the presence of cointegrated variables is considered. Moreover, a unique or identified parameterization based on the echelon form is presented. Model specification, estimation, model checking and forecasting are discussed. Special attention is paid to forecasting issues related to contemporaneously and temporally aggregated processes.