s2dverification (version 2.4.0)

Alpha: Estimates AutoCorrelation At Lag 1 following Guemas et al, BAMS, 2013b

Description

This function, relying on the FitAcfCoef() function, estimates the autocorrelation at lag 1 of the xdata array following the method described in Guemas V., Auger L., Doblas-Reyes F., JAMC, 2013. After applying a linear detrending and/or a filtering of any frequency peak if requested, the sample autocorrelation is estimated. Then the theoretical autocorrelation of an AR1 is fitted to the sample autocorrelation using the Cardano's formula (see FitAcfCoef()) to obtain the autocorrelation at lag 1. This method assumes xdata is an AR1 process.

Usage

Alpha(xdata, detrend = FALSE, filter = FALSE)

Value

Autocorrelation at lag 1

Examples

# Load sample data as in Load() example:
example(Load)
alpha <- Alpha(sampleData$mod[1, 1, , 1])
print(alpha)