logo

Step Function and Pulse Function 📂Statistical Analysis

Step Function and Pulse Function

Definition 1

  1. The function defined as follows St(T)S_{t}^{(T)} is called a step function. St(T):={1,tT0,otherwise S_{t}^{(T)} := \begin{cases} 1 & , t \le T \\ 0 & , \text{otherwise} \end{cases} 20190818_180235.png
  2. The function defined as follows Pt(T)P_{t}^{(T)} is called a pulse function. Pt(T):=St(T)=St(T)St1(T) \begin{align*} P_{t}^{(T)} :=& \nabla S_{t}^{(T)} \\ =& S_{t}^{(T)} - S_{t-1}^{(T)} \end{align*} 20190818_180247.png

Description

Step functions and pulse functions are useful for representing equations used in intervention analysis, and their properties per se don’t have significant meanings. The step function is named for its stair-like appearance on the graph, and the pulse function represents a short moment of impact graphically. [ NOTE: It’s interesting that such shapes and concepts also appear in mathematical physics. ]

In the form of Yt=mt+NtY_{t} = m_{t} + N_{t} in intervention analysis, the term intervening mtm_{t} can be represented by these functions. If the analysis drastically changes at some point TT, a step function can be used, or a pulse function can be used to handle just one exception. For example, it can be as follows: mt=ωSt(T) m_{t} = \omega S_{t}^{(T)}

mt=ωPt(T) m_{t} = \omega P_{t}^{(T)} Here, ω\omega is a coefficient. Since the function values of step and pulse functions are 00 or 11, such correction is needed. mtm_{t} can be used more freely than you might think. For instance, mtm_{t} itself can be assumed to follow some ARIMA model. The following resembles the form of ARMA(1,1)ARMA(1,1): mt=δmt1+ωPt1(T) m_{t} = \delta m_{t-1} + \omega P_{t-1}^{(T)} Similarly, δ\delta is a coefficient. This representation is interesting because by using the backshift BB, the following equation manipulation can be done: mt=δmt1+ωPt1(T)=δBmt+ωBPt(T) \begin{align*} m_{t} =& \delta m_{t-1} + \omega P_{t-1}^{(T)} \\ =& \delta B m_{t} + \omega B P_{t}^{(T)} \end{align*} If we move δBmt\delta B m_{t} to the left side, mtδBmt=ωBPt(T) m_{t} - \delta B m_{t} = \omega B P_{t}^{(T)} and divide both sides by (1δB)(1-\delta B), mt=ωB1δBPt(T) m_{t} = {{\omega B} \over {1-\delta B}} P_{t}^{(T)} that is, if mtm_{t} is not extremely complicated, it can be represented in a clean form like mt=ω(B)δ(B)Pt(T)m_{t} = {{\omega ( B) } \over { \delta (B) }} P_{t}^{(T)}. In this manner, we can also obtain the useful relationship, St(T)=11BPt(T) S_{t}^{(T)} = {{1} \over {1 - B}} P_{t}^{(T)} allowing us to develop equations more freely. Although in actual analysis there might not be a real occasion to use this, at least understanding that the mtm_{t} in intervention analysis is derived in this fashion and form is necessary.


  1. Cryer. (2008). Time Series Analysis: With Applications in R(2nd Edition): p250~251. ↩︎