README.md

pcdpca

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

Two examples in demo directory:

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


Try the pcdpca package in your browser

Any scripts or data that you put into this service are public.

pcdpca documentation built on Sept. 3, 2017, 9:02 a.m.