The SPECTRA Procedure

Kernels

Kernels are used to smooth the periodogram by using a weighted moving average of nearby points. A smoothed periodogram is defined by the following equation.

\begin{eqnarray*} \hat{J}_{i}(\textrm{l}(q)) = \sum _{{\tau } = -\textrm{l}(q)}^{\textrm{l}(q)} {\textrm{w}\left( \frac{{\tau }}{\textrm{l}(q)} \right) \tilde{J}_{i+{\tau }}} \nonumber \end{eqnarray*}

where ${\mr{w}(x)}$ is the kernel or weight function. At the endpoints, the moving average is computed cyclically; that is,

\begin{eqnarray*} \tilde{J}_{i+{\tau }} = \begin{cases} J_{i+{\tau }} & 0 <= i+\tau <= q \\ J_{-(i+{\tau })} & i+\tau < 0 \\ J_{q-(i+{\tau })} & i+\tau > q \end{cases}\nonumber \end{eqnarray*}

The SPECTRA procedure supports the following kernels. They are listed with their default bandwidth functions.

Bartlett: KERNEL BART

\begin{eqnarray*} \textrm{w}(x) & =& \begin{cases} 1-{|x|} & {|x|}{\le }1 \\ 0 & \mr{otherwise} \end{cases}\\ \mr{l}(q) & =& \frac{1}{2} q^{1 / 3} \nonumber \end{eqnarray*}

Parzen: KERNEL PARZEN

\begin{eqnarray*} \textrm{w}(x) & =& \begin{cases} 1-6{|x|}^{2} + 6{|x|}^{3} & {0{\le }{|x|}{\le }\frac{1}{2}} \\ 2(1-{|x|})^{3} & {\frac{1}{2}{\le }{|x|}{\le }1} \\ 0 & \mr{otherwise} \end{cases}\\ \mr{l}(q) & =& q^{1 / 5} \nonumber \end{eqnarray*}

Quadratic spectral: KERNEL QS

\begin{eqnarray*} \mr{w}(x) & =& \frac{25}{12{\pi }^{2} x^{2}} \left( \frac{{sin}(6{\pi }x/5)}{6{\pi }x/5} - {cos}(6{\pi }x/5) \right) \\ \textrm{l}(q) & =& \frac{1}{2} q^{1 / 5} \nonumber \end{eqnarray*}

Tukey-Hanning: KERNEL TUKEY

\begin{eqnarray*} \textrm{w}(x) & =& \begin{cases} (1+{cos}( {\pi } x))/2 & {|x|}{\le }1 \\ 0 & \mr{otherwise} \end{cases}\\ \mr{l}(q) & =& \frac{2}{3} q^{1 / 5} \nonumber \end{eqnarray*}

Truncated: KERNEL TRUNCAT

\begin{eqnarray*} \textrm{w}(x) & =& \begin{cases} 1 & {|x|}{\le }1 \\ 0 & \mr{otherwise} \end{cases}\\ \mr{l}(q) & =& \frac{1}{4} q^{1 / 5} \nonumber \end{eqnarray*}

A summary of the default values of the bandwidth parameters, c and e, associated with the kernel smoothers in PROC SPECTRA are listed below in Table 33.2:

Table 33.2: Bandwidth Parameters

Kernel

c

e

Bartlett

$1/2$

$1/3$

Parzen

1

$1/5$

quadratic

$1/2$

$1/5$

Tukey-Hanning

$2/3$

$1/5$

truncated

$1/4$

$1/5$


Figure 33.1: Kernels for Smoothing

Kernels for Smoothing


See Andrews (1991) for details about the properties of these kernels.