# Spec_compare_localize_freq_curvelength: Compare the spectral density operator of two Functional Time... In ftsspec: Spectral Density Estimation and Comparison for Functional Time Series

## Description

Compare the spectral density operator of two Functional Time Series and localize frequencies at which they differ, and (spatial) regions where they differ

## Usage

 1 2 Spec_compare_localize_freq_curvelength(X, Y, B.T = (dim(X)[1])^(-1/5), W, alpha = 0.05, accept = 0, reject = 1, verbose = 0, demean = FALSE) 

## Arguments

 X The T \times nbasis matrices of containing the coordinates, expressed in some functional basis, of the two FTS that to be compared. expressed in a basis. Y The T \times nbasis matrices of containing the coordinates, expressed in some functional basis, of the two FTS that to be compared. expressed in a basis. B.T The bandwidth of frequencies over which the periodogram operator is smoothed. If B.T=0, the periodogram operator is returned. W The weight function used to smooth the periodogram operator. Set by default to be the Epanechnikov kernel alpha level for the test accept,reject values for accepted, rejected regions verbose A variable to show the progress of the computations. By default, verbose=0. demean A logical variable to choose if the FTS is centered before computing its spectral density operator.

## Examples

  1 2 3 4 5 6 7 8 9 10 ma.scale2=ma.scale1=c(-1.4,2.3,-2) ma.scale2[3] = ma.scale1[3]+.4 a1=Generate_filterMA(10, 10, MA.len=3, ma.scale=ma.scale1) a2=Generate_filterMA(10, 10, MA.len=3, ma.scale=ma.scale2) X=Simulate_new_MA(a1, T.len=2^9, noise.type='wiener') Y=Simulate_new_MA(a2, T.len=2^9, noise.type='wiener') ans0=Spec_compare_localize_freq_curvelength(X, Y, W=Epanechnikov_kernel, alpha=.01, demean=TRUE) print(ans0) plot(ans0) rm(ma.scale1, ma.scale2, a1, a2, X, Y, ans0) 

ftsspec documentation built on May 1, 2019, 7:04 p.m.