Identification and Inference for Synthetic Control Methods with Spillover Effects: Estimating the Economic Cost of the Sudan Split

Shosei Sakaguchi ∗ and Hayato Tagawa †

October 8, 2024

## Abstract

The synthetic control method (SCM) is widely used for causal inference with panel data, particularly when there are few treated units. SCM assumes the stable unit treatment value assumption (SUTVA), which posits that potential outcomes are unaffected by the treatment status of other units. However, interventions often impact not only treated units but also untreated units, known as spillover effects. This study introduces a novel panel data method that extends SCM to allow for spillover effects and estimate both treatment and spillover effects. The method leverages a spatial autoregressive panel data model to account for spillover effects. The study also proposes a Bayesian inference method using Bayesian horseshoe priors for regularization. The proposed method is applied to two empirical studies: one evaluating the effect of the California tobacco tax on consumption, and the other estimating the economic impact of the 2011 division of Sudan on GDP per capita.

Keywords: Synthetic control; Spillover effect; Spatial data; Bayesian inference; Sudan split.

∗ Faculty of Economics, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan.

Email: sakaguchi@e.u-tokyo.ac.jp.

† Nowcast Inc, Sumitomo RD Kudan Bldg. 9F, 1-8-10 Kudankita, Chiyoda-ku, Tokyo 102-0073, Japan.

Email: h.tagawa8810@gmail.com.

## 1 Introduction

The synthetic control method (SCM) (Abadie and Gardeazabal, 2003; Abadie et al., 2010) is a causal inference approach used to estimate treatment effects when the number of units undergoing intervention is very small, such as a single unit. This approach identifies treatment effects from panel data by substituting the counterfactual outcomes of the treated unit in the absence of treatment with weighted averages of outcomes of untreated units. The method has been widely used in various fields, including political economy and marketing. Athey and Imbens (2017) describe SCM as 'arguably the most important innovation in the policy evaluation literature in the last 15 years.'

The SCM is typically applied to country- and district-level data. For instance, Abadie et al. (2015) apply SCM to country-level data to estimate the economic impact of the 1990 German reunification, while Bifulco et al. (2017) apply this method to school-district-level data to evaluate the effect of an educational program, among other examples. In such contexts, interventions may impact untreated units due to geographic or socioeconomic connections among countries or districts, which are known as spillovers. However, SCM relies on the stable unit treatment value assumption (SUTVA) (Rubin, 1978), which posits that one's outcome is not dependent on the treatment status of others, excluding spillovers. If spillovers occur to untreated units, SUTVA is violated, which can lead to biased estimates of treatment effects in the SCM.

This paper tackles this problem by extending the standard SCM to allow for spillover effects. We leverage the spatial autoregressive (SAR) model, a workhorse model in spatial data analysis, to address this issue with its capability to capture the correlation of outcomes among units. We propose a new approach to identify and estimate treatment effects by incorporating a SAR model into SCM, where we characterize the outcomes of untreated units with the SAR panel data model. This approach allows for spillover effects arising from the dependence of outcomes among treated and untreated units, thereby relaxing SUTVA in SCM. Furthermore, the extended SCM can identify and estimate spillover effects on untreated units, which are also parameters of interest in many empirical studies adopting SCM.

Building on the identification results, we propose a Bayesian inference method for the proposed SCM. Data targeted by SCM often involves panel data, but the number of units or the length of pretreatment periods is not always sufficiently large. In such cases, frequentist approaches may not provide accurate inference. Bayesian inference offers several advantages. First, it can yield more accurate estimates than frequentist methods when the sample size is small. Second, it facilitates sta-

tistical inference via the Markov Chain Monte Carlo (MCMC) procedure. Third, Bayesian modeling is flexible, allowing models to be less dependent on the assumptions about prior distributions.

We apply Bayesian regularization in the construction of synthetic controls. Following Kim et al. (2020), we use Bayesian horseshoe priors (Carvalho et al., 2010) to model the synthetic weights. The Bayesian horseshoe prior has strong regularization effects, and hence is effective in avoiding overfitting bias. This is particularly important in SCM, as it is often applied to data with a relatively large number of control units and short pretreatment periods, which can easily cause overfitting bias in the estimation of synthetic control outcomes. The simulation study in the paper examines the finite sample performance of the proposed Bayesian inference.

This paper presents two substantive empirical studies with the application of the proposed SCM. One study evaluates the impact of California's tobacco tax on its consumption (Abadie et al., 2010), and the other estimates the impact of Sudan's north-south split in 2011 on GDP per capita. Both studies involve spillover effects among US states and African countries. Regarding California's tobacco tax, our empirical results show a negative impact on tobacco consumption in California, supporting the findings of Abadie et al. (2010). We also estimate the spillover effects on other states, showing that California's tobacco tax reduced tobacco consumption in other states.

To the best of our knowledge, the second empirical study is the first to examine the economic impact of Sudan's north-south split on the Sudans (the region of the former united Sudan). In contrast, Mawejje and McSharry (2021) estimate the economic impact of the 2012 oil production halt in South Sudan on South Sudan itself using the standard SCM. Mawejje and McSharry (2021) intentionally exclude countries neighboring South Sudan from the control group to address spillover concerns. This is a common practice in dealing with spillover in SCM (Abadie, 2021); however, excluding some untreated units, particularly neighboring countries, can lead to poor fitting of synthetic controls and cause additional bias. The SCM proposed in this paper does not need to exclude any untreated units. The results of our empirical study show that Sudan's north-south split decreased GDP per capita in the Sudans by about 9.5%, with a cumulative reduction of 34% from 2011 to 2015. Our empirical results also show the negative economic impacts of Sudan's split on other African countries with strong economic connections to Sudan, such as Egypt and Kenya.

## Related Literature

Since Abadie and Gardeazabal (2003) and Abadie et al. (2010) pioneered SCM, along with its broad application in the social sciences, SCM has been theoretically and methodologically advanced by

many works. Abadie and L'hour (2021) present a penalized synthetic control method to reduce interpolation biases. Arkhangelsky et al. (2021) combine SCM with difference-in-differences (DID) and propose a new estimator with robustness properties. Ben-Michael et al. (2021) propose an augmented synthetic control method to correct bias resulting from imperfect pre-treatment fit. Li (2020) and Chernozhukov et al. (2021) propose inference methods for SCM. Ferman and Pinto (2021) point out a potential bias in the SCM when the perfect fit assumption is not satisfied.

Several works have proposed Bayesian estimation in SCM to facilitate statistical inference through MCMC sampling and/or to apply flexible modeling. Kim et al. (2020) introduce a Bayesian synthetic control method that uses the Bayesian horseshoe prior proposed by Carvalho et al. (2010) for the synthetic weights. Their simulations demonstrate more accurate predictions of counterfactual outcomes than the standard SCM, particularly when the number of units is relatively large compared to the length of the pretreatment periods. Brodersen et al. (2015) propose a method based on Bayesian structural time-series models. Pang et al. (2022) and Klinenberg (2023) propose approaches to correct bias in SCM induced by the dynamic characteristics of untreated units. Pang et al. (2022) propose a Bayesian posterior predictive approach with a latent factor term that incorporates unit-specific time trends. Klinenberg (2023) incorporates time-varying parameters by using a state space framework.

While the related literature mentioned so far assumes SUTVA in SCM, few studies have considered its relaxation. Cao and Dowd (2019) allow for spillover effects in SCM by assuming a specified structure of spillover effects. Our approach differs from theirs in two aspects: (i) our approach supposes a spatial correlation structure for the outcomes of units rather than a direct structure for spillover effects, and (ii) our approach supposes the existence of a synthetic control only for the treated unit, whereas their approach supposes synthetic controls for both the treated and untreated units. Menchetti and Bojinov (2022) assume that a unit's potential outcome depends on its own and neighboring treatment assignments. They propose a Bayesian structural time series method to identify and estimate treatment and spillover effects by estimating the predictive distribution of counterfactual outcomes. Grossi et al. (2020) focus on the average spillover effects within specified clusters and provide a method for constructing confidence intervals for treatment and average spillover effects. In the context of DID, Miguel and Kremer (2004) propose an approach for estimating spillover effects using RCT data, which utilizes non-compliance in treatment groups to identify spillover effects.

## Structure of the Paper

The remainder of the paper proceeds as follows. Section 2 describes the SCM setup with spillover effects. Section 3 introduces the SAR panel data model, and presents the identification results. Section 4 presents our proposal for Bayesian SCM. Section 5 shows simulation results to demonstrate the finite-sample performance of the proposed SCM. Section 6 presents the results of two empirical applications. Section 7 concludes the paper.

## 2 Setup

Consider units i = 0 1 , , · · · , N and time points t = 1 2 , , · · · , T . We suppose that i = 0 is the unit receiving the treatment, and i = 1 , · · · , N are the units in the control group. Let T 0 be the number of pretreatment periods with 1 ≤ T 0 &lt; T . The treatment status of a unit i at time t is denoted as d it ∈ { 0 1 , } , where d it = 1 means that unit i receives treatment at time t , and d it = 0 means otherwise. In our context, d it = 1 when i = 0 and t &gt; T 0 and d it = 0 otherwise. We define d t ≡ ( d 0 t , d 1 t , · · · , d Nt ) ∈ { 0 1 , } N +1 , the treatment status vector for time t .

We consider the potential outcome framework. While many studies assume that one's potential outcome is not influenced by the treatment status of others (SUTVA) (Rubin, 1978), this study allows for the potential outcome to be affected by the treatment status of others. Thus, for treatment status d t at time t , the potential outcome of unit i is denoted as Y it ( d t ). 1 The observed outcomes are

$$Y _ { i t } \equiv \begin{cases} Y _ { i t } ( 0, 0, \cdots, 0 ) & \text{if $t \leq T_{0}$}, \\ Y _ { i t } ( 1, 0, \cdots, 0 ) & \text{if $t > T_{0}$} \end{cases}$$

for all i = 0 1 , , · · · , N and t = 1 2 , , . . . , T .

We define the treatment effect ξ 0 t for the treated unit 0 at time t ( &gt; T 0 ) as follows:

$$\xi _ { 0 t } \equiv \underbrace { Y _ { 0 t } ( 1, 0, \cdots, 0 ) } _ { \text{observed outcome} } - \underbrace { Y _ { 0 t } ( 0, 0, \cdots, 0 ) } _ { \text{counterfactual outcome} }$$

where Y 0 t (1 0 , , · · · , 0) is the outcome when unit 0 is treated, and Y 0 t (0 , 0 , · · · , 0) is the outcome when the unit is not treated, a counterfactual outcome not observed in the data.

1 SUTVA means that the potential outcome of each unit i is represented by Y it ( d it ), which does not depend on the treatment status of others d t \ d it ,

The framework of this study can capture the spillover effects on units in the control group. Although no unit in the control group receives treatment, the influence of unit 0 receiving treatment may affect the outcomes of the control group units, which is referred to as spillover effects. We define the spillover effect ξ it for any untreated unit i ( ≥ 1) and time t ( &gt; T 0 ) as follows:

$$\xi _ { i t } \equiv \underbrace { Y _ { i t } ( 1, 0, \cdots, 0 ) } _ { \text{observed outcome} } - \underbrace { Y _ { i t } ( 0, 0, \cdots, 0 ) } _ { \text{counterfactual outcome} } \.$$

The spillover effect is the counterfactual difference in the outcomes of an untreated unit when unit i = 0 receives treatment versus when unit i = 0 does not.

For the outcomes of the treated unit, similar to Abadie et al. (2010), we suppose the following assumption holds.

Assumption 1 (Perfect Fit) . There exists a vector of weights α = ( α , α 1 2 , · · · , α N ) ⊤ ∈ R N that satisfies the following: For each t = 1 2 , , · · · , T ,

$$Y _ { 0 t } ( 0, 0, \cdots, 0 ) = \sum _ { i = 1 } ^ { N } \alpha _ { i } Y _ { i t } ( 0, 0, \cdots, 0 ) \text{ a.s.}$$

This assumption implies that, in the absence of treatment, the control outcome of the treated unit can be replicated by the weighted average of the control outcomes of the untreated units. As pointed out by Ferman and Pinto (2021), this is a fundamental assumption for the SCM and is explicitly or implicitly assumed in SCM studies. Because we allow outcomes to depend on the treatment status of others, Assumption 1 pertains to a state where all units are untreated.

Under Assumption 1, we can estimate the synthetic weights α by solving the following leastsquares problem:

$$\widehat { \alpha } = \underset { \alpha \in \mathbb { R } ^ { N } } { \arg \min } \underset { t = 1 } { \overset { T _ { 0 } } { \sum } } \left ( Y _ { 0 t } - \underset { i = 1 } { \overset { N } { \sum } } \alpha _ { i } Y _ { i t } \right ) ^ { 2 }.$$

Using ̂ α , the standard SCM (Abadie et al., 2010) estimates the treatment effect ξ 0 t ( t &gt; T 0 ) by Y 0 t -∑ N i =1 ˆ α Y i it . Under SUTVA and Assumption 1 (perfect fit), Abadie et al. (2010) show the identification and consistent estimation of treatment effects by SCM. However, when SUTVA is violated, the standard SCM can be biased in its estimation of treatment effects. The subsequent section clarifies the sources of this bias and proposes a new approach to identify treatment and spillover effects while allowing for spillover.

## 3 Identification

This section first clarifies the bias in standard SCM induced by spillover effects. Subsequently, we introduce the SAR panel data model and discuss its capability to capture spillover effects. We then discuss the identification of the treatment and spillover effects.

## 3.1 Bias for the Standard SCM

We first show that when the SUTVA is not satisfied, the standard SCM can be biased due to spillover. Assumption 1 implies that the treatment effects ξ 0 t can be expressed as Y 0 t (1 , 0 , · · · , 0) -∑ N i =1 α Y i it (0 0 , , · · · , 0). However, when the outcomes depend on the treatment statuses of other units, we cannot observe counterfactual control outcomes Y it (0 , 0 , · · · , 0) after the pre-treatment periods ( t &gt; T 0 ). This causes bias in the standard SCM as follows:

$$Y _ { 0 t } - \sum _ { i = 1 } ^ { N } \alpha _ { i } Y _ { i t } & = Y _ { 0 t } ( 1, 0, \cdots, 0 ) - \sum _ { i = 1 } ^ { N } \alpha _ { i } Y _ { i t } ( 1, 0, \cdots, 0 ) \\ & = \xi _ { 0 t } + \underbrace { \sum _ { i = 1 } ^ { N } \alpha _ { i } ( Y _ { i t } ( 0, 0, \cdots, 0 ) - Y _ { i t } ( 1, 0, \cdots, 0 ) ) } _ { \text{bias} },$$

where Y 0 t -∑ N i =1 α Y i it is the standard SCM estimator of ξ 0 t given knowledge of α . This result shows the existence of bias in the standard SCM when SUTVA is not satisfied.

## 3.2 Spatial Autoregressive Model

To address the issue of bias caused by spillover, we introduce a model that captures spillover effects between treated and untreated units. Specifically, for the outcomes of units in the control group, we assume the following SAR panel data model for each t ≥ 1:

$$Y _ { t } ^ { c } ( d _ { t } ) = \rho ( w Y _ { 0 t } ( d _ { t } ) + W Y _ { t } ^ { c } ( d _ { t } ) ) + X _ { t } \beta + u _ { t } ^ { d _ { t } },$$

where Y c t ( d t ) ≡ ( Y 1 t ( d t ) , Y 2 t ( d t ) , · · · , Y Nt ( d t )) ⊤ ∈ R N . The vector w ∈ R N and matrix W ∈ R N × N comprise spatial weights, which are to be specified a priori. For i = 1 2 , , . . . , N and j = 0 1 , , . . . , N , we denote by w ij the spatial weight between units i and j ; that is, the ( i, j + 1)-the element of ( w W , ) ∈ R N × (1+ N ) . Typical examples of spatial weights include adjacent weights, where w ij is 1 if units i and j are adjacent and 0 otherwise. Geographic distance (e.g., distance between the

capitals of two countries) and economic distance (e.g., trade amount between two countries) are also often employed as spatial weights in the literature on spatial data analysis (LeSage and Pace, 2009). X t ≡ ( X 1 t , . . . , X Nt ) ⊤ ∈ R N × k denotes the covariate matrix for control units, where X it denotes a vector of covariates for unit i at time t . u d t t ≡ ( u d t 1 t , · · · , u d t Nt ) ⊤ ∈ R N is a vector of error terms.

The SAR panel data model (3) captures the spillovers arising from the spatial dependence of the outcomes among treated and untreated units. The magnitude of the spillover effects depends on the value of ρ as well as those of the spatial weights w and W , which represent the strength of spatial correlation. Several estimation methods are proposed to estimate ρ in SAR panel data models (e.g., Fingleton, 2008; Lee and Yu, 2010; Su, 2012; Glass et al., 2016; Liang et al., 2022).

## 3.3 Identification

We discuss the identification of the treatment and spillover effects under Assumption 1 and the SAR panel data model (3). For simplicity, since the treatment state d it for all i ≥ 1 and t is always 0, we write

$$Y _ { i t } ( 1 ) = Y _ { i t } ( 1, 0, \cdots, 0 ) \quad ( \forall i, t ),$$

$$Y _ { i t } ( 0 ) = Y _ { i t } ( 0, 0, \cdots, 0 ) \quad ( \forall i, t ).$$

$$&, 0, \cdots, 0 ) \quad ( \forall i, t ), \\ &, 0, \cdots, 0 ) \quad ( \forall i, t ).$$

We assume the following condition holds.

Assumption 2. For t = 1 2 , , · · · , T and i = 1 2 , , · · · , N , the error term u d t it depends only on its own treatment state d it ; that is, u d t it = u d ′ t it a.s. for any d t = ( d 0 t , d 1 t , . . . , d 1 T ) ⊤ and d ′ t = ( d ′ 0 t , d ′ 1 t , . . . , d ′ 1 T ) ⊤ such that d it = d ′ it .

This assumption implies that, since d it = 0 for any control unit i at any time t , the error terms for the control units always satisfy u d t t = u 0 N +1 t , where 0 N +1 represents an ( N + 1)-dimensional vector of zeros. For simplicity, we denote u d t t (= u 0 N +1 t ) as u t .

We further assume that the model (3) and the synthetic weights α in Assumption 1 satisfy the following assumption.

Assumption 3. I N -ρ wα ⊤ -ρ W is full rank.

This assumption ensures that I N -ρ wα ⊤ -ρ W has an inverse matrix. Assumption 3 is testable given estimators of α and ρ .

Given the parameters α and ρ , the treatment effects and spillover effects can be identified as follows. For each t &gt; T 0 , under Assumption 1 and the SAR panel data model (3), we obtain

$$Y _ { t } ^ { c } ( 0 ) & = \rho w Y _ { 0 t } ( 0 ) + \rho W Y _ { t } ^ { c } ( 0 ) + X _ { t } \beta + u _ { t } \\ & = \rho w \alpha ^ { \top } Y _ { t } ^ { c } ( 0 ) + \rho W Y _ { t } ^ { c } ( 0 ) + X _ { t } \beta + u _ { t }.$$

Thus, because the inverse matrix of ( I N -ρ wα ⊤ -ρ W ) exists under Assumption 3, arranging the above equation yields:

$$Y _ { t } ^ { c } ( 0 ) = \left ( I _ { N } - \rho w \alpha ^ { \top } - \rho W \right ) ^ { - 1 } ( X _ { t } \beta + u _ { t } ),$$

which leads to

$$\xi _ { 0 t } & = Y _ { 0 t } ( 1 ) - \alpha ^ { \top } Y _ { t } ^ { c } ( 0 ) \\ & = Y _ { 0 t } ( 1 ) - \alpha ^ { \top } \left ( I _ { N } - \rho w \alpha ^ { \top } - \rho W \right ) ^ { - 1 } ( X _ { t } \beta + u _ { t } ).$$

Furthermore, under the model (3) and Assumption, 2, for each t &gt; T 0

$$( I _ { N } - \rho W ) Y _ { t } ^ { c } ( 1 ) - \rho w Y _ { 0 t } ( 1 ) = X _ { t } \beta + u _ { t }.$$

Therefore, for t &gt; T 0 , we obtain

$$\xi _ { 0 t } & = Y _ { 0 t } ( 1 ) - \alpha ^ { \top } ( I _ { N } - \rho w \alpha ^ { \top } - \rho W ) ^ { - 1 } ( ( I _ { N } - \rho W ) Y _ { t } ^ { c } ( 1 ) - \rho w Y _ { 0 t } ( 1 ) ) \\ & = Y _ { 0 t } - \alpha ^ { \top } ( I _ { N } - \rho w \alpha ^ { \top } - \rho W ) ^ { - 1 } ( ( I _ { N } - \rho W ) Y _ { t } ^ { c } - \rho w Y _ { 0 t } ).$$

If the parameters α and ρ are estimated from the sample of the pretreatment periods ( t ≤ T 0 ), ξ 0 t ( t &gt; T 0 ) can be estimated as equation (4). Note that α can be consistently estimated by the least-squares method (2) under Assumption 1, and several methods are applicable to estimate ρ in the literature of the SAR panel data model (e.g., Fingleton, 2008; Lee and Yu, 2010; Su, 2012; Glass et al., 2016; Liang et al., 2022).

The following theorem summarizes the identification result of ξ 0 t .

Theorem 1 (Identification of the Treatment Effect) . Suppose that Assumptions 1, 2, and 3 hold. Then, given ρ and α , the treatment effect ξ 0 t for t &gt; T 0 can be identified as follows:

$$\xi _ { 0 t } = Y _ { 0 t } - \alpha ^ { \top } \left ( I _ { N } - \rho w \alpha ^ { \top } - \rho W \right ) ^ { - 1 } \left ( ( I _ { N } - \rho W ) Y _ { t } ^ { c } - \rho w Y _ { 0 t } ).$$

Proof. The result follows from the discussion above.

The discussion above shows that the counterfactual outcome Y c t (0) for each untreated unit for t &gt; T 0 can be identified as ( I N -ρ wα ⊤ -ρ W ) -1 (( I N -ρ W Y ) c t -ρ w Y 0 t ) . Therefore, given the parameters α and ρ , the spillover effects on the untreated units are also identifiable.

Theorem 2 (Identification of the Spillover Effects) . Suppose that Assumptions 1, 2, and 3 hold. Then, given ρ and α , the spillover effects ξ c t = ( ξ 1 t , ξ 2 t , · · · , ξ Nt ) ⊤ ∈ R N on the N control units for t &gt; T 0 can be identified as follows:

$$\xi _ { t } ^ { c } = Y _ { t } ^ { c } - ( I _ { N } - \rho w \alpha ^ { \top } - \rho W ) ^ { - 1 } ( ( I _ { N } - \rho W ) Y _ { t } ^ { c } - \rho w Y _ { 0 t } ).$$

Proof. Note that Y c t = Y c t (1) for t &gt; T 0 and that ( I N -ρ wα ⊤ -ρ W ) -1 (( I N -ρ W Y ) c t -ρ w Y 0 t ) = Y c t (0). This leads to the result (6).

Theorems 1 and 2 show that by incorporating the SAR model into SCM, we can identify treatment and spillover effects without relying on SUTVA. The results in Theorems 1 and 2 also suggest that although the SAR model (3) includes the covariates X t and error term u t , it is unnecessary to estimate β or the distribution of u t for the estimation of treatment and spillover effects. Only ρ and the spatial weights ( w W , ) in the SAR model (3) are relevant for the estimation of treatment and spillover effects.

Remarkably, when there is no spillover (i.e., ρ = 0), our identification result (5) for ξ 0 t corresponds to the identification result of the standard SCM (Abadie et al., 2010). That is, when ρ = 0, the identification result (5) simplifies to

$$\xi _ { 0 t } & = Y _ { 0 t } - \alpha ^ { \top } ( I _ { N } - \rho w \alpha ^ { \top } - \rho W ) ^ { - 1 } ( ( I _ { N } - \rho W ) Y _ { t } ^ { c } - \rho w Y _ { 0 t } ) \\ & = Y _ { 0 t } - \alpha ^ { \top } I _ { N } ( Y _ { t } ^ { c } - 0 ) \\ & = Y _ { 0 t } - \alpha ^ { \top } Y _ { t } ^ { c },$$

where the final line is identical to the identification equation for the standard SCM (Abadie et al., 2010). This result indicates that the identification result obtained from the standard SCM is the special case of our identification result when ρ = 0 (no spillover).

## 4 Inference

In this section, we propose a Bayesian inference method for treatment effects ξ 0 t and spillover effects ξ it ( i = 1 2 , , . . . , N ), building on the identification results presented in Section 3. SCM is typically applied to data with small or no large pretreatment periods, where frequentist methods may not provide accurate estimations. On the other hand, Bayesian methods can provide precise statistical inference (e.g., credible intervals) through MCMC procedures even with short pretreatment periods. This section presents a Bayesian inference method for the model described in Sections 2 and 3.

## 4.1 Synthetic Weights

Regarding the estimation of the synthetic weights α , we adopt Kim et al.'s (2020) method, which utilizes the Bayesian horseshoe prior as the prior distribution for parameters. The Bayesian horseshoe prior places a hierarchical prior distribution on parameters, characterized by a strong shrinkage effect around zero. Park and Casella (2008) show that LASSO regression in linear models corresponds to estimation with a Laplace prior distribution on coefficients, while the Bayesian horseshoe prior proposed by Carvalho et al. (2010) suggests an even more selective coefficient choice. Thus, by placing this prior distribution on the synthetic weights, it is anticipated that the estimated weights will select the units in the control group that are more related to the treated unit.

The following hierarchical prior distributions are placed on α , α 1 2 , · · · , α N :

$$\alpha _ { i } \, | \, \lambda _ { i } \sim N ( 0, \lambda _ { i } ^ { 2 } ), \text{ \ for } i = 1, 2, \cdots, N, \\ \lambda _ { i } \, | \, \tau \sim Half-Cauchy ( 0, \tau ), \text{ \ for } i = 1, 2, \cdots, N, \\ \tau \, | \, \sigma _ { 1 } \sim Half-Cauchy ( 0, \sigma _ { 1 } ), \\ \sigma _ { 1 } \sim Half-Cauchy ( 0, 1 0 ).$$

When these prior distributions are set, as shown by Makalic and Schmidt (2015), the full conditionals of each parameter can be derived analytically, enabling parameter sampling using the Gibbs sampler. The Appendix shows the full conditionals for each parameter.

## 4.2 Spatial Autoregressive Panel Data Model

Regarding the model (3), we set the Bayesian horseshoe prior for β as in estimating the synthetic weights, that is,

$$\beta _ { k } \, | \, \kappa _ { k } \sim N ( 0, \kappa _ { i } ^ { 2 } ), \ \text{ for } k = 1, 2, \cdots, K, \\ \kappa _ { k } \, | \, \psi \sim \text{ Half-Cauchy} ( 0, \psi ), \ \text{ for } k = 1, 2, \cdots, K, \\ \psi \, | \, \sigma _ { 2 } \sim \text{ Half-Cauchy} ( 0, \sigma _ { 2 } ), \\ \sigma _ { 2 } \sim \text{ Half-Cauchy} ( 0, 1 0 ).$$

We assume that u t follows a multilevel latent factor model u t = ηγ t + e t , where γ t ∈ R p denote factors at time t , η ∈ R N × p are factor loadings, and e t ∈ R N is an error vector. We assume γ t follow the AR(1) model γ t = ϕ γ γ t -1 + v t , where v t ∼ N ( 0 p , σ 2 γ I p ). For each control unit i (= 1 , · · · , N ), η i follows p -dimensional normal distribution N p ( 0 p , σ 2 η Σ η ), where Σ η = diag( ω , 2 1 · · · , ω 2 p ). The priors of σ η and ω , 1 · · · , ω p are half-Cauchy(0 10). , These settings are the same as those used by Pang et al. (2022).

We adopt the method proposed by LeSage (1997) and LeSage and Pace (2009) for modeling ρ . Given all parameters except ρ , the conditional distribution of ρ is

$$p ( \rho \, | \, \text{rest} ) \in \prod _ { t = 1 } ^ { T _ { 0 } } | I _ { N } - \rho W | \exp \left ( - \frac { 1 } { 2 \sigma ^ { 2 } } \tilde { u } _ { t } ^ { \top } \tilde { u } _ { t } \right ),$$

where ˜ u t = AY c t -ρ w Y 0 t -X β t -γ ⊤ t η and A = I N -ρ W . However, since this is not a standard form, it is impossible to analytically derive the full conditional distribution. Therefore, sampling is conducted using the Metropolis algorithm. If the value of ρ in the m -th sampling is ρ ( m ) , then the proposal value ρ ∗ for the ( m +1)-th sampling is determined by ρ ∗ = ρ ( m ) + k · N (0 1), where , k is a tuning parameter adjusted to achieve an acceptance rate of approximately 40% to 60%. The MCMC procedure for all parameters is described in the Appendix.

Using the sampled parameters from the posterior distributions, we can estimate the treatment and spillover effects through the identification results (5) and (6) in Theorems 1 and 2. Specifically, if M samplings are performed, the estimated values of the treatment and spillover effects for the m -th sample are as follows:

$$\xi _ { 0 t } ^ { ( m ) } = Y _ { 0 t } - \alpha ^ { \top ( m ) } ( I _ { N } - \rho ^ { ( m ) } w \alpha ^ { \top ( m ) } - \rho ^ { ( m ) } W ) ^ { - 1 } ( ( I _ { N } - \rho ^ { ( m ) } W ) Y _ { t } ^ { c } - \rho ^ { ( m ) } w Y _ { 0 t } ),$$

$$\xi _ { t } ^ { ( m ) } = Y _ { t } ^ { c } - ( I _ { N } - \rho ^ { ( m ) } w \alpha ^ { \top ( m ) } - \rho ^ { ( m ) } W ) ^ { - 1 } ( ( I _ { N } - \rho ^ { ( m ) } W ) Y _ { t } ^ { c } - \rho ^ { ( m ) } w Y _ { 0 t } ).$$

Calculating this for each m (= 1 2 , , · · · , M ) allows us to construct the distribution of the treatment and spillover effects. Then the posterior means of ξ 0 t and ξ t can be obtained by (1 /M ) ∑ M m =1 ξ ( m ) 0 t and (1 /M ) ∑ M m =1 ξ ( m ) t , respectively. Given the estimators of α and ρ , the estimation of the treatment effect ξ 0 t and spillover effect ξ t depends only on the observed outcomes Y 0 t and Y c t over time. This suggests that time-varying components in the SAR panel data model, such as factors γ t , are not required to be estimated in the post-treatment periods t ( &gt; T 0 ).

## 5 Simulation Study

## 5.1 Simulation Design

We conduct a simulation study to examine the finite sample performance of the proposed Bayesian inference method. Three scenarios are considered for the number of untreated units N : N = 16, 36, and 64. For the length of time periods, we consider two scenarios: ( T, T 0 ) = (30 20) and (60 50), , , where each scenario has a post-treatment length of ten.

In each scenario of ( N,T,T 0 ), we set up the DGPs as follows. We consider a network of N untreated units represented by a rook matrix. 2 Subsequently, we set the spatial weight ω ij of any two untreated units to be 1 if units i and j are connected in the network and 0 otherwise. Each element of the adjacency vector w takes the value of 1 if i ∈ { 1 2 3 4 , , , } and 0 otherwise. The synthetic weights α are set to provide large weights only to control units adjacent to the treatment unit, as follows:

$$\alpha _ { i } = \begin{cases} 0. 5 & \text{ if } i = 1, \\ - 0. 2 & \text{ if } i = 2, \\ 0. 4 & \text{ if } i \in \{ 3, 4 \}, \\ 0. 1 / 6 & \text{ if } i \in \{ 5, 6, \dots, 10 \}, \\ 0 & \text{ otherwise.} \end{cases}$$

For each t (= 1 2 , , · · · , T ), the control outcomes of untreated units Y c t (0) are distributed from

2 We use a rook matrix based on an r board (so that N = r 2 ) to represent the network of N untreated units. The rook matrix represents a square tessellation with connectivity of four for the inner fields on the chessboard, and two and three for the corner and border fields, respectively.

a SAR panel data model as follows: Y c t (0) = ρ w Y 0 t (0) + ρ WY c t (0) + X β t + u t , where X t = ( X ,...,X 1 t NT ) ⊤ with each element being i.i.d as N (0 , 1), u t = ( u 1 t , . . . , u NT ) ⊤ with each element being i.i.d as N (0 , 1), and β = 1 0. . The control outcome of the treated unit Y 0 t (0) is generated as Y 0 t (0) = ∑ N i =1 α Y i i,t (0). We consider seven scenarios of ρ : ρ = -0 8 . , -0 3 . , -0 1 . , 0 0 . , 0 1 . , 0 3 . , 0 8. . A larger absolute value of ρ implies a stronger spatial correlation among the units.

We generate the treatment outcomes of the treated unit Y 0 t (1) as Y 0 t (1) = Y 0 t (0) + N (1 , 1). We also set Y c t (1) = ρ w Y 0 t (1) + ρ WY c t (1) + X β t + u t . The observed outcome for each i and t is Y it = Y it (0) · 1 { t ≤ T 0 } + Y it (1) · 1 { t &gt; T 0 } .

Following the DGPs, we conduct 1000 Monte Carlo simulations. For each simulation r (= 1 , · · · , 1000), we draw M (= 5000) samples by MCMC and compute bias and root mean squared error (RMSE) of the posterior mean of treatment effect as follows:

$$B i a s & = \frac { 1 } { 1 0 0 0 } \sum _ { r = 1 } ^ { 1 0 0 0 } \frac { 1 } { T - T _ { 0 } } \sum _ { t = T _ { 0 } + 1 } ^ { T } \left ( \xi _ { 0 t } ^ { ( r ) } - \widehat { \xi } _ { 0 t } ^ { ( r ) } \right ), \\ R M S E & = \sqrt { \frac { 1 } { 1 0 0 0 } \sum _ { r = 1 } ^ { 1 0 0 0 } \frac { 1 } { T - T _ { 0 } } \sum _ { t = T _ { 0 } + 1 } ^ { T } \left ( \xi ^ { ( r ) } - \widehat { \xi } _ { 0 t } ^ { ( r ) } \right ) ^ { 2 } },$$

$$\begin{matrix} 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 \\ \end{matrix}$$

where ξ ( r ) 0 t denotes a true treatment effect in the r -th simulation, ̂ ξ ( r ) 0 t = ∑ M m =1 ξ ( r,m ) 0 t /M is the estimated treatment effect with ξ ( r,m ) 0 t being the estimate of the treatment effect in the m -th iteration for the r -th simulation.

We compare the performance of the proposed method (labeled 'Proposed') with those of the standard SCM (labeled 'SCM') (Abadie et al., 2010) and the Bayesian SCM of Kim et al. (2020) (labeled 'BSCM'). 3 We also calculate the 95% coverage rate of treatment effect for the proposed method.

## 5.2 Results

Table 1 presents the simulation results for bias and RMSE. The key finding is that the error in estimating the treatment effect using SCM and BSCM becomes large as the absolute value of ρ increases. Particularly in the case of a strong positive spatial correlation ( ρ = 0 8), SCM and BSCM . exhibit substantial biases and large RMSEs. These results indicate that, when spillovers exist, the treatment effects estimated by SCM and BSCM are biased. Conversely, the proposed method exhibits a small bias and RMSE in each simulation scenario. The bias and RMSE of the proposed

3 SCM does not involve MCMC sampling.

method do not increase with the magnitude of the spatial correlation. These results indicate that the proposed method is robust to the spillover effects arising from the spatial correlation of the outcomes.

$$[ Table 1 about here.]$$

Table 2 presents the coverage rate of 95% credible interval of the treatment effect for the proposed method. In each scenario, the coverage rate is close to 95%, indicating that the proposed inference method performs adequately. The inference performs well even when the length of pretreatment periods is not very large ( T 0 = 20) and/or the spatial correlation is strong (i.e., | ρ | is large).

$$[ Table 2 about here.]$$

## 6 Empirical Application

We conduct two empirical studies applying the proposed method: one estimates the impact of the California tobacco tax on consumption (Abadie et al., 2010), and the other estimates the economic impact of the 2011 division of Sudan on GDP per capita. This section presents the results of the two empirical studies.

## 6.1 Application I: California Tobacco Tax

In this section, we apply the proposed method to estimate the effect of California cigarette tax (Proposition 99) on consumption (Abadie et al., 2010). Proposition 99 is an anti-tobacco law issued in California in 1988 to promote awareness of the health risks associated with tobacco. It increased the tobacco excise tax on cigarettes in California by 25 cents per pack. Abadie et al. (2010) estimate the treatment effect of Proposition 99 on cigarette sales by SCM without accounting for spillovers. Applying the proposed method, we estimate both the treatment and spillover effects of Proposition 99 while allowing for spillovers.

## 6.1.1 Data

We use annual cigarette sales data (per state) for the U.S. from 1970 to 2000, as used in Abadie et al. (2010). 4 The outcome of interest is state-level annual per capita cigarette consumption. We

4 This dataset is available from Cunningham (2021) on GitHub. It has undergone preprocessing as described in Abadie et al. (2010).

compare the proposed method with the standard SCM (Abadie et al., 2010) (labeled as 'SCM'), which is not robust to the existence of spillover effects. For SCM, we use the synthetic weights estimated by Abadie et al. (2010). 5

In the estimation using the proposed method, the average retail cigarette prices for each year t are included as covariates X t . As for the spatial weights w ij in the SAR model (3), we use the adjacency weights of the states in the dataset. We also normalize the spatial weights matrix ( w W , ) such that each row of ( w W , ) sums to one.

## 6.1.2 Results

Figure 1 shows the estimation results for the synthetic control outcomes for California and the treatment effects of the California tobacco tax on consumption, where the posterior means for the proposed method are shown. Figure 1 (a) illustrates that both the SCM and the proposed method fit well with per-capita cigarette sales in California during the pretreatment periods. As shown in Figure 1 (b), each estimation method suggests that the California cigarette tax (Proposition 99) decreased cigarette consumption in California. The proposed method exhibits higher estimates of treatment effects compared to the SCM. The 90% credible interval for the proposed method is negative overall for any year following the intervention, which began in 1988. This finding suggests that Proposition 99 negatively impacted cigarette sales for a decade, supporting the findings of Abadie et al. (2010).

[Figure 1 about here]

Figure 2 shows the estimates of spillover effects for all states in the control group. The results suggest that the California tobacco tax reduced tobacco consumption in many other states. The particularly affected states are Louisiana, Texas, Mississippi, and Oklahoma, all of which are located in the Southern U.S.

[Figure 2 about here]

## 6.2 Application II: The Economic Cost of the 2011 Sudan Split

In this section, we assess the impact of Sudan's north-south split in 2011 on GDP per capita in the Sudans (the region of the former united Sudan) and other African countries.

5 The estimated values are listed in Table 2 of Abadie et al. (2010).

## 6.2.1 Background

Sudan has long been divided along ethnic and religious lines, with Arabs (primarily Muslims) predominantly in the north and Africans (primarily Christians) in the south, leading to many conflicts. Particularly, the Darfur conflict, driven by the Arab versus non-Arab ethnic divide, has persisted for many years in western Sudan. This conflict has been marked by large-scale atrocities, including mass killings carried out by Arab militias known as 'Janjaweed.'

Amid these unending conflicts, the Sudanese government and the Sudan People's Liberation Army (SPLA), the main rebel force in Southern Sudan, signed the Comprehensive Peace Agreement (CPA) in 2005. This agreement, aimed at ending Sudan's civil wars, allowed South Sudan to establish its own government, achieve autonomy, and pursue independence through a referendum. Following the CPA, South Sudan voted for independence in January 2011 and was officially recognized as a nation on July 9, 2011.

Since 2011, South Sudan's independence has caused several economic disturbances in both Sudan and South Sudan. In particular, South Sudan's oil production shutdown in 2012 and its relapse into conflict in 2013 provoked a severe macroeconomic crisis in the region (Mawejje, 2020). This study estimates the economic impact of Sudan's south-north split, with a focus on GDP per capita in the Sudans (the region of the former united Sudan) and other African countries. 6

## 6.2.2 Data

We use data from African countries obtained from the World Bank DataBank. Our outcome of interest is 'GDP per capita (constant 2015 US $ )' post-division. The covariates we use include 'exports of goods and services (% of GDP)', 'merchandise trade (% of GDP)', 'access to electricity (% of population)', 'inflation measured by the consumer price index (annual %)', 'net migration', and 'trade (% of GDP)'. Countries with missing values for the outcome or covariates were excluded. We focus on GDP per capita in the Sudans (the region of the former united Sudan) as the outcome Y 0 t of the treated unit. 7 The control group comprises N = 29 African countries with complete

6 Using the standard SCM, Mawejje and McSharry (2021) estimate the economic losses in South Sudan, owing to the oil production halt in 2012, finding a nearly 70% loss in per capita real GDP from 2012 to 2018. They excluded neighboring countries of South Sudan from their SCM analysis to avoid bias caused by spillovers; however, this practice can result in a poorly fitting synthetic control, making the perfect-fit assumption (Assumption 1) less plausible and potentially causing additional bias.

7 The GDP per capita in the Sudans after the Sudan split is calculated by dividing the sum of GDP in North and South Sudan by the sum of their populations. Before the split, this measure corresponds to the GDP per capita in (the united) Sudan.

observations from 2000 to 2015. 8 Since South Sudan's independence occurred in July 2011, we set pre-treatment periods to be 2000 to 2010 and post-treatment periods to be 2011 to 2015. We cannot compute GDP per capita for the Sudans in 2011 because of incomplete data caused by the Sudan split; however, this does not affect the estimation of synthetic control outcomes in the pre-treatment periods or the estimation of treatment effects from 2012 onward.

Regarding the spatial weights in the model (3), we specify w ij as the average amount of international trade between countries i and j as follows:

$$w _ { i j } = \frac { \text{the average amount of trade between countries $i$ and $j$} } { \sum _ { j } \text{the average amount of trade between countries $i$ and $j$} },$$

where the trade data come from the IMF and the average amounts of trade are computed from the data of the pre-intervention periods. Each weight w ij , measured by the amount of trade, reflects the strength of the economic connections between the two countries. Countries with stronger economic connections to the former united Sudan are expected to experience stronger spillover effects from the Sudan split.

## 6.2.3 Results

Figure 3 illustrates the estimation results for the synthetic outcomes of the Sudans and the treatment effects of the Sudan split. For each method, the estimates indicate a negative impact of southern independence on GDP per capita in the Sudans. In particular, the estimation results from the proposed method show that GDP per capita decreased by about 100 USD in 2012, which represents approximately 7 8% of the GDP per capita for that year. . In 2015, the estimated GDP per capita in the synthetic Sudan is about 9 5% higher than the actual GDP per capita. . In addition, the cumulative losses in the Sudans are estimated to have reached 34% from 2011 to 2015. 9

$$[ F i g u r e 3 \text{ about here} ]$$

Figure 4 shows the estimated spillover effects of the Sudan split on other African countries. Countries with significant trade volumes with the former united Sudan, such as Egypt and Kenya, experienced substantial negative spillover effects from the split.

8 The list of countries used in the control group is as follows: Algeria, Angola, Benin, Botswana, Burundi, Cameroon, Central African Republic, Chad, Egypt, Gabon, Ghana, Ivory Coast, Kenya, Madagascar, Mali, Mauritania, Mauritius, Morocco, Niger, Nigeria, Republic of the Congo, Rwanda, Senegal, South Africa, Tanzania, Togo, Tunisia, Uganda, Zambia.

9 The losses in the Sudans are defined by 100 × ( ξ 0 t /Y 0 t (0)) (%) for each t &gt; T 0 and those in control countries are defined by 100 × ξ it /Y it (0) (%) for each i = 1 2 , , · · · , N and t &gt; T 0 . The cumulative losses are computed by 100 × ∑ 2015 t =2011 ξ it /Y it (0) for each i = 0 1 , , · · · , 29. We estimate these by the posterior means of these losses.

The political instability and north-south split in Sudan significantly altered its industrial structure, which, in turn, profoundly impacted trade. For example, while South Sudan is rich in oil and other natural resources, the conflict between the North and South Sudans halted oil production in South Sudan in 2012, resulting in the loss of crucial export commodities. This disruption extended to countries with close economic ties through trade channels. Thus, the north-south split in Sudan inflicted negative economic impacts on the Sudans themselves and caused negative spillover effects on other countries. This empirical study illustrates that political and economic changes in one country can potentially affect other countries with close economic ties.

[Figure 4 about here]

## 7 Conclusion

This study extends SCM to allow for spillover effects. While SCM is often applied to spatial data that may involve spillovers, conventional SCM relies on SUTVA and can produce biased estimates of treatment effects when spillovers are present. We propose a novel SCM that leverages the SAR panel data model to address these spillover effects. We also introduce a Bayesian inference method for estimating both treatment and spillover effects, while employing horseshoe priors for regularization. We apply the method to two empirical studies: evaluating the impact of the California tobacco tax on its consumption (Abadie et al., 2010) and assessing the economic impact of the 2011 Sudan division on GDP per capita. The first study demonstrates the negative impact of the tax on tobacco consumption in California and other US states, while the second study shows that the Sudan split reduced GDP per capita in the Sudans and negatively affected other African countries with strong economic connections to the former united Sudan.

## Acknowledgments

We thank Ryo Okui, Yasuyuki Sawada, Kaoru Irie, and participants in various seminars and workshops for valuable comments. The authors gratefully acknowledge the financial support from JSPS KAKENHI Grant (number 24K16342).

## References

Abadie, Alberto. 2021. 'Using synthetic controls: Feasibility, data requirements, and methodological aspects.' Journal of Economic Literature 59 (2): 391-425.

Abadie, Alberto, Alexis Diamond, and Jens Hainmueller. 2010. 'Synthetic control methods for comparative case studies: Estimating the effect of California's tobacco control program.' Journal of the American Statistical Association 105 (490): 493-505.

Abadie, Alberto, Alexis Diamond, and Jens Hainmueller. 2015. 'Comparative politics and the synthetic control method.' American Journal of Political Science 59 (2): 495-510.

Abadie, Alberto, and Javier Gardeazabal. 2003. 'The economic costs of conflict: A case study of the Basque Country.' American Economic Review 93 (1): 113-132.

Abadie, Alberto, and J´r´my L'hour. e e 2021. 'A penalized synthetic control estimator for disaggregated data.' Journal of the American Statistical Association 116 (536): 1817-1834.

Arkhangelsky, Dmitry, Susan Athey, David A Hirshberg, Guido W Imbens, and Stefan Wager. 2021. 'Synthetic difference-in-differences.' American Economic Review 111 (12): 40884118.

Athey, Susan, and Guido W Imbens. 2017. 'The state of applied econometrics: Causality and policy evaluation.' Journal of Economic Perspectives 31 (2): 3-32.

Ben-Michael, Eli, Avi Feller, and Jesse Rothstein. 2021. 'The augmented synthetic control method.' Journal of the American Statistical Association 116 (536): 1789-1803.

Bifulco, Robert, Ross Rubenstein, and Hosung Sohn. 2017. 'Using synthetic controls to evaluate the effect of unique interventions: The case of say yes to education.' Evaluation Review 41 (6): 593-619.

Brodersen, Kay H., Fabian Gallusser, Jim Koehler, Nicolas Remy, and Steven L. Scott. 2015. 'Inferring causal impact using Bayesian structural time-series models.' The Annals of Applied Statistics 9 (1): 247 - 274.

Cao, Jianfei, and Connor Dowd. 2019. 'Estimation and inference for synthetic control methods with spillover effects.' arXiv preprint arXiv:1902.07343 .

Carvalho, Carlos M, Nicholas G Polson, and James G Scott. 2010. 'The horseshoe estimator for sparse signals.' Biometrika 97 (2): 465-480.

Chernozhukov, Victor, Kaspar W¨thrich, and Yinchu Zhu. u 2021. 'An exact and robust conformal inference method for counterfactual and synthetic controls.' Journal of the American Statistical Association 116 (536): 1849-1864.

Cunningham, Scott. 2021. Causal Inference: The Mixtape . Yale University Press.

Ferman, Bruno, and Cristine Pinto. 2021. 'Synthetic controls with imperfect pretreatment fit.' Quantitative Economics 12 (4): 1197-1221.

Fingleton, Bernard. 2008. 'A generalized method of moments estimator for a spatial panel model with an endogenous spatial lag and spatial moving average errors.' Spatial Economic Analysis 3 (1): 27-44.

Glass, Anthony J, Karligash Kenjegalieva, and Robin C Sickles. 2016. 'A spatial autoregressive stochastic frontier model for panel data with asymmetric efficiency spillovers.' Journal of Econometrics 190 (2): 289-300.

Grossi, Giulio, Patrizia Lattarulo, Marco Mariani, Alessandra Mattei, and O Oner. 2020. 'Synthetic control group methods in the presence of interference: The direct and spillover effects of light rail on neighborhood retail activity.' arXiv preprint arXiv:2004.05027 .

Kim, Sungjin, Clarence Lee, and Sachin Gupta. 2020. 'Bayesian synthetic control methods.' Journal of Marketing Research 57 (5): 831-852.

Klinenberg, Danny. 2023. 'Synthetic control with time varying coefficients: A state space approach with Bayesian shrinkage.' Journal of Business &amp; Economic Statistics 41 (4): 1065-1076.

Lee, Lung-fei, and Jihai Yu. 2010. 'Estimation of spatial autoregressive panel data models with fixed effects.' Journal of Econometrics 154 (2): 165-185.

LeSage, James P. 1997. 'Bayesian estimation of spatial autoregressive models.' International Regional Science Review 20 (1-2): 113-129.

LeSage, James, and Robert Kelley Pace. 2009. Introduction to Spatial Econometrics . Chapman and Hall/CRC.

Li, Kathleen T. 2020. 'Statistical inference for average treatment effects estimated by synthetic control methods.' Journal of the American Statistical Association 115 (532): 2068-2083.

Liang, Xuan, Jiti Gao, and Xiaodong Gong. 2022. 'Semiparametric spatial autoregressive panel data model with fixed effects and time-varying coefficients.' Journal of Business &amp; Economic Statistics 40 (4): 1784-1802.

Makalic, Enes, and Daniel F Schmidt. 2015. 'A simple sampler for the horseshoe estimator.' IEEE Signal Processing Letters 23 (1): 179-182.

Mawejje, Joseph. 2020. The Macroeconomic Environment for Jobs in South Sudan: Jobs, Recovery, and Peacebuilding in Urban South Sudan-Technical Report II . World Bank.

Mawejje, Joseph, and Patrick McSharry. 2021. 'The economic cost of conflict: Evidence from South Sudan.' Review of Development Economics 25 (4): 1969-1990.

Menchetti, Fiammetta, and Iavor Bojinov. 2022. 'Estimating the effectiveness of permanent price reductions for competing products using multivariate Bayesian structural time series models.' The Annals of Applied Statistics 16 (1): 414-435.

Miguel, Edward, and Michael Kremer. 2004. 'Worms: identifying impacts on education and health in the presence of treatment externalities.' Econometrica 72 (1): 159-217.

Pang, Xun, Licheng Liu, and Yiqing Xu. 2022. 'A Bayesian alternative to synthetic control for comparative case studies.' Political Analysis 30 (2): 269-288.

Park, Trevor, and George Casella. 2008. 'The Bayesian lasso.' Journal of the American Statistical Association 103 (482): 681-686.

Rubin, Donald B. 1978. 'Bayesian inference for causal effects: The role of randomization.' The Annals of Statistics 34-58.

Su, Liangjun. 2012. 'Semiparametric GMM estimation of spatial autoregressive models.' Journal of Econometrics 167 (2): 543-560.

Wand, Matthew P., John T. Ormerod, Simone A. Padoan, and Rudolf Fr¨hwirth. u 2011. 'Mean field variational Bayes for elaborate distributions.' Bayesian Analysis 6 (4): 847 - 900.

## Tables

Table 1: Simulation Results for Bias and RMSE

Panel (a): T 0 = 20

|      |      | N = 16   | N = 16   | N = 16   | N = 36   | N = 36   | N = 36   | N = 64   | N = 64   | N = 64   |
|------|------|----------|----------|----------|----------|----------|----------|----------|----------|----------|
|      | ρ    | Proposed | SCM      | BSCM     | Proposed | SCM      | BSCM     | Proposed | SCM      | BSCM     |
|      | -0.8 | 0.003    | 0.875    | 1.022    | 0.002    | 0.904    | 1.422    | 0.002    | 0.752    | 1.410    |
|      | -0.3 | 0.004    | 0.404    | 0.421    | 0.000    | 0.366    | 0.490    | -0.002   | 0.307    | 0.484    |
|      | -0.1 | 0.007    | 0.144    | 0.150    | 0.003    | 0.113    | 0.166    | 0.003    | 0.101    | 0.165    |
| Bias | 0.0  | 0.009    | -0.006   | 0.000    | 0.004    | -0.004   | 0.000    | 0.004    | -0.010   | 0.000    |
|      | 0.1  | 0.011    | -0.156   | -0.163   | 0.003    | -0.145   | -0.173   | 0.002    | -0.139   | -0.171   |
|      | 0.3  | 0.011    | -0.557   | -0.550   | 0.009    | -0.482   | -0.556   | 0.002    | -0.449   | -0.549   |
|      | 0.8  | 0.061    | -2.642   | -2.613   | 0.018    | -2.051   | -2.109   | 0.007    | -2.002   | -2.100   |
|      | -0.8 | 0.350    | 1.274    | 1.043    | 0.053    | 1.354    | 1.449    | 0.057    | 1.330    | 1.440    |
|      | -0.3 | 0.067    | 0.871    | 0.429    | 0.066    | 0.923    | 0.501    | 0.064    | 0.987    | 0.497    |
|      | -0.1 | 0.085    | 0.738    | 0.153    | 0.069    | 0.806    | 0.174    | 0.075    | 0.869    | 0.179    |
| RMSE | 0.0  | 0.091    | 0.713    | 0.000    | 0.076    | 0.783    | 0.038    | 0.070    | 0.832    | 0.052    |
|      | 0.1  | 0.098    | 0.714    | 0.166    | 0.079    | 0.758    | 0.180    | 0.075    | 0.822    | 0.183    |
|      | 0.3  | 0.132    | 0.878    | 0.561    | 0.094    | 0.879    | 0.568    | 0.086    | 0.895    | 0.562    |
|      | 0.8  | 0.350    | 2.777    | 2.664    | 0.180    | 2.188    | 2.150    | 0.167    | 2.150    | 2.142    |

Panel (b): T 0 = 50

|      |      | N = 16   | N = 16   | N = 16   | N = 36   | N = 36   | N = 36   | N = 64   | N = 64   | N = 64   |
|------|------|----------|----------|----------|----------|----------|----------|----------|----------|----------|
|      | ρ    | Proposed | SCM      | BSCM     | Proposed | SCM      | BSCM     | Proposed | SCM      | BSCM     |
|      | -0.8 | -0.002   | 0.962    | 1.023    | -0.002   | 1.109    | 1.456    | -0.001   | 1.036    | 1.461    |
|      | -0.3 | 0.001    | 0.443    | 0.423    | -0.001   | 0.444    | 0.500    | -0.001   | 0.427    | 0.499    |
|      | -0.1 | 0.002    | 0.158    | 0.149    | 0.001    | 0.152    | 0.170    | 0.001    | 0.146    | 0.170    |
| Bias | 0.0  | 0.002    | 0.002    | 0.000    | 0.003    | 0.000    | 0.000    | -0.001   | 0.001    | 0.000    |
|      | 0.1  | 0.009    | -0.175   | -0.162   | 0.004    | -0.163   | -0.175   | 0.002    | -0.163   | -0.176   |
|      | 0.3  | 0.008    | -0.578   | -0.551   | 0.003    | -0.531   | -0.559   | 0.002    | -0.526   | -0.558   |
|      | 0.8  | 0.041    | -2.665   | -2.611   | 0.015    | -2.116   | -2.118   | 0.009    | -2.111   | -2.118   |
|      | -0.8 | 0.021    | 1.274    | 1.044    | 0.026    | 1.401    | 1.485    | 0.020    | 1.358    | 1.491    |
|      | -0.3 | 0.041    | 0.826    | 0.431    | 0.035    | 0.837    | 0.510    | 0.026    | 0.838    | 0.509    |
|      | -0.1 | 0.052    | 0.671    | 0.152    | 0.039    | 0.692    | 0.173    | 0.028    | 0.693    | 0.173    |
| RMSE | 0.0  | 0.057    | 0.638    | 0.000    | 0.041    | 0.648    | 0.000    | 0.028    | 0.665    | 0.000    |
|      | 0.1  | 0.062    | 0.655    | 0.166    | 0.044    | 0.671    | 0.179    | 0.033    | 0.670    | 0.179    |
|      | 0.3  | 0.081    | 0.854    | 0.562    | 0.052    | 0.821    | 0.570    | 0.038    | 0.829    | 0.569    |
|      | 0.8  | 0.215    | 2.777    | 2.663    | 0.112    | 2.233    | 2.160    | 0.086    | 2.230    | 2.160    |

Notes: Panels (a) and (b) show the simulation results for the bias and RMSE for T 0 = 20 and 50, respectively. Each panel shows the simulation results for each of the proposed method, SCM, and BSCM, and each of ρ ∈ {-0 8 . , -0 3 . , -0 1 . , 0 0 . , 0 1 . , 0 3 . , 0 8 . } and N ∈ { 16 36 64 , , } .

Table 2: Coverage Rate of 95% Credible Interval for the proposed method

|      | T 0 = 20   | T 0 = 20   | T 0 = 20   | T 0 = 50   | T 0 = 50   | T 0 = 50   |
|------|------------|------------|------------|------------|------------|------------|
| ρ    | N = 16     | N = 36     | N = 64     | N = 16     | N = 36     | N = 64     |
| -0.8 | 0.955      | 0.963      | 0.961      | 0.960      | 0.951      | 0.939      |
| -0.3 | 0.936      | 0.962      | 0.967      | 0.956      | 0.943      | 0.950      |
| -0.1 | 0.938      | 0.975      | 0.975      | 0.947      | 0.951      | 0.951      |
| 0.0  | 0.951      | 0.970      | 0.975      | 0.950      | 0.943      | 0.968      |
| 0.1  | 0.954      | 0.968      | 0.975      | 0.952      | 0.950      | 0.944      |
| 0.3  | 0.949      | 0.965      | 0.975      | 0.961      | 0.953      | 0.943      |
| 0.8  | 0.947      | 0.963      | 0.979      | 0.970      | 0.941      | 0.952      |

Notes: This table shows the coverage rate of 95% credible interval for the proposed method for each scenario of ρ , T 0 , and N . The coverage rate is computed over 1000 simulations.

## Figures

Figure 1: Estimates of the Counterfactual Outcomes and Treatment Effects of California Tobacco Tax

## (a) Observed Outcomes and Synthetic Control Outcomes

## (b) Treatment Effect Estimation

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00291v2-figure-78a718dc-5ec4-495b-a63a-b0197fbe02cb.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00291v2-figure-2987eeb4-bf64-42d7-a42c-2c3d85826c9c.jpg)

Note: Panel (a) shows the observed outcomes (per-capita cigarette sales in packs) for California (black dotted line) and estimates of the counterfactual control outcomes for California using the proposed method (red solid line) and SCM (blue dashed line). Panel (b) shows the estimates of the treatment effects of Proposition 99 using the proposed method (red solid line) and the SCM (black dashed line). Each panel also illustrates a 95% credible interval by the proposed method. Post-mean of ρ is 0 477 with a 95% posterior credible . interval of [0 427 . , 0 528]. .

Figure 2: Estimates of the Spillover Effects of California Tobacco Tax

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00291v2-figure-0577a53d-7a6f-4870-9370-4a69a0a6c45a.jpg)

Notes: This figure shows the spillover effects of Proposition 99 estimated by the proposed method for each state in the control group. For ease of reading, only three states with particularly strong spillover effects are indicated.

Figure 3: Estimates of the Counterfactual Outcomes and Treatment Effects of Sudan Split

- (a) Observed Outcomes and Synthetic Control Outcomes

## (b) Estimates of Treatment Effects

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00291v2-figure-b8d251d5-2e71-4fe0-8b45-4a2644d14fc3.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00291v2-figure-188f320b-b2dd-4707-85de-eae6e6d34813.jpg)

Notes: Panel (a) shows the observed outcomes (GDP per capita) of the Sudans (black dotted-dashed line) and estimates of the counterfactual control outcomes for the Sudans using the proposed method (red solid line) and the SCM (blue dashed line). Panel (b) shows the estimates of the treatment effects of the Sudan split using the proposed method (red solid line) and the SCM (blue dashed line). In each panel, missing values are illustrated with dotted lines. Each panel also illustrates a 95% credible interval by the proposed method. Post-mean of ρ is about 0 1875 (95% posterior credible interval is [0 1352 . . , 0 2401]). .

Figure 4: Estimates of the Spillover Effects of Sudan Split

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00291v2-figure-bec72b4a-3bb9-4630-be57-de15ef4f6830.jpg)

Notes: This figure shows the spillover effects of the Sudan split, estimated by the proposed method for each country in the control group. Dotted lines illustrate missing values. To ensure readability, only four countries with strong spillover effects are indicated.

## Appendix

To derive full conditional distributions, we use the following proposition:

Proposition 1 (Wand et al. (2011)) .

$$X \ | \ a \sim \text{Half-Cauchy} ( 0, a ) \Longleftrightarrow \begin{cases} X ^ { 2 } \ | \ b \sim \text{Inverse-Gamma} \begin{pmatrix} 1 & 1 \\ 2 \end{pmatrix} \\ b \ | \ a \sim \text{Inverse-Gamma} \begin{pmatrix} 1 & 1 \\ 2 \end{pmatrix} \end{cases}$$

This proposition implies that a half-Cauchy distribution is equivalent to a hierarchical form of an inverse gamma distribution with an auxiliary variable.

## A Synthetic Weights

By using auxiliary variables, we can derive the full conditional distributions as follows:

$$\begin{array} { l l l } \text{By using auxiliary variables, we can derive the full conditional distributions as follows:} \\ \\ \alpha | rest \sim \mathcal { N } _ { N } ( A ^ { - 1 } ( Y ^ { c } ( 0 ) ) ^ { T } Y _ { 0 } ( 0 ), \, \sigma _ { 1 } ^ { 2 } A ^ { - 1 } ), \\ \text{where } A = ( Y ^ { c } ( 0 ) ) ^ { T } Y ^ { c } ( 0 ) + \sigma _ { 1 } ^ { 2 } \text{diag} ( 1 / \lambda _ { 1 } ^ { 2 }, 1 / \lambda _ { 2 } ^ { 2 }, \cdots, 1 / \lambda _ { N } ^ { 2 } ) \in \mathbb { R } ^ { N \times N } \\ \\ \lambda _ { 1 } ^ { 2 } | rest \sim \text{Inverse-Gamma} ( 1, \, \sum _ { i = 1 } ^ { N } \frac { 1 } { \lambda _ { 1 } ^ { 2 } } + \frac { 1 } { \nu _ { \lambda } } ), \, \text{for } i = 1, 2, \cdots, N \\ \\ \nu _ { \lambda } | rest \sim \text{Inverse-Gamma} ( 1, \, \sum _ { i = 1 } ^ { N } \frac { 1 } { \lambda _ { 1 } ^ { 2 } } + \frac { 1 } { \nu _ { \nu } } ) \\ \\ \tau ^ { 2 } | rest \sim \text{Inverse-Gamma} ( 1, \, \frac { 1 } { \nu _ { \lambda } } + \frac { 1 } { \nu _ { \nu } } ) \\ \\ \nu _ { \tau } | rest \sim \text{Inverse-Gamma} ( 1, \, \frac { 1 } { \tau ^ { 2 } } + \frac { 1 } { \tau ^ { 2 } } ) \\ \\ \sigma _ { 1 } ^ { 2 } | rest \sim \text{Inverse-Gamma} ( 1 + \tau _ { 2 } ^ { 2 }, \, \frac { 1 } { \nu _ { \tau } } + \frac { 1 } { \nu _ { \sigma _ { 1 } } } + \frac { 1 } { ( Y _ { 0 } ( 0 ) - \alpha ^ { T } Y ^ { c } ( 0 ) ) ^ { T } ( Y _ { 0 } ( 0 ) - \alpha ^ { T } Y ^ { c } ( 0 ) ) ) \\ \\ \nu _ { \sigma _ { 1 } } | rest \sim \text{Inverse-Gamma} ( 1, \, \frac { 1 } { \sigma _ { 1 } ^ { 2 } } + \frac { 1 } { 1 0 ^ { 2 } } ) \\ \\ \text{B} & \text{Spatial Autoregressive Panel Data Model} \end{array}$$

## B Spatial Autoregressive Panel Data Model

For the remaining parameters, we specify priors as follows.

- 1. Prior of β :

$$\beta \, | \, \text{rest} \sim \mathcal { N } _ { k } ( A _ { \beta } ^ { - 1 } \tilde { u }, \, \sigma _ { 2 } ^ { 2 } A _ { \beta } ^ { - 1 } )$$

$$\text{where}$$

$$\text{ where } \ A _ { \beta } & = X ^ { \top } X + \sigma _ { 2 } ^ { 2 } \text{diag} ( 1 / \lambda _ { \beta _ { 1 } } ^ { 2 }, 1 / \lambda _ { \beta _ { 2 } } ^ { 2 }, \cdots, 1 / \lambda _ { \beta _ { k } } ^ { 2 } ) \\ \lambda _ { \beta _ { j } } ^ { 2 } \ | \, \text{rest} \sim \text{Inverse-Gamma} ( 1, \, \frac { \beta _ { j } ^ { 2 } } { 2 } + \frac { 1 } { \nu _ { \lambda _ { \beta _ { 0 } } } } ), \, \text{ for } j = 1, 2, \cdots, k \\ \nu _ { \lambda _ { \beta } } \ | \, \text{rest} \sim \text{Inverse-Gamma} ( 1, \, \frac { 1 } { \tau _ { \beta } ^ { 2 } } + \sum _ { j = 1 } ^ { k } \frac { 1 } { \lambda _ { \beta _ { j } } ^ { 2 } } ) \\ \tau _ { \beta } ^ { 2 } \ | \, \text{rest} \sim \text{Inverse-Gamma} ( 1, \, \frac { 1 } { \nu _ { \lambda _ { \beta } } } + \frac { 1 } { \nu _ { \tau _ { \beta } } } ) \\ \nu _ { \tau _ { \beta } } \ | \, \text{rest} \sim \text{Inverse-Gamma} ( 1, \, \frac { 1 } { \tau _ { \beta } ^ { 2 } } + \frac { 1 } { \sigma _ { 2 } ^ { 2 } } ) ; \\ \text{ups of $ok$ and $\sigma^{2}$.}$$

## 2. Priors of ϕ γ and σ 2 γ :

$$& \phi _ { \gamma } \, | \, \text{rest} \sim N \left ( \frac { \sum _ { t = 1 } ^ { T _ { 0 } } \gamma _ { t - 1 } ^ { T } \gamma _ { t } } { \sum _ { t = 1 } ^ { T _ { 0 } } \gamma _ { t - 1 } ^ { T } \gamma _ { t - 1 } }, \, \frac { \sigma _ { \gamma } ^ { 2 } } { \sum _ { t = 1 } ^ { T _ { 0 } } \gamma _ { t - 1 } ^ { T } \gamma _ { t - 1 } } \right ), \\ & \sigma _ { \gamma } ^ { 2 } \, | \, \text{rest} \sim \text{Inverse-Gamma} \left ( \frac { 1 } { 2 } + \frac { T _ { 0 } } { 2 }, \, \frac { 1 } { \nu _ { \sigma _ { \gamma } } } + \frac { 1 } { 2 } \sum _ { t = 1 } ^ { T _ { 0 } } ( \gamma _ { t } - \phi _ { \gamma } \gamma _ { t - 1 } ) ^ { T } ( \gamma _ { t } - \phi _ { \gamma } \gamma _ { t - 1 } ) \right ), \\ & \nu _ { \sigma } \, | \, \text{rest} \sim \text{Inverse-Gamma} \left ( 1, \, \frac { 1 } { \sigma _ { \gamma } ^ { 2 } } + \frac { 1 } { 1 0 ^ { 2 } } \right ). \\.. & \. & \. & \. & \. & \. & \.. & \.. & \.. & \.. & \..$$

After drawing samples from the above full conditional distributions, we update γ t = ϕ γ γ t -1 + v t .

## 3. Priors η 0 , σ 2 η and ω , 2 1 · · · , ω 2 p :

$$\underset { p \times 1 } { \eta _ { i } } \ | \ \text{rest} \sim \mathcal { N } _ { p } \left ( C ^ { - 1 } \sum _ { t = 1 } ^ { T _ { 0 } } \gamma _ { t } ^ { \top } ( Y ^ { c } ( 0 ) - u _ { i t } ), \ \sigma _ { 2 } ^ { 2 } C ^ { - 1 } \right ), \ \text{ for } i = 1, 2, \cdots, N$$

$$& \underset { p \times 1 } { \eta _ { i } } \ | \ \text{rest} \sim \mathcal { N } _ { p } \left ( \mathcal { C } \ \, \nu \sum _ { t = 1 } ^ { \ } \gamma _ { t } \left ( \mathbf Y \ \kappa ( 0 ) - u _ { i t } \right ) \\ & \text{where } \ C = \sum _ { t = 1 } ^ { T _ { 0 } } \gamma _ { t } \gamma _ { t } ^ { \top } + \sigma _ { 2 } ^ { 2 } \sigma _ { \eta } ^ { - 2 } \mathbf S _ { \eta } ^ { - 1 },$$

$$\text{and} \ u _ { i t } = \rho w _ { i } Y _ { 0 t } ( 0 ) + \rho \sum _ { j = 1 } ^ { N } W _ { i j } Y _ { j t } ( 0 ) + X _ { i t } \beta _ { 0 }$$

$$\sigma _ { \eta } ^ { 2 } \, | \, \text{rest} \sim \text{Inverse-Gamma} \left ( \frac { 1 } { 2 } + \frac { p N } { 2 }, \, \frac { 1 } { \nu _ { \sigma _ { \eta } } } + \frac { 1 } { 2 } \sum _ { i = 1 } ^ { N } \eta _ { i } ^ { 0 \top } \Sigma _ { \eta } ^ { - 1 } \eta _ { i } ^ { 0 } \right )$$

$$\omega _ { j } ^ { 2 } \, | \, \text{rest} \sim \text{Inverse-Gamma} { \left ( \frac { 1 } { 2 } + \frac { N } { 2 }, \, \frac { 1 } { \nu _ { \omega _ { k } } } + \frac { 1 } { 2 } \sum _ { i = 1 } ^ { N } \frac { ( \eta _ { i j } ^ { 0 } ) ^ { 2 } } { \sigma _ { \eta } ^ { 2 } } \right ) }, \, \text{ for } j = 1, 2, \cdots, p$$

$$\nu _ { \omega _ { j } } \ | \ \text{rest} \sim \text{Inverse-Gamma} \| \text{$(1,\ \frac{1}{\omega_{j}^{2}}+\frac{1}{10^{2}}$} \| ), \ \text{for $j=1,2,\cdots,p$.}$$

$$& \overleftarrow { j = 1 } \\ \i a m m a & \left ( \frac { 1 } { 2 } + \frac { p N } { 2 }, \ \frac { 1 } { \nu _ { \sigma _ { \eta } } } + \frac { 1 } { 2 } \sum _ { i = 1 } ^ { N } \eta _ { i } ^ { 0 \top } \Sigma _ { \eta } ^ { - 1 } \eta _ { i } ^ { 0 } \right ) \\ \i a m m a & \left ( \frac { 1 } { 2 } + \frac { N } { 2 }, \ \frac { 1 } { \nu _ { \omega _ { k } } } + \frac { 1 } { 2 } \sum _ { i = 1 } ^ { N } \frac { ( \eta _ { i j } ^ { 0 } ) ^ { 2 } } { \sigma _ { \eta } ^ { 2 } } \right ), \ \text{ for } j = 1, \colon \\ \i a m m a & \left ( 1, \ \frac { 1 } { \omega _ { j } ^ { 2 } } + \frac { 1 } { 1 0 ^ { 2 } } \right ), \ \text{ for } j = 1, 2, \cdots, p.$$

## 4. Prior of σ 2 2 :

$$\sigma _ { 2 } ^ { 2 } \, | \, \text{rest} \, \sim \, \text{Inverse-Gamma} \, & \left ( 1 + \frac { N T _ { 0 } } { 2 }, \, \frac { 1 } { \nu _ { \beta } } + \frac { 1 } { \nu _ { \sigma _ { 2 } } } + \sum _ { t = 1 } ^ { T _ { 0 } } u _ { t } ^ { \top } u _ { t } \right ) \\ \text{where} \, \ u _ { t } = Y _ { t } ^ { c } ( 0 ) - \rho w Y _ { 0 t } ( 0 ) - \rho W Y _ { t } ^ { c } ( 0 ) - X _ { t } \beta _ { 0 } - \eta ^ { 0 } \gamma _ { t } \\ \nu _ { \sigma _ { 2 } } \, | \, \text{rest} \, \sim \, \text{Inverse-Gamma} \, & \left ( 1, \, \frac { 1 } { \sigma _ { 2 } ^ { 2 } } + \frac { 1 } { 1 0 ^ { 2 } } \right ).$$

- 5. Regarding ρ , because we cannot analytically derive the full conditional distribution, we use the Metropolis algorithm:
- (a) Let ρ ( c ) be the current value. Generate ρ ∗ from N ρ ( ( c ) , k 2 ).
- (b) Using the following, we compute likelihood p ρ ( c | rest) and p ρ ( ∗ | rest), respectively:

p ρ ( | rest)

$$= \prod _ { t = 1 } ^ { T _ { 0 } } | I _ { N } - \rho W | \exp \left ( - \frac { 1 } { 2 \sigma ^ { 2 } } ( A Y _ { t } ^ { c } - \rho w Y _ { 0 t } - X _ { t } \beta - \eta ^ { 0 } \gamma _ { t } ) ^ { \top } ( A Y _ { t } ^ { c } - \rho w Y _ { 0 t } - X _ { t } \beta - \eta ^ { 0 } \gamma _ { t } ) \right ),$$

where A = I N -ρ W .

- (c) Let r = min 1 { , p ( ρ ∗ | rest) /p ρ ( ( c ) | rest) } be the acceptance probability. Subsequently, we generate u ∼ U (0 , 1) and update ρ ( c ) as ρ ∗ if r &gt; u , otherwise stay ρ ( c ) .

This algorithm is referred from LeSage and Pace (2009).
Shosei Sakaguchi Hayato Tagawa
Identification and Inference for Synthetic Control Methods with Spillover Effects: Estimating the Economic Cost of the Sudan Split
The synthetic control method (SCM) is widely used for causal inference with
panel data, particularly when there are few treated units. SCM assumes the
stable unit treatment value assumption (SUTVA), which posits that potential
outcomes are unaffected by the treatment status of other units. However,
interventions often impact not only treated units but also untreated units,
known as spillover effects. This study introduces a novel panel data method
that extends SCM to allow for spillover effects and estimate both treatment and
spillover effects. This method leverages a spatial autoregressive panel data
model to account for spillover effects. We also propose Bayesian inference
methods using Bayesian horseshoe priors for regularization. We apply the
proposed method to two empirical studies: evaluating the effect of the
California tobacco tax on consumption and estimating the economic impact of the
2011 division of Sudan on GDP per capita.