## Heterogeneous Endogenous Effects in Networks ∗

Sida Peng †

August 5, 2019

## Abstract

This paper proposes a new method to identify leaders and followers in a network. Prior works use spatial autoregression models (SARs) which implicitly assume that each individual in the network has the same peer effects on others. Mechanically, they conclude the key player in the network to be the one with the highest centrality. However, when some individuals are more influential than others, centrality may fail to be a good measure. I develop a model that allows for individual-specific endogenous effects and propose a two-stage LASSO procedure to identify influential individuals in a network. Under an assumption of sparsity: only a subset of individuals (which can increase with sample size n ) is influential, I show that my 2SLSS estimator for individual-specific endogenous effects is consistent and achieves asymptotic normality. I also develop robust inference including uniformly valid confidence intervals. These results also carry through to scenarios where the influential individuals are not sparse. I extend the analysis to allow for multiple types of connections (multiple networks), and I show how to use the sparse group LASSO to detect which of the multiple connection types is more influential. Simulation evidence shows that my estimator has good finite sample performance. I further apply my method to the data in Banerjee et al. (2013) and my proposed procedure is able to identify leaders and effective networks.

key words: key players, network, endogenous effects, spillovers, high-dimensional models, LASSO, model selection, robust inference

∗ I would like to thank Francesca Molinari, Matthew Backus, David Easley, Marten Wegkamp, Donald Kenkel, Zhuan Pei, Douglas Miller, Joris Pinkse, Peter Hull, Yanlei Ma and participants in seminars and conferences at which this paper was presented. All remaining errors are mine.

† Microsoft Research, sidpeng@microsoft.com

## 1 Introduction

How an individual's behavior is affected by the behavior of her neighbors in an exogenously given network is an important research question in applied economics. With the increasing availability of detailed data documenting connections among individuals, spatial autoregression models (SARs) have been widely applied in the empirical networks literature to estimate endogenous effects.

In SARs, an individual's behavior depends on the weighted average of other individuals' behaviors (see Anselin, 1988; Kelejian and Prucha, 1998). Standard SARs assume that the peer effects/endogenous effects are the same across individuals in a network. Each individual influences her neighbors at the same rate regardless of who she is. However, in many contexts, some individuals are clearly more influential than others. For example, Mas and Moretti (2009) finds that the magnitude of spillovers varies dramatically among workers with different skill levels. Clark and Loheac (2007) also notes that popular teenagers in a school have much stronger influence on their classmates' smoking decisions than their less popular peers.

I propose a novel SAR model which allows for heterogeneous endogenous effects. Each individual in a network simultaneously generates an outcome that takes into account all her neighbors' behaviors. Unlike standard SARs, each individual has an individual-specific effect on her neighbors. As a result, there are as many coefficients for individual-specific endogenous effects as there are individuals in the network. To achieve identification, I assume that 'truly-influential' individuals only constitute a small fraction of the total population. In other words, individual-specific coefficients are assumed to be sparse. This assumption allows me to estimate the model via the least absolute shrinkage and selection operator (LASSO). The LASSO procedure penalizes the l 1 norm for the coefficients of heterogeneous endogenous effects. The geometry of the l 1 norm enforces the sparsity in the LASSO estimators. If a coefficient is selected by LASSO (i.e. the estimated coefficient is nonzero), the individual associated with this coefficient can influence all her neighbors at her specific rate. Otherwise the LASSO estimator will indicate that the individual has no influence on her neighbors. With some restrictions on the network structures, I show that the LASSO estimates for heterogeneous endogenous effects have near oracle performance (see B¨hlmann and van de Geer, u 2011). In other words, the selection of influential individuals is consistent and the convergence rate of non-zero LASSO estimates is the same as the convergence rate that would have been achieved if the truly influential individuals were known.

One challenge in my estimation process is the presence of endogeneity in the spatial lag and error term. As with standard SARs, the dependent variable in my model is used to construct spatial lags as an independent variable. As a result, the regressors are correlated with the error term and

estimates would be biased if we were to apply LASSO directly.

First I propose a set of novel instruments to address the endogeneity. Following Kelejian and Prucha (1998), I express the dependent variable as an infinite sum of functions consisting of exogenous characteristics and an adjacency matrix. I show that the exogenous characteristics of influential individuals can be used as instruments for their neighbors. Then I design a two-stage estimation process for heterogeneous endogenous effects using LASSO at each stage. In the first stage , I use LASSO to estimate the coefficients for the instruments. These estimated coefficients and instruments are then used to create a synthetic dependent variable. In the second stage , I replace the dependent variable in the spatial lags with the synthetic variable to perform the LASSO estimation.

The next challenge is to construct robust confidence intervals for my LASSO type two-stage estimator. As pointed out in Leeb and Potscher (2008), it is impossible to estimate the distribution for post-model-selection estimators. Consistent model selection by LASSO is only guaranteed when all non-zero coefficients are large enough to be distinguished from zero in a finite sample (i.e. a condition usually named 'beta-min'). LASSO may fail to select regressors with very small coefficients, resulting in omitted variable bias in the post LASSO inference.

I propose a bias correction for my two-stage estimator following the recent LASSO inference literature (see Chernozhukov et al., 2018; Belloni et al., 2015; van de Geer et al., 2014; Javanmard and Montanari, 2014; Zhang and Zhang, 2011; Zhu, 2018). The idea is to correct the first order bias and make the estimators independent from the model selection. Heuristically, shrinkage bias due to the l 1 penalty in LASSO can be expressed as a function of the LASSO estimators. Normality can still be achieved after adjusting for this bias. I show that this strategy also works in my two-stage estimation process under the presence of spatial errors. I derive the asymptotic normality for my 'de-bias' two-stage LASSO estimator and conduct robust inference including confidence intervals.

My model can be extended to allow for more flexible structures of influencing. One real world scenario is a network which consists of many local leaders. Local leaders may only influence individuals within their own cliques/groups but have no influence on individuals outside their cliques/groups. Influence from local leaders can be represented by a homogenous endogenous effects. Exogenous effect and homophily within cliques may also be present.

To solve the problem in this scenario, I modify my model by bringing back the classical SAR model. I assume that there are both local leaders and global leaders in the network. In contrast to local leaders, global leaders can influence individuals across different cliques. I show that homogeneous endogenous effects across local leaders, exogenous effects, and correlated effects can be identified under assumptions similar to those Bramoull´ et al. (2009). e Under a sparsity assumption, the en-

dogenous effects of global leaders, whose influence remain individual-specific, can also be identified under the same set of instruments proposed in my main model.

Another real world scenario is the existence of multiple types of connections among individuals. For example, connections among individuals can be classified as social (e.g. friendship, kinship) or economic (e.g. lending, employment). In epidemiology, infectious disease can be spread through air, insects, or direct contact. It is important to identify which networks are more efficient at transmitting the endogenous effects.

I model different types of connections as multiple networks. I propose the use of sparse group LASSO to estimate a heterogeneous endogenous effects model with multiple networks. The sparse group LASSO penalizes both the l 1 norm and the l 2 norm for each coefficient in each type of connection and thus selects both the influential individuals and effective types of connections that generate spillover. I derive the convergence rate and prove the consistency of selection for this estimator. To the best of my knowledge, my paper is the first to show statistical properties for sparse group LASSO.

I provide simulation evidence for networks of different sizes and different data generating algorithms. The empirical coverage of my proposed estimators is close to the nominal level in all scenarios. Similar results are also found in models with multiple networks and with cliques.

I apply my method to study villagers' decisions to participate in micro-finance programs in rural areas of India as in Banerjee et al. (2013). Instead of simulating the contagion process and reconstructing the data into panel format as in Banerjee et al. (2013), my method allows researchers to directly analyze cross-sectional data under the standard equilibrium assumptions. Among different social and economic networks, my method shows that some networks such as 'visit go-come' and 'borrow money', are much more effective at influencing villagers' decisions than other networks such as 'go to temple together' and 'medical help'. I further show that individuals in certain careers such as agricultural workers, Anganwadi teachers and small business owners are more likely to influence villagers.

## 1.1 Literature Review

This paper brings together literature on spatial autoregression model, LASSO and networks.

## SARs :

SARs have been widely applied in empirical studies. For instance, they have been used to study

peer effects in labor productivity (see Mas and Moretti, 2009; Guryan et al., 2009; Bandiera et al., 2009), smoking behavior among teenagers (see Krauth, 2005; Clark and Loheac, 2007; Nakajima, 2007), educational achievements among different student groups (see Sacerdote, 2001; Neidell and Waldfogel, 2010), systemic risk in finance (see Bonaldi et al., 2015; Denbee et al., 2015), and the adoption of new agricultural technologies (see Coelli et al., 2002; Conley and Udry, 2010). My paper proposes a novel extension of standard SARs that could be used to identify influential individuals in a given network. My methodology for estimating such a model could easily be adopted in existing empirical SARs analyses to identify influential individuals who influence their peers productivity, smoking decisions, or financial holdings.

More specifically, my model extends existing SARs literature by introducing heterogeneous endogenous effects. Until very recently, SARs always assume a constant rate of dependence for endogenous effects across different individuals. Moreover, row-sum normalization is widely adopted in the estimation, which may further result misspecification under the existence of heterogeneous endogenous effects. (see Cliff and Ord (1973), the first monograph on the topic, and the later studies, Upton and Fingleton (1985); Anselin (1988); Cressie (1993); Lee and Liu (2010); Lee and Yu (2010); Jin and Lee (2018)). Recent developments in social interaction literature incorporate individual characteristics into SARs, essentially allows for a limited degrees of pre-stipulated heterogeneity, for example, (see Pinkse et al., 2002). In contrast, this paper shows that heterogeneous endogenous effects can be identified from individuals' outcomes instead of being pre-specified through individuals' characteristics. Another article studying individual heterogeneity from random coefficient aspect is Masten (2018). Masten (2018)'s model assumes each individual has different rate of receiving the influence while my model assumes each individual has different rate of influencing her neighbors. The two different models yield two different identification strategies. This paper contributes to the literature by showing heterogeneous endogenous effects can be directly identified using cross-sectional data.

To estimate the heterogeneous endogenous effects in my model, I propose a methodology that is different from standard SARs literature. In classic SARs, there is only one endogenous variable and hence it is sufficient to identify the model through only one instrument. In my model, the number of potentially endogenous variables increases as the number of observations increases. As a result, I propose a set of instruments that contain the same number of instruments as the total number of individuals. Each instrument is essentially a decomposition from the standard SARs instruments as in Kelejian and Prucha (1998), Lee (2002), Lee (2003) and Lee (2004).

I show that my model can be combined with SAR model under homogeneous endogenous effects, exogenous effects and correlated effects. To solve the classical 'reflection problem' as in Manski

(1993), I adopt the same strategy as in Bramoull´ et al. (2009) and show that 'neighbors' neighbor' e instruments can be combined with my set of instruments to identify additional structures in the spillovers.

This paper also contributes to literature that models multiple networks through SARs. In standard SARs, multiple networks are modeled as higher order spatial lags (see Lee and Liu, 2010). My model allows each individual to have her own specific endogenous effects in each network. This design enables the selection at the network level which implies some networks can be classified as completely irrelevant to decision-making. Furthermore, as the asymptotic allows the number of networks q to go to infinity, my model can handle those cases when researchers observe a large number of networks ( q &gt; n ).

## LASSO :

My paper contributes to the growing literature on endogenous regressors in LASSO estimators. For instance, Belloni et al. (2014) proposes the double selection mechanism to study confounded treatment effects. Fan and Liao (2014) proposes a GMM type estimator to deal with many endogenous regressors. Gautier and TsyBakov (2014) proposes a Self Tuning Instrumental Variables (STIV) estimator. The paper that is closest to mine is Zhu (2018), which studies the statistical properties of two-stage least square procedure with high-dimensional endogenous regressors. The two-stage estimator proposed in Zhu (2018) assumes an i.i.d error term for the first stage. However, this assumption is incompatible with SAR model as the structure assumptions in SAR lead to a first stage with correlated unobservables. In this paper, I derive the rate of convergence and consistency of selection for a two-stage LASSO estimator under the presence of spatial errors. I show that a modified 'de-bias' LASSO estimator accounting for spatial errors can be constructed for my estimator in a manner similar to Zhang and Zhang (2011), B¨hlmann (2013), van de Geer et al. u (2014), and Zhu (2018). I derive its asymptotic distributions and show how to perform inference.

This paper also extends LASSO literature by deriving statistical bounds and consistency of selection for sparse group LASSO estimator. Yuan and Lin (2006) proposes the group LASSO, in which explanatory variables are represented by different groups. The group LASSO assumes that sparsity exists only among groups, i.e. some groups of variables are relevant while other groups are not. Simon et al. (2013) proposes the sparse group LASSO, which further allows sparsity within each group, i.e. some regressors within the relevant groups can also be irrelevant. Bunea et al. (2014) derives statistical properties for the square-root group LASSO, which combines group LASSO and square-root LASSO. When estimating a heterogeneous endogenous effects model with multiple networks, I establish both statistical bounds and consistency of selection for the sparse group LASSO estimator. Sparse group LASSO differs from group LASSO by allowing the number of

regressors in each group to also goes to infinity as the number of groups goes to infinity. To the best of my knowledge, this paper is the first to show asymptotic statistical properties for the sparse group LASSO estimator.

## Networks :

My paper shares similar microfoundations with SARs as discussed in Blume et al. (2015), where the individual utility function can be written as a linear summation of the private and social components. The private component is a quadratic loss function on individual's efforts. The social components depend on the network structure as well as the efforts of one's neighbors. While the marginal rate of substitution between the private and social components of utility is assumed fixed in SARs, I assume this rate is individual-specific and depends on one's neighbors. My paper applies and extends LASSO approaches to deal with a high-dimensional problem in networks. The total number of possible edges in a network is n 2 , however, the social interaction networks we often observe are far more sparse. This is an ideal setting where penalized estimators like LASSO could be applied. Manresa (2013) studies the heterogeneous exogenous effects in a network using LASSO. de Paula et al. (2015) explores the use of LASSO to recover network structures. Both these two papers consider panel data and rely on repeated observations of the same network to identify their models. My model considers cross-sectional data. To identify an individual's endogenous effects, I rely on the variations in her neighbors' outcome.

My paper also relates to the literature on identifying the key players in the network following Ballester et al. (2006), Calv´-Armengol et al. (2009), and Horrace et al. (2016). o Under the framework of SARs, every individual is assumed to have the same endogenous effects. As a result, individuals who are well-connected in the network (with high centrality measure) are considered as the key players. However, well connected individuals may effectively have zero effects on their neighbors under heterogeneous endogenous effects. Indeed, as shown in the empirical application, well connected villagers such as tailors, hotel workers, veterans, and barbers are not influential in other villagers' decisions to join the micro-finance program.

The rest of this paper is organized as follows: in Section 2, I introduce the model; in Section 3, I discuss assumptions; in Section 4, I design estimation procedures and derive consistency and asymptotic properties; in Section 5, I show finite sample performance using Monte Carlo simulations; in Section 6, I apply my proposed model to study influential individuals and effective networks in promoting micro-finance programs in rural India; and in Section 7, I conclude.

## 2 Models

## 2.1 Benchmark Endogenous Effects Model

Let n denotes the total number of observed individuals in a network. The outcome of individual i is denoted as d i and is the variable of interest. Here d i can represent outcome variable of interest associated with individual i , such as whether to join a program or i 's labor productivity. In standard SAR model, it is assumed that the outcome of each neighbor of individual i impacts her outcome homogeneously through a constant rate λ 0 :

$$d _ { i } = \lambda _ { 0 } \sum _ { j \in N _ { i } } d _ { j } + x _ { i } \beta _ { 0 } + \epsilon _ { i },$$

where the set N i is defined as individual i 's neighbors. The matrix form of this model is expressed as follows:

$$D _ { n } = \lambda _ { 0 } M _ { n } D _ { n } + X _ { n } \beta _ { 0 } + \epsilon _ { n },$$

where D n = ( d , d 1 2 , · · · , d n ) ′ is the n -dimensional vector of observable outcomes. The n by k matrix X n represents the observable exogenous characteristics of individuals. When glyph[epsilon1] n is specified as an n -dimensional vector of independent and identically distributed disturbances with zero mean and a constant variance σ 2 , equation (2) is also called a mixed regression model.

The spatial weight matrix M n is of size n by n , where the ( i, j )th entry represents the connection between individual i and individual j . In empirical studies, the spatial weight matrix is often replaced by the adjacency matrix (see Ammermuller and Pischke, 2009; Acemoglu et al., 2012; Banerjee et al., 2013): the ( i, j )-th entry of the matrix M n takes value 1 if individual i and individual j are connected and takes value 0 otherwise; the diagonal entries of the matrix M n are always 0s. In the SAR literature, spatial weight matrix or adjacency matrix is taken as exogenous. My method follows this assumption and is designed for cross-sectional data. However, it is important to recognize that social networks do change over time and it is then important to take network formation into modeling (see Christakis et al. (2010), Sheng (2016)).

In mixed regression model, endogenous effects (see Manski, 1993) or network effects (see Bramoull´ e et al., 2009) are captured by the scalar λ 0 . An implicit assumption in equation (2) is that λ 0 , the rate of endogenous effects, is identical across all individuals in the network. Although a limited degree of heterogeneity can be built into the adjacency matrix via pre-specified structure assumptions, the identification potential for heterogeneous endogenous effects has not been fully explored. This limitation has been noted in various studies (see Ammermuller and Pischke, 2009; de Paula

et al., 2015). I relax this assumption by proposing a more flexible model that allows and identifies individual-specific endogenous effects as discussed below.

## 2.2 Heterogeneous Endogenous Effects Model

I propose the following model:

$$d _ { i } = \sum _ { j \in N _ { i } } d _ { j } \eta _ { j } + x _ { i } \beta + \epsilon _ { i }$$

where N i represents the set of individual i 's neighbors and η j represents the endogenous effects of individual j on the outcome of all her neighbors i ∈ N j . the model can be rewritten in matrix form as:

$$D _ { n } = \left ( M _ { n } \circ D _ { n } \right ) \eta _ { 0 } + X _ { n } \beta _ { 0 } + \epsilon _ { n },$$

glyph[negationslash]

where η 0 = ( η , η 1 2 , · · · , η n ) ′ is a vector of parameter of size n by 1. The i th entry in η 0 represents the endogenous effects of individual i on her neighbors. This model allows for individual heterogeneity to interact with endogenous effects so that every individual is allowed to have her own coefficient η i . My model allows some η j = 0. In other word, there are individuals that impose no endogenous effects on their neighbors. I define those individuals with η j = 0 as influential.

The operator ◦ is defined between a n by n matrix M n and a n by 1 vector D n as

$$M _ { n } \circ D _ { n } = M _ { n } \cdot \text{ diag} ( D _ { n } ) = C,$$

where diag ( ) is the diagonalization operator and · C i,j = M d i,j j .

Note that in contrast to fixed rate λ 0 specified in equation (2), even though each neighbor of individual j is assumed to receive the same influence d η j j from her, each individual is allowed to influence her neighbors at her own rate η j .

A more generalized form of the model is to replace η j with η ij . This generalization can further capture the different perceiving rate of endogenous effects on the same leader but among different followers. However, the number of parameters increase from n to n 2 and the model becomes too saturated to estimate.

Another direction of modeling is to assume heterogeneous perceiving rate of endogenous effects but maintaining the assumption of homogeneous influencing rate, for example

$$d _ { i } = \lambda _ { i } ^ { * } \sum _ { j \in N _ { i } } d _ { j } + x _ { i } \beta + \epsilon _ { i }$$

This model has been studied using random coefficients approach in Masten (2018). This specification differs from equation (3) as it prevents any sparsity pattern on λ ∗ i to be passed to a first stage regression as shown in proposition 1. As a result, the identification strategies are fundamentally different.

Equation (4) can be derived from a bayesian Nash Equilibrium. Let ( x , glyph[epsilon1] i i ) denotes an individual's type, where x i is publicly observed characteristics and glyph[epsilon1] i is private characteristics only observable by i . Individual i 's utility depends on her own action and characteristics as well as her neighbors' actions. Individual i chooses action d i to maximize the following utility:

$$U _ { i } ( d _ { i }, d _ { - i } ) = ( x _ { i } \beta + \epsilon _ { i } ) d _ { i } - \frac { 1 } { 2 } d _ { i } ^ { 2 } + \sum _ { j \in N _ { i } } d _ { j } d _ { i } \eta _ { j }$$

The first order condition yields equation (4). The micro-foundations derived above is similar to the one for SARs as discussed in Blume et al. (2015).

## 2.3 Motivating Examples

## Peer Effects in Labor Productivity:

Understanding the mechanism and magnitude of the dependence of labor productivity on coworkers is an important question for economists and policy makers. As found in Mas and Moretti (2009), workers respond more to the presence of coworkers with whom they frequently interact. Another modern example is on code sharing platforms like GitHub. Programmers provide efforts to an open source project depending on how much other programers are contributing. In such cases, the influence level of each individual to hers coworkers is not the same. Equation (4) can be used to incorporate such differences.

$$y _ { i } = \sum _ { j \in N _ { i } } y _ { j } \eta _ { j } + x _ { i } \beta + \epsilon _ { i },$$

where y i is individual i 's productivity, and η j represents the size of influence of coworker j - all else being equal, the additional effect on individual i 's productivity if individual j becomes her coworker.

Note that if we restrict the parameters η j to be the same across different workers, then we are back to the classical SARs setting as laid out in equation (2). Thus, λ = 1 n ∑ n j =1 η j can be interpreted as the averaged spillover effects in the canonical sense. Notice that λ converges to 0 when the spillover effects are sparse and it may lead to the conclusion that there is no spillover in the network under such scenario.

## Define

glyph[negationslash]

glyph[negationslash]

$$\lambda ^ { * } = \frac { 1 } { \sum { \mathbf 1 } _ { \eta _ { j } \neq 0 } } \sum _ { j = 1 } ^ { n } \eta _ { j } { \mathbf 1 } _ { \eta _ { j } \neq 0 }$$

as the averaged endogenous effects for influential workers. The estimand λ ∗ is arguably better than λ for average endogenous effects as it excludes individuals who are not influential to their neighbors.

## Online Opinion Leaders:

A decision can represent whether to 'tweet' a news story seen online. When individuals make such decisions, they are often influenced by several online opinion leaders - whether those people also 'tweet' the news or not. Political figures may have stronger influence on people's tweet for political news than celebrities. And vice versa for entertainment news. Assume a binary decision (0 1) is made from a bayesian Nash Equilibrium, such that ,

$$d _ { i } ^ { * } = \sum _ { j \in N _ { i } } d _ { j } ^ { * } \eta _ { j } + x _ { i } \beta + \epsilon _ { i },$$

where d ∗ i is the probability of individual i playing action 1, and ∑ j ∈ N j d η ∗ j j is the expected endogenous effects from i 's neighbors N i . Define glyph[negationslash]

$$S = \{ j \, \colon \eta _ { j } \neq 0 \}$$

as truly influential opinion leaders. Mymethod provides a way to estimate ˆ which is asymptotically S consistent with S . This is an important metric to policy makers or private sectors as targeting or nudging the influential individuals is usually more efficient than targeting the entire population.

## 2.4 Two Extensions of the Model

## Heterogeneous Endogenous Effects Model with Cliques:

There are two important assumptions for equation (4). First is the sparsity assumption, which requires non-influential individuals to have completely zero influence on their neighbors. Second it assumes away the exogenous and correlated effects. Consider a network composed of many cliques (small groups of connected individuals). Each clique has its local leaders who only influence individuals within their own clique. Figure 1 provides an example of such a network structure. Note that in Figure 1, node S 2 , S 3 and S 4 represent local leaders who only influence individuals within their own cliques. On the contrary, node S 1 represents a global leader who can influence individuals across different cliques. For example, one can think about the local leaders S 2 , S 3 and

S 4 as local news channels while S 1 is the national news channel. Furthermore, within each clique, there might exist exogenous and correlated effects. I propose an extension to my heterogeneous endogenous effects model which could address such concerns. I assume that all local leaders will influence their followers at a small but similar rate while global leaders can have different effects on their audience.

Figure 1: Local Leader

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00663v1-figure-5338502f-1200-4f95-982d-0a9988897935.jpg)

First, I consider the following extension to introduce only homogeneous endogenous effects:

$$d _ { i } = \sum _ { j \in N _ { i } } d _ { j } \eta _ { j } + \gamma _ { 0 } \sum _ { j \in N _ { i } } d _ { j } + x _ { i } \beta _ { 0 } + \epsilon _ { i },$$

which can be represented in matrix form as:

$$D _ { n } = \left ( M _ { n } \circ D _ { n } \right ) \eta _ { 0 } + M _ { n } D _ { n } \gamma _ { 0 } + X _ { n } \beta _ { 0 } + \epsilon _ { n },$$

where η ′ 0 = ( η , η 1 2 , · · · , η n ) . ′ The new term γ 0 ∑ j ∈ N i d j captures influence from the local level. Note that this is the same term as the spatial lag in the benchmark spatial autoregression model. The vector η 0 captures the heterogeneous endogenous effects of global leaders.

To further include exogenous and correlated effects, consider the following model:

$$d _ { i } = \sum _ { j \in N _ { i } } d _ { j } \eta _ { j } + \gamma _ { 0 } ^ { e n d } \sum _ { j \in N _ { i } } d _ { j } + \gamma _ { 0 } ^ { e x o } \sum _ { j \in N _ { i } } x _ { j } + x _ { i } \beta _ { 0 } + \mu _ { c } + \epsilon _ { i },$$

Besides the first term representing heterogeneous endogenous effects, equation (6) is the same as the model in Manski (1993).

## Heterogeneous Endogenous Effects Model with Multiple Networks:

Individuals are often connected with each other through more than one type of network. For example, one's financial network (for borrowing/lending) may be different from hers relative network

or even friendship network. A common strategy in empirical application is to 'pool' the networks by mixing all types of connections into one network. This approach may increase the noise in the network measurement when the outcome variables only depend on certain types of networks.

To capture different types of connections among the same set of individuals, we can incorporate multiple networks in the heterogeneous endogenous model. More specifically, a separate adjacency matrix can be constructed for each type of network. For instance, the ( i, j )-th entry of the adjacency matrix representing friendship takes value 1 if individual i and individual j are friends and takes value 0 otherwise; that representing the borrowing/lending network takes value 1 if individual i and individual j lend money to each other and takes value 0 otherwise.

Let q be the total number of different types of networks. Define M l n as the adjacency matrix for the l th network. The heterogeneous endogenous effects model with multiple networks is defined as

$$d _ { i } = \sum _ { l = 1 } ^ { q } \sum _ { k \in N _ { i } } d _ { k } ^ { l } \eta _ { k } ^ { l } + x _ { i } \beta _ { 0 } + \epsilon _ { i }$$

Note that in this model, different networks could potentially bear different endogenous effects for the same individual. In equation (7), coefficient η l k represents the rate of endogenous effect of individual k through network l . As a result, we have nq + k coefficients for endogenous effects. In addition, I assume endogenous effects from different types of networks are linearly additive. The model can also be rewritten in matrix form as:

$$D _ { n } = \sum _ { l = 1 } ^ { q } \left ( M _ { n } ^ { l } \circ D _ { n } \right ) \eta _ { 0 } ^ { l } + X _ { n } \beta _ { 0 } + \epsilon _ { n },$$

glyph[negationslash]

where M l n is the adjacency matrix for network l . η l = ( η , η l 1 l 2 , · · · , η l n ) ′ is an n by 1 vector for l = 1 2 , , · · · , q . Define a network l as efficient network if η l i = 0 for at least one individual i = 1 2 , , · · · , n . In this specification, leaders can only influence their neighbors through efficient networks and nonefficient networks are completely independent from the outcome variable.

## 3 Assumptions

The assumptions discussed in this section combine both standard SARs type assumptions and LASSO type assumptions. SARs type assumptions ensure the existence of valid instruments. LASSO type assumptions provide sufficient and necessary conditions for valid inference. The identification for influential individuals is achieved through its coincidence with the sparsity pattern in the reduced form estimator.

First recall the benchmark SAR model:

$$D _ { n } = \lambda _ { 0 } M _ { n } D _ { n } + X _ { n } \beta _ { 0 } + \epsilon _ { n },$$

By rearranging the above equation, we can express endogenous variable M D n n solely as a function of X n and M n , since:

$$D _ { n } = J _ { n } ^ { - 1 } X _ { n } \beta _ { 0 } + J _ { n } ^ { - 1 } \epsilon _ { n }$$

where I n is the n by n identity matrix and J n = I n -λ M 0 n . It is straightforward that J -1 n X n can serve as valid instruments for M D n n . As a result, the identification and estimation of equation (9) can be achieved through either 2SLS or GMM as proposed in papers such as Kelejian and Prucha (1995), Kelejian and Prucha (1998) Lee (2002), Lee (2003), and Lee (2004). Following the same strategy, I derive a set of instruments by solving D n as a function of exogenous variables and the adjacency matrix.

## 3.1 Assumptions for the Heterogeneous Endogenous Effects Model

Without additional restrictions, equation (4) could not be estimated through canonical method as the number of parameters n + k is greater than the number of observations n .

glyph[negationslash]

Assumption 1 (Sparsity) . Let S n ⊂ { 1 2 , , · · · , n } denotes the set of influential individuals (i.e. η j = 0 ). Let s n = | S n | be the number of elements in S n .

$$s _ { n } = o \left ( \frac { \sqrt { n } } { \log n } \right ), \quad a s \ n \to \infty$$

Assumption 1 is usually referred to as 'sparsity' assumption. The assumption that most individuals in a network are not influential is plausible under many circumstances. For example, opinion leaders on social media only constitute a very small fraction of internet users. Star programmers who can encourage other programmers to work with them on Github are also a small portion of all Github users. On the other hand, assumption 1 can be easily violated when there exists many local leaders. For example, when studying the peer effects in obesity among school children, each class/grade can be treated as a clique and as the number of cliques increases, the number of leaders may increase at a rate of O n ( ). In next subsection, I propose an extension of the current model to partially address this problem.

Assumption 2 (SAR restrictions) .

- -There exists an η max &lt; 1 such that ‖ η 0 ‖ ∞ ≤ η max

- -The glyph[epsilon1] j are i.i.d sub-Gaussian random variable with 0 mean and variance σ 2
- -The regressors x i in X n are non-stochastic and uniformly bounded for all n. lim n →∞ X X /n ′ n n exists and is nonsingular
- -The minimum eigenvalue of ( I n -M n ◦ η 0 ) , Λ min , is uniformly bounded away from 0.

The restriction on η 0 excludes the unit root process and ensures the uniqueness of equilibrium. It guarantees the invertibility of ( I n -M n ◦ η 0 ) . This is a fundamental assumption required in all SAR literatures see Upton and Fingleton (1985); Anselin (1988); Jin and Lee (2018)).

The assumption on the error term currently excludes exogenous effects and correlated effects from my model. The SAR model with this exclusion restriction is known as mixed regression model as in (see Lee, 2002). The main challenge to relax this assumption is known as 'reflection problem' as pointed out in Manski (1993). I adopt a strategy similar to Bramoull´ et al. (2009) to incorporate e both exogenous and correlated effects in next subsection. I also require the error term to be sub-Gaussian so that concentration inequalities can be derived to bound the empirical process. Sub-Gaussian process is known to have 'almost' bounded support due to the fast decay of its tails. This assumption is usually required for high dimensional estimators as in Belloni et al. (2014), Zhu (2018) and etc.

The assumption on the regressors follows the convention in the SAR literatures see Anselin (1988); Jin and Lee (2018). This deterministic design assumption can be extended to the random design cases under my model. In the following context, I focus on the case where X n is an n by 1 vector and study identification as in Bramoull´ et al. (2009). e

I also require the minimum eigenvalue of ( I n -M n ◦ η 0 ) to be uniformly bounded away from 0. This is to prevent the spatial errors from accumulating too fast. The spatial error term after the matrix inversion is ( I n -M n ◦ η 0 ) -1 glyph[epsilon1] . A similar assumption can be found in Kelejian and Prucha (1995) and Kelejian and Prucha (1998), which requires the uniformly boundedness of ( I n -M n ◦ η 0 ) -1 .

To proceed, recall the definition of the operator ' ◦ ' as M n ◦ D n = M n · diag( D n ), where diag( ) is · the diagonalization operator. Note the following property of the ' ◦ ':

$$\left ( M _ { n } \circ D _ { n } \right ) \eta _ { 0 } = \left ( M _ { n } \circ \eta _ { 0 } \right ) D _ { n }$$

If the invertibility of ( I n -M n ◦ η 0 ) is guaranteed, then

$$D _ { n } = \left ( M _ { n } \circ D _ { n } \right ) \eta _ { 0 } + X _ { n } \beta _ { 0 } + \epsilon _ { n } \Leftrightarrow D _ { n } = \sum _ { i = 0 } ^ { \infty } \left ( M _ { n } \circ \eta _ { 0 } \right ) ^ { i } ( X _ { n } \beta _ { 0 } + \epsilon _ { n } )$$

Since ( M n ◦ D n ) η 0 is correlated with glyph[epsilon1] n and η 0 is sparse (i.e. having at most s n non-zero elements), we need at least s n instruments to deal with the endogeneity in the model. Using equation (10), we can express the expectation of D n as follows:

$$E ( D _ { n } ) = X _ { n } \beta _ { 0 } + \left ( M _ { n } \circ X _ { n } \right ) ( \beta \eta _ { 0 } ) + \sum _ { i = 2 } ^ { \infty } \left ( M _ { n } \circ \eta _ { 0 } \right ) ^ { i } X _ { n } \beta _ { 0 },$$

Let ( ) · S denote the operator such that ( M n S ) is a sub matrix of M n with its columns restricted to columns corresponding to the elements of S . The first and second terms of equation (11) suggest that X n and ( M n ◦ X n S ) can serve as valid instruments to point identify β 0 and η 0 .

## Proposition 1 (First Stage Equivalence) . Under assumption 1 and 2

$$E ( D _ { n } ) = X _ { n } \beta _ { 0 } + \left ( M _ { n } \circ X _ { n } \right ) \tilde { \eta },$$

where ˜ η j = η g η , β j ˜( 0 0 , X n , M n ) for some function ˜ g depends on η 0 , β 0 , X n , and M n .

Proposition 1 shows the sparsity pattern is preserved when solving the simultaneous equations since ˜ η j = 0 as long as η j = 0 1 . As a result, the sparsity assumption is also satisfied in equation (12), and I can thus estimate equation (12) as the first stage using a LASSO type estimator.

Define W n as the projection matrix on to the orthogonal space of X n : W n = I n -X n ( X X ′ n n ) -1 X ′ n

Assumption 3 (Independence) . W M n ( n ◦ X n S ) has full column rank.

The linear independence among ( M n ◦ X n S ) requires the assumption that any two influential individuals may not necessarily connect with identical neighbors. Moreover, assumption 3 also requires that neighbors of an influential individual cannot be a linear combination of neighbors of several other influential individuals, which rules out network structures as depicted in Figure 2.

1 de Paula et al. (2015) noticed that sparsity pattern will generally not be preserved during matrix inversion when there is no further restriction on the adjacency matrix. Proposition 1 shows that with pre-existing network structure and homogeneous influence assumption for a given influential individual, sparsity pattern can still be preserved after matrix inversion.

Figure 3: Examples of networks which satisfies assumption 3

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00663v1-figure-887f3e32-65a6-4d31-8572-f63f4e7727e8.jpg)

The influence of S 1 can be identified by comparing red (right shaded) and grey groups (plain), while the influence of S 2 can be identified by comparing blue (left shaded) and grey (plain) groups. Or the influence of S 1 can be identified by comparing green (dotted) and blue (left shaded) groups, while the influence of S 2 can be identified by comparing red (right shaded) and green (dotted) groups.

S

1

. . .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00663v1-figure-de6dc611-8ca7-4ace-8b50-c67155754277.jpg)

S

2

(1)

Figure 2: Examples of networks which violate assumption 3

- (1) Two influential individual S 1 and S 2 share the exact same neighbors. (2) The neighbors of an influential individual S 1 is a linear combination of S 2 and S 3 's neighbors.

In other words, as long as each influential individual has a neighbor that is not connected with any other influential individuals, assumption 3 is satisfied. One can think of the identification here as estimating fixed effects from influential individuals as illustrated in Figure 3. Collinearity arises when the fixed effects of two influential individuals are imposed on exactly the same observations.

Assumption 3 is essentially a restriction on the topology of network structures. It rules out cases like complete network or cases (1) and (2) as in Figure 2. To achieve identification under a crosssectional network data, one has to rely on certain network structures. Similar assumptions can be found in Bramoull´ et al. (2009) and de Paula et al. (2015). e

At this point, if the truly influential individuals set S n were available to us, we would be able

to estimate the model using 2SLS method or GMM. However, in most cases, S n is not known beforehand. Notice that the identification for the set S n in the structure model (4) coincides with the sparsity pattern in the reduced form (12). I propose to use a LASSO type estimator to both recover the set of influential individuals and estimate the model. For LASSO to achieve correct recovery, I need the following assumptions:

## Assumption 4. [LASSO restrictions]

(Irrepresentable Condition) There is a ϑ ∈ (0 , 1) such that

$$\max _ { \| u \| _ { \infty } \leq 1 } \left \| d i a g ( f _ { S ^ { c } } ) \cdot \Sigma _ { 2, 1 } ^ { M } ( \Sigma _ { 1, 1 } ^ { M } ) ^ { - 1 } \cdot d i a g ( f _ { S } ) ^ { - 1 } \cdot u \right \| _ { \infty } < \vartheta$$

where and

$$f = ( I - M _ { n } \circ \eta _ { 0 } ) ^ { - 1 } X _ { n } \beta _ { 0 }$$

(Beta Min Condition) There exists N ∈ N and a m&gt; 0 such that ∀ n ≥ N ,

$$\min ( | \eta _ { 0 } | ) _ { S } \geq m / \sqrt { n },$$

Here define the operator ( ) · S as the sub-matrix/vector restricted to the columns/entries corresponding to influential individuals. Similarly, ( ) · S c represents the sub-matrix/vector restricted by the columns/entries corresponding to non-influential individuals. Also notice that the invertibility of Σ M , 1 1 is guaranteed by assumption 3.

We prove in theorem 2 that assumption 4 is a sufficient condition for the LASSO estimator to achieve a consistent selection for the set S n in the second stage. While assumption 3 restricts how influential individual may connect with their neighbors, assumption 4 restricts how non-influential individual may connect. For example, Irrepresentable Condition prevents the neighbors of a noninfluential individual to be exactly the same as the neighbors of any influential individual. This is because when two individuals connect with exactly the same neighbors, we cannot distinguish which individual is the true source of influence. This is illustrated in Figure 4 (1). On the other hand, assumption 4 does not require full independence between influential individuals and non-influential individuals. This is illustrated in Figure 4 (2).

$$\Sigma \coloneqq \frac { 1 } { n } M _ { n } ^ { \prime } W _ { n } M _ { n } = \left ( \begin{array} { c c } \Sigma _ { 1, 1 } ^ { M } & \Sigma _ { 1, 2 } ^ { M } \\ \Sigma _ { 2, 1 } ^ { M } & \Sigma _ { 2, 2 } ^ { M } \end{array} \right )$$

Figure 4: Examples of networks which violate and satisfy assumption 4

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00663v1-figure-fd649850-6506-4e00-9cbd-5b8df5484b9b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00663v1-figure-b76e0477-7ad1-42af-a69f-b38a8d8f9641.jpg)

In (1), the influence of S 1 can not be separately identified from S 4 although S 1 is influential and S 4 is non-influential. In (2), the influence of S 1 can be identified by comparing red (right shaded) and blue (left shaded) groups. And then we can find out S 4 is non-influential by comparing blue (left shaded) and green (dotted) groups.

The Beta Min Condition requires the magnitude of the endogenous effects to be sufficiently strong in order to be detected by LASSO. For example, there does not exist a sequence of individuals whose influence decay to 0 faster than the rate of 1 / √ n . Beta Min Condition is restrictive such that it rules out uniform inference and creates additional problems when constructing confidence intervals.

Equation (4) can still be consistently estimated under weaker conditions than Irrepresentable Condition and Beta Min Condition. The stronger version is assumed above to ensure selection consistency. As shown in Zhao and Yu (2006), the Irrepresentable Condition together with the Beta Min Condition are necessary and sufficient conditions for LASSO to achieve consistent model selection.

The following compatibility condition will guarantee the valid inference using 'de-bias' estimator proposed in the next section.

Assumption 4-1 (Compatibility Condition) . For some φ 0 &gt; 0 (independent of n) and for all η satisfying ‖ η S c ‖ 1 ≤ 3 ‖ η S ‖ 1 , it holds that

$$\| \eta _ { S } \| _ { 1 } ^ { 2 } \leq ( \eta ^ { \prime } ( M _ { n } ) ^ { \prime } _ { S } W _ { n } ( M _ { n } ) _ { S } \eta ) s _ { 0 } / \phi _ { 0 } ^ { 2 },$$

Assumption 4-1 is referred as compatibility condition as in van de Geer et al. (2014). It also restricts the network topology (i.e complete network is still ruled out) but is weaker than assumption 3 + assumption 4. Under assumption 4-1, we can still construct uniformly valid inference for de-bias

estimator proposed in the next section even though consistent model selection is not guaranteed.

## 3.2 Assumptions for Two Extensions

## Heterogeneous Endogenous Effects Model with Cliques:

The heterogeneous endogenous effects model with cliques can be written in matrix form as:

$$D _ { n } = \left ( M _ { n } \circ D _ { n } \right ) \eta _ { 0 } + M _ { n } D _ { n } \gamma _ { 0 } + X _ { n } \beta _ { 0 } + \epsilon _ { n }$$

Assumption 1' (Sparsity with Cliques) . Among n individuals in the network, let S n ⊂ { 1 2 , , · · · , n } be the set of global leaders. Let s n = | S n | be the number of elements in S n . Assume:

$$\dots \lambda _ { n } \quad _ { 1 } \lambda _ { n _ { 1 } } \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \\ s _ { n } = o \left ( \frac { \sqrt { n } } { \log n } \right ), \quad a s \ n \to \infty$$

Assumption 1' relaxes the exact sparsity in assumption 1 without imposing any restriction on the number of local leaders. For example, it does not rule out situations where everyone is (locally) influential. Local leaders' influence will be captured by the γ 0 , coefficient of classical spatial lag. The number of global leaders is restricted to be sparse and we can identify these leaders similar as in previous case.

To ensure invertibility of the matrix ( I n -M n ◦ η 0 -M γ n 0 ) , I modify assumption 2 as:

Assumption 2' (SAR Restrictions with Cliques) . In addition to assumption 2, there exists an η max &lt; 1 such that ‖ η 0 + γ 0 ‖ ∞ ≤ η max

Similar to assumption 2, this assumption excludes unit root processes. Since there exists a local level influence γ 0 in the network, global level influence η 0 need to be further bounded away from 1. As a result, equation (5) can be transformed into the following:

$$E ( D _ { n } ) = X _ { n } \beta _ { 0 } + \left ( M _ { n } \circ X _ { n } \right ) ( \beta _ { 0 } \eta _ { 0 } ) + M _ { n } X _ { n } ( \beta _ { 0 } \gamma _ { 0 } ) + \sum _ { i = 2 } ^ { \infty } \left ( M _ { n } \circ \eta _ { 0 } + \gamma _ { 0 } M _ { n } \right ) ^ { i } \beta _ { 0 } X _ { n } \quad \quad ( 1 3 )$$

Proposition 2 (First Stage Equivalence with Cliques) . Under assumption 1' and 2'

$$E ( D _ { n } ) = X _ { n } \beta + ( M _ { n } \circ X _ { n } ) \tilde { \eta } ^ { *, \infty } + \sum _ { i = 1 } ^ { \infty } \gamma _ { 0 } ^ { i } M _ { n } ^ { i } X _ { n } \beta _ { 0 } + \sum _ { i = 1 } ^ { \infty } M _ { n } ^ { i } ( M _ { n } \circ X _ { n } ) \tilde { \tilde { \eta } } ^ { *, ( i, \infty ) },$$

where

$$\tilde { \eta } _ { j } ^ { *, \infty } = \eta _ { 0, j } \tilde { g } _ { k, j } ^ { \infty } ( \eta _ { 0 }, \gamma _ { 0 }, \beta _ { 0 }, X _ { n }, M _ { n } ), \quad \tilde { \bar { \eta } } _ { j } ^ { *, ( i, \infty ) } = \eta _ { 0, j } \tilde { h } _ { k, j } ^ { ( i, \infty ) } ( \eta _ { 0 }, \gamma _ { 0 }, \beta _ { 0 }, X _ { n }, M _ { n } )$$

for some function ˜ g ∞ and ˜ h depend on β 0 , γ 0 , η 0 , X n , and M n

) k,j ( i, ∞ ) k,j .

From proposition 2, in order for the coefficient γ 0 introduced in equation (5) to be identified. An extra instrument M X n n need to be introduced. And thus we need to assume additional independence:

Assumption 3' (Independence with Cliques) . [ W M n ( n ◦ X n S ) , W n M X n n ] has full column rank.

To further incorporate exogenous and correlated effects, recall equation (6) can be written in matrix form as

$$D _ { n } = \left ( M _ { n } ^ { * } \circ D _ { n } \right ) \eta _ { 0 } + M _ { n } ^ { * } D _ { n } \gamma _ { 0 } ^ { e n d } + M _ { n } ^ { * } X _ { n } \gamma _ { 0 } ^ { e x o } + X _ { n } \beta _ { 0 } + \mu _ { c } + \epsilon _ { n }$$

Here, I define M ∗ n as the row sum normalized version of M n . The main challenge for the identification is known as 'reflection problem' (Manski (1993)). Bramoull´ e et al. (2009) proposed the 'neighbor's neighbor' instruments as a solution. By taking the local difference, we can obtain the following form

$$( I - M _ { n } ^ { * } ) D _ { n } = ( I - M _ { n } ^ { * } ) \Big { ( } M _ { n } ^ { * } \circ D _ { n } \Big { ) } \eta _ { 0 } + ( I - M _ { n } ^ { * } ) M _ { n } ^ { * } D _ { n } \gamma _ { 0 } ^ { e n d } + ( I - M _ { n } ^ { * } ) M _ { n } ^ { * } X _ { n } \gamma _ { 0 } ^ { e z o } + ( I - M _ { n } ^ { * } ) X _ { n } \beta _ { 0 } + ( I - M _ { n } ^ { * } ) \epsilon _ { n }$$

and invert the simultaneous equation, we have

$$E ( ( I - M _ { n } ^ { * } ) D _ { n } ) = \left ( I - M _ { n } ^ { * } \circ \eta _ { 0 } - M _ { n } ^ { * } \gamma _ { 0 } ^ { e n d } \right ) ^ { - 1 } \left ( \beta _ { 0 } \cdot I + \gamma _ { 0 } ^ { e x o } M _ { n } ^ { * } \right ) ( I - M _ { n } ^ { * } ) X _ { n },$$

Proposition 3 (First Stage Equivalence with Cliques + exogenous and correlated effects) . Under assumption 1' and 2' and assume M ∗ n is row sum normalized

$$E ( ( I - M _ { n } ^ { * } ) D _ { n } ) & = ( I - M _ { n } ^ { * } ) X _ { n } \beta _ { 0 } + ( \gamma _ { 0 } ^ { e x o } + \gamma _ { 0 } ^ { e n d } \beta _ { 0 } ) M _ { n } ^ { * } ( I - M _ { n } ^ { * } ) X _ { n } + ( M _ { n } ^ { * } \circ X _ { n } ) \tilde { \eta } ^ { * 2, \infty } \\ & + \sum _ { i = 1 } ^ { \infty } ( \gamma _ { 0 } ^ { e n d } ) ^ { i } ( \gamma _ { 0 } ^ { e x o } + \gamma _ { 0 } ^ { e n d } ) M _ { n } ^ { * ( i + 1 ) } ( I - M _ { n } ^ { * } ) X _ { n } \beta _ { 0 } + \sum _ { i = 1 } ^ { \infty } M _ { n } ^ { * i } ( M _ { n } \circ X _ { n } ) \tilde { \tilde { \eta } } ^ { * 2, ( i, \infty ) }$$

where

$$\tilde { \eta } _ { j } ^ { * 2, \infty } = \eta _ { 0, j } \tilde { g } _ { 2, k, j } ^ { \infty } ( \eta _ { 0 }, \gamma _ { 0 } ^ { e x o }, \gamma _ { 0 } ^ { e n d }, \beta _ { 0 }, X _ { n }, M _ { n } ), \quad \tilde { \tilde { \eta } } _ { j } ^ { * 2, ( i, \infty ) } = \eta _ { 0, j } \tilde { h } _ { 2, k, j } ^ { ( i, \infty ) } ( \eta _ { 0 }, \gamma _ { 0 } ^ { e x o }, \gamma _ { 0 } ^ { e n d }, \beta _ { 0 }, X _ { n }, M _ { n } )$$

for some function ˜ g ∞ 2 ,k,j and ˜ h ( i, ∞ ) 2 ,k,j depend on β 0 , γ exo 0 , γ end 0 , η 0 , X n , and M n .

glyph[negationslash]

Assumption 3'-1. Assume γ exo 0 + β γ 0 end 0 = 0 and [ W M n ( ∗ n ◦ X n S ) , W n M X ,W M ∗ n n n ∗ 2 n X ,W M n n ∗ 3 n X n ] is full rank and M ∗ is row sum normalizable.

Similar to proposition 4 in Bramoull´ e et al. (2009), assumption 3'-1 requires the independence between M X ∗ n n and M X ∗ 2 n n in order to use the 'neighbor's neighbors' as instruments to address the 'reflection' problem. Furthermore, the independence of M X ∗ 3 n n will pin down the identification for correlated effects.

## Heterogeneous Endogenous Effects Model with Multiple Networks:

The heterogeneous endogenous effects model with multiple networks can be represented in matrix form as follows:

$$D _ { n } = \sum _ { j = 1 } ^ { q } \left ( M _ { n } ^ { j } \circ D _ { n } \right ) \eta _ { 0 } ^ { j } + X _ { n } \beta _ { 0 } + \epsilon _ { n }$$

glyph[negationslash]

The number of coefficients in this model becomes nq + k and the number of observed networks q may also increase as the number of observations n increases. As a result, the sparsity assumption will be imposed on both the influential individuals and the effective networks. I assume that some of the networks are completely irrelevant (i.e. η j 0 = 0) and that relevant networks are not necessarily passing influence for everyone (i.e. η j 0 = 0 but η j 0 ,i = 0 for some i ).

Second, to ensure invertibility, for any matrix norm ‖ ‖ . :

$$\left \| \sum _ { j = 1 } ^ { q } \left ( M _ { n } ^ { j } \circ \eta _ { 0 } ^ { j } \right ) \right \| \leq \sum _ { j = 1 } ^ { q } \left \| \left ( M _ { n } ^ { j } \circ \eta _ { 0 } ^ { j } \right ) \right \| \leq \sum _ { j = 1 } ^ { q } \| \eta _ { 0 } ^ { j } \| _ { \infty } \right \| \left ( M _ { n } ^ { j } \right ) \right \|$$

Because M j n is the adjacency matrix such that each entry is 0 or 1, ∑ q j =1 ‖ η j 0 ‖ ∞ &lt; 1 guarantees the invertibility of I -∑ q j =1 ( M j n ◦ η j 0 ) .

Proposition 4 (First Stage Equivalence - Multiple Networks) . Under assumption assumption 1* and assumption 2*

$$E ( D _ { n } ) = X _ { n } \beta _ { 0 } + \sum _ { j = 1 } ^ { q } \left ( M _ { n } ^ { j } \circ X _ { n } \right ) ( \tilde { \eta } _ { 0 } ^ { j } )$$

where ˜ η j k = η g j k ˜ ( j η , β 0 0 , X n , M n ) for some function ˜ g j depends on η j 0 , β 0 , X n , and M j n .

Third, I require [ X , n ( M 1 n ◦ X n ) S , ( M 2 n ◦ X n ) S , · · · , ( M q n ◦ X n ) ] S to be full rank. Compared with the standard model, this assumption requires the independence condition to hold across different networks. The four assumptions required for multiple networks are listed formally in the appendix as assumption assumption 1*, 2*, 3*, and 4*.

## 4 Estimators and Asymptotics

The proposed estimator is similar to the two-stage least square method but use LASSO in both stages.

## Two-Stage LASSO Estimator:

## -First Stage:

$$\lambda ( \tilde { \beta }, \tilde { \eta } ) = \arg \min _ { \beta, \eta } \| D _ { n } - X _ { n } \beta - \left ( M _ { n } \circ X _ { n } \right ) \eta \| _ { 2 } ^ { 2 } + \lambda _ { 1 } | \eta | _ { 1 }$$

Obtain a LASSO fitting ˆ D n

$$\hat { D } _ { n } = X _ { n } \tilde { \beta } + \left ( M _ { n } \circ X _ { n } \right ) \tilde { \eta }$$

-Second Stage:

$$\lambda ( \hat { \beta }, \hat { \eta } ) = \arg \min _ { \beta, \eta } \| D _ { n } - \left ( M _ { n } \circ \hat { D } _ { n } \right ) \eta - X _ { n } \beta \| _ { 2 } ^ { 2 } + \lambda | \eta | _ { 1 }$$

As shown in section 3, ( M n ◦ D n ) is correlated with glyph[epsilon1] n . Thus equation (4), equation (5) and equation (8) cannot be estimated directly using LASSO or sparse group LASSO. The instruments proposed in section 3 are [ X , n ( M n ◦ X n S ) ]. We do not observe the set S but [ X , n ( M n ◦ X n )] is a set of regressors that contains the valid instruments.

## 4.1 de-bias 2SLSS Estimator

glyph[negationslash]

The estimator ˆ and ˆ suffer from LASSO shrinkage bias. Moreover, post model selection inference β η conditioning on the selected model ˆ S n = { | i η ˆ = 0 } suffers from the omitted variable bias and thus is not uniformly valid (see Leeb and Potscher, 2005, 2008, 2009). I construct a 'de-bias' estimator under my setting and derive the asymptotic distribution for it. I propose the following de-bias LASSO estimator:

## de-bias 2SLSS Estimator:

$$\hat { e } = \hat { \eta } + \hat { \Theta } \tilde { \mathcal { X } } _ { n } ^ { \prime } ( M _ { n } \circ \hat { D } _ { n } ) ^ { \prime } W _ { n } ( D _ { n } - ( M _ { n } \circ D _ { n } ) \hat { \eta } ) / n$$

$$\hat { b } = \hat { \beta } - \frac { ( X _ { n } - ( M _ { n } \circ D _ { n } ) \hat { \gamma } _ { \beta } ) ^ { \prime } \Big ( ( M _ { n } \circ ( \hat { D } _ { n } - D _ { n } ) ) \hat { \eta } \Big ) } { ( X _ { n } - ( M _ { n } \circ D _ { n } ) \hat { \gamma } _ { \beta } ) ^ { \prime } X _ { n } }$$

ˆ β and ˆ are η estimators from the 2SLSS. Define W n = ( I -X n ( X X ′ n n ) -X ′ n ) . ˜ X n = 1 n ( M n ◦ D n ) ′ W M n ( n ◦ ˆ D n ), Θ are defined by the nodewise regression as in Meinshausen and B¨hlmann (2006) ˆ u on ˜ X n and ˆ γ β are again defined by the nodewise regression on between X n and ( M n ◦ D n ). Nodewise regression explores the correlation between the columns of the design matrix ˜ X n by regressing each column on all the rest of the columns while penalizing the coefficients. An approximation of the

glyph[negationslash]

inverse of the matrix 1 n ˜ X X ′ n ˜ n can be constructed based on nodewise regression. Further, define ˆ S n = { | i η ˆ = 0 } , which represents the LASSO selected active set. The estimators (ˆ ˆ ) are adjusted e, b for the LASSO shrinkage bias and are a consistent estimator for η and β . They are similar to the estimators proposed in van de Geer et al. (2014), but are constructed through a two-stage process. The de-bias estimator also differs from the two-stage estimators proposed in Zhu (2018) due to spatial correlation.

Theorem 1. Under assumption 1, assumption 2 and assumption 4-1. There exist constant c 1 , c 2 , c 3 and ι : ‖ ‖ ι 0 &lt; ∞ such that for the first stage tuning parameter λ 1 ≥ 2 √ 4 σ c 2 1 Λ -2 min log n n and second stage tuning parameter λ ≥ 2 √ 4 σ c 2 3 Λ -2 min log n n +2 c 2 λ 1 +2 c 3 λ 1 , the de-bias estimator

$$\sqrt { n } ( \iota ^ { \prime } \hat { e } - \iota ^ { \prime } \eta _ { 0 } ) = \frac { 1 } { \sqrt { n } } \iota ^ { \prime } \hat { \Theta } \tilde { \mathcal { X } } _ { n } ^ { \prime } \epsilon - \iota ^ { \prime } \Delta \sim N ( 0, \sigma ^ { 2 } \iota ^ { \prime } \hat { \Theta } \tilde { \mathcal { X } } _ { n } ^ { \prime } \Omega _ { n } \tilde { \mathcal { X } } _ { n } \hat { \Theta } ^ { \prime } \iota )$$

$$\sqrt { n } ( \hat { b } - \beta _ { 0 } ) = \frac { ( X _ { n } - ( M _ { n } \circ D _ { n } ) \hat { \gamma } _ { \beta } ) ^ { \prime } \epsilon } { ( X _ { n } - ( M _ { n } \circ D _ { n } ) \hat { \gamma } _ { \beta } ) ^ { \prime } X _ { n } } + \Delta _ { \beta } \sim N \left ( 0, \sigma ^ { 2 } \frac { ( X _ { n } - ( M _ { n } \circ D _ { n } ) \hat { \gamma } _ { \beta } ) ^ { \prime } ( X _ { n } - ( M _ { n } \circ D _ { n } ) \hat { \gamma } _ { \beta } ) } { \left ( ( X _ { n } - ( M _ { n } \circ D _ { n } ) \hat { \gamma } _ { \beta } ) ^ { \prime } X _ { n } \right ) ^ { 2 } } \right )$$

where ‖ ∆ ‖ ∞ = o p (1) , ‖ ∆ β ‖ ∞ = o p (1) , ˜ X n = 1 n ( M n ◦ D n ) ′ W M n ( n ◦ ˆ D n ) , Ω n = 1 n ( M n ◦ ˆ D n ) ′ W M n ( n ◦ ˆ D n ) , ˆ Θ are defined by the nodewise regression on ˜ X n and ˆ γ β are defined by the nodewise regression on between X n and ( M n ◦ D n ) .

Theorem 1 shows that the 2SLSS estimator achieves normality at the standard rate √ n . The shifts ∆ and ∆ β represent the bias from using nodewise regression and they are shown to be o p (1) with the proper choice of tuning parameters.

Theorem 2. Under assumption 1-4, there exist ˜ γ &gt; 1 and ˜ = ϑ 1+ ϑ 1 -ϑ . For λ ≥ 4(˜ γ ∨ ˜ ) ϑ √ σ c 2 3 Λ -2 min log n n + 2 c 2 λ 1 +2 c 3 λ 1 ,

$$\lim _ { n \to \infty } \mathbb { P } ( \hat { S } _ { n } = S ) = 1$$

;

## 4.2 2SLSS for Two Extensions

## Two-Stage LASSO Estimator with Homogenous Effects:

- -First Stage: for a pre-specified constant k ,

$$( \tilde { \beta }, \tilde { \gamma }, \tilde { \eta } ) = \arg \min _ { \beta, \gamma, \eta } \left \| D _ { n } - X _ { n } \beta - \sum _ { i = 1 } ^ { k } M _ { n } ^ { i } X _ { n } \gamma ^ { 1, i } - \sum _ { i = 0 } ^ { k } M _ { n } ^ { i } \left ( M _ { n } \circ X _ { n } \right ) \eta ^ { 1, i } \right \| _ { 2 } ^ { 2 } + \lambda \left ( \sum _ { i = 1 } ^ { k } | \eta ^ { 1, i } | _ { 1 } + | \gamma ^ { 1, i } | \right )$$

Obtain a LASSO fitting ˆ D n

$$\hat { D } _ { n } = X _ { n } \tilde { \beta } + M _ { n } X _ { n } \tilde { \gamma } + \left ( M _ { n } \circ X _ { n } \right ) \tilde { \eta }$$

- -Second Stage:

$$( \hat { \beta }, \hat { \gamma }, \hat { \eta } ) = \arg \min _ { \beta, \gamma, \eta } \| D _ { n } - M _ { n } \hat { D } _ { n } \gamma - \left ( M _ { n } \circ \hat { D } _ { n } \right ) \eta - X _ { n } \beta \| _ { 2 } ^ { 2 } + \lambda ( | \eta | _ { 1 } + | \gamma | )$$

From proposition 2, the sparsity pattern can not be fully preserved by ( M n ◦ X n ) so additional structures like M M i n ( n ◦ X n ) need to be included in the first stage. Those terms represent the influence from global leaders but passing i th times through local leaders. By assumption 2', the influence represented by M M i n ( n ◦ X n ) decreases as i increases.

The second stage of 2SLSS with Cliques case can be viewed as a special case of the standard 2SLSS estimator. For example, one can rewrite the second stage as

$$( \tilde { \beta }, \tilde { \gamma }, \tilde { \eta } ) = \arg \min _ { \beta, \gamma, \eta } \left \| D _ { n } - X _ { n } \beta - \left ( \ M _ { n } \hat { D } _ { n } \ \left ( M _ { n } \circ \hat { D } _ { n } \right ) \ \right ) \cdot \left ( \begin{array} { c } \gamma \\ \eta \end{array} \right ) \right \| _ { 2 } ^ { 2 } + \lambda _ { 1 } | \eta | _ { 1 } + \lambda _ { 1 } | \gamma | _ { 1 }$$

As a result, asymptotics as theorem 1 and theorem 2 follows.

To incorporate the structure of multiple networks, I propose the use of sparse group LASSO.

## Two-Stage LASSO Estimator with Multiple Networks:

## -First Stage:

$$( \tilde { \beta }, \tilde { \eta } ) = \arg \min _ { \beta, \eta } \left \{ \left \| D _ { n } - X _ { n } \beta - \sum _ { j = 1 } ^ { q } ( M _ { n } ^ { j } \circ X _ { n } ) \eta ^ { j } \right \| _ { 2 } ^ { 2 } + \left ( \sum _ { j = 1 } ^ { q } \left ( \lambda _ { 1 } \| \eta ^ { j } \| _ { 2 } + \lambda _ { 2 } \| \eta ^ { j } \| _ { 1 } \right ) \right ) \right \}$$

Obtain a LASSO fitting ˆ D n

$$\hat { D } _ { n } = X _ { n } \tilde { \beta } + \sum _ { j = 1 } ^ { q } ( M _ { n } ^ { j } \circ X _ { n } ) \tilde { \eta } ^ { j }$$

## -Second Stage:

$$( \hat { \beta }, \hat { \eta } ) = \arg \min _ { \beta, \eta } \left \{ \left \| D _ { n } - X _ { n } \beta - \sum _ { j = 1 } ^ { q } ( M _ { n } ^ { j } \circ \hat { D } _ { n } ) \eta ^ { j } \right \| _ { 2 } ^ { 2 } + \left ( \sum _ { j = 1 } ^ { q } \left ( \lambda _ { 1 } \| \eta ^ { j } \| _ { 2 } + \lambda _ { 2 } \| \eta ^ { j } \| _ { 1 } \right ) \right ) \right \}$$

glyph[negationslash]

The sparse group LASSO introduces two tuning parameters, λ 1 and λ 2 , to penalize both the l 1 and the l 2 norm in each network. Similar to the LASSO estimator, the geometric shape of the penalties allows the sparse group LASSO to identify sparsity not only within each network (group) but also among networks (groups). In other words, some networks could be completely irrelevant (i.e. η j = 0) and within relevant networks, some individuals can have no influence on their neighbors (i.e. η j = 0 but η j i = 0 for some i ).

## de-bias 2SLSS Estimator under Multiple Networks:

$$\hat { e } ^ { M } = \hat { \eta } ^ { M } + \hat { \Theta } _ { Z } \tilde { \mathcal { Z } } _ { n } ^ { \prime } \hat { Z } _ { n } ^ { \prime } W _ { n } ( D _ { n } - Z _ { n } \hat { \eta } ^ { M } ) / n$$

$$\hat { b } ^ { M } = \hat { \beta } - \frac { ( X _ { n } - Z _ { n } \hat { \gamma } _ { \beta } ) ^ { \prime } ( \hat { Z } _ { n } - Z _ { n } \hat { \eta } ^ { M } ) } { ( X _ { n } - Z _ { n } \hat { \gamma } _ { \beta } ) ^ { \prime } X _ { n } }$$

ˆ and ˆ β η M 0 = (ˆ η 1 ′ , η ˆ 2 ′ , · · · , η ˆ q ′ ) ′ are estimators from the 2SLSS estimator with multiple networks. Let ˆ Z n = ( [ M 1 n ◦ ˆ D n ) , ( M 2 n ◦ ˆ D n ) , · · · ( M q n ◦ ˆ D n ) ] , Z n = ( [ M 1 n ◦ D n ) , ( M 2 n ◦ D n ) , · · · ( M q n ◦ D n ) ] . Define ˜ Z n = 1 n Z W Z ′ n n ˆ n , ˆ Θ Z are defined by the nodewise regression as in Meinshausen and B¨hlmann u (2006) on ˜ Z n and ˆ γ β are again defined by the nodewise regression on between X n and Z n .

Theorem 3. Under assumption 1*, assumption 2* and assumption 4-1*. There exists constant

d , d 1 2 , d 3 and ι : ‖ ‖ ι 0 &lt; ∞ such that for the first stage tuning parameter λ 1 1 , ≥ 2 √ 4 σ d 2 1 Λ -2 min (log n +log q ) n and second stage tuning parameter λ 2 1 , ≥ 2 √ 4 σ d 2 3 Λ -2 min (log n +log q ) n +2 d λ 2 1 1 , +2 d λ 3 1 1 , , the de-bias estimators

$$\sqrt { n } ( \iota ^ { \prime } \hat { e } ^ { M } - \iota ^ { \prime } \eta _ { 0 } ^ { M } ) = \frac { 1 } { \sqrt { n } } \iota ^ { \prime } \hat { \Theta } _ { Z } \tilde { Z } _ { n } ^ { \prime } \epsilon - \iota ^ { \prime } \Delta _ { M } \sim N ( 0, \sigma ^ { 2 } \iota ^ { \prime } \hat { \Theta } _ { Z } \tilde { Z } _ { n } ^ { \prime } \Omega _ { n } ^ { m } \tilde { Z } _ { n } \hat { \Theta } _ { Z } ^ { \prime } \iota )$$

$$\sqrt { n } ( \hat { b } ^ { M } - \beta _ { 0 } ) = \frac { ( X _ { n } - Z _ { n } \hat { \gamma } _ { \beta } ) ^ { \prime } \epsilon } { ( X _ { n } - Z _ { n } \hat { \gamma } _ { \beta } ) ^ { \prime } X _ { n } } + \Delta _ { M, \beta } \sim N \left ( 0, \sigma ^ { 2 } \frac { ( X _ { n } - Z _ { n } \hat { \gamma } _ { \beta } ) ^ { \prime } ( X _ { n } - Z _ { n } \hat { \gamma } _ { \beta } ) } { \left ( ( X _ { n } - Z _ { n } \hat { \gamma } _ { \beta } ) ^ { \prime } X _ { n } \right ) ^ { 2 } } \right )$$

where ‖ ∆ M ‖ ∞ = o p (1) , ‖ ∆ M,β ‖ ∞ = o p (1) and ˜ Z n = 1 n Z W Z ′ n n ˆ n and Ω m n = 1 n ˆ Z W Z ′ n n ˆ n . ˆ Θ Z are defined by the nodewise regression on ˜ Z n and ˆ γ β are defined by the nodewise regression on between X n and Z n .

Theorem 3 requires the rate of convergence for sparse group LASSO at the first stage. This is proved in Lemma ?? in the appendix. In group LASSO, λ 2 1 , controls the convergence of the l 2norm and thus it might be tempting to set it larger than λ 1 1 , to achieve l 2 convergence. However, group LASSO requires the number of regressors in each group to be finite. In my estimator, the number of regressors in each group equals to n . In the next theorem, I show that λ 2 1 , will need to be chosen in the same order as λ 1 1 , in order to achieve consistent selection.

Theorem 4. Let ˜ = c λ 2 1 , λ 2 2 , and ˜ ϑ mul 1 = 1+ ϑ mul 1 1 -ϑ mul 1 . Choose λ 2 1 , and λ 2 2 , such that ˜ c &gt; ϑ mul 1 1 -ϑ mul 1 . Under assumption 1*-4* and for λ 2 1 , ≥ 4(˜ γ ∨ ˜ ϑ mul 1 ) √ σ d 2 3 Λ -2 min (log n +log q ) n +2 d λ 2 1 1 , +2 d λ 3 1 1 , , there exists ˜ γ &gt; 1 ,

$$\lim _ { n \to \infty } \mathbb { P } ( \hat { S } _ { n } = S ) = 1$$

It is worth pointing out that assumption ?? is a weaker assumption than assume the full irrepresentable condition on the design matrix

$$\Sigma ^ { m u l } \coloneqq \frac { 1 } { n } [ M _ { n } ^ { 1 }, M _ { n } ^ { 2 }, \cdots, M _ { n } ^ { q } ] ^ { \prime } W _ { n } [ M _ { n } ^ { 1 }, M _ { n } ^ { 2 }, \cdots, M _ { n } ^ { q } ]$$

As the number of regressors in each group also goes to infinity as the number of groups goes to infinity, the multicollinearity in Σ mul can be severe. Instead, assumption ?? decomposes the multicollinearity into between-group and within-group and thus allows a more flexible dependence structure. As a result, sparse group LASSO can be used to recover the influential individuals as well as the effective networks that generate spillovers.

## 5 Simulations

In this section, I report Monte Carlo simulation results for the models proposed above. My results are robust when applied to networks generated by different algorithms and to networks of different sizes.

## 5.1 Heterogeneous Endogenous Effects Model

First, I use the Erdos-Renyi algorithm to simulate a network of size n . Individuals are added into the graph one at a time. When one individual is added to the network, she has probability p of generating a link with all existing individuals independently. I choose p = 0 1 and . p = 0 2 in the . simulation. I avoid a large p because collinearity among regressors may arise when links become very dense, violating assumption 3 or 4.

I set the first 5 individuals to be influential by letting their coefficients η j be non-zero. To guarantee the existence of endogenous effects, I arbitrarily specify the connections among these five individuals. The adjacency matrix M n for the five influential individuals is given in the appendix. If the connections among these five individuals are not fixed, there is a possibility that no connections are formed among these five and thus there is no endogeneity in the network. The results will be too good in such a case. The true parameters are fixed as β 0 = 3, η 0 1 , = η 0 2 , = η 0 3 , = η 0 4 , = η 0 5 , = 0 5, . and η 0 ,j = 0 for j &gt; 5. Individual characteristics X n are generated from a standard normal distribution. Individual outcomes Y n are then generated as Y n = ( I -M n ◦ η 0 ) -1 ( X β n 0 + glyph[epsilon1] n ) where glyph[epsilon1] n is drawn independently from a standard normal distribution.

I use ( M ,X ,Y n n n ) as observations and apply my two-stage LASSO estimator. I construct the debias 2SLSS estimator and repeat the above process 200 times in a manner similar to van de Geer et al. (2014). I report the average coverage probability (Avgcov) and average length (Avglength) of confidence intervals for the coefficients for influential individuals, { η , 1 · · · , η 5 } , the coefficient for individual characteristics, β 0 , and the coefficients for non-influential individuals, the η j s ( j &gt; 5). For example:

$$\ A v g { \cos } \ S _ { 0 } = s _ { 0 } ^ { - 1 } \sum _ { j \in S _ { 0 } } \mathbb { P } [ \eta _ { 0, j } \in C I _ { j } ]$$

$$\ A v g { \text{length $S_{0}=s_{0}^{-1}\sum_{j\in S_{0}} l e n g t h ( C I_{j})$} }$$

I separately report the average coverage and average length for each of the five influential individuals. As shown in appendix table A1, the coverage is around the nominal 95% level and the length

of the confidence intervals decreases as the sample size grows.

Since we can construct confidence intervals for all n coefficients, joint inference can be performed under the control of False Discover Rate (FDR). As shown in equation (18), the power reported in appendix table A1 represents the average percentage in the active set (i.e. { 1 2 3 4 5 , , , , } ) that is significant after controlling for the False Discover Rate (FDR) at 5% using the BenjaminiHochberg method. The FDR reported in appendix table A1 represents the average percentage of the non-active set (i.e. { 6 7 , , · · · , n } ) that is significant after controlling the FDR at 5% using the Benjamini-Hochberg method. The exact definition is as in equation (19).

$$\text{Power } = s _ { 0 } ^ { - 1 } \sum _ { j \in S _ { 0 } } \mathbb { P } [ H _ { 0, j } \text{ is rejected} ]$$

$$F D R \, = \sum _ { j \in S _ { 0 } ^ { c } } \mathbb { P } [ H _ { 0, j } \text{ is rejected} ] / \sum _ { j = 1 } ^ { n } \mathbb { P } [ H _ { 0, j } \text{ is rejected} ]$$

The power varies because the networks change when the sample size increases. It is strictly increasing when the network is sparse (i.e. p = 0 1). . The power decreases in the p = 0 2 case as . the problem of endogeneity increases when the network is dense. The empirical FDR is controlled well, which all under the 5% rate. Notice that the confidence interval's length is large when the sample size equals 50. This is because when the number of individuals is small, some individuals might only connect to 1 or 2 other individuals. This means that the regressors that represent this individual are all 0s except for a small numbers of non-zero terms, which leads to a large standard error.

The two-stage LASSO estimator requires the choice of two tuning parameters (i.e. the two λ s from both stages as in section 4.1). Moreover, when calculating Θ in the de-bias 2SLSS estimator ˆ (section 4.2) and using the nodewise regression, one also need to choose a tuning parameter. I use a benchmark choice of λ nodewise glyph[greaterorsimilar] √ log( n /n ) and λ 1 glyph[greaterorsimilar] √ log( n /n ) in the nodewise regression and first stage. Then I use cross-validation to pick λ in the second stage.

I further increase the number of influential individuals to 10 and report the results in appendix table A2. Again, to guarantee the existence of endogeneity, the adjacency matrix for these ten individuals is set as shown in the appendix. All average coverages and average confidence interval lengths are separately reported for these ten individuals. The choice of the tuning parameters is similar to those used to generate appendix table A1 for networks with 50 and 200 individuals. For networks with 500 individuals, I use benchmark λ to replace cross validation in the second stage.

As shown in appendix table A2, all coverages are very close to the nominal levels. The average lengths of confidence intervals is slightly larger compared with appendix table A1. This is due to the increase in influential individuals; it is more difficult to differentiate them from those irrelevant individuals.

Appendix table A3 presents the result when a network is generated using the Watts-Strogatz mechanism or the 'small world' network. Define the pN (even number) as the mean degree for each node and a special parameter ω = 0 4. . The WattsStrogatz mechanism works as follows:

- -construct a graph with N nodes each connected to pN neighbors, which pN 2 on each side.
- -For each node n i , take every edge ( n , n i j ) with i &lt; j and rewrite it with probability ω . Rewrite means replace ( n , n i j ) with ( n , n i k ) where k is choosing uniformly among all nodes that are not currently connected with n i

The influential individuals are chosen as the 1st, 5th, 15th, 40th and 50th individuals in the network. As shown in appendix table A3, my estimator is robust under a 'small world' algorithm. Nominal level is reached as the size of the network grows and the length of confidence intervals is slightly smaller than in the standard case.

## 5.2 Heterogeneous Endogenous Effects Model with Cliques

Appendix table A5 presents results for the heterogeneous endogenous effects model with cliques. The outcome variable Y n is now generated as Y n = ( I -M n ◦ η 0 -M γ n 0 ) -1 ( X β n 0 + glyph[epsilon1] n ). The coefficient of the homogeneity effects γ 0 is set at 0.05.

The choice of the tuning parameters is similar to that used to generate appendix table A1 for networks with 50 and 200 individuals. For networks with 500 individuals, I use benchmark λ (i.e. λ glyph[greaterorsimilar] √ log( n /n ) ) to replace cross validation in the second stage.

The coverage is above the 95% nominal level in all cases. I also report the mean coverage and average length of the confidence interval for the coefficient of the homogeneous effects. My model gives above 95% coverage in all cases. I also report the empirical probability of rejecting a null hypothesis of zeros effects at 95% nominal level. The probability of rejecting the test converges to 1 when the sample size grows to 500.

## 5.3 Heterogeneous Endogenous Effects Model with Multiple Networks

In this Monte Carlo exercise, I include two different networks generated by the Erdos-Renyi algorithm, where one is influential and the other is not. I use the two-stage LASSO estimator with multiple networks to estimate the parameters. The sparse group LASSO requires two tuning parameters, one for the l 2 norm and the other for the l 1 norm. I set the two parameters to be equal to each other as the correlations among the columns of the adjacency matrices are very small. The choice of tuning parameters is similar to that used to generate appendix table A1 for networks with 50 and 200 individuals. For networks with 500 individuals, I use benchmark λ instead of cross-validation in the second stage. Appendix table A4 summarizes the results. As in previous results, all coverages exceed the nominal 95% level.

I report the empirical probabilities such that at least one individual is detected in a given network controlling for the FDR at 5% using the Benjamini-Hochberg method. I also report the average number of detections conditioning on at least one individual who is detected in a given network. Appendix table A4 shows that network 1, which is the relevant network, is more likely to be detected in all cases than network 2, the irrelevant network. The average number of identified individuals for network 1 is also more than that of network 2.

## 6 Empirical Application

I use the proposed estimator to study the importance of different networks in spreading the participation in a micro finance program within rural Indian villages. I show that different kinds of networks have different effects on individuals decisions. I identify the influential individuals in each village. My analysis shows that leaders among agricultural laborers, Anganavadi teachers, construction workers, small business owners and mechanics are very likely to be influential in the villages.

## 6.1 Background

A non-profit organization named Bharatha Swamukti Samsthe (BSS) has been running micro finance programs in rural southern Karnataka, India since 2007. It provides small loan products to poor women and, through them, to their families. The villages covered by the program are geographically isolated and heterogeneous in terms of caste.

When BSS initially introduces a micro finance program to a village, the credit officers of BSS first approached a number of 'predefined leaders', such as teachers, shopkeepers and village elders. BSS held a private meeting with these leaders and explained the program. Then these predefined leaders passed the information onto other villagers. Those who were interested in the program and contacted BSS were trained and assigned to groups to receive credit. Each group consisted of 5 borrowers and group members were jointly liable for loans. Loans were around 10,000 rupees (approximately $200) at an annualized rate of approximately 28%. Note that 74.5 percent of the households in rural area said the monthly income of their highest earning member is less than 5,000 rupees (source: Socio-Economic Caste Census-2011). This loan had to be repaid within 50 weeks.

In 2006, 75 villages in Karnataka were surveyed 6 months before the initiation of the BSS micro finance program. This survey consisted of a village questionnaire and a detailed follow-up survey conducted among a subsample of villagers. The village questionnaire gathered demographic information on all households in a village including GPS coordinates, age, gender, number of rooms, whether the house had electricity, and whether the house had a latrine. The data set also contains information on the 'pre-defined leaders' set who helped spread the information to the entire village. The follow-up survey collected data from a villager sample stratified according to age, education level, caste, occupancy, etc. It also asked questions about social network structures along 12 dimensions, including:

- -Friends: Name the 4 non-relatives whom you speak to the most.
- -Visit-go: In your free time, whose house do you visit?
- -Visit-come: Who visits your house in his or her free time?
- -Borrow-kerorice: If you needed to borrow kerosene or rice, to whom would you go?
- -Lend-kerorice: Who would come to you if he/she needed to borrow kerosene or rice?
- -Borrow-money: If you suddenly needed to borrow Rs. 50 for a day, whom would you ask?
- -Lend-money: Who do you trust enough that if he/she needed to borrow Rs. 50 for a day you would lend it to him/her?
- -Advice-come: Who comes to you for advice?
- -Advice-go: If you had to make a difficult personal decision, whom would you ask for advice?
- -Medical-help: If you had a medical emergency and were alone at home whom would you ask for help in getting to a hospital?

- -Relatives: Name any close relatives, aside from those in this household, who also live in this village.
- -Temple-company: Do you visit a temple/mosque/church? Do you go with anyone else? What are the names of these people?

For the 43 villages where micro finance was introduced by the time of 2011, BSS also collects information on which villagers have joined the program. These survey questions reveal the underlying structures for connections among any two individuals in the network. Figure 5 presents all those connections at the household-level in a graph. Each node in the graph represents a household. A black node indicates that the household joined the micro finance program, while a white node indicates that it did not. Bigger nodes represent those households in which at least one family member has been chosen as being among the 'pre-defined leaders'. An edge between two nodes signifies that the two nodes are connected in at least one of the 12 networks. The darker the color of the edge, the more connections it represents.

This dataset provides an ideal framework for application of the heterogeneous endogenous effects model. First, it allows me to model endogenous effects. An individual may decide to join the micro finance program if her neighbors or friends plan to join. Second, the endogenous effects are individual specific. Given the diversity of the villagers, it is possible that some villagers are more influential than others. Third, it allows me to implement the heterogeneous endogenous effects model with multiple networks. The questions asked regarding multiple dimensions of the network structure allow me to explore which network is most influential.

## 6.2 Data

In this empirical study, I focus on the 38 villages that have been introduced to the micro finance programs by BSS. For each village, I can observe both its social network structure and the villagers' decisions about joining the program. I drop the data for one village (Village 46) that contains incorrect entries on the index of households. Appendix table ?? summarizes the descriptive statistics for each village.

Among the 12 questions about the social network structure, 4 pairs essentially capture the same connections among the villagers 2 . Therefore, I consolidate each pair of questions into one dimension:

2 Assuming every villager truthfully answers a pair of questions, the adjacency matrices associated with each question are the same. It is also plausible to treat villagers' answers to each question as a separate directed graph.

Figure 5: Network in Village 1

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00663v1-figure-a5b2bbb9-8048-4f56-8369-2e0704b0506f.jpg)

| Visit-go-come        |    In your free time, whose house do you visit? Who visits your house in his or her free time?                                                                                    |
|----------------------|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Borrow-Lend-kerorice |    If you needed to borrow kerosene or rice, to whom would you go? Who would come to you if he/she needed to borrow kerosene or rice?                                             |
| Borrow-Lend-money    |          If you needed to borrow Rs.50 for a day, whom would you ask? Who do you trust enough that if he/she needed to borrow Rs.50 for a day you would lend it to him/her? |
| Help decision        |         Who comes to you for advice? If you had to make a difficult personal decision, whom would you ask for advice?                                                        |

I restructure all the data at the household level as only women are allowed to apply for the micro finance program because the goal of BSS is to support families through the women in them. As a result, a woman's decision to join or not join the micro finance program becomes her family's decision. A connection between two villagers becomes a connection between two families. A 'predefined leader' is a villager selected by BSS to help spread information about the micro finance program to the other villagers. At the household level, I use the term 'predefined leader' for a household that contains at least one such villager.

## 6.3 Sparsity and Equilibrium

To demonstrate how my method identifies influential households, I model families' decisions regarding joining the micro finance program as a network game with Bayesian Nash Equilibrium. For household i , let d ∗ i be the expected probability that i chooses to join the micro finance program. The decision of household i depends on its neighbors' decisions as well as the types of connections between them. The decision also depends on its characteristics X i and on unobserved information glyph[epsilon1] i . Formally, it can be written as:

$$d _ { i } ^ { * } = \sum _ { l \in N _ { i } } d _ { l } ^ { * } ( \sum _ { j = 1 } ^ { q } \eta _ { l } ^ { j } ) + x _ { i } \beta + \epsilon _ { i }$$

However, these questions do not allow for clear determination of the directions. For example, if villager A visits villager B 's house, it is not clear whether villager A influences villager B or vice versa

Rewritten in matrix form:

$$\colon \\ D _ { n } ^ { * } = \sum _ { j = 1 } ^ { q } \left ( M _ { n } ^ { j } \circ D _ { n } ^ { * } \right ) \eta ^ { j } + X _ { n } \beta + \epsilon _ { n }$$

I assume that only a small number of households are influential over their neighbors. Leaders and followers are usually observed in those rural villages. Big decisions are often made by the village elders or by the more educated among the villagers. BSS recognized the importance of leaders and gathered a group of predefined leaders, asking them to inform the rest of the villagers about their program. I do not consider the local level influence in these villages given the size and how complicated the network structures are. Households are closely connected by these 8 networks as shown in Figure 5 and there is no form of clique visible.

Because the villages are considered geographically isolated, I apply my estimator separately to each of the 38 villages. I use the number of rooms per person in a household as the independent variable X n . Number of rooms per person is a proxy for the wealth in the family. As shown in table 1, it is negatively correlated with the decision to join the micro finance program. The richer the family, the less likely the family is to participate in the micro-finance program. On the other hand, it is arguably an exogenous variable to other factors (career, education, etc) that might affect the decision to join the micro-finance program. I further check the robustness of my independent variable by including additional controls. The adjacency matrix M j n is constructed from the questions in the survey. Households i and k are connected in network j if either i or k reported the other in question j . Finally, d ∗ i is replaced with the household's decision.

The instruments are constructed as ( M j n ◦ X n ) for j = 1 2 , , · · · , 8. I use the heterogeneous endogenous effects model with multiple networks to: 1) Identify the effective networks affecting a household's decision and 2) Identify that households that are leaders in the village and study the association between observable characteristics and leader status. If a new program is going to try to recruit these households, the organizers can target those influential households and try to persuade them to join first.

## 6.4 Results

## 6.4.1 Identifying Effective Networks

First, I study how LASSO selects networks. I define a coefficient for a household's endogenous effect in a network as significant according to two different criteria. The first criterion, 'Cross-Validation',

Table 1: Predictive Power of Characteristics X n

|                           | (1) Participate   | (2) Participate   | (3) Participate   |
|---------------------------|-------------------|-------------------|-------------------|
| Average Num. rooms x100   | -8.19 ∗∗∗         | -7.12 ∗∗∗         | -3.36 ∗∗∗         |
| Household Size x100       | (1.30)            | (1.30) 0.45 ∗∗    | (1.12) 0.40 ∗∗    |
| Electricity x100          |                   | (0.21)            | (0.20)            |
|                           |                   |                   | 1.29              |
| Latrine x100              |                   |                   | -5.66 ∗∗∗         |
| Average Num. workers x100 |                   |                   | (1.38) 0.64 ∗     |
| Average age x100          |                   |                   | (0.50) -0.27 ∗∗∗  |
| Village Fixed Effects     | Y                 | Y                 | Y 8,375           |
| n                         | 8,375             | 8,375             |                   |
| # of villages             | 37                | 37                | 37                |
| R 2                       | 0.05              | 0.05              | 0.06              |

Table 1 provides a robustness check for variable X n : Average Num. rooms, which used to construct instruments. Standard errors in parentheses * p &lt; 0 1, ** . p &lt; 0 05, *** . p &lt; 0 01. . Standard deviation clustered at village level. Dependent variable is households' decision on whether to join the micro finance program or not. All design control village fixed effects.

determines a coefficient to be significant if LASSO predicts the coefficient to be non-zero after cross-validation. The second criterion, 'de-bias', first constructs a bias-adjusted coefficient and calculates its standard error. It then determines a coefficient to be significant if the BenjaminiHochberg method rejects the null hypothesis of zero effect at the 5% false discovery rate. A network is defined as significant if at least one coefficient for heterogeneous endogenous effects in this network is significant.

Table 2 presents the empirical probability of the 8 networks being significant among the 37 villages. Note that certain types of networks (such as visit go-come) are more likely to pass influence then others (such as temple company). For example, by cross-validation criterion, the visit go-come network is detected as significant in 26 out of the 37 villages (i.e. 70%) whereas temple company is detected as significant in none of the 37 villages (i.e. 0%). I also present the average number of households associated with significant endogenous effects in each significant network. For example, according to the cross-validation criterion, an average of 101 households are detected having significant coefficients associated with the visit go-come network among the 26 village. The average of the absolute magnitude of the heterogeneous endogenous effects is also reported.

In terms of variable selection, if Assumption 4 holds, the cross-validation criterion may consistently select the truly influential households with high probability even in a finite sample. On the other hand, the de-bias criterion is likely to be conservative because of its use of the false discovery control process. In terms of coefficients estimated, de-bias estimators are asymptotically consistent. On the other hand, estimates based on the LASSO estimator suffer from shrinkage bias and are not consistent.

Table 2: Influential Networks

|            |               | visit go-come   | friendship   | borrow-lend keroric   | borrow-lend money   | relatives   | help decision   | medical help   | temple company   |
|------------|---------------|-----------------|--------------|-----------------------|---------------------|-------------|-----------------|----------------|------------------|
| Cross 1    | probability 3 | 70%             | 54%          | 38%                   | 49%                 | 22%         | 30%             | 14%            | 0%               |
| Validation | identified 4  | 101             | 89           | 88                    | 86                  | 69          | 70              | 80             | 0                |
|            | probability 3 | 54%             | 35%          | 30%                   | 19%                 | 16%         | 8%              | 5%             | 0%               |
| De-bias 2  | identified 4  | 11              | 8            | 8                     | 9                   | 7           | 14              | 8              | 0                |
| De-bias 2  | magnitudes 5  | 0.0166          | 0.0095       | 0.0070                | 0.0061              | 0.0050      | 0.0040          | 0.0028         | 0.0018           |

Table 2 reports the probability of detection for different networks among the 38 villages. A network is detected as influential if at least one leader is detected within this network. 1. Cross Validation represents those networks detected by lasso using cross validation. 2. De-bias represents those networks detected by significant de-bias estimators under FDR control. 3. Probability reports the empirical probability that at least one coefficient ˆ e j i is significant in network j . 4. Identified reports the averaged number of significant ˆ e j i in the network j conditioning on the network being detected. 5. Magnitudes is the mean of | ˆ e j i | and represents the average endogenous effects through network j

The results in Table 2 suggest villagers are more likely to discuss the micro-finance program when they visit each other, chat with friends, and meet with people to whom they are economically connected. Villagers are not likely to talk about the micro finance program when they go to the temple or need medical help.

## 6.4.2 Identifying Influential Households

Second, I focus on how LASSO selects households. I compare the LASSO selected influential households with the BSS selected 'predefined leaders'. It is important to point out that these 'predefined leaders' are not necessarily influential villagers in a network. Recall that predefined leaders are a set of villagers that BSS select to help spread the information about the micro finance program. The fact that a villager is selected as a 'predefined leader' to pass information about the micro finance program does not a priori guarantee her or her family's influence - her decision to join the micro finance program may not lead to her neighbors' decisions to join. In the analyses below, I will examine how influential villagers are associated with 'predefined leaders' and explore their potential differences.

## 1. Influential Predefined Households

In table 3, I report results indicating that influential households selected by LASSO partly overlap with 'predefined leaders'. This is intuitive because some 'predefined leaders' such as school headmasters and village elders are highly respected figures in a village. Therefore, their decisions are likely to be followed by others in the village. On average, BSS selected 27 villagers as 'predefined leaders' in each village. In comparison, Cross-Validation criterion selects around 136 villagers and de-bias criterion selects around 19. Furthermore, on average, 19 out of 136 influential villagers (i.e. 14%) selected by Cross-Validation criterion are also BSS 'predefined leaders'; 3 out of 19 influential villagers (i.e. 17%) selected by de-bias criterion are also BSS 'predefined leaders'. The likelihood of selected by the two methods are both higher than the percentage of predefined leaders in the entire village (11%). Comparing with a random guess of influential individuals, table 3 suggests the LASSO detected influential individuals are more likely to overlap with the predefined leaders. In Table 6 below, I show that small business owners are more likely to be both influential and selected as 'predefined leaders'.

## 2. Influential Non-Predefined Households

In this and the following section, I focus on understanding the differences between the influential households selected by LASSO and the 'predefined households' selected by BSS. I investigate the likelihood that a household being selected by LASSO or by BSS, as associated with the careers of

Table 3: Coverage of predefined leaders

|                    | % of predefined leaders among:   | % of predefined leaders among:   | Average number of discovery 3   |
|--------------------|----------------------------------|----------------------------------|---------------------------------|
|                    | LASSO detected 1                 | entire village 2                 | Average number of discovery 3   |
| Cross Validation 4 | 14%                              | 11%                              | 136                             |
| De-bias 5          | 17%                              | 11%                              | 19                              |

Table 3 depicts the overlapping between influential households selected by LASSO and 'predefined leaders'. Predefined leaders are a set of villagers defined by BSS, who helped spread the information about the micro-finance program. 1. LASSO detected reports the percentage of households detected by LASSO and also selected as 'predefined leaders' in total LASSO detected households. 2. Entire village reports the percentage of 'predefined leaders' among the entire village. 3. Average number of discovery reports the total number of individuals discovered by LASSO using each method. 4. Cross Validation represents those individuals identified from lasso using cross validation. 5. De-bias represents those individuals identified from significant de-biased estimators controlling FDR. 6. The average number of predefined leaders in one village is 27.

its family members. Table 4 and 5 present linear regression results using career dummy variables of family members to explain whether a household is selected as 'predefined leader' (Column 1 in table 4), whether a household is selected by LASSO as influential using cross-validation (Column 2 in table 4), and whether a household is selected by LASSO as influential using de-bias estimator (Column 3 in table 4). The full results of these regressions are reported in appendix table A7.

Table 4 summarizes all careers that have a significant impact ( p &lt; 0 001) on the likelihood of a . household being selected by LASSO as influential using de-bias method. Agriculture is the backbone of those rural villages and over 67% of the villagers are agricultural laborers. The prior that there exists influential individuals among agricultural laborers coincides with the LASSO detection. Anganwadi teacher is a set of groups that provides pre-school education to the children. They are part of the government's health care system in the rural areas. Note that except for small business owners and Anganavadi teacher all the other careers in this table are either not significantly or negatively associated with the likelihood of a household being selected as 'predefined leaders'. Compare with Anganavadi teacher and small business owners, regular teachers and animal skin business owners are not targeted by BSS as 'predefined leaders' but LASSO estimators indicate they may also become influential. Other careers that are correlated with LASSO selection include construction/mud worker, truck/tractor driver, factory worker, daily labor and wood cutter. Industrial workers are likely to form circles of collaboration due to the natural of their jobs and thus it seems compelling that they are selected as influential individuals.

Table 5 summarizes all careers that have a significant impact ( p &lt; 0 01) on the likelihood of a . household being selected by BSS as being among the 'predefined leaders'. Poojari are Indian priests in those villages and they are very likely to be included as 'predefined leaders'. However, they are not likely to influence people to join the micro finance program. Other careers as tailor,

hotel workers, veteran and barber are included as 'predefined leaders' because individuals doing these jobs can spread information quickly in the village. However, LASSO does not find these individuals to be influential.

Table 6 reports the counter factual study when selected leaders all decide to join the micro-finance program. The participation rate for non-leaders in the data is 16%. When all 'predefined leaders' decide to join, the participation rate for non-leaders will increase to 20%. And when all LASSO selected leaders decide to join, the participation rate for non-leaders will further increase to 33%.

Table 4: Second Stage: LASSO selected leaders' careers

|                                     | Predefined   | Selected by LASSO   | Selected by LASSO   |
|-------------------------------------|--------------|---------------------|---------------------|
|                                     | leaders      | Cross-validation    | De-bias             |
| Agriculture labour                  | 0.00         | 0.31 ∗∗∗            | 0.06 ∗∗∗            |
|                                     | (0.01)       | (0.01)              | (0.00)              |
| Anganavadi Teacher                  | 0.14 ∗       | 0.11                | 0.13 ∗∗∗            |
|                                     | (0.06)       | (0.07)              | (0.04)              |
| Construction/mud worker             | 0.00         | 0.17 ∗∗∗            | 0.15 ∗∗∗            |
|                                     | (0.02)       | (0.03)              | (0.02)              |
| Truck/Tractor Driver                | -0.03        | 0.16 ∗∗∗            | 0.08 ∗∗∗            |
|                                     | (0.03)       | (0.03)              | (0.02)              |
| Factory worker (bricks/stones/mill) | -0.00        | 0.17 ∗∗∗            | 0.07 ∗∗∗            |
|                                     | (0.02)       | (0.03)              | (0.01)              |
| Small business                      | 0.22 ∗∗∗     | 0.18 ∗∗∗            | 0.05 ∗∗∗            |
|                                     | (0.02)       | (0.03)              | (0.01)              |
| Teacher                             | 0.05         | 0.22 ∗∗∗            | 0.09 ∗∗∗            |
|                                     | (0.04)       | (0.05)              | (0.03)              |
| Daily labourer                      | -0.05 ∗      | 0.16 ∗∗∗            | 0.08 ∗∗∗            |
|                                     | (0.03)       | (0.03)              | (0.02)              |
| Wood cutter                         | -0.03        | 0.15 ∗              | 0.12 ∗∗∗            |
|                                     | (0.06)       | (0.07)              | (0.04)              |
| Animal skin business                | 0.36         | 0.62 ∗              | 0.52 ∗∗∗            |
|                                     | (0.23)       | (0.28)              | (0.15)              |
| Control other careers               | Y            | Y                   | Y                   |
| Control village fix effect          | Y            | Y                   | Y                   |

Table 4 summarizes all careers that have a significant impact ( p &lt; 0 001) on the likelihood of a household being selected by . LASSO de-bias method from appendix table ?? . The first column uses whether one is predefined leaders as response variable, the second column uses whether one joins the micro-finance program as response variable and the third column uses whether one is selected by lasso as response variable. Standard errors in parentheses * p &lt; 0 1, ** . p &lt; 0 01, *** . p &lt; 0 001. .

Table 5: Second Stage: predefined leaders' careers

|                            | Predefined leaders   | Selected by LASSO   | Selected by LASSO   |
|----------------------------|----------------------|---------------------|---------------------|
|                            |                      | Cross-validation    | De-bias             |
| Small business             | 0.22 ∗∗∗             | 0.18 ∗∗∗            | 0.05 ∗∗∗            |
|                            | (0.02)               | (0.03)              | (0.01)              |
| Tailor Garment worker      | 0.08 ∗∗              | 0.13 ∗∗∗            | 0.03                |
|                            | (0.03)               | (0.04)              | (0.02)              |
| Hotel worker               | 0.32 ∗∗∗             | 0.24 ∗∗             | -0.02               |
|                            | (0.07)               | (0.09)              | (0.05)              |
| Poojari                    | 0.53 ∗∗∗             | 0.15                | -0.05               |
|                            | (0.13)               | (0.16)              | (0.09)              |
| Veterinary clinic          | 0.87 ∗∗              | -0.03               | 0.00                |
|                            | (0.32)               | (0.39)              | (0.21)              |
| Barber/saloon              | 0.41 ∗∗∗             | -0.01               | -0.02               |
|                            | (0.10)               | (0.12)              | (0.06)              |
| Control other careers      | Y                    | Y                   | Y                   |
| Control village fix effect | Y                    | Y                   | Y                   |

Table 5 summarizes all careers that have a significant impact (( p &lt; 0 01)) on the likelihood of a household being selected by . BSS as being among the 'predefined leaders' from appendix table ?? . The first column uses whether one is predefined leaders as response variable, the second column uses whether one joins the micro-finance program as response variable and the third column uses whether one is selected by lasso as response variable. Standard errors in parentheses * p &lt; 0 1, ** . p &lt; 0 01, *** . p &lt; 0 001. .

Table 6: Participation Rate when Targeting Different Leaders

|                                  | In data   | Predefined Leaders   | LASSO Leaders   |
|----------------------------------|-----------|----------------------|-----------------|
| Participation Rate (non-leaders) | 16%       | 20%                  | 33%             |

Table 6 reports the participation rate of non-leaders when all targeted leaders decided to join. The true participation rate in data is 16%. If all predefined leaders decide to join, the participation rate will increase to 20%. If all LASSO detected leaders decide to join, the participation rate will increase to 33%.

## 7 Conclusions

In this paper, I propose a novel spatial autoregression model which allows for heterogeneous endogenous effects. Specifically, each individual has an individual-specific endogenous effect on her neighbors. My approach is useful for modeling a network with leaders and followers.

I propose a set of instruments as well as a two stage LASSO (2SLSS) method to estimate my model. The instruments are constructed as a function of the independent variables and an adjacency matrix. I use a LASSO type estimator to select the valid instruments in the first stage and the influential individuals in the second stage. I propose a bias correction for my two-stage estimator following van de Geer et al. (2014). I derive the asymptotic normality for my 'de-bias' two-stage LASSO estimator and conduct robust inference including confidence intervals.

My model can be extended to allow for more flexible structures. To apply LASSO, I assume that the number of influential individuals is sparse. I propose heterogeneous endogenous effects model with cliques to incorporate locally influential individuals, where the sparsity assumption is only applied to globally influential individuals. My model can also be extended to situations where there are multiple networks. I propose the use of the sparse group LASSO in my 2SLSS process. I derive the convergence rate and prove the consistency of selection for the sparse group LASSO estimator.

I apply my method to study villagers' decisions to participate in micro-finance programs in rural areas of Indian. I show that leaders in those villages have significant influence over their neighbors' decision to join the micro-finance program, and I provide rankings for the different social and economic networks among villagers. Based on how effectively each network spreads the impact of influential individuals' decisions, my method shows that some networks such as 'visit go-come' and 'borrow money' are much more effective in influencing villagers' decisions than other networks such as 'temple company' and 'medical help'. I further show that individuals from certain careers such as agricultural workers, Anganwadi teachers and small business owners are more likely to influence other villagers and the 'predefined leaders' selected by BSS are different than the LASSO detected influential individuals.

## References

Acemoglu, D., Garc´ ıa-Jimeno, C., and Robinson, J. A. (2012). Finding eldorado: Slavery and long-run development in colombia. NBER WORKING PAPER SERIES.

| Ammermuller, A. and Pischke, J.-S. (2009). Peer effects in european primary schools: Evidence from pirls. Journal of Labor Economics , 27(3):315-348.                                                            |
|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Anselin, L. (1988). Spatial Econometrics: Methods and Models . Boston: Kluwer.                                                                                                                                   |
| Ballester, C., Calv´-Armengol, o A., and Zenou, Y. (2006). Who's who in networks. wanted: The key player. Econometrica , 74(5):1403-1417.                                                                        |
| Bandiera, O., Barankay, I., and Rasul, I. (2009). Social connections and incentives in the workplace: Evidence from personnel data. Econometrica , 77(4):1047-1094.                                              |
| Banerjee, A., Chandrasekhar, A., Duflo, E., and Jackson, M. (2013). The diffusion of microfinance. Science , 341(6144).                                                                                          |
| Belloni, A., Chernozhukov, V., and Hansen, C. (2014). Inference on treatment effects after selection amongst high-dimensional controls. The Review of Economic Studies , 81(2):608-650.                          |
| Belloni, A., Chernozhukov, V., and Kato, K. (2015). Uniform post selection inference for lad regression and other z-estimation problems. Biometrika , 102:77-94.                                                 |
| Blume, L. E., Brock, W. A., Durlauf, S. N., and Jayaraman, R. (2015). Linear social interactions models. Journal of Political Economy , 123(2):444-496.                                                          |
| Bonaldi, P., Hortacsu, A., and Kastl, J. (2015). An empirical analysis of funding costs spillovers in the euro-zone with application to systemic risk. NBER Working Paper.                                       |
| Bramoull´, e Y., Djebbari, H., and Fortin, B. (2009). Identification of peer effects through social networks. Journal of Econometrics , 150(1):41-55.                                                            |
| B¨hlmann, u P. (2013). Statistical significance in high-dimensional linear models. Bernoulli , 41(2):802-837.                                                                                                    |
| B¨hlmann, u P. and van de Geer, S. (2011). Statistics for High-Dimensional Data . Springer.                                                                                                                      |
| Bunea, F., Lederer, J., and She, Y. (2014). The square root group lasso: theoretical properties and fast algorithms. IEEE-Information Theory , 60:1313-1325.                                                     |
| Calv´-Armengol, o A., Patacchini, E., and Zenou, Y. (2009). Peer effects and social networks in education. Review of Economic Studies , 76(4):1239-1267.                                                         |
| Chernozhukov, V., Chetverikov, D., Demirer, M., Duflo, E., Hansen, C., Newey, W., and Robins, J. (2018). Double/debiased machine learning for treatment and causal parameters. Econometrics Journal , 21:C1-C68. |

Christakis, N. A., Fowler, J. H., Imbens, G. W., and Kalyanaraman, K. (2010). An empirical model for strategic network formation. NNBER working paper.

Clark, A. E. and Loheac, Y. (2007). 'it wasn't me, it was them!' social influence in risky behavior by adolescents. Journal of Health Economics , 26:763-784.

Cliff, A. and Ord, J. K. (1973). Spatial autocorrelation . London: Pion.

- Coelli, T., Rahman, S., and Thirtle, C. (2002). Technical, allocative, cost and scale efficiencies in bangladesh rice cultivation: A nonparametric approach. Journal of Agricultural Economics , 53(3):607-626.
- Conley, T. G. and Udry, C. R. (2010). Learning about a new technology: Pineapple in ghana. AMERICAN ECONOMIC REVIEW , 100(1):35-69.

Cressie, N. A. C. (1993). Statistics for Spatial Data . John Wiley &amp; Sons, Inc.

- de Paula, A., Rasul, I., and Souza, P. C. (2015). Recovering social networks from panel data: identification, simulations and an application.
- Denbee, E., Julliard, C., Li, Y., and Yuan, K. (2015). Network risk and key players: A structural analysis of interbank liquidity.

Fan, J. and Liao, Y. (2014). Endogeneity in high dimensions. The Annals of Statistics , 42(3):872917.

- Gautier, E. and TsyBakov, A. B. (2014). High-dimensional instrumental variables regression and confidence sets. TSE Working Paper.
- Guryan, J., Kroft, K., and Notowidigdo, M. J. (2009). Peer effects in the workplace: Evidence from random groupings in professional golf tournaments. American Economic Journal: Applied Economics , 1(4):34-68.
- Horrace, W. C., Liu, X., and Patacchini, E. (2016). Endogenous network production functions with selectivity. Journal of Econometrics , 190(2):222-232.
- Javanmard, A. and Montanari, A. (2014). Confidence intervals and hypothesis testing for highdimensional regression. Journal of Machine Learning Research , 15(1):2869-2909.
- Jin, F. and Lee, L.-F. (2018). Lasso maximum likelihood estimation of parametric models with singular information matrices. Econometrics , 6(1):1-24.

Kelejian, H. H. and Prucha, I. R. (1995). A generalized moments estimator for the autoregressive parameter in a spatial model. INTERNATIONAL ECONOMIC REVIEW , 40.

Kelejian, H. H. and Prucha, I. R. (1998). A generalized spatial two-stage least squares procedure for estimating a spatial autoregressive model with autoregressive disturbances. Journal of Real Estate Finance and Economics , 17(1):99-121.

Krauth, B. V. (2005). Peer effects and selection effects on smoking among canadian youth. Canadian Journal of Economics , 38(3):735-757.

Lee, L. (2002). Consistency and efficiency of least squares estimation for mixed regressive, spatial. Econometric Theory , 18(2):252-277.

Lee, L. (2003). Best spatial two-stage least squares estimators for a spatial autoregressive model with autoregressive. Econometric Reviews , 22(4):305-335.

Lee, L. (2004). Asymptotic distributions of quasi-maximum likelihood estimators for spatial econometric models. Econometrica , 72:1899-1926.

Lee, L. and Liu, X. (2010). Efficient gmm estimation of high order spatial autoregrssive models with autoregressive disturbances. Econometric Theory , 26:187-230.

Lee, L.-f. and Yu, J. (2010). A spatial dynamic panel data model with both time and individual effects. Econometric Theory , 26:564-597.

Leeb, H. and Potscher, B. M. (2005). Model selection and inference: facts and fiction. Econometric Theory , 21(1):21-59.

Leeb, H. and Potscher, B. M. (2008). Can one estimate the unconditional distribution of postmodelselection estimators? Econometric Theory , 24(2):38-376.

Leeb, H. and Potscher, B. M. (2009). Model selection. Handbook of Financial Time Series , pages 889-925.

Manresa, E. (2013). Estimating the structure of social interactions using panel data.

Manski, C. (1993). Identification of endogenous social effects: The reflection problem. The Review of Economic Studies , 60(3):531-542.

Mas, A. and Moretti, E. (2009). Peers at work. American Economic Review , 99(1):112-145.

Masten, M. A. (2018). Random coefficients on endogenous variables in simultaneous equations models. The Review of Economic Studies , 85(2):1193-1250.

Meinshausen, N. and B¨hlmann, P. (2006). u High-dimensional graphs and variable selection with the lasso. The Annals of Statistics , 34(1436-1462).

| Nakajima, R. (2007). Measuring peer effects on youth smoking behaviour. The Review of Economic Studies , 74(3):897-935.                                                                            |
|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Neidell, M. and Waldfogel, J. (2010). Cognitive and noncognitive peer effects in early education. Review of Economics and Statistics , 92(3):562-576.                                              |
| Pinkse, J., Slade, M., and Brett, C. (2002). Spatial price competition: a semiparametric approach. Econometrica , 70(3):1111-1153.                                                                 |
| Sacerdote, B. (2001). Peer effects with random assignment: Results for dartmouth roommates. The Quarterly Journal of Economics , 116(2):681-704.                                                   |
| Sheng, S. (2016). A structural econometric analysis of network formation games through subnet- works. Conditional Acceptance at Econometrica.                                                      |
| Simon, N., Friedman, J., Hastie, T., and Tibshirani, R. (2013). The sparse group lasso. Journal of Computational and Graphical Statistics , 22(2):231-245.                                         |
| Upton, G. and Fingleton, B. (1985). Spatial data analysis by example. Volume 1: Point pattern and quantitative data. John Wiley and Sons Ltd.                                                      |
| van de Geer, S., Buhlmann, P., Ritov, Y., and Dezeure, R. (2014). On asymptotically optimal confidence regions and tests for high-dimensional models. The Annals of Statistics , 42(3):1166- 1202. |
| Yuan, M. and Lin, Y. (2006). Model selection and estimation in regression with grouped variables. Journal of the Royal Statistical Society , B(68):49-67.                                          |
| Zhang, C.-H. and Zhang, S. S. (2011). Confidence intervals for low-dimensional parameters in high-dimensional linear models. Journal of the Royal Statistical Society , 76(1):217-242.             |
| Zhao, P. and Yu, B. (2006). On model selection consistency of lasso. Journal of Machine Learning Research , 7:2541-2563.                                                                           |
| Zhu, Y. (2018). Sparse linear models and l1regularized 2sls with high-dimensional endogenous regressors and instruments. Journal of Econometrics , 202:196-213.                                    |
Sida Peng
Heterogeneous Endogenous Effects in Networks
This paper proposes a new method to identify leaders and followers in a
network. Prior works use spatial autoregression models (SARs) which implicitly
assume that each individual in the network has the same peer effects on others.
Mechanically, they conclude the key player in the network to be the one with
the highest centrality. However, when some individuals are more influential
than others, centrality may fail to be a good measure. I develop a model that
allows for individual-specific endogenous effects and propose a two-stage LASSO
procedure to identify influential individuals in a network. Under an assumption
of sparsity: only a subset of individuals (which can increase with sample size
n) is influential, I show that my 2SLSS estimator for individual-specific
endogenous effects is consistent and achieves asymptotic normality. I also
develop robust inference including uniformly valid confidence intervals. These
results also carry through to scenarios where the influential individuals are
not sparse. I extend the analysis to allow for multiple types of connections
(multiple networks), and I show how to use the sparse group LASSO to detect
which of the multiple connection types is more influential. Simulation evidence
shows that my estimator has good finite sample performance. I further apply my
method to the data in Banerjee et al. (2013) and my proposed procedure is able
to identify leaders and effective networks.