log_lin_ic_plot {PNAR}

Scatter plot of information criteria versus the number of lags in the log-linear Poisson NAR(p) model with p lags and q covariates (log-PNAR(p))

Description

Scatter plot of information criteria versus the number of lags in log-linear Poisson Network Autoregressive model of order p with q covariates (log-PNAR(p)).

Usage

log_lin_ic_plot(y, W, p = 1:10, Z = NULL, uncons = FALSE, ic = "QIC")

Arguments

y

A TT x N time series object or a TT x N numerical matrix with the N multivariate count time series over TT time periods.

W

The N x N row-normalized non-negative adjacency matrix describing the network. The main diagonal entries of the matrix should be zeros, all the other entries should be non-negative and the maximum sum of elements over the rows should equal one. The function row-normalizes the matrix if a non-normalized adjacency matrix is provided.

p

A vector with integer numbers, the range of lags in the model, for which the AIC, BIC and QIC will be computed.

Z

An N x q matrix of covariates (one for each column), where q is the number of covariates in the model. Note that they must be non-negative.

uncons

Logical, if TRUE an unconstrained optimization without stationarity constraints is performed (default is FALSE).

ic

The information criterion you want to plot, "QIC" (default value), "AIC" or "BIC".

Details

The function computes the AIC, BIC or QIC for a range of lag orders of the log-linear Poisson Network Autoregressive model of order p with q covariates (PNAR(p)).

Value

A scatter plot with the lag order versus either QIC (default), AIC or BIC, and a vector with their values, for each lag order.

Author(s)

Mirko Armillotta, Michail Tsagris and Konstantinos Fokianos.

References

Armillotta, M. and K. Fokianos (2022). Poisson network autoregression. https://arxiv.org/abs/2104.06296

See Also

log_lin_estimnarpq, lin_ic_plot

Examples

data(crime)
data(crime_W)
log_lin_ic_plot(crime, crime_W, p = 1:3)