- class statsmodels.sandbox.tsa.fftarma.ArmaFft(ar, ma, n)¶
fft tools for arma processes
This class contains several methods that are providing the same or similar returns to try out and test different implementations.
TODO: check whether we do not want to fix maxlags, and create new instance if maxlag changes. usage for different lengths of timeseries ? or fix frequency and length for fft
check default frequencies w, terminology norw n_or_w
some ffts are currently done without padding with zeros
returns for spectral density methods needs checking, is it always the power spectrum hw*hw.conj()
normalization of the power spectrum, spectral density: not checked yet, for example no variance of underlying process is used
Theoretical autocorrelation function of an ARMA process.
acf2spdfreq(acovf[, nfreq, w])
not really a method just for comparison, not efficient for large n or long acf
Theoretical autocovariances of stationary ARMA processes
A finite-lag AR approximation of an ARMA process.
A finite-lag approximate MA representation of an ARMA process.
Fourier transform of AR polynomial, zero-padded at end to n
Fourier transform of ARMA polynomial, zero-padded at end to n
Fourier transform of MA polynomial, zero-padded at end to n
filter a timeseries with the ARMA filter
filter a time series using fftconvolve3 with ARMA filter
from_coeffs([arcoefs, macoefs, nobs])
Create ArmaProcess from an ARMA representation.
Create an ArmaProcess from the results of an ARIMA estimation.
generate_sample([nsample, scale, distrvs, ...])
Simulate data from an ARMA.
Compute the impulse response function (MA representation) for ARMA process.
Make MA polynomial invertible by inverting roots inside unit circle.
autocovariance from spectral density
Theoretical partial autocorrelation function of an ARMA process.
construct AR and MA polynomials that are zero-padded to a common length
padarr(arr, maxlag[, atend])
pad 1d array with zeros at end to have length maxlag function that is a method, no self used
Periodogram for ARMA process given by lag-polynomials ar and ma.
plot4([fig, nobs, nacf, nfreq])
raw spectral density, returns Fourier transform
power spectral density using padding to length n done by fft
ma only, need division for ar, use LagPolynomial
spectral density from MA polynomial representation for ARMA process
spectral density for frequency using polynomial roots
power spectral density using fftshift
Roots of autoregressive lag-polynomial
Arma process is invertible if MA roots are outside unit circle.
Arma process is stationary if AR roots are outside unit circle.
Roots of moving average lag-polynomial