Plot results of a evaluate.pm
analysis
evaluate.pm
analysisCreate 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=" ")
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
.No return value.
Louca, S., Doebeli, M. (2015) Detecting cyclicity in ecological time series, Ecology 96: 1724--1732
Stilianos Louca
evaluate.pm
# 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);
Useful links