variation such that it is possible to generalize to a given target counterfactual treatment. Further,
NSI obtains an estimate within error εwith high probability under the proposed experiment design
when Opoly(d)/ε4training samples per unit are available. This is a significant improvement over
the O(exp(d)/ε2) training samples that a naive procedure would require.
We conclude with simulations showing that NSI is robust to spillovers under which existing
estimators are biased.
1.1 Related work
The literature on causal inference with network interference or spillover effect has mostly considered
the setting of a single measurement per unit, whether in the setting of a randomized experiment or
an observational study. Under fully arbitrary interference, it has been shown that it is impossible
to estimate any desired causal estimands as the model is not identifiable (Manski 2013,Aronow
et al. 2017,Basse and Airoldi 2018a,Karwa and Airoldi 2018). Subsequently, various models
have been proposed in the literature that impose restrictions on the exposure functions (Manski
2013,Aronow et al. 2017,Viviano 2020,Auerbach and Tabord-Meehan 2021,Li et al. 2021),
interference neighborhoods (Ugander et al. 2013,Bargagli-Stoffi et al. 2020,Sussman and Airoldi
2017a,Bhattacharya et al. 2020), parametric structure (Toulis and Kao 2013,Basse and Airoldi
2018b,Cai et al. 2015,Gui et al. 2015,Eckles et al. 2017), two-sided platforms (Johari et al.
2022,Bajari et al. 2021) or a combination of these, each leading to a different solution concept. A
comprehensive review on network interference models is given by De Paula (2017). In this work,
we focus on network interference that is additive across the neighbors, referred to in the literature
as the joint assumptions of neighborhood interference, additivity of main effects, or additivity of
interference effects (Sussman and Airoldi 2017a,Yu et al. 2022,Cortez et al. 2022a,b).
Distinct to our work is that we consider a panel data setting in which there are multiple mea-
surements (e.g., a time series) for each unit. The potential outcomes function is thus also dependent
on both the unit and the measurement. Additionally, we allow for the estimation of unit-specific
counterfactuals under multiple treatments, whereas the existing literature has largely focused on
binary treatments. Key to our approach is a novel latent factor model that takes into account net-
work interference and is a generalization of the factor models typically used in panel data settings.
Previous work has focused on causal estimands that capture population-level effects, such as the
average direct treatment effect (the average difference in outcomes if only one unit and none of its
neighbors get treated (Basse and Airoldi 2018b,Jagadeesan et al. 2020,S¨avje et al. 2021,Sussman
and Airoldi 2017a,Leung 2019,Ma and Tresp 2021)) and the average total treatment effect (the
average difference in outcomes if all units get treated versus if they do not (Ugander et al. 2013,
Eckles et al. 2017,Chin 2019,Yu et al. 2022,Cortez et al. 2022a,b)). Alternately there has been
some literature that focuses on hypothesis testing for the presence of network interference (Aronow
2012,Bowers et al. 2013,Athey et al. 2018,Pouget-Abadie et al. 2017,Saveski et al. 2017); these
results do not immediately extend to estimation as they are based on randomization inference with
a fixed network size, and focus on testing the sharp null hypotheses.
While a majority of the literature focuses on randomized experiments, there is a growing in-
terest in the literature to account for network interference when analyzing observational studies.
The existing literature generally assumes partial interference, where the network consists of many
disconnected sub-communities (Tchetgen and VanderWeele 2012,Perez-Heydrich et al. 2014,Liu
et al. 2016,DiTraglia et al. 2020,Vazquez-Bare 2022). Without this strong clustering condition,
other works impose strong parametric assumptions on the potential outcomes function, assuming
that the potential outcomes only depend on a known statistic of the neighborhood treatment, e.g.
the number or fraction of treated (Verbitsky-Savitz and Raudenbush 2012,Chin 2019,Ogburn et al.
3