- •Distribution Overview
- •Discrete Distributions
- •Continuous Distributions
- •Probability Theory
- •Random Variables
- •Transformations
- •Expectation
- •Variance
- •Inequalities
- •Distribution Relationships
- •Probability and Moment Generating Functions
- •Multivariate Distributions
- •Standard Bivariate Normal
- •Bivariate Normal
- •Multivariate Normal
- •Convergence
- •Statistical Inference
- •Point Estimation
- •Empirical distribution
- •Statistical Functionals
- •Parametric Inference
- •Method of Moments
- •Maximum Likelihood
- •Delta Method
- •Multiparameter Models
- •Multiparameter delta method
- •Parametric Bootstrap
- •Hypothesis Testing
- •Bayesian Inference
- •Credible Intervals
- •Function of parameters
- •Priors
- •Conjugate Priors
- •Bayesian Testing
- •Exponential Family
- •Sampling Methods
- •The Bootstrap
- •Rejection Sampling
- •Importance Sampling
- •Decision Theory
- •Risk
- •Admissibility
- •Bayes Rule
- •Minimax Rules
- •Linear Regression
- •Simple Linear Regression
- •Prediction
- •Multiple Regression
- •Model Selection
- •Non-parametric Function Estimation
- •Density Estimation
- •Histograms
- •Kernel Density Estimator (KDE)
- •Smoothing Using Orthogonal Functions
- •Stochastic Processes
- •Markov Chains
- •Poisson Processes
- •Time Series
- •Stationary Time Series
- •Estimation of Correlation
- •Detrending
- •ARIMA models
- •Causality and Invertibility
- •Spectral Analysis
- •Math
- •Gamma Function
- •Beta Function
- •Series
- •Combinatorics
Moving average polynomial |
|
|
|
|
|
|
|
|
(z) = 1 + 1z + + qzq |
|
z 2 C ^ q 6= 0 |
||||||
Moving average operator |
|
|
|
|
|
|
|
|
(B) = 1 + 1B + + pBp |
||||||||
MA (q) (moving average model order q) |
|
|
|
|
||||
xt = wt + 1wt 1 + + qwt q () xt = (B)wt |
||||||||
|
|
|
q |
|
|
|
|
|
|
Xj |
|
|
|
|
|||
E [xt] = jE [wt j] = 0 |
||||||||
|
|
=0 |
P |
|
|
|||
t+h |
t |
|
(0 |
q h |
h > q |
|||
|
|
|
2 |
|
|
j j+h 0 h q |
||
(h) = Cov [x |
; x ] = |
|
w |
j=0 |
||||
MA (1) |
|
|
|
|
|
|
|
|
|
xt = wt + wt 1 |
h = 1 |
||||||
(h) = 8 w2 |
|
|
|
|||||
|
> |
(1 + 2) 2 |
h = 0 |
|||||
|
|
|
|
|
|
w |
|
|
|
<0 |
|
|
|
|
h > 1 |
||
|
> |
|
|
|
|
|
|
|
|
: |
|
|
|
|
h = 1 |
||
|
|
|
(1+ 2) |
|||||
(h) = (0 |
|
|
h > 1 |
|||||
ARMA (p; q) |
|
|
|
|
|
|
|
|
xt = 1xt 1 + + pxt p + wt + 1wt 1 + + qwt q |
||||||||
|
(B)xt = (B)wt |
|||||||
Partial autocorrelation function (PACF) |
|
|
|
|
||||
xih 1 , regression of xi on fxh 1; xh 2; : : : ; x1g |
||||||||
hh = corr(xh xhh 1; x0 x0h 1) |
h 2 |
|
||||||
E.g., 11 = corr(x1; x0) = (1) |
|
|
|
|
|
|
ARIMA (p; d; q)
rdxt = (1 B)dxt is ARMA (p; q)(B)(1 B)dxt = (B)wt
Exponentially Weighted Moving Average (EWMA)
xt = xt 1 + wt wt 1
1
X
xt = (1 ) j 1xt j + wt when j j < 1
j=1
x~n+1 = (1 )xn + x~n
Seasonal ARIMA
Denoted by ARIMA (p; d; q) (P; D; Q)s
P (Bs) (B)rDs rdxt = + Q(Bs) (B)wt
21.4.1Causality and Invertibility
ARMA (p; q) is causal (future-independent) () 9f |
|
1 |
|
||||||||
jg : Pj=0 j < 1 such that |
|||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
Xj |
wt j = (B)wt |
|
|
|
||||
|
|
|
xt = |
|
|
|
|||||
|
|
|
=0 |
|
|
|
|
|
|
|
|
|
|
|
() 9f jg |
1 |
|
|
|
|
|
|
|
ARMA (p; q) is invertible |
: Pj=0 j < 1 such that |
||||||||||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
Xj |
|
|
|
||||
|
|
|
(B)xt = Xt j = wt |
|
|
|
|||||
|
|
|
|
=0 |
|
|
|
|
|
|
|
Properties |
|
|
|
|
|
|
|
|
|
|
|
ARMA (p; q) causal |
() roots of (z) lie outside the unit circle |
||||||||||
|
|
|
1 |
|
|
(z) |
|
|
|
||
|
|
|
Xj |
|
|
|
jzj 1 |
||||
|
|
|
|
|
|||||||
|
|
|
(z) = |
jzj = (z) |
|||||||
|
|
|
=0 |
|
|
|
|
|
|
|
|
ARMA (p; q) invertible () roots of (z) lie outside the unit circle |
|||||||||||
|
|
|
1 |
(z) |
|
|
|
||||
|
|
|
Xj |
|
|
jzj 1 |
|
|
|||
|
|
|
|
|
|
||||||
|
|
|
(z) = |
jzj = (z) |
|
|
|||||
|
|
|
=0 |
|
|
|
|
|
|
|
|
Behavior of the ACF and PACF for causal and invertible ARMA models |
|||||||||||
|
|
|
|
|
|
|
|||||
|
|
AR (p) |
MA (q) |
|
ARMA (p; q) |
|
|||||
|
ACF |
tails o |
cuts o after lag q |
tails o |
|
||||||
|
PACF |
cuts o after lag p |
tails o q |
|
tails o |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
21.5Spectral Analysis
Periodic process
xt = A cos(2 !t + )
=U1 cos(2 !t) + U2 sin(2 !t)
Frequency index ! (cycles per unit time), period 1=!
25