bayespecon.models.flow_panel.PoissonFlowPanel¶
- class bayespecon.models.flow_panel.PoissonFlowPanel(y, G, X, T, **kwargs)[source]¶
Non-spatial Bayesian OD-flow Poisson gravity model for balanced panel data.
Panel analogue of
PoissonFlowand count analogue ofOLSFlowPanel. Models stacked panel flow counts with a log-linear gravity mean and no spatial-lag terms,\[y_{ij,t} \sim \operatorname{Poisson}(\lambda_{ij,t}), \qquad \log \boldsymbol{\lambda}_{t} = X_{t}\beta,\]on a balanced panel of \(T\) periods. Provided as the canonical aspatial count baseline for panel flow data.
- Parameters:¶
- y : array-like of int¶
Stacked panel flow counts in shape
(T, n, n),(T, n^2), or(n^2 * T,). Non-integer values raise an error if they do not round exactly.- G : libpysal.graph.Graph¶
Row-standardised graph on
nunits. Required for API symmetry but not used in estimation.- X¶
Stacked panel design matrix in time-first order.
- T : int¶
Number of panel periods.
- col_names : list of str, optional
Feature names for
X. Inferred from a DataFrame if omitted.- k : int, optional
Number of destination/origin covariate pairs used by flow effects; inferred from columns prefixed
dest_if omitted.- priors : dict, optional
Override default priors. Supported keys:
beta_mu(float, default 0.0): Normal prior mean for \(\beta\).beta_sigma(float, default 10.0): Normal prior std for \(\beta\).
Spatial keys (
rho_*),sigma_sigma, and therobustflag are ignored (Poisson has no scale parameter to robustify).- symmetric_xo_xd : bool, optional
Whether to constrain origin and destination covariate effects to be equal. Forwarded to
FlowPanelModel.
Notes
Currently supports pooled panels only (
model=0). Within-transformed fixed-effects panels are not valid for Poisson counts (they break the non-negative integer support), matching the restriction enforced byPoissonSARFlowPanel.Methods
__init__(y, G, X, T, **kwargs)fit([draws, tune, chains, target_accept, ...])Draw samples from the posterior.
fit_approx([draws, n, method, random_seed, ...])Fit a variational approximation and return posterior draws.
posterior_predictive([n_draws, random_seed, ...])Draw posterior-predictive flow counts for the panel Poisson gravity model.
Run Bayesian LM specification tests for flow panel models.
spatial_diagnostics_decision([alpha, ...])Return a model-selection decision from Bayesian LM test results.
spatial_effects([draws, ...])Summarise posterior origin/destination/intra/network/total effects.
summary([var_names])Return posterior summary table via ArviZ.
Attributes
Return the most recent PyMC variational approximation, if any.
Return posterior draws from the most recent fit.
Return the most recently built PyMC model.
-
fit(draws=
2000, tune=1000, chains=4, target_accept=0.9, random_seed=None, store_lambda=False, idata_kwargs=None, **sample_kwargs)[source]¶ Draw samples from the posterior.
-
fit_approx(draws=
2000, n=10000, method='advi', random_seed=None, store_lambda=False, compute_log_likelihood=True, **fit_kwargs)[source]¶ Fit a variational approximation and return posterior draws.
- property inference_data : arviz.data.inference_data.InferenceData | None[source]¶
Return posterior draws from the most recent fit.
-
posterior_predictive(n_draws=
None, random_seed=None, parallel=-1)[source]¶ Draw posterior-predictive flow counts for the panel Poisson gravity model.
- property pymc_model : pymc.model.core.Model | None[source]¶
Return the most recently built PyMC model.
- spatial_diagnostics()[source]¶
Run Bayesian LM specification tests for flow panel models.
See
bayespecon.models.base.SpatialModel.spatial_diagnostics()for the column schema.- Raises:¶
RuntimeError – If the model has not been fit yet.
-
spatial_diagnostics_decision(alpha=
0.05, format='graphviz', theme='default')[source]¶ Return a model-selection decision from Bayesian LM test results.
Walks the panel-flow decision tree using Bayesian p-values from
spatial_diagnostics()and recommends eitherOLSFlowPanel(no spatial dependence detected) orSARFlowPanel(at least one direction is significant).- Parameters:¶
- alpha : float, default 0.05¶
Significance level for the Bayesian p-values.
- format : {"graphviz", "ascii", "model"}, default "graphviz"¶
Output format.
"model"returns the recommended model name string."ascii"returns an indented box-drawing tree."graphviz"returns agraphviz.Digraph(with ASCII fallback if graphviz is not installed).
- Return type:¶
str or graphviz.Digraph
-
spatial_effects(draws=
None, return_posterior_samples=False, ci=0.95, mode='auto', parallel=-1)[source]¶ Summarise posterior origin/destination/intra/network/total effects.
See
bayespecon.models.flow.FlowModel.spatial_effects()for themodesemantics (auto / combined / separate destination-origin sides per Thomas-Agnan & LeSage 2014, §83.5.2) and theparallelkwarg.