pybhatlib 0.2.0__tar.gz → 0.3.1__tar.gz

{pybhatlib-0.2.0 → pybhatlib-0.3.1}/.gitignore

@@ -30,6 +30,14 @@ htmlcov/
 *.pdf
 results/
 
+# Original GAUSS source
+chandra_code/
+
+# Project-local docs (not for repo)
+CLAUDE.md
+GAUSS_INTEGRATION_PLAN.md
+IMPLEMENTATION_PLAN.md
+
 # OS
 .DS_Store
 Thumbs.db

{pybhatlib-0.2.0 → pybhatlib-0.3.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pybhatlib
-Version: 0.2.0
+Version: 0.3.1
 Summary: Python implementation of BHATLIB: matrix-based inference for advanced econometric models
 Project-URL: Homepage, https://github.com/UMN-Choi-Lab/pybhatlib
 Project-URL: Repository, https://github.com/UMN-Choi-Lab/pybhatlib
@@ -96,6 +96,22 @@ results.summary()
 - **Dual backend**: NumPy (default) + optional PyTorch with GPU support
 - **Post-estimation**: Average Treatment Effects (ATE), forecasting
 
+## Verification
+
+pybhatlib reproduces Table 1 from the BHATLIB paper (Bhat 2018) using synthetic
+TRAVELMODE data with 3 travel modes (DA, SR, TR) and 210 observations:
+
+| Model | Specification | Target LL | Achieved LL | Status |
+|-------|--------------|-----------|-------------|--------|
+| (a)(i) | IID errors | -670.956 | -670.956 | Exact match |
+| (a)(ii) | Flexible covariance | -661.111 | -661.111 | Exact match |
+| (b) | + AGE45 demographics | -659.285 | -659.284 | Exact match |
+| (c) | + Random coeff. OVTT | -635.871 | -635.871 | Exact match |
+| (d) | 2-segment mixture | -634.975 | -632.912 | Close (synthetic data) |
+
+Models (a)–(c) match the published results exactly. Model (d) achieves a slightly
+better log-likelihood due to multiple local optima in mixture models with synthetic data.
+
 ## References
 
 1. Bhat, C. R. (2018). New Matrix-Based Methods for the Analytic Evaluation of the

{pybhatlib-0.2.0 → pybhatlib-0.3.1}/README.md

@@ -58,6 +58,22 @@ results.summary()
 - **Dual backend**: NumPy (default) + optional PyTorch with GPU support
 - **Post-estimation**: Average Treatment Effects (ATE), forecasting
 
+## Verification
+
+pybhatlib reproduces Table 1 from the BHATLIB paper (Bhat 2018) using synthetic
+TRAVELMODE data with 3 travel modes (DA, SR, TR) and 210 observations:
+
+| Model | Specification | Target LL | Achieved LL | Status |
+|-------|--------------|-----------|-------------|--------|
+| (a)(i) | IID errors | -670.956 | -670.956 | Exact match |
+| (a)(ii) | Flexible covariance | -661.111 | -661.111 | Exact match |
+| (b) | + AGE45 demographics | -659.285 | -659.284 | Exact match |
+| (c) | + Random coeff. OVTT | -635.871 | -635.871 | Exact match |
+| (d) | 2-segment mixture | -634.975 | -632.912 | Close (synthetic data) |
+
+Models (a)–(c) match the published results exactly. Model (d) achieves a slightly
+better log-likelihood due to multiple local optima in mixture models with synthetic data.
+
 ## References
 
 1. Bhat, C. R. (2018). New Matrix-Based Methods for the Analytic Evaluation of the