Plotting varmaPredict Forecasts

/*
** Plot differenced data and forecast
*/
 
// Declare 'myPlot' to be a plotControl structure
// and fill with default values
struct plotControl myPlot;
myPlot = plotGetDefaults("xy");
 
// Turn grid off
plotSetGrid(&myPlot, "off");
 
// Font settings
font_name = "Helvetica Neue";
font_color = "black";
 
plotSetTitle(&myPlot, "Differenced data and VAR(2) model forecast", 
             font_name, 18, font_color);
 
// Add a tick label every 10 years, starting at 1860
plotSetXTicInterval(&myPlot, 10, 1860);
 
// Legend
plotSetLegend(&myPlot, "Y1 Obs." $| "Y2 Obs.");
plotSetLegendFont(&myPlot, font_name, 14, font_color);
 
// Plot differenced 'y' data
plotTS(myPlot, 1850, 1, y);
 
// Set line style to dots for forecasted data
plotSetLineStyle(&myPlot, 3);
 
// Legend
plotSetLegend(&myPlot, "Y1 Forecasts"$|"Y2 Forecasts");
plotSetLegendFont(&myPlot, font_name, 14, font_color);
 
// Add forecasts to plot
plotAddTS(myPlot, 1910, 1, f[.,2:3]);