bayespecon.dgp.simulate_panel_sar_fe

bayespecon.dgp.simulate_panel_sar_fe(N, T, rho=0.4, beta=None, sigma=1.0, sigma_alpha=0.5, err_hetero=False, rng=None, seed=None, W=None, gdf=None, n=None, contiguity='queen', create_gdf=False, geometry_type='polygon', wide=False)[source]

Simulate SAR panel data in time-first stacking order.

DGP

For each period t: y_t = (I-rho W)^(-1) (X_t beta + alpha + eps_t).

param N:

Number of units and time periods.

type N:

int

param T:

Number of units and time periods.

type T:

int

param rho:

Spatial lag coefficient.

type rho:

float, default=0.4

param beta:

Coefficients including intercept.

type beta:

np.ndarray, optional

param sigma:

Idiosyncratic noise scale.

type sigma:

float, default=1.0

param sigma_alpha:

Unit effect scale.

type sigma_alpha:

float, default=0.5

param err_hetero:

If True, generate heteroskedastic innovations with observation-specific standard deviations \(\sigma_i = \sigma \sqrt{1 + \|x_{it}\|^2}\) per period.

type err_hetero:

bool, default=False

param rng:

Random state controls.

param seed:

Random state controls.

param W:

Spatial structure input and GeoDataFrame neighbor rule.

param gdf:

Spatial structure input and GeoDataFrame neighbor rule.

param contiguity:

Spatial structure input and GeoDataFrame neighbor rule.

returns:

Keys: y, X, unit, time, W_dense, W_graph, params_true.

rtype:

dict