lin_narpq_init {PNAR}R Documentation

Starting values for the linear Poisson NAR(p) model model with p lags and q covariates (PNAR(p))

Description

Starting values for the linear Poisson Network Autoregressive model of order pp with qq covariates (PNAR(pp)).

Usage

lin_narpq_init(y, W, p, Z = NULL)

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

The number of lags in the model.

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.

Details

The function computes starting values to be used in the function lin_estimnarpq. These are simply the ordinary least squares estimators with a correction. If any of the the resulting coefficients is negative they become equal to 0.001

Value

A vector with the initial values.

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

Examples

data(crime)
data(crime_W)
x0 <- lin_narpq_init(crime, crime_W, p = 2)

[Package PNAR version 1.6 Index]