Inference for Network Count Time Series with the R Package PNAR

We introduce a new R package useful for inference about network count time series. Such data are frequently encountered in statistics and they are usually treated as multivariate time series. Their statistical analysis is based on linear or log-linear models. Nonlinear models, which have been applied successfully in several research areas, have been neglected from such applications mainly because of their computational complexity. We provide R users the flexibility to fit and study nonlinear network count time series models which include either a drift in the intercept or a regime switching mechanism. We develop several computational tools including estimation of various count Network Autoregressive models and fast computational algorithms for testing linearity in standard cases and when non-identifiable parameters hamper the analysis. Finally, we introduce a copula Poisson algorithm for simulating multivariate network count time series. We illustrate the methodology by modeling weekly number of influenza cases in Germany.

Mirko Armillotta (Vrije Universiteit Amsterdam) , Michail Tsagris (University of Crete) , Konstantinos Fokianos (University of Cyprus)

0.1 Supplementary materials

Supplementary materials are available in addition to this article. It can be downloaded at

M. Armillotta and K. Fokianos. Count network autoregression. Journal of Time Series Analysis, jtsa.12728: 2023a. URL
M. Armillotta and K. Fokianos. Nonlinear Network Autoregression. The Annals of Statistics, 51(6): 2526–2552, 2023b. URL
J. Bracher and L. Held. Endemic-epidemic models with discrete-time serial interval distributions for infectious disease prediction. International Journal of Forecasting, 38: 1221–1233, 2020.
R. Brent. Algorithms for minimization without derivatives. Prentice-Hall, 1973.
G. Csardi and T. Nepusz. The igraph software package for complex network research. InterJournal, Complex Systems: 1695, 2006. URL
R. Davies. Hypothesis testing when a nuisance parameter is present only under the alternative. Biometrika, 74: 33–43, 1987.
P. Erdös and A. Rényi. On random graphs I. Publicationes Mathematicae Debrecen, 6: 1959.
K. Fokianos. Multivariate count time series modelling. To appear in Econometrics and Statistics, 2022.
K. Fokianos, B. Støve, D. Tjøstheim and P. Doukhan. Multivariate count autoregression. Bernoulli, 26: 471–499, 2020.
C. Gourieroux, A. Monfort and A. Trognon. Pseudo maximum likelihood methods: theory. Econometrica, 52: 681–700, 1984.
B. Hansen. Inference when a nuisance parameter is not identified under the null hypothesis. Econometrica, 64: 413–430, 1996.
L. Held, M. Höhle and M. Hofmann. A statistical framework for the analysis of multivariate infectious disease surveillance counts. Statistical Modelling, 5: 187–199, 2005.
L. Held and M. Paul. Modeling seasonality in space-time infectious disease surveillance data. Biometrical Journal, 54: 824–843, 2012.
D. Kraft. Algorithm 733: TOMPFortran modules for optimal control calculations. ACM Transactions on Mathematical Software (TOMS), 20: 262–281, 1994.
K. Leeming, G. Nason, M. Nunes and J. Wei. GNAR: Methods for fitting network time series models. 2023. URL R package version 1.1.2.
S. Meyer and L. Held. Power-law models for infectious disease spread. The Annals of Applied Statistics, 8: 1612–1639, 2014.
S. Meyer, L. Held and M. Höhle. Spatio-temporal analysis of epidemic phenomena using the R package surveillance. Journal of Statistical Software, 77: 1–55, 2017.
W. Pan. Akaike’s information criterion in generalized estimating equations. Biometrics, 57: 120–125, 2001.
M. Papadakis, M. Tsagris, M. Dimitriadis, S. Fafalios, I. Tsamardinos, M. Fasiolo, G. Borboudakis, J. Burkardt, C. Zou, K. Lakiotaki, et al. Rfast: A collection of efficient and extremely fast r functions. 2023a. URL R package version 2.0.8.
M. Papadakis, M. Tsagris, S. Fafalios and M. Dimitriadis. Rfast2: A collection of efficient and extremely fast r functions II. 2023b. URL R package version
M. Paul and L. Held. Predictive assessment of a non-linear random effects model for multivariate time series of infectious disease counts. Statistics in Medicine, 30: 1118–1136, 2011.
M. Paul, L. Held and A. M. Toschke. Multivariate modelling of infectious disease surveillance data. Statistics in Medicine, 27: 6250–6267, 2008.
M. Tsagris, M. Armillotta and K. Fokianos. PNAR: Poisson Network Autoregressive Models. 2023. URL R package version 1.6.
Y. J. Wang and G. Y. Wong. Stochastic blockmodels for directed graphs. Journal of the American Statistical Association, 82: 8–19, 1987.
R. W. Wedderburn. Quasi-likelihood functions, generalized linear models, and the Gauss-Newton method. Biometrika, 61: 439–447, 1974.
J. Ypma and S. G. Johnson. Nloptr: R interface to NLopt. 2022. URL R package version 2.0.3.
X. Zhu, R. Pan, G. Li, Y. Liu and H. Wang. Network vector autoregression. The Annals of Statistics, 45: 1096–1123, 2017.



Text and figures are licensed under Creative Commons Attribution CC BY 4.0. The figures that have been reused from other sources don't fall under this license and can be recognized by a note in their caption: "Figure from ...".


For attribution, please cite this work as

Armillotta, et al., "Inference for Network Count Time Series with the R Package PNAR", The R Journal, 2024

BibTeX citation

  author = {Armillotta, Mirko and Tsagris, Michail and Fokianos, Konstantinos},
  title = {Inference for Network Count Time Series with the R Package PNAR},
  journal = {The R Journal},
  year = {2024},
  note = {},
  doi = {10.32614/RJ-2023-094},
  volume = {15},
  issue = {4},
  issn = {2073-4859},
  pages = {255-269}