plotReport function

Plot results of a evaluate.pm analysis

Create a simple plot of a time series and the results of a evaluate.pm analysis (including the periodogram and the fitted OUSS power spectrum).

plotReport(name="", times=NULL, signal=NULL, 
           report=NULL, plotFile=NULL, dataFile=NULL, sep=" ")

Arguments

• name: Character. A short name for the time series to be used for the plots (e.g. 'long-term study' or 'hare population').
• times: Numeric vector. The time points of the time series used for the analysis. Set to NULL to not plot the time series.
• signal: Numeric vector. The time series values (signal) used for the analysis. Set to NULL to not plot the time series.
• report: The value returned by evaluate.pm. This is a list with the periodogram analysis results. Set to NULL to not plot any analysis results.
• plotFile: An optional path to a PDF file to be generated with the plot.
• dataFile: An optional path to a data file for storing the time series and the results of the analysis.
• sep: Separator to be used for the data file. Only relevant if dataFile is not NULL.

Returns

No return value.

References

Louca, S., Doebeli, M. (2015) Detecting cyclicity in ecological time series, Ecology 96: 1724--1732

Author(s)

Stilianos Louca

See Also

evaluate.pm

Examples

# generate cyclic time series by adding a periodic signal to an OUSS process
times  = seq(0,20,0.25);
signal = 0.6 * cos(2*pi*times) + generate_ouss(times, mu=0, sigma=1, lambda=1, epsilon=0.5);

# find periodogram peak and estimate statistical significance
report = evaluate.pm( times=times, 
                      signal=signal, 
                      minPeakFreq=0.3, 
                      minFitFreq=0.3,
                      startRadius=2);

# plot overview of periodogram peak analysis
plotReport(sprintf("Example"), times=times, signal=signal, report=report);