REPPERIODOGRAM procedure

Gives periodogram-based analyses for replicated time series (R.P. Littlejohn).


Options

PRINT = string
What to print (pair, randomization, glm); default * i.e. none

PLOT = string
What graphs to plot (group, mean, logmean, cumulative, cv, pair); default mean, logm

TITLE = text
Title for each page of graphs

REPRESENTATION = string
Form of data in SERIES (timeseries, meanperiodogram); default time

LENGTH = scalar or variate
Scalar specifying that the first N units of the series are to be used, or a variate specifying the first and last units of the series to be used

SEED = scalar
Seed for randomization; default 0

NRANDOMIZATIONS = scalar
Number of randomizations; default 99

TREATMENTS = factor
Contains ordered classification of SERIES

PAIR = variates
Treatment pair levels for pairwise comparisons

COLOUR = scalars
One colour for each level of TREATMENTS; default 1...nl where nl = number of levels of TREATMENTS

MEANPERIODOGRAM = pointer
Saves mean periodograms according if REPRESENTATION=timeseries

REPLICATION = scalar or variate
Inputs or saves number of replicate series if REPRESENTATION=timeseries; scalar can be used for equal replication


Parameter

SERIES = variates
Specify the time series to be analysed


Description

REPPERIODOGRAM gives periodogram-based analyses of replicated time series. The data are supplied in a list of variates using the SERIES parameter, either as the original time series (option REPRESENTATION set to timeseries) with the level for each series given by the factor specified by the TREATMENTS option, or as the mean periodograms for each treatment level (option REPRESENTATION set to meanperiodogram), with levels and labels optionally given by the TREATMENTS factor and the multiplicity of each treatment defined by the REPLICATION option. In the former case the LENGTH option can specify that only part of each series is to be used, using either a scalar N to indicate that the first N values are to be used, or a variate of length two, holding the values of the first and last units of the required subseries. This may be used to eliminate missing values, which are otherwise not permitted. Further, when REPRESENTATION=timeseries, periodogram means and the replication variate can be saved using the MEANPERIODOGRAM and REPLICATION options, respectively.

   Graphical output is controlled by the PLOT option. For the group (REPRESENTATION=timeseries only), mean, logmean and cumulative periodogram and cv graphs, the COLOUR option can be used to code for treatments. The cv plot (REPRESENTATION=timeseries only) gives a scatterplot of coefficients of variation for each treatment group at each frequency, together with lines for the means of these cvs at each frequency for those treatments with replication greater than one, and cv=1, the theoretical value if there is no subject-specific variation. For these graphs a title can be supplied using the TITLE option. Graphs are also given for the differences between pairs of log periodograms as defined by PAIR (see below), with 95% confidence intervals on the sample and null (equal periodograms) distribution.

   Output of various test statistics for pairwise comparison of treatment levels described by Diggle (1990) and Diggle & Fisher (1991) is controlled by the PRINT and PAIR options. PAIR is a list of 2-unit variates representing treatment levels, e.g.

PAIR=!(1,2),!(3,4)

gives tests comparing treatment levels 1 and 2, followed by tests for levels 3 and 4. With PRINT=pair, the maximum absolute value and range of the difference of log periodograms give (weak) tests against the null hypotheses of equal and proportional spectra, respectively. With PRINT=random, a randomization test is given for the equality of cumulative spectra, which is insensitive to the alternative of proportional spectra. The seed for the ramdomizations can be set using the SEED option, and the number of randomizations is specified by NRANDOMIZATIONS (default 99). This is available only if the treatments in the pair have equal replication.

   When PRINT=glm, a generalized linear model is fitted to the mean periodograms for all treatments, adjusting for frequency, and testing for differences with treatment in constant (proportional spectra), linear (power shift) and quadratic (power spread) contrasts with frequency (Diggle 1990). Results are presented in the accumulated analysis of deviance table and tables of parameter estimates, within which the Intercept-Difference, Slope-Difference and Curve-Difference estimates relate to the above hyptheses.

 

Options: PRINT, PLOT, REPRESENTATION, LENGTH, TREATMENTS, PAIR, SEED, NRANDOMIZATIONS, COLOUR, TITLE, MEANPERIODOGRAM, REPLICATION.

Parameter: SERIES.


Method

The series are mean-corrected, but not trend corrected, before transformation, and are not smoothed. Critical values for the Range test are obtained from tables in Potscher & Reschenhofer (1988) and Coates & Diggle (1986). Random numbers are generated using URAND. The analysis for PRINT=glm is obtained from fitting a generalized linear model with DISTRIBUTION=gamma, LINK=log and DISPERSION=1/nr, where nr is number of replicates of the treatments.


Action with RESTRICT

The SERIES may not be restricted; restriction of the input series to a contiguous set of units may be achieved by use of the LENGTH parameter.


References

Coates, D.S. & Diggle, P.J. (1986) Tests for comparing two estimated spectral densities. Journal of Time Series Analysis, 7, 7-20.

Diggle, P.J. (1990). Time Series: A Biostatistical Introduction. Oxford, Clarendon Press.

Diggle, P.J. & Fisher, N.I. (1991). Nonparametric comparisons of cumulative periodograms. Applied Statistics, 40, 423-434.

Potscher, B.M. & Reschenhofer, E. (1988). Discriminating between two spectral densities in case of replicated observations. Journal of Time Series Analysis, 9, 221-224.