PCDPCA
Switch branches/tags
Nothing to show
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Permalink
Failed to load latest commit information.
R
demo
man
tests
.Rbuildignore
.gitignore
DESCRIPTION
NAMESPACE
README.md
cran-comments.md
pcdpca.Rproj

README.md

pcdpca

Implementation of "Dynamic principal components of periodically correlated functional time series".

Two examples in demo directory:

  • pm10 data from Graz (comparizon with DFPCA paper)
  • simplation with parametrized periodicity

Installation

library("devtools")
install_github("kidzik/pcdpca")

Running a demo

library("pcdpca")
demo("simulation")
demo("pcdpca.pm10")

Usage

Let X be a multivariate time series, a matrix with n observations and d covariates, periodic with period = 2. Then

FF = pcdpca(X, period=2)  # finds the optimal filter
Yhat = pcdpca.scores(X, FF)  # applies the filter
Yhat[,-1] = 0 # forces the use of only one component
Xhat = pcdpca.inverse(Yhat, FF)  # deconvolution
cat(sum((X-Xhat)^2) / sum(X^2)) # variance explained