networkinference 0.0.4

Links

A tutorial can be found at

https://nbviewer.org/github/mpleung/networkinference/blob/main/docs/tutorial/tutorial.ipynb

Online documentation, which contains minimal working examples for all functions and classes, is available at

https://networkinference.readthedocs.io/en/latest/

The latest development version is hosted at

https://github.com/mpleung/networkinference

Installation

Install via command line using pip3, which is included with Python 3.4+:

$ pip3 install networkinference

Overview

The package contains three main classes:

• OLS: OLS estimator.

• TSLS: 2SLS estimator.

• IPW: Horvitz-Thompson estimator (inverse-probability weighting estimator with known propensity scores).

Each class contains five methods for constructing confidence intervals and two methods implementing scalar equality tests. These methods are based on three classes of inference procedures. The first uses a network HAC variance estimator [1] [2]. The second constructs network clusters using spectral clustering and applies a cluster-robust inference method [4]. Both require network data and account for a general form of network dependence in which observations are less correlated when further apart in the network (in the sense of shortest path distance) [1]. The third involves resampling and can be applied to weakly dependent network, clustered, time series, or spatial data (or combination or these) without knowledge of the particular type of dependence [3] [5]. However, it is more computationally intensive and requires a larger amount of data to ensure adequate size control and power.

The methods are also available as independent functions through the core class for use with custom estimators.

The core class and utils subpackage contain various utilities for computing summary statistics, plotting the spectrum of the Laplacian, constructing friends-of-friends instruments for estimating linear-in-means models, and simulating data.

Example Usage

>>> import networkinference as ni               # main package
>>> from networkinference.utils import FakeData # utilities for generating fake data
>>> Y, X, W, A = FakeData.tsls()                # simulate data from linear-in-means model
>>> tsls = ni.TSLS(Y, X, W, A)                  # load data into tsls object
>>> tsls.network_se()                           # print estimates with network-robust SEs
>>> help(tsls)                                  # displays documentation for tsls object

References

Dependencies

• Matplotlib v3.7+

• NetworkX v3.0+

• NumPy v1.24+

• Scikit-learn v1.2+

• SciPy v1.10+

• Seaborn v0.12+

• Tabulate v0.9+