1 Introduction
Social actors are often embedded in webs of relationships that profoundly shape political and economic outcomes (Franzese and Hays, 2008; Ward, Stovel and Sacks, 2011). One challenge in analyzing networks arises in situations where an analyst cannot fully observe the nature of relational ties. In many dyadic interactions — treaties, marriages — outsiders can observe ties only if both agents agree, that is, the payoff for forming a tie exceeds its cost for both members of a dyad (Jackson and Wolinsky, 1996). The observed network is therefore composed of symmetric (“undirected”) ties even though the social process at work contains important relational asymmetries. The pursuit of a tie by one party may not be reciprocated to the same extent by another.
As an illustration, suppose , , and are three warring factions deciding whether to sign bilateral peace agreements. We observe a network in which dyads and have signed agreements, but the dyad continues fighting. This observed network of ‘peaceful’ ties could be generated from any of three unobserved sets of relations: it may be that failed to reciprocate ’s pursuit of peace, or vice versa, or neither nor pursued peace. To identify conditions that drive factions to sign peace agreements we must account for these unobserved asymmetries. Here, the observed symmetric graph is an incomplete representation of the underlying, asymmetric network, which is frequently the object of scientific interest. We refer to this situation as “partial observability.”
We present the partial observability generalized bilinear mixed effects model (PGBME) to address this challenge. The model is a synthesis of the generalized bilinear mixed effects (GBME) model (Hoff, 2005) and the bivariate or “partial observability” probit model (Poirier, 1980; Przeworski and Vreeland, 2002). The model can probabilistically reconstruct the directed network from which the observed, undirected graph emerged. The model enables the study of network ties in a regression framework by accounting for interdependencies as well as unobserved asymmetries in network relations. The stochastic actororiented model (SAOM) for networks (Snijders and Pickup, 2017) also allows for partial observability. However, SAOM was designed to assess how specific network features (e.g.,
star triangles) give rise to an observed network. The latent network approach is not used to study the role of specific network statistics. Rather, latent network models aim to account for broad patterns of network interdependence using a variance decomposition regression framework.
^{1}^{1}1See Minhas, Hoff and Ward (2016)for detailed discussion. Other approaches based on generalized spatiotemporal dependence can also recover directed predicted probabilities
Franzese, Hays and Kachi (2012).We illustrate the PGBME model by applying it to the bilateral investment treaties (BIT) network for each year in 19902012. The model substantially improves predictive accuracy relative to both conventional logit and standard GBME. As important, PGBME extracts new information about the factors that drive treaty preferences, identifies important structural changes in the network, and highlights possible “hidden” agreements that are easily overlooked when latent network asymmetries are ignored.
2 The Model
Building on the random utility framework (McFadden, 1980), we model an actor as having net utility from forming a tie with another : with representing the systematic component (that depends on observables) and representing the stochastic error. To account for interdependencies in actors’ utilities from having ties, we use the “latent space” approach (Hoff, 2005) and model these utilities as follows:
(1) 
The correlation, , captures the “reciprocity” between the utilities that actors derive from tie formation. Parameters and
are sender and receiverspecific random effects, respectively, and they capture secondorder network dependencies. The vectors
and represent the location of actor in the latent space of ‘senders’ and ‘receivers,’ respectively. These random effects capture higherorder dependencies in network ties: derives a large utility from forming a tie with , if ’s location in the latent space of ‘senders’ is close to ’s location in the latent space of ‘receivers’ (so that the crossproduct is large).^{2}^{2}2As in previous literature, these random effects are modeled as , , and , where , , and are unknown parameters. The choice of , dimensionality of the latent space, is discussed supplementary materials. We express the systematic components of actors’ utilities as linear functions of predictors:(2)  
(3) 
A researcher cannot directly observe net utility (the ’s). We only observe agents’ behaviors, in this case undirected bilateral ties. A directional tie is formed if and only if ’s net gain from doing so is strictly positive, . Accordingly, the bilateral tie is formed if and only if both actors derive a net positive payoff from having a tie so that and . A researcher observes an undirected (bilateral) tie arising from the following data generating process:
(4) 
Under a standard identifying restriction , the model is a partially observable probit regression (Poirier, 1980), augmented with randomeffects to capture unobserved heterogeneity and interdyadic dependencies. Vectors and represent the senderspecific and receiverspecific covariates, respectively. A model for the directional link (eq. 11), uses variables as senderspecific predictors, but these same predictors become receiverspecific in the model for the directional link (eq. 12). Vector contains dyadspecific variables. These dyadspecific variables might be symmetric, (e.g., distance between countries), or not (e.g., exportimport).
A partially observed probit model requires at least one of the following identifying restrictions: (1) regression equations 11 and 12 must have the same parameters and/or (2) one equation contains a predictor not included in another equation (Poirier, 1980). If the dyadic predictors are asymmetric, , then condition (2) is satisfied. Furthermore, regression equations 11 and 12 have the same parameters, and so condition (1) holds as well; thus, the above model is parametrically identified. However, we impose an additional restriction that . While this restriction is not required for parametric identification, Rajbhandari (2014)
showed that finite sample estimates of
are sensitive to the starting values and generally cannot be treated as reliable.^{3}^{3}3The estimation algorithm provided with this paper allows to be estimated, but caution should be used when utilizing this option.We estimate the model in a Bayesian framework using Markov chain Monte Carlo. In the supplemental materials give a more detailed exposition of the model, prior assumptions, the sampling algorithm. We also provide results from a simulation study demonstrating that the model successfully recovers known parameter values.
3 Application: Bilateral Investment Treaties
We apply the PGBME model to the network of bilateral investment treaties (BITs) from 1990 to 2012 using the standard United Nations BIT database. There is a vibrant debate on whether BITs boost FDI (Jandhyala, Henisz and Mansfield, 2011; Simmons, 2014; Minhas, 2016), but a proper resolution of this debate requires a convincing empirical model of treaty formation (Rosendorff and Shin, 2012). Partial observability is one key challenge in building such model: the observed network of signed bilateral treaties is symmetric, while the underlying preferences for these treaties are asymmetric.
We fit the PGBME model separately for each year of data using a suite of covariates that closely follows the existing empirical literature (see supplementary materials). Our model improves on the previous literature by accounting for both network interdependence and partial observability.
3.1 Predictive Performance
In Table 1 we compare the insample and outofsample predictive performance of the PGBME to that of pooled probit, which assumes dyadic independence and ignores partial observability, and GBME, which models dyadic interdependencies, but not partial observability. The predictive accuracy of the GBME model in this case is similar to the pooled probit. Adding the partial observability component to the GBME model, however, produces an additional substantial improvement in the predictive accuracy as shown by all metrics for the PGBME model.
Insample  Outofsample  

ROC  PR  ROC  PR  
Pooled probit  0.75  0.48  0.73  0.44 
GBME  0.76  0.47  0.77  0.48 
PGBME  0.90  0.71  0.91  0.78 
Predictive performance in BIT data: area under the receiver operating characteristic curve (ROC) and area under the precision recall curve (PR).
3.2 Regression Parameters
Existing models, including GBME, estimate a single coefficient for each predictor. This assumes away the possibility that the same factor differentially “affects” ’s demand for a treaty with and ’s attractiveness to . The PGBME recovers directed sender and receivereffects for nodelevel covariates. For instance, our estimates suggest that, countries faster growth in GDP per capita were no more inclined to sign BITs with others (sender effect). But highgrowth countries were more attractive BIT partners to others (receiver effect). Supplementary materials describes regression parameter estimates in detail.
3.3 The Structure of Latent Treaty Preferences
The PGBME model allows us to extract “latent preferences” for treaty formation – the estimated probability that country demands a treaty from , and vice versa. Figure 1
displays the mean posterior predicted probabilities relevant to China in 1995 and 2010. The horizontal axis represents a country’s attractiveness to China as a BIT partner and the vertical axis is China’s attractiveness to that country. Countries above the diagonal line find China a more attractive BIT partner than China finds them, and vice versa. Color identifies observed BITs.
The plots reveal how China’s position in the BIT network changed over time. In 1995, China was moving aggressively to demand BITs around the world, forming ties that the model views as relatively unlikely. By 2010, China had many more BITs in place and is more likely to be a treaty target by the remaining countries, an indication of China’s expanded role in the global economy.
1995  2010 

Figure 2 illustrates a different use of the PGBME model. It shows the predicted probabilities that the USA is demanded (left) and demands (right) a BIT in 2010. Observe that the model predicts Peru, Mexico, and Chile—the top 3 nonBITs—demand a BIT with the USA at probabilities close to one, and are also likely BIT targets of the USA with probabilities exceeding . Closer inspection of UN treaty data reveals that, by 2010, all three had signed other agreements with the US that contain provisions functionally equivalent to BITs; these agreements do not appear in the BIT dataset commonly used in the literature.^{4}^{4}4The treaties are 2006 PeruUSA Free Trade Agreement (FTA), NAFTA in 1992, and the 2003 ChileUSA FTA. These agreements included investment provisions that mirror the terms of a BIT almost exactly (see http://investmentpolicyhub.unctad.org/Download/TreatyFile/5454). The PGBME model nevertheless highlights these “hidden” agreements as dyads likely to have a BIT. This suggests that researchers studying BITs need to carefully examine the dataset they employ and perhaps expand the set of treaties considered relevant.
4 Conclusion
Partial observability occurs whenever the observed graph is undirected yet the underlying process implies directed relationships. We introduced a model that can reconstruct the latent directed network ties, and illustrated its advantages on an example of bilateral trade agreements. The future work in this area could focus on several extensions. First, we accounted for network dependencies using bilinear mixed effects framework (GBME), which could be generalized using the recently developed additive and multiplicative effects network model (Hoff, 2015; Minhas, Hoff and Ward, 2016). Second, estimating this type of network model in a fully dynamic setting (as opposed to slicing data by time, as we did here) remains a challenge, especially when the set of nodes changes in time.
Supplemental Materials
Appendix A Details on the PGBME
As detailed in the paper, the partial observability generalized bilinear mixed effects (PGBME) framework treats the observed symmetric outcome, , as resulting from a joint decision taken by a pair of actors. We formalize the joint decision making process using a bivariate probit model with a standard normal link function:
(5)  
(6)  
(7)  
(8)  
(9) 
and represent sender and receiver random effects that account for first order dependence patterns that often arise in relational data, while captures the likelihood of a pair of actors interacting with one another based on third order dependence patterns such as transitivity, balance, and clustering. For identification purposes, we fix and . The former is a standard restriction in probit frameworks with a binary outcome. We undertake the latter restriction because Rajbhandari (2014) shows that in this framework it is difficult to recover reliable estimates for as the parameter is highly sensitive to the initial value.
The sender and receiver random effects ( and
) are drawn from a multivariate normal distribution centered at zero with a covariance matrix,
, parameterized as follows:(10) 
The nodal effects are modeled in this way to account for the fact that in many relational datasets we often find that actors who send a lot of ties are also more likely to receive a lot of ties. Heterogeneity in the the sender and receiver effects is captured by and , respectively, and describes the covariance between these two effects.
represents the systematic component of actors’ utilities and is expressed as a linear function of sender , receiver , and dyadic covariates:
(11)  
(12) 
This formulation allows us to incorporate exogenous actor and dyad level characteristics into how actors make decisions within the partial probit framework. Following Hoff (2005), to enable a more efficient estimation, we reparameterize the model to implement hierarchical centering of the random effects (Gelfand, Sahu and Carlin, 1995):
(13)  
(14)  
(15) 
a.1 Parameters and Priors
To estimate the parameters discussed in the previous section, we utilize conjugate priors and a Monte Carlo Markov Chain (MCMC) algorithm. Prior distributions for the parameters are specified as follows:
^{5}^{5}5For details on the full conditional distributions of each of the parameters see Hoff (2005).
, , and are each drawn from multivariate normals with mean zero and a covariance matrix in which the covariances are set to zero and variances to 10


, and are each drawn from an i.i.d. inverse gamma(1,1).
Starting values for each of the parameters are determined using maximum likelihood estimation.
a.2 The MCMC algorithm
To estimate this model a Gibbs sampler is used. This sampler follows the procedure laid out in Hoff (2005, 2009) with the exception of the first step in which we extend the GBME by accounting for the possibility that seemingly symmetric events are the result of a joint decision between a pair of actors. This first step involves sampling from a truncated normal distribution, we show the full conditional distribution below.

Modeling partially observable outcome. Conditional on there being an observed link between and , and conditional on other parameters, we draw the latent variables and from the bivariate normal distribution such that both latent variables are positive:
Conditional (there is no observed link between and ), we sample the latent variables from the bivariate normal distribution where at least one of the latent variables, or , is constrained to be negative:

Additive effects

Sample

Sample (linear regression)

Sample from full conditional distribution


Multiplicative effects^{6}^{6}6See Hoff (2009) for further details on how multiplicative effects are estimated in a directed context within the GBME framework.

For :

Sample (linear regression)

Sample (linear regression)


a.3 Simulation Exercise
To test the capabilities of the PGBME framework in representing the data generating process for a partially observable outcome we conduct a simulation exercise. In each simulation, we randomly construct a directed network from a pair of dyadic covariates, nodal covariates, and the random effects structure detailed in the previous section. The regression parameters for the dyadic covariates are set at 1 and 1/2, and the parameters for the nodal covariates are set at 0 and 1/2. At this stage, the network simulated from this data generating process is directed. We modify the simulated network so that a link between a dyad only appears in the network if both the dyad and the dyad both have a link in the simulated network, thus making the network appear undirected.
Next, we examine whether the PGBME model can recover the data generating process underlying the partially observed simulated network. We run the PGBME model in every simulation for 20,000 iterations with a 10,000 burnin period. We repeat this simulation process 100 times.
With the simulation results our first step is to examine whether the PGBME accurately recovers the regression parameter estimates for the dyadic and nodal covariates. To test whether this is the case we calculate the mean regression parameter estimate from the MCMC results for each simulation, and we summarize these results in Figure 1. For each parameter we indicate its true value by a colored horizontal line and summarize the distribution of the mean regression values estimated from the PGBME using a boxplot. Given that for each of the parameters the true value almost exactly crosses the median value indicated in the box plot, this simulation shows that the PGBME is quite effective in estimating the true parameter values underlying a partially observed outcome.
We also examine the proportion of times that the true value falls within the 95% credible interval of the estimated regression parameter. In over ninety percent of the simulations, the true value falls within the 95% credible interval of each of the estimated regression parameters from the PGBME. Specifically, for
the coverage rate is 0.85, for 0.93, for 0.93, and for 0.97.Appendix B Estimation and Application
b.1 Data
Table B1 provides a description for each of the variables used in the analysis.
Variable  Level  Definition  Source 

BIT  dyadic  if & Signed a BIT by year  UNCTAD^{7}^{7}7http://investmentpolicyhub.unctad.org/IIA 
UDS (median)  dyadic  UDSUDS  Pemstein et al. (2010) 
Law & Order  dyadic  LOLO  ICRG^{8}^{8}8http://epub.prsgroup.com/products/internationalcountryriskguideicrg 
Log(GDP per capita)  dyadic  Log(GDPcap)Log(GDPcap)  WDI 
OECD  dyadic  if & Both OECD members by year  OECD 
Distance  dyadic  Minimum distance between &  Gleditsch & Ward (2001) ^{9}^{9}9Also see Weidmann and Gleditsch (2015). 
FDI/GDP  node  Net FDI inflow as % GDP in year  WDI 
ICSID Disputes  node  Cumulative number of disputes by year  ICSID^{10}^{10}10https://icsid.worldbank.org/en/Pages/cases/AdvancedSearch.aspx 
GDP per capita growth  node  Level of GDP per capita growth by year  WDI 
PTAs  node  Cumulative number of PTAs signed by year  DESTA^{11}^{11}11https://www.designoftradeagreements.org/ 
A shortcoming of the existing GBME framework is its inability to account for applications where there is missingness in the set of exogenous covariates used in the model. For our application, a number of the nodal covariates had varying levels of missingness. Additionally, most of the dyadic covariates that we construct from nodal variables, such as the unified democracy scores, also have varying levels of missingness. The table below shows how much missingness we had for the variables included in our analysis:
Variable  Proportion of Cases Missing 

Law & Order  16.6% 
FDI/GDP  2.8% 
GDP per capita  1.4% 
GDP per capita growth  1.4% 
Unified Democracy Scores (UDS)  0.7% 
ICSID Disputes  0% 
PTAs  0% 
OECD  0% 
In general, the level of missingness is not high. The only exception here is with the Law & Order variable from the ICRG dataset, for this variable we had approximately 17% of countryyear observations missing from 1990 to 2012. The only true dyadic variable we include in our analysis is a calculation of the minimum distance between countries, and this variable has no missingness. Additionally, our dependent variable measuring whether or not two countries had signed a BIT by year also has no missingness.
A number of works have noted the issues that can arise when simply using listwise deletion,^{12}^{12}12See, for example, King et al. (2001).
thus before running the PGBME sampler we impute missingness among the covariates used in our model with a Bayesian, semiparametric copula imputation scheme.
^{13}^{13}13See Hoff (2007); Hollenbach et al. (2016) for details on this imputation scheme and how it differs from other approaches frequently utilized in political science. We generate 1,000 imputed datasets from this imputation scheme and save 10 for use in the PGBME MCMC sampler.To account for missingness within the PGBME, at the beginning of every iteration of the MCMC for model, we draw a randomly sampled imputed dataset from the posterior of the Copula, calculate the parameters associated with the PGBME using the imputed dataset, and repeat this process for every iteration of the sampler for the model. This approach directly incorporates imputation uncertainty into our posterior distributions of the PGBME parameters without having to run and combine separate models.
b.2 Estimation details
In our application we estimate the PGBME separately for each year from 1990 to 2012 using the prior distributions and MCMC algorithm described above. For each year, we ran the PGBME MCMC sampler for 20,000 iterations, discarding the first 10,000 iterations as burnin. We thinned the chain by saving only every 10 value.
The following trace plot describes MCMC convergence for all parameters in the 2012 PGBME model.
b.3 Regression Parameters
Figure 3 displays the posterior mean and the 90% (thicker) and 95% (thinner) credible intervals (CIs). The first column contains the dyadic covariates; the second, senderlevel covariates; and, the third, receiver covariates. In each panel, we show the parameter estimates for that variable from 1990 to 2012. The dotted horizontal line is and the thicker grey line is the posterior mean, pooling the posterior draws across all years.
The PGBME recovers directed sender and receivereffects for nodelevel covariates from an observed undirected network, something that other approaches, including the GBME, are unable to do. Our estimation in this application indicates substantial instability in these estimates over time, both relative to a baseline of and relative to the pooled posterior mean. This instability is consistent with substantive arguments that the incentives to sign BITs have changed over time (Jandhyala, Henisz and Mansfield, 2011).
Dyadic covariates tend to be more stable across years in this application. But they, too, show that the BIT formation process has changed over time. Economic and political “distance”, for example, has become less important as the network evolved and more lowerincome countries have signed BITs with each other. Geographic distance, on the other hand, continues to be strongly related to the formation of BITs.
PGBME covariate estimates shed light on the evolving processes producing the observed BIT network in the 19902012 period. The changing values of covariate parameters over time also indicates that common practice of pooling dyads and assuming the existence of temporally stable parameters may be dangerous.
b.4 Choosing dimension of the multiplicative effects,
One of the parameters that users are able to set within the PGBME to account for third order dependence patterns is – see the MCMC algorithm section above for more details on this parameter and its relation to the model. In the results reported in the paper, we set = 2. To understand whether or not a higher value of is necessary users of this approach can compare the insample fit of the model with varying values for . In our application exercise, we varied from 1 to 3 to settle on an appropriate value of that can represent the data generating process of the network. Results are shown in Table B3 below.
AUC (ROC)  AUC (PR)  

=1  0.87  0.63 
=2  0.90  0.71 
=3  0.92  0.71 
As you can see after =2, the subsequent insample performance improvement notably declines. There is a slight increase in performance from =2 to =3, however, every time one increases we are also adding more parameters to the PGBME model. Adding this many more parameters can easily lead one to overfit the data in an outofsample context.
References
 (1)
 Franzese and Hays (2008) Franzese, Robert J. and Jude C. Hays. 2008. “Interdependence in Comparative Politics: Substance, Theory, Empirics, Substance.” Comparative Political Studies 41(4/5):742–780.
 Franzese, Hays and Kachi (2012) Franzese, Robert J., Jude C. Hays and Aya Kachi. 2012. “Modeling History Dependence in NetworkBehavior Coevolution.” Political Analysis 20(2):175–190.

Gelfand, Sahu and Carlin (1995)
Gelfand, Alan E., Sujit K. Sahu and Bradley P. Carlin. 1995.
“Efficient parametrisations for normal linear mixed models.”
Biometrika 82(3):479–488.  Gleditsch and Ward (2001) Gleditsch, Kristian S. and Michael D. Ward. 2001. “Measuring Space: A Minimum Distance Database and Applications to International Studies.” Journal of Peace Research 38(6):749–768.
 Hoff (2009) Hoff, Peter. 2009. “Multiplicative Latent Factor Models for Description and Prediction of Social Networks.” Computational and Mathematical Organization Theory 15(4):261–272.
 Hoff (2005) Hoff, Peter D. 2005. “Bilinear MixedEffects Models for Dyadic Data.” Journal of the American Statistical Association 100(4690):286–295.
 Hoff (2007) Hoff, Peter D. 2007. “sbgcop: Semiparametric Bayesian Gaussian copula estimation and imputation.” R Package.
 Hoff (2015) Hoff, Peter D. 2015. “Dyadic Data Analysis with amen.” arxiv arxiv:1506.08237:1–48.
 Hollenbach et al. (2016) Hollenbach, Florian M., Iavor Bojinov, Shahryar Minhas, Nils W. Metternich, Michael D. Ward and Alexander Volfovsky. 2016. “Principled Imputation Made Simple: Multiple Imputation Using Gaussian Copulas.” Working Paper .
 Jackson and Wolinsky (1996) Jackson, Matthew O. and Asher Wolinsky. 1996. “A strategic model of social and economic networks.” Journal of Economic Theory 71(1):44–74.
 Jandhyala, Henisz and Mansfield (2011) Jandhyala, Srividya, Witold J. Henisz and Edward D. Mansfield. 2011. “Three waves of BITs: The global diffusion of foreign investment policy.” Journal of Conflict Resolution 55(6):1047–1073.
 King et al. (2001) King, Gary, James Honaker, Anne Joseph and Kenneth Scheve. 2001. “Analyzing Incomplete Political Science Data: An Alternative Algorithm for Multiple Imputation.” American Political Science Review 95(1):49–69.
 McFadden (1980) McFadden, Daniel. 1980. “Econometric models for probabilistic choice among products.” Journal of Business pp. S13–S29.
 Minhas (2016) Minhas, Shahryar. 2016. Learning from Incredible Commitments: Evolution and Impact of Bilateral Investment Treaties Ph. D. Thesis Duke University Durham, North Carolina: .

Minhas, Hoff and Ward (2016)
Minhas, Shahryar, Peter D. Hoff and Michael D. Ward. 2016.
“Let’s Say Amen for Latent Factor Models.” Working paper.
https://arxiv.org/abs/1611.00460  Pemstein, Meserve and Melton (2010) Pemstein, Daniel, Steven A. Meserve and James Melton. 2010. “Democratic Compromise: A Latent Variable Analysis of Ten Measures of Regime Type.” Political Analysis 18(4):426–449.
 Poirier (1980) Poirier, Dale J. 1980. “Partial Observability in Bivariate Probit Models.” Journal of Econometrics 12(2):209–0217.
 Przeworski and Vreeland (2002) Przeworski, Adam and James Raymond Vreeland. 2002. “A Statistical Model of Bilateral Cooperation.” Political Analysis 10(2):101–112.
 Rajbhandari (2014) Rajbhandari, Ashish. 2014. Identification and MCMC Estimation of Bivariate Probit Models with Partial Observability. In Bayesian Inference in the Social Sciences, ed. Ivan Jeliazkov and XinShe Yang. Wiley.
 Rosendorff and Shin (2012) Rosendorff, B. Peter and Kongjoo Shin. 2012. “Importing transparency: The political economy of BITs and FDI Flows.” Manuscript, New York University Political Science Department. New York, NY: NYU. .
 Simmons (2014) Simmons, Beth A. 2014. “Bargaining over BITs, Arbitrating Awards: The Regime for Protection and Promotion of International Investment.” World Politics 66(1):12–46.
 Snijders and Pickup (2017) Snijders, Tom A. B. and Mark Pickup. 2017. Stochastic Actor Oriented Models for Network Dynamics. In Oxford Handbook of Politicial Networks, ed. Jennifer Nicoll Victor, Alexander H. Montgomery and Mark Lubell. Oxford University Press.
 Ward, Stovel and Sacks (2011) Ward, Michael D, Katherine Stovel and Audrey Sacks. 2011. “Network analysis and political science.” Annual Review of Political Science 14:245–264.

Weidmann and Gleditsch (2015)
Weidmann, Nils B. and Kristian Skrede Gleditsch. 2015.
cshapes: CShapes Dataset and Utilities.
package version 0.51.
https://CRAN.Rproject.org/package=cshapes