This package provides functions for estimation, testing and diagnostic checking of generalized lin-ear autoregressive moving average (GLARMA) models for discrete valued time series with regres-sion variables. These are a class of observation driven non-linear non-Gaussian state space mod-els. The state vector consists of a linear regression component plus an observation driven compo-nent consisting of an autoregressive-moving average (ARMA) filter of past predictive residu-als. Currently three distributions (Poisson, negative binomial and binomial) can be used for the re-sponse series. Three options (Pearson, score-type and unscaled) for the residuals in the observa-tion driven component are available. Estimation is via maximum likelihood (...