pysteps.timeseries.autoregression.estimate_var_params_yw#
- pysteps.timeseries.autoregression.estimate_var_params_yw(gamma, d=0, check_stationarity=True)#
Estimate the parameters of a VAR(p) model
\(\mathbf{x}_{k+1}=\mathbf{\Phi}_1\mathbf{x}_k+ \mathbf{\Phi}_2\mathbf{x}_{k-1}+\dots+\mathbf{\Phi}_p\mathbf{x}_{k-p}+ \mathbf{\Phi}_{p+1}\mathbf{\epsilon}\)
from the Yule-Walker equations using the given correlation matrices \(\mathbf{\Gamma}_0,\mathbf{\Gamma}_1,\dots,\mathbf{\Gamma}_n\), where n=p.
- Parameters:
gamma (list) – List of correlation matrices \(\mathbf{\Gamma}_0,\mathbf{\Gamma}_1,\dots,\mathbf{\Gamma}_n\). To obtain these matrices, use
pysteps.timeseries.correlation.temporal_autocorrelation_multivariate()
with window_radius=np.inf.d ({0,1}) – The order of differencing. If d=1, the correlation coefficients gamma are assumed to be computed from the differenced time series, which is also done for the resulting parameter estimates.
check_stationarity (bool) – If True, the stationarity of the resulting VAR(p) process is tested. An exception is thrown if the process is not stationary.
- Returns:
out – List of VAR(p) coefficient matrices \(\mathbf{\Phi}_1, \mathbf{\Phi}_2,\dots\mathbf{\Phi}_{p+1}\), where the last matrix corresponds to the innovation term.
- Return type:
list
Notes
To estimate the parameters of an integrated VARI(p,d) model, compute the correlation coefficients gamma by calling
pysteps.timeseries.correlation.temporal_autocorrelation_multivariate()
with d>0. Estimation of the innovation parameter \(\mathbf{\Phi}_{p+1}\) is not currently implemented, and it is set to a zero matrix.