Good day everyone, I am using the data below to fit Intervention Time Series model for two policies introduced in 2002 (54th data point) and 2003 (55th data point) respectively. Please can anyone give me a complete R code for modeling the two step functions? I have already modeled the pre-intervention period (1948 - 2001) but I need codes for the full intervention model that estimates the impact parameters(omega) and the decay parameters(delta) as well. Thank you all for your assistance. data(1948 - 2011) 1 207559 2 278372 3 247834 4 262223 5 210663 6 246982 7 211016 8 223317 9 232448 10 259788 11 209765 12 259572 13 322223 14 439159 15 415186 16 428018 17 427782 18 591031 19 415753 20 381353 21 430665 22 355588 23 417457 24 427894 25 469864 26 421843 27 354634 28 378759 29 400321 30 324111 31 271339 32 265076 33 296419 34 257974 35 224882 36 178626 37 158956 38 174809 39 219044 40 227765 41 188177 42 300101 43 296051 44 293352 45 242817 46 312123 47 254653 48 309454 49 403872 50 322488 51 409383 52 397675 53 436947 54 389772 55 340562 56 496846 57 736975 58 599318 59 740458 60 614532 61 680781 62 710642 63 632037 64 1024553 [[alternative HTML version deleted]]