Rfssa: Rfssa: A Package for Functional Singular Spectrum Analysis...

RfssaR Documentation

Rfssa: A Package for Functional Singular Spectrum Analysis and Related Methods.


The Rfssa package provides the collection of necessary functions to implement functional singular spectrum analysis (FSSA)-based methods for analyzing univariate and multivariate functional time series (FTS). Univariate and multivariate FSSA are novel, non-parametric methods used to perform decomposition and reconstruction of univariate and multivariate FTS respectively. In addition, the FSSA-based routines may be performed on FTS whose variables are observed over a one or two-dimensional domain. Finally, one may perform FSSA recurrent or FSSA vector forecasting of univariate or multivariate FTS observed over one-dimensional domains. Forecasting of FTS whose variables are observed over domains of dimension greater than one is under development.


The use of the package begins by defining an fts object by providing the constructor with the raw data, basis specifications, and grid specifications. We note that the FTS object may be univariate or multivariate and variables may be observed over one (curves) or two-dimensional (images) domains. Validity checking of the S4 object constructor inputs is included to help guide the user. The user may leverage the plot.fts method to visualize the fts object. A variety of plotting options are available for variables observed over a one-dimensional domain and a visuanimation is offered for variables observed over a two-dimensional domain. Next, the user provides the fts object and a chosen lag parameter to the FSSA routine (fssa) to obtain the decomposition. We note that the decomposition function leverages the RSpectra and RcppEigen R packages, and the Eigen C++ package to speed up the routine. The plot.fssa method may be used to visualize the results of the decomposition and to choose an appropriate grouping of the eigentriples for reconstruction (freconstruct) or forecasting (fforecast). The freconstruct routine can be used to reconstruct a list of fts objects specified by the grouping while the fforecast function returns a list of fts objects that contain predictions of the signals specified by the grouping. We note that when forecasting is performed, usually the user specifies one group that captures the assumed deterministic, extracted signal that is found within the FTS and all other modes of variation are excluded. We also note that currently, forecasting only supports FTS whose variables are observed over a one-dimensional domain with two-dimensional domain forecasting to be added in the future.

Other functionalities offered by the package include:

  • FTS arithmetic - Allows the user to perform FTS-FTS arithmetic and FTS-scalar arithmetic (such as addition, subtraction, etc.).

  • eval.fts - Allows the user to evaluate the FTS object over a new specified grid.

  • load_github_data - Allows the user to load any .RData file hosted on GitHub including the Callcenter, Jambi, and Montana datasets.

  • fwcor - Returns the weighted correlation matrix corresponding to the decomposition of an FTS.

  • cor.fts - Returns the correlation between two fts objects.

  • launchApp - Launches the built-in R Shiny app that can be used to interactively explore the FSSA-based routines on various datasets.

The first update we include in this version of the Rfssa R package, is the eval.fts method used to evaluate an fts object over a new, specified grid. We updated the plot.fts method to allow for custom tick labels and new choices in visuanimation colors (for variables observed over two-dimensional domains) that are offered by the ggplot2 package. We have also updated the plot.fssa method to allow for new choices in visuanimation colors offered by the ggplot2 package when plotting left singular functions that correspond with variables observed over two-dimensional domains. A user may now specify a character vector that contains the time when each observation is made when building an fts object and we improved various plot fonts for readability. Finally, we include many other small updates that further improve plotting quality, code readability, documentation improvements, and other details that add to the professionalism of the package.


Haghbin, H., Morteza Najibi, S., Mahmoudvand, R., Trinka, J., and Maadooliat, M. (2021). Functional singular spectrum analysis. Stat. e330 STAT-20-0240.R1.

Trinka J., Haghbin H., Maadooliat M. (Accepted) Multivariate Functional Singular Spectrum Analysis: A Nonparametric Approach for Analyzing Multivariate Functional Time Series. In: Bekker A., Ferreira, J., Arashi M., Chen D. (eds) Innovations in Multivariate Statistical Modeling: Navigating Theoretical and Multidisciplinary Domains. Emerging Topics in Statistics and Biostatistics. Springer, Cham.

Trinka J. (2021) Functional Singular Spectrum Analysis: Nonparametric Decomposition and Forecasting Approaches for Functional Time Series [Doctoral dissertation, Marquette University]. ProQuest Dissertations Publishing.

Trinka, J., Haghbin, H., and Maadooliat, M. (2021). Functional time series forecasting: Functional singular spectrum analysis approaches. Version 4 retrieved from https://arxiv.org/abs/2011. 13077.

See Also

fssa, freconstruct, fforecast fwcor, wplot, fts, plot.fts, plot.fssa, cor.fts, launchApp

Rfssa documentation built on Sept. 9, 2022, 1:07 a.m.