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]]

