J. Indones. Math. Soc. Vol. No. , pp. 1Ae11. Double Intervention Analysis on The Arima Model of Covid-19 Cases in Bali Nurfitri ImroAoah1 and NurAoainul Miftahul Huda2O Department of Statistics. Universitas Tanjungpura. Indonesia. nurfitriimroah@math. Department of Mathematics. Universitas Tanjungpura. Indonesia. nurAoainul@fmipa. Abstract. The time series process is not only influenced by previous observations, but some phenomena result in drastic changes to observations in the time series process so that there is a change in the average or only a temporary change in For example, there is a policy from the government towards handling a case. This is referred to as an intervention. Therefore, it is necessary to do time series modeling with intervention factors. One form of intervention in the current pandemic era is a policy issued by the government. In this study, the time series model used is ARIMA. This study aimed to analyze the effect of an intervention on the ARIMA model on Covid-19 cases in Bali. This study uses data on the number of new Covid-19 cases in Bali from 24 April 2020 to 31 May 2021. There are two interventions used in this study, namely restrictions on activities for the Panca Yadnya ceremony and crowds in Bali and restrictions on traveling outside the area and/or going home and/or leaving for employees of the State Civil Apparatus during the Covid-19 pandemic. The results of this study show that two policies issued by the Bali provincial government can handle the addition of new cases of Covid-19. It can be seen from the decline in the number of new Covid-19 cases in Bali until the end of May 2021. Key words and Phrases: Intervention, step function, pulse function. MSE, restriction INTRODUCTION Time series analysis is a method for modeling a time series process, namely a stochastic process, a series of random variables with a time index . Equations are generated as a result of time series analysis that reflect the progression of data across time. Therefore, the fundamental goal of time series analysis is to forecast Corresponding author 2020 Mathematics Subject Classification: Received: 05-01-2024. Accepted: 13-08-2024. ImroAoah and Huda how a process will evolve throughout subsequent periods . Time series analysis can be used to simulate a variety of topics that are encountered in everyday life . Some examples include anticipating the movement of stock prices and the number of tourists who will visit in the next months. Therefore, time series analysis is a particularly applicable science in the sense that it helps solve a great deal of everyday issues (. , . The time series analysis only applies excellent form patterns to some observations . An illustration of this would be a decision made by the government regarding how to handle a particular instance, which leads to an abrupt rise in the The spread of the Covid-19 virus is tied to one of the measures that the government has in place . The government has created a policy in response to the Covid-19 outbreak that emphasizes reducing the risk of transmission within the population . In general, the governmentAos strategy for avoiding the spread of infectious diseases is broken up into three distinct categories: in and around the place of residence, when traveling, and when engaging in activities that take place outside the home. Researchers such as . , . , . , and . have successfully applied time series analysis to the COVID-19 data. Aside from that, there are also a large number of studies that integrate time series analysis with intervention One such study was undertaken by . , . , and . , who investigated the effects of the COVID-19 pandemic using intervention analysis. The data on the spread of Covid-19 in several areas, including in Bali, has increased again with a slowing recovery rate . The emergence of clusters of Covid-19 cases stems from community interactions, so Parisadha Hindu Dharma Bali Province and the Traditional Village Council The Province of Bali issued circulars Number 081/PHDI-Bali/IX/2020 and Number: 007/SE/MDA-Bali Province/IX/2020 concerning Restrictions on the Activities of Panca Yadnya Ceremonies and Crowds in Bali in the Gering Agung situation. The ceremony is attempted to be carried out with provisions, a maximum of 1 . day unless other provisions require more than 1 . day while still implementing strict health In addition, the Bali provincial government has also issued a circular regarding Restrictions on Activities for Traveling Outside the Region and/or Homecoming and/or Leave for State Civil Apparatus Employees during the Corona Virus Disease 2019 (Covid-. Pandemic. This is done in the context of preventing and overcoming the Corona Virus Disease 2019 (Covid-. , which has the potential to increase due to people traveling during Christmas Day 2021 and New Year 2022. The two policies issued by the Bali provincial government can be said to be a form of intervention to suppress the addition of positive cases of Covid-19 in Bali. This study analyzes the time series model while considering intervention considerations. The autoregressive integrated moving average, sometimes known as the ARIMA model, is a popular choice for modeling time series. The number of new Covid-19 cases in Bali from 24 April 2020 to 31 May 2021 was used as the dependent variable in this study. Meanwhile, there are two intervention factors, namely: Restrictions on the Activities of Panca Yadnya Ceremonies and Crowds in Bali in the Gering Agung situation . tated as the first interventio. and Restrictions on Double Intervention Analysis on The ARIMA Model of Covid-19 Cases in Bali Activities for Traveling Outside the Region and/or Homecoming and/or Leave for State Civil Apparatus Employees during the Corona Virus Disease 2019 (Covid-. Pandemic . tated as the second interventio. Because there are two intervention factors, this study analyzes the time series model with double intervention. The first part of this study describes the background. The second and third sections explain an overview of the ARIMA time series model and the ARIMA model with double intervention factors. The data analysis and conclusions can be seen in the fourth and fifth sections. ARIMA TIME SERIES MODEL Assuming that the sequence of random variables Zt adheres to the ARIMA process, the definition of Zt can be expressed as . ip (B). Oe B)d Zt = 0 q (B)at . where ip (B) = . Oei1 B OeA A AOeip B ), q (B) = . Oe1 B OeA A AOeq B ). B is backshift operator and at is error at time t. The model in Equation . is called the ARIMA model with the order . , d, . and is written as the ARIMA model . , d, . If p = 0, then ARIMA . , d, . model is also called IMA model with order . , . Likewise, if q = 0, then the ARIMA . , d, . model is called an ARI model with order . , . ARIMA model is a model used for data stationery with the following three main stages . Order identification, can be accomplished by utilizing Auto-Correlation Function (ACF) plots and Partial Auto-Correlation Function (PACF) plots. Parameter estimation, minimum-squares and maximum-likelihood estimation. Diagnostic test, independent and normal residual tests. After the three primary steps have been completed, the best model will have been obtained, and it will now serve as a point of reference for any subsequent forecasts, h . ARIMA MODEL WITH INTERVENTION FACTOR Step and pulse functions are the two varieties found in the time series model with intervention elements . Step function, which takes place at time T and always leaves the effect of the intervention after it, symbolizes the intervention at that moment. One definition of an intervention with a step function type is . (T ) for t < T, for t Ou T. ImroAoah and Huda . Pulse function, refers to an intervention that will only have an effect within a specific period. alternatively, the effect will end at t = T or some point after t = T . The following is a definition of the intervention type known as the pulse function . (T ) for t = T, for t = T. Both a step function and a pulse function are viable options for accurately representing models incorporating intervention factors . Because the findings of the step function type differentiation can be used to derive the pulse function type intervention, this is the case . (T ) (T ) (T ) Oe StOe1 = . Oe B)St (T ) = Pt Numerous other answers could be taken in response to either form of intervention. The following are some of the various intervention responses that frequently come up . The constant impact of the intervention was felt b the period after the (T ) OB b St (T ) OB b Pt . For the step function type, the intervention will constantly impact O. for the pulse function type, the intervention will only impact O at the time of the T event. The perceived impact of the intervention b period after the intervention, but the resulting response is gradual. For the step function, write: OB b (T ) S . Oe B) t Furthermore, for the pulse function, written: OB b (T ) Oe B) t with 0 O O 1. In this instance, an intervention based on a step function has an impact that gradually increases or decreases until it reaches a constant state. Meanwhile, the intervention based on the pulse function has a gradual impact, but it will run out of steam at a certain point in the future. The following is a textual representation of the generic ARIMA model of the time series, which incorporates the intervention components . Zt = Nt ft Oc (B) b (T ) B It r (B) Double Intervention Analysis on The ARIMA Model of Covid-19 Cases in Bali and Nt is the time series model before adding the intervention factor. Oc (B) = (T ) 1 Oe O1 B Oe A A A Oe Oc B c , r (B) = 1 Oe 1 B Oe A A A Oe r B r . It is the intervention variable . tep function or pulse functio. , and B is the backshift operator . The time series modeling stage with the addition of intervention factors is presented in Figure Figure 1. Flowchart of time series modelling with intervention DATA ANALYSIS The number of daily new cases of COVID-19 that occurred in Bali between April 24, 2020, and May 31, 2021, was used as a data source for this investigation. The data came from the National Disaster Management Agency and included multiple observations of 403 cases. Meanwhile, there are two intervention factors, namely: Restrictions on the Activities of Panca Yadnya Ceremonies and Crowds in Bali in the Gering Agung situation . tated as the first interventio. and Restrictions on Activities for Traveling Outside the Region and/or Homecoming and/or Leave for State Civil Apparatus Employees during the Corona Virus Disease 2019 (Covid. Pandemic . tated as second interventio. The first step in modeling time series with intervention factors is to classify data from before or after the intervention. This studyAos pre-intervention data started from April 24, 2020, to September 13. The Figure 2 presents the number of new daily cases of Covid-19 in Bali in ImroAoah and Huda Figure 2. Time-series plot of the number of new daily cases of Covid-19 in Bali The intervention implemented by the Bali provincial government is expected to prevent more comprehensive transmission of Covid-19. However, the average number of new cases increased from 49. 02 to 162. ee Figure . This indicates that the intervention did not impact the number of new cases. Even after the first intervention was implemented, new Covid-19 cases in Bali reached 532 cases on January 26, 2021. This number was the highest number during the observation So the government implemented a second intervention to prevent more comprehensive transmission on April 8, 2021. As seen in Figure 2, the number of new cases decreased after this second interventionAos implementation. The first thing that needs to be done in the data modeling process is to apply the ARIMA . , d, . model to the data collected before the intervention. Figure 3 presents a plot of the pre-intervention data. Based on Figure 3(A), it can be seen that the pre-intervention data was not stationary. Therefore, differentiation is carried out to eliminate trend patterns in the data to fulfill a constant average. After doing the first differencing, it can be seen that the pre-intervention data is stationary because it has a constant mean and variance . ee Figure 3(B)). Furthermore, the ACF and PACF plots will be reviewed from the stationary differencing results to identify the order in the ARIMA model. Based on Figure 3(C), several candidate models were obtained, namely ARIMA . , 1, . ARIMA . , 1, . and ARIMA . , 1, . The parameter estimation. AIC values, and plot of residual diagnostic tests for each model are presented in Table 1. In determining the model for the pre-intervention, the model must fulfill the assumption of white noise. namely, the residuals are independent and normally Based on Table 1, it can be seen that the ARIMA . , 1, . model is the model with the smallest AIC value . ed tex. and fulfill the assumption of white noise . hird column in Table . Therefore, it can be concluded that the ARIMA . , 1, . model is the best model for pre-intervention data, namely: Yt = Oe0. 4559 et where Yt is the number of new daily cases at time t, and et is an error at time t with t = 1, 2, . , 143. Double Intervention Analysis on The ARIMA Model of Covid-19 Cases in Bali . Figure 3. Plot of the pre-intervention data. (A) actual observation. (B) first differencing. (C) ACF and PACF Plot from (B) In time series modeling with intervention factors on the data on the number of new Covid-19 cases in Bali, the intervention model used is as follows: f144 (O, ) = O0 Pt 1 Oe B t . f350 (O, ) = OBSt Equation . is a model for the first intervention on September 14, 2020. The model is a pulse function type to describe that the intervention applied has an effect that will disappear gradually. At the same time. Equation . states a model with a step function type for the second intervention that occurs on April 8, 2021. The second intervention model shows that the intervention applied can bring a permanent effect. The following is a complete model of time series modeling with intervention factors in this study. Zt = . Oe B)et O1 O0 Pt O2 BSt . Oe B) 1 Oe B t so that . Zt Oe ZtOe1 = (Oe Oe )etOe1 etOe2 (O0 O1 )Pt . (OeO0 Oe O0 Oe O1 )PtOe1 . O0 PtOe2 O2 StOe1 (OeO2 Oe O2 )StOe2 O2 StOe3 et . ImroAoah and Huda Table 1. The result of parameter estimation and AIC value from the ARIMA . , d, . on preintervention data Model Order . , d, . , 1, . , 1, . , 1, . y 10Oe2 Oe39. Oe45. Oe49. AIC Residuals plots After obtaining a time series model with intervention factors, the next step is to estimate the parameters contained in the model . n this case, the model in Equation . The blue text in Equation . is the parameters to be estimated. this study, the Least Square method is used to estimate the parameters. Equation . is the ARIMA . , 1, . model with the addition of two intervention factors so that the fitting results can be seen in Figure 4. Wt = Oe0. 575etOe1 Oe 0. 003etOe2 Oe 43. Oe 45. 272PtOe1 Oe 49. 285PtOe2 Oe 45. 876StOe1 Oe 28. 392StOe2 68. 029StOe3 . Figure 4 presents the plot of fitted ARIMA . , 1, . model with two intervention factors and the plot of actual observations. Based on the Figure 4, it can be seen that the model used is quite good in modeling the data on the number of new Covid-19 cases in Bali. This is also supported by the MSE value generated from the model, which is 1,970. This value is smaller than the MSE generated by the ARIMA . , 1, . model without adding an intervention factor of 2,041. CONCLUDING REMARKS Two policies issued by the Bali provincial government from April 24, 2020, to May 31, 2021, impacted the number of new cases of Covid-19 in Bali. These policies Double Intervention Analysis on The ARIMA Model of Covid-19 Cases in Bali Figure 4. Plot of fitted values vs observations are Restriction of Panca Yadnya Ceremony and Crowd Activities in Bali in the Gering Agung Covid-19 Situation and Restrictions on Traveling Outside the Region and/or Homecoming and/or Leave for State Civil Apparatus Employees During the Corona Virus Disease 2019 (Covid-. Pandemic. Prior to the issuance of the policy, the average number of new Covid-19 cases had increased. Even after the first policy was issued, the number of new cases had reached 532 per day. This value is the highest number during the observation period. Therefore, the Bali provincial government issued a second policy to suppress the number of new cases added. The data before the implementation of the policy . re-intervention dat. followed the ARIMA . , 1, . The addition of two intervention factors . olicies issued by the provincial government of Bal. to the ARIMA . , 1, . model shows better This can be seen from the resulting MSE value of 1,970. This value is smaller than the MSE of the ARIMA . , 1, . model without intervention factors, which is 2,041. In addition, it can also be seen from the decreasing number of positive cases of Covid-19 in Bali until the end of May 2021. The impact after the implementation of the policy occurs sustainably. So, it can be said that issuing a single policy does not guarantee that the number of new Covid-19 cases will However, the provincial government of Bali must issue further policies to get optimal results . decrease in the number of new case. REFERENCES