Jurnal Ekonomi Indonesia, 14 . , 2025: 1-23 p-ISSN: 0854-1507 . e-ISSN: 2721-222X Investigating inflation dynamics in Indonesia: Identifying the inflation spillover for enhancing regional inflation control Taly Purwaa. Ni Made Inna Dariwardanib*. Diyang Gita Cendekia c aBPS-Statistics of Bali Province. Denpasar. Indonesia b Faculty of Tourism. Universitas Udayana. Denpasar. Indonesia cBPS- Statistics of South Kalimantan Province. Banjarbaru. Indonesia Submitted: 29 September 2024 Ai Revised: 07 March 2025 Ai Accepted: 12 April 2025 Abstract This study employs a hybrid New Keynesian Phillips curve (NKPC) model, comprising the dynamic panel, spatial dynamic panel, and semiparametric spatial dynamic panel models, to examine the dynamics of regional inflation in Indonesia from Q1 1997 to Q2 2023. The latter model, which employed the inverse distance weight (IDW), was identified as the most optimal The findings indicated that the lagged inflation and output gap variables had a notable impact on regional inflation in Indonesia. Furthermore, the significance of the spatial spillover parameters indicates the existence of spatial regional inflation spillover in Indonesia. Jakarta. Central Java, and East Java represent Indonesia's three most significant focal points of inflation dynamics, exerting the most considerable influence on other regions. South Kalimantan. Jakarta, and East Java were the top three regions experiencing the most significant impact from inflationary pressures originating elsewhere. Sensitivity analysis found that a flattened Phillips curve was present throughout the period under study. It is therefore recommended that policies be intensified to control prices, strengthen policy synergies, and improve the efficiency of goods distribution channels, both for the focal point and for regions experiencing inflationary pressures due to spillovers. Keywords: direct and indirect effect. non-linear model. semiparametric model. spatial spillover JEL Classification: C31. E31. F36. Recommended Citation: Purwa. Taly. Dariwardani. Ni Made. , and Cendekia. Gita . Investigating inflation dynamics in Indonesia: Identifying the inflation spillover for enhancing regional inflation control. Jurnal Ekonomi Indonesia. , 2025: 1-23. DOI: https://doi. org/10. 52813/jei. Available at: https://jurnal. id/index. php/isei Copyright A2025 ISEI. This article is distributed under a Creative Commons Attribution-Share Alike 4. 0 International Jurnal Ekonomi Indonesia is published by the Indonesian Economic Association. OCorresponding Author: dariwardani. 2290711008@students. Jurnal Ekonomi Indonesia, 14 . , 2025: 1-23 Taly Purwa. Ni Made Dariwardani & Diyang CendekiaA Introduction Following the COVID-19 pandemic, the global economy faced high inflation pressures in response to . tightened monetary policies in various countries, . tighter policies related to the COVID-19 pandemic in China, and . increased geopolitical pressures as a result of the Russian-Ukrainian war (Bank Indonesia. Inflation in developed countries such as the United States and the European Union has reached its highest point in recent decades (Bank Indonesia, 2. Countries across various regions also experience the phenomenon of increasing post-pandemic inflation in developed countries. Despite similar pressures. Indonesian inflation is not as high as in the United States and the European Union. Indonesian inflation in 2022 was 5. 51 percent, lower than that of other Asian countries (Bank Indonesia, 2023. The interregional inflation correlation, known as the spatial spillover of inflation, is driven by the increased interstate dependence on global supply chains and the liberalization of international trade. Empirical studies found the existence of inflation spillover between regions and countries. (Osorio and Unsal, 2013. Halka and Szafranek, 2017. Istiak et al. , 2021. Al-Nassar and Albahouth, 2. In addition to the global economy, the existence of regional inflation spillover in a country was also empirically found by Aginta . Nagayasu . Yang et al. , and Yesilyurt & Elhorst . In the field of inflation control, studies related to inflation spillover in the regional economy of Indonesia, as an archipelagic country, are Thus, an in-depth analysis of regional inflation correlations is expected to increase the effectiveness of inflation control by understanding the mechanisms of interregional relations and estimating the potential impact on other regions. (Aginta, 2. Inflation control policy in Indonesia is placed within the policy framework of achieving inflation targets as the primary focus of monetary policy. However, the complexity of inflation determinants in Indonesia requires a policy mix, both fiscal and other structural policies, as an attempt to control the inflation rate along with monetary policy (Bank Indonesia, 2. As a response, the government has established the Regional Inflation Control Team or Tim Pengendalian Inflasi Daerah (TPID), which is comprised of related stakeholders involved in commodity price control efforts at the regional level. Price control at the regional level is one of the keys to price control on the national In 2023, the government provided fiscal incentives to districts that were consistently able to maintain inflation levels (Ministry of Finance Republic of Indonesia, 2. This attempt was made to motivate the regional government to control the price volatility in its region. Furthermore, with an interregional linkage or inflation spillover in Indonesia, it is also necessary to consider providing fiscal incentives to regions that potentially contribute to inflation in other regions. other words, price volatility control in inflation focal point regions will influence price stability in other regions. As the largest archipelagic country. Indonesia's regional growth pattern offers immense opportunities for spatial study and policy considerations (Aginta, 2. However, relatively few studies have been conducted regarding spatial spillover Jurnal Ekonomi Indonesia,14 . , 2025: 1-23 Investigating inflation dynamic in IndonesiaA inflation at the regional economic level, even though it has great potential as a reference for effective price control policies. Therefore, this study aims to carry out spatial inflation modeling to see the relationship of inflation between regions in Indonesia through regional level data . during the period Quarter I 1997Ae Quarter II 2023. The panel data will be analyzed further using a dynamic panel model and a dynamic panel spatial model based on the hybrid New Keynesian Phillips curve (NKPC) model. It was developed by Galy and Gertler . by relating inflation to lagged inflation, expected inflation, and real marginal cost. This study employs the output gap as a proxy of real marginal cost . redictor variabl. as discussed by Basarac, krabic, and Soric . This research will also explore the possibility of a non-linear relationship between predictor and response variables and then accommodate them in a dynamic panel spatial semiparametric model. addition, by utilizing the spatial model, this study attempts to find the existence of flattened Phillips curve phenomena during the period of the study through a sensitivity analysis. This research aims to provide several novelties. First, this research complements previous research related to the spatial spillover of regional inflation in Indonesia, such as that carried out by Aginta . , by utilizing an extended of a longer period of data. Second, to the best of the author's knowledge, this study will be the first to use data from the COVID-19 pandemic period while also covering the Monetary Crisis period as a result of the Asian Financial Crisis (AFC) in 1997Ae1998 and the Global Financial Crisis (GFC), which occurred around 2007Ae2009, in the field of spatial spillover of inflation in Indonesia. The results of this research are expected to provide input for policy-making related to regional inflation control to strengthen the mix and synergy of inflation control between regions in Indonesia. The Phillips curve model explains the relationship between inflation and the output gap, where a decrease in unemployment raises inflation. However, in recent decades, this relationship has weakened, known as the flattening of the Phillips Curve (Aginta, 2. As a result, central banks have difficulty determining the output gap that is in line with the inflation target (Galy and Gertler, 1. The flattening Phillips Curve phenomenon has led to the New Keynesian Phillips Curve (NKPC) emergence, which explains the relationship between inflation and the output gap. When this relationship weakens, inflation expectations and prices become increasingly crucial in determining inflation. With a forward-looking approach, this theory allows for explaining long-term inflation. However, the NKPC model often fails to describe the dynamics of inflation empirically. thus, it is developed by adding inflation lags . ackward-lookin. to account for past inflation behavior (Yesilyurt and Elhorst, 2. Galy & Gertler . introduced a hybrid NKPC model that combines the effects of inflation lags as a rule of thumb in inflation Various factors influence the formation of inflation in several regions using the hybrid NKPC model. One focus of the study is the inter-regional linkages or spatial spillover, which shows the influence of spatial interactions on inflation dynamics Jurnal Ekonomi Indonesia, 14 . , 2025: 1-23 Taly Purwa. Ni Made Dariwardani & Diyang CendekiaA (Halka and Szafranek, 2017. Aginta, 2. Inflation spillovers usually occur through supply chains and production networks (Auer. Levchenko, and Saure, 2021. Bilgin, 2. Production networks accelerate the transmission of inflation between sectors and countries when there are supply shocks, such as oil prices, but do not play a significant role in demand shocks (Bilgin, 2. Developed countries are more vulnerable to global shocks and have stronger spillover impacts, while developing countries are more influenced by domestic factors (Osorio and Unsal. Al-Nassar and Albahouth, 2. Empirical studies conducted by Al-Nassar & Albahouth . Halka & Szafranek . Istiak et al. , and Osorio & Unsal . show that inflation spillover has a significant impact on the formation of inflation between regions and between countries. Al-Nassar & Albahouth . and Istiak et al. conclude that Japan and the United States are the two countries with the largest spillover impacts in the G-7 and G-20. Inflation spillovers between the Euro Area and its five economic partners are wider outside the Eurozone (Halka and Szafranek, 2. , while China has the largest impact in ASEAN and India (Osorio and Unsal, 2. Osorio & Unsal's . research shows that more than 60 percent of inflation dynamics in Asia are influenced by domestic factors, especially in countries with a strong domestic demand base, such as China. India. Indonesia, and developed countries such as Japan and Korea. For Indonesia, which is an archipelagic country with significant price variations between regions, it is important to understand inflation dynamics at the regional level and their impact on the national economy (Aginta, 2. Several studies of inflation in Indonesia generally examine inflation convergence, spatial, and temporal interactions at the provincial and district/city levels (Tirtosuharto and Adiwilaga, 2014. Ridhwan, 2016. Purwono. Yasin, and Mubin, 2. Spillover studies in regional economies are still rare, although they are very important for understanding inflation dynamics in heterogeneous regions such as Indonesia. This research helps evaluate the impact of economic and monetary policies, as well as develop more specific regional policy recommendations (Aginta. Several studies on inflation spillovers to see price dynamics at the regional level include those conducted by Aginta . Nagayasu . Yang et al. and Yesilyurt & Elhorst . Yang et al. concluded that there was a significant spillover effect in the Chinese property market. Yesilyurt & Elhorst . conducted a study in Turkey using the output gap as an independent variable and found that the inflation lag in the province and its surrounding areas had a significant impact on inflation, while the output gap did not show a significant effect. On the other hand. Nagayasu . in Japan found a real effect of inflation spillover between provinces, with labor productivity as an independent variable. Aginta . found that spatial spillover of regional inflation in Indonesia has a significant effect on national inflation dynamics. This study uses spatial weighting in the Phillips curve model and finds differences in spillover effects between general inflation and core inflation. Significant spillover effects only occur in provinces with non-oil and gas sectors. In contrast to Yesilyurt & Elhorst . , this study shows that the output gap has significant regional spillover, although it does not yet cover the post-COVID-19 period and an in-depth analysis of the provinces that dominate inflation spillover. Jurnal Ekonomi Indonesia,14 . , 2025: 1-23 Investigating inflation dynamic in IndonesiaA Methodology This study utilizes panel data from 34 provinces in Indonesia covering the period from Quarter I 1997 to Quarter II 2023 . total of 106 quarter. , amounting to 3,604 The observation period includes the Asian Financial Crisis (AFC), which had an impact on the Monetary Crisis in Indonesia around 1997Ae1998, the Global Financial Crisis (GFC), which occurred around 2007Ae2009, as well as the COVID-19 pandemic during the 2020Ae2022 period. The main response and predictor variables used, respectively, are year-on-year . -on-. inflation and the output gap from Gross Regional Domestic Product (GRDP) at constant prices, obtained from the BPS-Statistics Indonesia. The following are the stages of analysis carried out in this research: Data preparation. y-on-y quarterly inflation data for 34 provinces is obtained by the following process Acquisition of initial monthly consumption value (CV) data for the period January 2020AeMarch 2023 as well as monthly month-to-month . -to-. inflation for the period January 1996AeJune 2023 from each city in Indonesia that is obtained from BPS-Statistics Indonesia. From these two data sets, the monthly CV for the period January 1996AeJune 2023 for each city in Indonesia is obtained. Hence quarterly inflation data can be obtained for the period Quarter I 1997Ae Quarter II 2023 for each province. For example, y-on-y inflation for the first quarter of 2023 for province A is calculated using the following formula: 2023Oeya = ( ya Ocycuycn=1 Oc2023ycAycayc yc=2023yaycaycu yaycOycnyc ycu ya Ocycn=1 Oc2022ycAycayc yc=2022yaycaycu yaycOycnyc Oe . y 100% . with ycuya being the number of city at province A and yaycOycnyc is the consumption value in the ycn-th city of CPI at the yc-th period. Before calculating the GDP output gap for Quarter I 1997AeQuarter II 2023 for each province, it is first necessary to obtain quarterly GRDP at constant prices data from 34 provinces in Indonesia. However, quarterly data has only been available since the 2010 period. Meanwhile, before that period. BPS-Statistics Indonesia only calculated provincial GDP in annual aggregates. Therefore, the following steps are taken: Initial data acquisition of Annual GRDP at constant prices 2010Ae2022 from 34 provinces with base year 2010 . 0 = . Annual GRDP at constant prices 2000Ae2013 . 0 = . , and 1996Ae2001 . 3 = . Then proceed by converting the base year of GRDP at constant prices to 2010 = 100. There are several provinces for which data is not yet available because the province has not yet been formed, i. Bangka Belitung in 1996-1999 . division of South Sumater. Riau Islands in 1996-2002 . division of Ria. Banten in 1996-1999 . division of West Jav. North Kalimantan in 1996-2009 . ivision from East Kalimanta. Gorontalo in 1996-1999 . ivision from North Sulawes. West Sulawesi in 1996-2004 . ivision from South Sulawes. North Maluku in 1996-1999 . ivision from Maluk. , and West Papua in 1996-2001 . ivision of Jurnal Ekonomi Indonesia, 14 . , 2025: 1-23 Taly Purwa. Ni Made Dariwardani & Diyang CendekiaA Papu. Therefore, annual GRDP at constant prices interpolation is carried out using the development pattern of the parent or initial province. The process of disaggregating annual GRDP at constant prices data for the 1996Ae 2022 period into quarterly data using the Denton-Cholette (Dagum & Cholette, 2. and Chow-Lin (Chow and Lin. , 1. temporal disaggregation methods, both without and with representative quarterly data indicators that describe development patterns of GRDP at constant prices. The indicators used in this study are Indonesia's quarterly GRDP at constant prices and the quarterly average IDX Composite Stock Price Index (JKSE) from Yahoo Finance. The best disaggregation results for GRDP at constant prices for Quarter I 1996AeQuarter IV 2009 are determined based on the smallest root mean square error (RMSE) value produced by the data disaggregated with actual data for the period Quarter I 2010AeQuarter IV 2022. A summary of the temporal disaggregation results is presented in Appendix 1. The seasonal adjustment process by removing or extracting seasonal patterns from the quarterly GRDP at constant prices series for each province, uses the X13ARIMA-SEATS method developed by the US Census Bureau (Monsell, 2. The process of calculating potential output, i. , the maximum potential of a region's GRDP in a certain period, by extracting long-term trend patterns from the quarterly GRDP at constant prices series, that has undergone seasonal The method used is the Hodrick-Prescott (HP) filter method (Hodrick and Prescott, 1. Finally, calculate the quarterly output gap for each province. For example, the output gap calculation formula for province A at the first quarter of 2023 is as ycuycycycyycycyciycaycyya. 2023Oeya = ycoycuyci ( ycaycaycycycayco yaycIyaycEya. 2023Oeya ycyycuycyceycuycycnycayco yaycIyaycEya. 2023Oeya Conduct data exploration to see developments in inflation and output gaps during the research period, both under normal and crisis conditions, as well as visually see the relationship between variables. Estimation of the inflation model. This study adopts the inflation model with the Phillip curve specification used by Aginta . , where in its first stage the following dynamic panel regression model is used: ycycnyc = yuycycnycOe1 yuycuycnyc yuAycnyc where ycycnyc and ycycnycOe1 are the inflation of the ycn-th province at yc and yc Oe 1 periods ycuycnyc is the output gap predictor variable of the ycn-th province at yc period, yu and yu are the parameters that will be estimated, and yuAycnyc is an error Next, include spatial interaction elements in the model through the Spatial Durbin Model (SDM) as follows: ycycnyc = yuUycOycycnyc yuycycnycOe1 yuycOycycnycOe1 yuycuycnyc yuycOycuycnyc yuAycnyc Jurnal Ekonomi Indonesia,14 . , 2025: 1-23 Investigating inflation dynamic in IndonesiaA where ycO is a weight matrix with size ycA y ycA, which shows the neighborhood structure of each province in Indonesia, with N = 34 . umber of province. , and yuU, yu, yu, yu, and yu are the parameters to be estimated. This research also uses a spatial autoregressive (SAR) model to accommodate if the parameters yu and yu are not significant with the following specifications: ycycnyc = yuUycOycycnyc yuycycnycOe1 yuycuycnyc yuAycnyc From the two models in Equations . , each spatial weighting matrix . cO) will be used, which represents the interaction or neighborhood of each analysis unit . based on the concepts of inverse distance weighting (IDW), contiguity of Thiessen polygons (Anselin. Syabri and Kho, 2010. Aginta, 2. , and K-nearest neighbor . NN). In this study, the ycO matrix is assumed to not change during the research period (Elhorst, 2. This research also adds elements of non-linear relationships to the SDM and SAR models above based on the scatter plot results between predictor variables and response variables. Hence a dynamic panel spatial semiparametric model is obtained (Mynguez. Basile and Durbyn, 2020, 2. as follows: ycycnyc = yuUycOycycnyc yuyc. cycnycOe1 ) yuycOyc. cycnycOe1 ) yuyc. cuycnyc ) yuycOyc. cuycnyc ) yuAycnyc or called as the PS-SDM model, and . ycycnyc = yuUycOycycnyc yuyc. cycnycOe1 ) yuyc. cuycnyc ) yuAycnyc as the PS-SAR model. Where yc. cycnycOe1 ) and yc. cuycnyc ) are P-splines functions for variables ycycnycOe1 and ycuycnyc , respectively. The selection of the best model from the various models explained in the third step above is carried out based on the AIC and RMSE values. Interpretation of the best model parameter estimation results. If the best model is the SDM. SAR. PS-SDM, or PS-SAR model, then the indirect impact . pillover effec. will be interpreted, i. the impact on inflation in a region caused by an increase or decrease in the value of a variable in the surrounding regions. Meanwhile, the direct impact is the impact on inflation in a region that is caused by changes in a variable in the same region. For example, to get the marginal effect matrix with direct and indirect effects in the PS-SDM model. Equation . first needs to be changed to: ycycnyc = . a Oe yuUycO)Oe1 . cycnycOe1 ) yuycOyc. cycnycOe1 ) yuyc. cuycnyc ) yuycOyc. cuycnyc ) yuAycnyc ) . Hence the marginal effect matrix for ycycnycOe1 is obtained as follows: cycnyc ) yccya. cycnyc ) ] = . a Oe yuUycO)Oe1 . uyaycA yuycO) A yccyc1yc yccycycAyc Jurnal Ekonomi Indonesia, 14 . , 2025: 1-23 Taly Purwa. Ni Made Dariwardani & Diyang CendekiaA The main diagonal and non-main diagonal of the matrix in Equation . respectively, show direct effects and indirect effects (LeSage and Pace, 2. Last, this study also tested the sensitivity of the parameter estimation results in the model by estimating the best model in several different ranges during the research period, i. the 1997-1999, 2000-2009, and 2010-2023 periods. The data analysis process in this study was carried out using Rstudio statistical The temporal disaggregation stage uses the temdisagg package (Sax and Steiner. , 2. , the seasonal adjustment using the X-13ARIMA-SEATS method uses the seasonal package (Sax and Eddelbuettel, 2. , and the HP filtering stage uses the hpfilter package (Monahov, 2. The Thiessen polygons are created using GeoDa software. Results and Discussion Inflation and Output Gap in Time of Various Crises Regional inflation rates in the nearly four decades of the Indonesian economy currently show spikes and declines in several periods, particularly during times of The Asian Financial Crisis (AFC), also known as the "monetary crisis," has exerted the greatest pressure on Indonesian inflation in the past four decades, when the inflation reached hyperinflationary conditions. Almost all provinces in Indonesia experienced high inflation during this period. even at the national level, inflation reached 77. 63 percent in 1998. The 2008 Global Financial Crisis (GFC) also brought inflationary pressures to almost all provinces in Indonesia, with national inflation reaching double digits . 06 percen. After the GFC. Indonesia's inflation was relatively controlled at single-digit levels as an implementation of the Inflation Targeting Framework (ITF) policy by the central bank (Bank Indonesi. However, during the COVID-19 pandemic that started in 2020. Indonesia's inflation experienced a decline and even reached the lowest record as calculated by BPS-Statistics Indonesia, i. 68 percent. This pattern is clearly shown in Figure 1 that inflation is higher in crisis conditions, both AFC and GFC, than in normal conditions. While inflation tends to be lower in COVID19 conditions. Jurnal Ekonomi Indonesia,14 . , 2025: 1-23 Investigating inflation dynamic in IndonesiaA Figure 1. Inflation . ear on yea. of 34 Provinces in Indonesia, 1997-2023 Source: Processed by Author Figure 2. Output Gap (Base on Regional GDP at Constant Pric. in 34 Provinces in Indonesia, 1997-2023 Source: Processed by Author In line with inflation movements, the regional output gap value also shows strong pressure during times of crisis, as presented in Figure 2. During the AFC period, the regional output gap showed a value moving away from zero, indicating a shock to the economy. In 1998, a positive output gap value was seen as an indication that the Indonesian economy was experiencing excess demand, causing the prices of goods and services . to rise. However, a year later, a recovery was observed, with the output gap value becoming negative, accompanied by relatively low inflation Meanwhile, during the GFC, the shock to the output gap did not appear to be as large as during the AFC. Meanwhile, during the COVID-19 pandemic, the output gap value tended to be positive, especially during the pandemic outbreak . , which implied that the Indonesian economy was operating below its potential capacity, as indicated by the low economic growth and inflation during this period. Jurnal Ekonomi Indonesia, 14 . , 2025: 1-23 Taly Purwa. Ni Made Dariwardani & Diyang CendekiaA Next, we will see the pattern of the relationship between the variables used in this research (Figure . A unidirectional linear relationship pattern can be seen between the variables inflation and its lag. While the output gap variable has a nonlinear relationship pattern visually. Based on these results, the output gap variable is used with a non-linear function (P-splin. in the modeling phase. In addition, a scatterplot between lagged inflation from lags 2 to 5 and actual inflation is presented, as well as a scatterplot between lagged output gap from lags 1 to 4 and actual inflation in Appendix 2. Figure 3. Scatter Plot of lag of Inflation . and lag of Output Gap . Actual Inflation . Source: Processed by Author Formation of Spatial Weighting Matrix (Spatial Weigh. The Inverse Distance Weight (IDW) spatial weighting matrix is calculated based on the distance of each centroid between provinces. The contiguity matrix of the polygon map of Indonesia cannot be formed because several provinces do not have direct borders. Therefore it needs to be converted into Thiessen polygon before a contiguity-based matrix can be constructed. The theoretical source of Thiessen polygon construction can be found in Yamada . Determining the number of neighbors (K) in the K-nearest neighbor . NN) based weighting matrix is done by exploratory model estimation with the number K=2,3,. ,11 and then selecting K based on the AIC value. In this study, the optimal number of K in the SDM. SAR. PSSDM, and PS-SAR models is 8. The results of the determination of the optimal number of K are presented in Appendix 3. Addition of Spatial and Non-Linear Components: Does It Improve Model Performance? This study uses the hybrid NKPC model specification (Galy and Gertler, 1. with and without the addition of spatial variables in the model estimation. There is one dynamic panel model . ithout adding spatial variable. and four spatial model specifications with a combination of three types of spatial weighting matrices, namely Inverse Distance Weight (IDW). Thiessen polygon, and K-nearest neighbor . NN). Thus, 13 model specifications were estimated, which were used to find the best estimation model based on the criteria for the lowest AIC and RMSE values. Jurnal Ekonomi Indonesia,14 . , 2025: 1-23 Investigating inflation dynamic in IndonesiaA summary of the estimation results from the dynamic panel model (PD), the dynamic panel spatial (SDM and SAR), and the dynamic panel spatial semiparametric models (PS-SDM and PS-SAR) is presented in Table 1. The results of the non-spatial model estimation demonstrate that the inflation lag and output gap exert a considerable influence on inflation, a finding that aligns with the predictions of the standard Phillips curve model. In addition, the spatial spillover parameters for inflation, inflation lag, and output gap variables in all spatial model specifications demonstrate statistical significance, indicating that there is a spillover effect among Indonesian provinces for these three variables. This result confirms the findings of the previous study conducted by Aginta . , but differs from the findings of the study conducted by Yesilyurt & Elhorst . , which found insignificant regional spillovers for the output gap. Therefore, with regard to regional inflation spillovers, this study validates the findings of previous studies conducted by Aginta . Nagayasu . Yang et al. , and Yesilyurt & Elhorst . Table 1. Summary of Estimation Results Dynamic Panel Intercept W inflation . Inflation (-. nflation (-. ) Output Gap W(Output Ga. SDM IDW SAR Thiessen IDW Thiessen s(Output Ga. (Output Ga. ) K kNN AIC 19,673 10,508 11,879 11,195 12,266 12,551 12,301 RMSE Source: Processed by Author A **) smallest AIC and RMSE values Note: *) significant at the significance level yu= 0. Taly Purwa. Ni Made Dariwardani & Diyang CendekiaA Jurnal Ekonomi Indonesia, 14 . , 2025: 1-23 Table 2. Summary of Estimation Results (Continue. PS-SDM IDW Thiessen PS-SAR IDW Thiessen Intercept W inflation . Inflation (-. nflation (-. ) s(Output Ga. (Output Ga. ) AIC 10,477** 11,793 11,099 12,247 12,513 12,277 RMSE Output Gap W(Output Ga. K kNN Source: Processed by Author Note: *) significant at the significance level yu= 0. 05, **) smallest AIC and RMSE values In general, the parameter estimation results for the whole model are not significantly different and have the same direction. Different directions were found only for the estimated intercept parameters. The addition of spatial elements to the SDM and SAR models can improve the goodness of fit of the model based on AIC values compared to the dynamic panel model. However, this is not the case when looking at the ability of the model to make predictions based on the RMSE values. Moreover, the use of the spatially lagged inflation and output gap variables in the SDM and PS-SDM models is able to improve the model's goodness of fit based on AIC values compared to the use of only the SAR and PS-SAR models. The use of the output gap variable in a non-linear form, namely the P-spline, is able to improve model performance in terms of both AIC and RMSE values, although it is not optimal as this variable is not significant in the PS-SDM model with an IDW weight matrix. The best model based on the minimum AIC and RMSE criteria is the semiparametric spatial dynamic panel model (PS-SDM) specification with the IDW spatial weighting matrix. The optimal model, as determined by the minimum AIC and RMSE criteria, is the semiparametric spatial dynamic panel model (PS-SDM) specification with an IDW spatial weighting matrix. These findings will enhance the existing body of literature on the subject, particularly in relation to the use of the Spatial Durbin Model (SDM) to examine regional inflation dynamics and spatial spillovers, as demonstrated by Aginta . and Nagayasu . The incorporation of semiparametric models has been demonstrated to enhance the performance of the resulting models, as evidenced by the findings of this study. this model specification, the inflation lag, the spatial aspects of inflation, and the spatial aspects of the inflation lag have a significant impact on the formation of Jurnal Ekonomi Indonesia,14 . , 2025: 1-23 Investigating inflation dynamic in IndonesiaA regional inflation in Indonesia. Meanwhile, the semiparametric element of the output gap and its spatial form have no significant effect on the inflation rate. Dynamics of Inflation and Output Gap in Provinces in Indonesia Based on the best model, i. the PS-SDM model, interpretation will then be carried out on the direct and indirect effects . pillover effect. of the lagged of inflation and the output gap on the actual inflation. The marginal effect of these two variables in the PS-SDM model can be obtained from the matrix in Equation 5. Table 3. Summary of Marginal Effect on the PS-SDM Model Variable Average Direct Impact Average Indirect Impact (Spatial Spillove. Inflation (-. Output Gap Source: Processed by Author Interesting findings regarding the impact of lagged inflation and the output gap on the actual inflation at the provincial level in Indonesia are presented in Table 2. Looking at the magnitude of the impact, both direct and indirect, it can be seen that the largest impact is transmitted by the output gap variable. This means that the actual inflation in a particular province is influenced by the development of its economic level in the same period. The magnitude of the direct impact of the output gap is -3. in other words, if there is an increase . in the output gap of a province by 1 percent, it will cause a decrease . in inflation in the same province by 3. 469 percent on average. Meanwhile, the magnitude of the output gap spillover is -0. 013, meaning that an increase . in the output gap of 1 percent in neighboring provinces will cause a decrease . in inflation in another province by 0. 013 percent on average. In contrast to the output gap, which has an impact in the same direction . oth marginal effects are negativ. , the lagged inflation has a marginal effect in a different The direct impact has a positive direction, i. 972, while the spillover impact is -0. In other words, the lagged inflation in a certain province will have a positive effect on the actual inflation in the same province. Meanwhile, on average, the laged inflation in neighboring provinces has a negative impact on the actual inflation in a certain province. Based on the magnitude of these two impacts, the magnitude of the direct impact from the province itself is much greater than of the spillover impact from neighboring provinces. Focal Points of Inflation in Indonesia This study will diverge from previous studies, such as by Aginta . Nagayasu . Yang et al. , and Yesilyurt & Elhorst . , which concentrate more on general regional inflation spillovers. Instead, this study will endeavor to Taly Purwa. Ni Made Dariwardani & Diyang CendekiaA Jurnal Ekonomi Indonesia, 14 . , 2025: 1-23 investigate which regions are both contributors and recipients of inflation To elaborate, the matrix of marginal effects generated in the preceding discussion can be analyzed at the provincial and island levels. In a spatial context, a change in the lagged inflation variable and the output gap in a specific province will undoubtedly give rise to alterations in actual inflation in other neighboring provinces . he spillover effec. Table 4. Average Impact Created by Each Island on Actual Inflation in Other Provinces Island Inflation (-. Rank Output Gap Rank Sumatera Java Bali Nusra Kalimantan Sulawesi Maluku & Papua Source: Processed by Author Table 4 summarises the average impact received by surrounding provinces when there is a change in a variable in one province on each island in Indonesia. It can be seen that, on average, the provinces in Sumatera have the largest role . in influencing the actual level of inflation in other provinces. The second largest influence comes from the island of Java, followed by Sulawesi. Kalimantan. Maluku, and Papua. This order holds for both the inflation lag and the output gap variables. The magnitude of the impact created by each province is presented in Appendix 4. Looking more closely at the national level, the provinces with the largest impact on the inflation levels are the DKI Jakarta. Central Java, and East Java. Knowing the role of each region, both islands and provinces, we then see how vulnerable the actual inflation on an island or province is to changes in a variable in other provinces. The impact summary results in Table 4 show the same ranking as the results in Table 3, with Sumatera being the most vulnerable island to inflation when there is a change in the inflation lag and output gap in other islands. The average impact of each province is shown in more detail in Appendix 5. When there is a change in the inflation lag or output gap in the surrounding provinces, the provinces most affected are South Kalimantan. DKI Jakarta, and East Java. Jurnal Ekonomi Indonesia,14 . , 2025: 1-23 Investigating inflation dynamic in IndonesiaA Table 5. Average Impact Received by Actual Inflation for Each Island Due to Changes in Variables in Other Provinces Island Inflation (-. Rank Output Gap Rank Sumatera Java Bali Nusra Kalimantan Sulawesi Maluku & Papua Source: Processed by Author 6 Sensitivity of Estimation Results in Different Periods: Is There a Flattened Phillips Curve Phenomenon? Aginta . used the idea of a flattening Phillips curve in his research on regional inflation spillovers in Indonesia. However, because of the study's limited data coverage, the flattening Phillips curve phenomenon on regional inflation in Indonesia cannot be empirically demonstrated. As a result, the study does not objectively demonstrate this occurrence on regional inflation data in Indonesia. However using a larger dataset, this research proved the empirical validity of this phenomenon, especially when considering the three major crises of the AFC. GFC, and COVID-19 pandemic. The estimation results of the best model, i. the PS-SDM, in different periods, with each period containing crisis conditions (AFC. GFC, and COVID-. , show different The degree of actual inflation spillover between provinces has a positive value with an increasing trend between periods, from 0. 562 in the 1997Ae1999 period to 0. 761 in 2010Ae2023. Jurnal Ekonomi Indonesia, 14 . , 2025: 1-23 Taly Purwa. Ni Made Dariwardani & Diyang CendekiaA Table 6. Summary of Sensitivity Test Results PS-SDM Intercept W inflation . Inflation (-. nflation (-. ) s(Output Ga. (Output Ga. ) Source: Processed by Author Note: *) significant at the significance level yu= 0. The inflation lag variable has a positive influence with a pattern of decreasing In other words, the influence of lagged inflation in a certain province is decreasing relative to the actual inflation of the same province. The spatial lag inflation variable from other provinces has a negative influence on the actual inflation in a given province. The magnitude decreases over time, from -0. 543 in the 1997-1999 period to -0. 601 in the 2010-2023 period. Finally, the estimated parameters of the semiparametric variable output gap show that there is a decrease in the Phillips curve, from 5. 842 in the 1997-1999 period to 0. 040 in the 2010-2023 period, which means that a flattened Phillips curve has occurred in Indonesia during the this period. Conclusion This study has empirically examined the dynamics of regional inflation in Indonesia based on the Phillips curve theory . he hybrid New Keynesian Phillips curv. by using data from 34 provinces and a quarterly data series for the period 1997-2023, which covers three major crisis conditions in Indonesia . he Asian financial crisis, the global financial crisis, and the COVID-19 pandemi. The addition of spatial variables in this modeling was done considering that Indonesia is the largest archipelagic country in the world, so the inter-regional linkage factor is an important aspect to consider in explaining regional inflation dynamics. Inflation lags and output gaps have been shown to have a significant impact on regional inflation in Indonesia, both in modeling with and without the addition of spatial elements. However, addition of spatial elements to the model specifications has been shown to increase the performance of the model and prove the importance of inflation spillovers in regional inflation dynamics in Indonesia. The dynamic panel spatial semiparametric model (PS-SDM) with the IDW weighting matrix is found to be the best model that can be used to examine regional inflation dynamics in Indonesia. This model specification concludes that the spatial Jurnal Ekonomi Indonesia,14 . , 2025: 1-23 Investigating inflation dynamic in IndonesiaA aspect of inflation, the lagged inflation, and the spatial aspects of the lag of inflation have a significant impact on regional inflation in Indonesia. The level of inflation in a province is more dominantly influenced by the performance of its economic level . hich is reflected in the value of the marginal effect of the output ga. and the regional lag of inflation. Meanwhile, the spillover effects of both the output gap and inflation lag from surrounding areas are still worth considering, although the effects are relatively small. The island of Sumatra is shown to have the largest inflation spillover effect compared to other islands in Indonesia. The three provinces that play the largest role in influencing the level of inflation are DKI Jakarta. Central Java, and East Java. This has implications for the need for more intensive price control policies in these focus areas, such as the provision of fiscal incentives. With these fiscal incentives, it is hoped that regional governments, especially the TPID, will have more room to design price control programs in their regions. This is important given that price control in focal areas will have an impact on other areas that receive spillover effects from that region. Meanwhile, the provinces that have received the highest inflation spillovers from the fluctuation of inflation to their surroundings are the provinces of South Kalimantan. DKI Jakarta, and East Java. These regions need to synergize price control policy with the regions that have an impact on inflation spillovers. The proportional fiscal stimulus is also important for these provinces to maintain price stability from the impact of inflation fluctuations from other provinces as inflation focal points. In addition, the disruptions of distribution channels and transportation costs betweenregions must also be anticipated to reduce price fluctuations. In addition, the results of the sensitivity test show that there is a shift in the Phillips curve over the period 1997-2023, with the magnitude of the estimated parameters for the output gap variables tending to decrease . he flattened Phillips For future studies, it would be interesting to examine inflation spillover studies using core inflation data, so that temporary shock factors in headline inflation data can be eliminated. When constructing the output gap variable, it is also necessary to consider separating the regional economies dominated by the oil and gas sector from the non-oil and gas sector, as they generally have different economic References