lin_ic_plot {PNAR}R Documentation

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

Description

Scatter plot of information criteria versus the number of lags in the linear Poisson Network Autoregressive model of order pp with qq covariates (PNAR(pp)).

Usage

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

Arguments

y

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

W

The NN x NN 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 NN x qq matrix of covariates (one for each column), where qq 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 linear Poisson Network Autoregressive model of order pp with qq covariates (PNAR(pp)).

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

lin_estimnarpq, log_lin_ic_plot

Examples

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

[Package PNAR version 1.6 Index]