naginterfaces.library.tsa.uni_spectrum_lag¶
- naginterfaces.library.tsa.uni_spectrum_lag(nx, mtx, px, iw, mw, ic, c, kc, l, lg, xg)[source]¶
uni_spectrum_lag
calculates the smoothed sample spectrum of a univariate time series using one of four lag windows – rectangular, Bartlett, Tukey or Parzen window.For full information please refer to the NAG Library document for g13ca
https://support.nag.com/numeric/nl/nagdoc_30.2/flhtml/g13/g13caf.html
- Parameters
- nxint
, the length of the time series.
- mtxint
If covariances are to be calculated by the function (), must specify whether the data are to be initially mean or trend corrected.
For no correction.
For mean correction.
For trend correction.
- pxfloat
If covariances are to be calculated by the function (), must specify the proportion of the data (totalled over both ends) to be initially tapered by the split cosine bell taper.
If covariances are supplied , must specify the proportion of data tapered before the supplied covariances were calculated and after any mean or trend correction. is required for the calculation of output statistics.
A value of implies no tapering.
- iwint
The choice of lag window.
Rectangular.
Bartlett.
Tukey.
Parzen.
- mwint
, the ‘cut-off’ point of the lag window. Windowed covariances at lag or greater are zero.
- icint
Indicates whether covariances are to be calculated in the function or supplied in the call to the function.
Covariances are to be calculated.
Covariances are to be supplied.
- cfloat, array-like, shape
If , must contain the covariances for lags from to , otherwise need not be set.
- kcint
If , must specify the order of the fast Fourier transform (FFT) used to calculate the covariances.
If , that is covariances are supplied, is not used.
- lint
, the frequency division of the spectral estimates as . Therefore, it is also the order of the FFT used to construct the sample spectrum from the covariances.
- lgint
Indicates whether unlogged or logged spectral estimates and confidence limits are required.
Unlogged.
Logged.
- xgfloat, array-like, shape
If the covariances are to be calculated, then must contain the data points. If covariances are supplied, may contain any values.
- Returns
- cfloat, ndarray, shape
If , will contain the calculated covariances.
If , the contents of will be unchanged.
- xgfloat, ndarray, shape
Contains the spectral estimates, , for in to respectively (logged if ). The elements , for contain .
- ngint
The number of spectral estimates, , in .
- statsfloat, ndarray, shape
Four associated statistics. These are the degrees of freedom in , the lower and upper confidence limit factors in and respectively (logged if ), and the bandwidth in .
- Raises
- NagValueError
- (errno )
On entry, and .
Constraint: .
- (errno )
On entry, and .
Constraint: .
- (errno )
On entry, and .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: , , or .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, and .
Constraint: if , .
- (errno )
On entry, , and .
Constraint: if , .
- (errno )
On entry, and : .
- (errno )
On entry, and : .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, , and .
Constraint: if , .
- (errno )
On entry, and .
Constraint: .
- Warns
- NagAlgorithmicWarning
- (errno )
One or more spectral estimates are negative.
Unlogged spectral estimates are returned in , and the degrees of freedom, unloged confidence limit factors and bandwidth in .
- (errno )
The calculation of confidence limit factors has failed.
Spectral estimates (logged if requested) are returned in , and degrees of freedom and bandwidth in .
- Notes
The smoothed sample spectrum is defined as
where is the window width, and is calculated for frequency values
where denotes the integer part.
The autocovariances may be supplied by you, or constructed from a time series , as
the fast Fourier transform (FFT) being used to carry out the convolution in this formula.
The time series may be mean or trend corrected (by classical least squares), and tapered before calculation of the covariances, the tapering factors being those of the split cosine bell:
where and is the tapering proportion.
The smoothing window is defined by
which for the various windows is defined over by
rectangular:
Bartlett:
Tukey:
Parzen:
The sampling distribution of is approximately that of a scaled variate, whose degrees of freedom is provided by the function, together with multiplying limits , from which approximate confidence intervals for the true spectrum may be constructed as . Alternatively, log may be returned, with additive limits.
The bandwidth of the corresponding smoothing window in the frequency domain is also provided. Spectrum estimates separated by (angular) frequencies much greater than may be assumed to be independent.
- References
Bloomfield, P, 1976, Fourier Analysis of Time Series: An Introduction, Wiley
Jenkins, G M and Watts, D G, 1968, Spectral Analysis and its Applications, Holden–Day