Integrating Antonio's code (untested)
edited the code to convert the multidimensional time series to an autoregressive form. checked that the resulting matrix is the same as original data
Rewrote the fitting code and added a plotting method
Updated some more of Antonio's code and added some model test plots
Added diagnostic plots and examples