bartz 0.4.0__py3-none-any.whl → 0.4.1__py3-none-any.whl

bartz/__init__.py

@@ -10,10 +10,10 @@
 # to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
 # copies of the Software, and to permit persons to whom the Software is
 # furnished to do so, subject to the following conditions:
-# 
+#
 # The above copyright notice and this permission notice shall be included in all
 # copies or substantial portions of the Software.
-# 
+#
 # THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
 # IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
 # FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
@@ -23,7 +23,7 @@
 # SOFTWARE.
 
 """
-A jax implementation of BART
+Super-fast BART (Bayesian Additive Regression Trees) in Python
 
 See the manual at https://gattocrucco.github.io/bartz/docs
 """

bartz/_version.py

@@ -1 +1 @@
-__version__ = '0.4.0'
+__version__ = '0.4.1'

bartz/jaxext.py

@@ -1,6 +1,6 @@
 # bartz/src/bartz/jaxext.py
 #
-# Copyright (c) 2024, Giacomo Petrillo
+# Copyright (c) 2024-2025, Giacomo Petrillo
 #
 # This file is part of bartz.
 #
@@ -45,27 +45,6 @@ def castto(func, type):
         return func(*args, **kw).astype(type)
     return newfunc
 
-def pure_callback_ufunc(callback, dtype, *args, excluded=None, **kwargs):
-    """ version of `jax.pure_callback` that deals correctly with ufuncs,
-    see `<https://github.com/google/jax/issues/17187>`_ """
-    if excluded is None:
-        excluded = ()
-    shape = jnp.broadcast_shapes(*(
-        a.shape
-        for i, a in enumerate(args)
-        if i not in excluded
-    ))
-    ndim = len(shape)
-    padded_args = [
-        a if i in excluded
-        else jnp.expand_dims(a, tuple(range(ndim - a.ndim)))
-        for i, a in enumerate(args)
-    ]
-    result = jax.ShapeDtypeStruct(shape, dtype)
-    return jax.pure_callback(callback, result, *padded_args, vectorized=True, **kwargs)
-
-    # TODO when jax solves this, check version and piggyback on original if new
-
 class scipy:
 
     class special:
@@ -74,9 +53,11 @@ class scipy:
         def gammainccinv(a, y):
             a = jnp.asarray(a)
             y = jnp.asarray(y)
+            shape = jnp.broadcast_shapes(a.shape, y.shape)
             dtype = float_type(a.dtype, y.dtype)
+            dummy = jax.ShapeDtypeStruct(shape, dtype)
             ufunc = castto(special.gammainccinv, dtype)
-            return pure_callback_ufunc(ufunc, dtype, a, y)
+            return jax.pure_callback(ufunc, dummy, a, y, vmap_method='expand_dims')
 
     class stats:
 

bartz-0.4.1.dist-info/METADATA

@@ -0,0 +1,52 @@
+Metadata-Version: 2.3
+Name: bartz
+Version: 0.4.1
+Summary: Super-fast BART (Bayesian Additive Regression Trees) in Python
+License: MIT
+Author: Giacomo Petrillo
+Author-email: info@giacomopetrillo.com
+Requires-Python: >=3.10
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Programming Language :: Python :: 3
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Programming Language :: Python :: 3.13
+Requires-Dist: jax (>=0.4.35,<1)
+Requires-Dist: jaxlib (>=0.4.35,<1)
+Requires-Dist: numpy (>=1.25.2,<3)
+Requires-Dist: scipy (>=1.11.4,<2)
+Project-URL: Bug Tracker, https://github.com/Gattocrucco/bartz/issues
+Description-Content-Type: text/markdown
+
+[![PyPI](https://img.shields.io/pypi/v/bartz)](https://pypi.org/project/bartz/)
+[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.13931477.svg)](https://doi.org/10.5281/zenodo.13931477)
+
+# BART vectoriZed
+
+An implementation of Bayesian Additive Regression Trees (BART) in JAX.
+
+If you don't know what BART is, but know XGBoost, consider BART as a sort of Bayesian XGBoost. bartz makes BART run as fast as XGBoost.
+
+BART is a nonparametric Bayesian regression technique. Given training predictors $X$ and responses $y$, BART finds a function to predict $y$ given $X$. The result of the inference is a sample of possible functions, representing the uncertainty over the determination of the function.
+
+This Python module provides an implementation of BART that runs on GPU, to process large datasets faster. It is also good on CPU. Most other implementations of BART are for R, and run on CPU only.
+
+On CPU, bartz runs at the speed of dbarts (the fastest implementation I know of) if n > 20,000, but using 1/20 of the memory. On GPU, the speed premium depends on sample size; it is convenient over CPU only for n > 10,000. The maximum speedup is currently 200x, on an Nvidia A100 and with at least 2,000,000 observations.
+
+[This Colab notebook](https://colab.research.google.com/github/Gattocrucco/bartz/blob/main/docs/examples/basic_simdata.ipynb) runs bartz with n = 100,000 observations, p = 1000 predictors, 10,000 trees, for 1000 MCMC iterations, in 5 minutes.
+
+## Links
+
+- [Documentation (latest release)](https://gattocrucco.github.io/bartz/docs)
+- [Documentation (development version)](https://gattocrucco.github.io/bartz/docs-dev)
+- [Repository](https://github.com/Gattocrucco/bartz)
+- [Code coverage](https://gattocrucco.github.io/bartz/coverage)
+- [List of BART packages](https://gattocrucco.github.io/bartz/docs-dev/pkglist.html)
+
+## Citing bartz
+
+Article: Petrillo (2024), "Very fast Bayesian Additive Regression Trees on GPU", [arXiv:2410.23244](https://arxiv.org/abs/2410.23244).
+
+To cite the software directly, including the specific version, use [zenodo](https://doi.org/10.5281/zenodo.13931477).
+

bartz-0.4.1.dist-info/RECORD

@@ -0,0 +1,13 @@
+bartz/BART.py,sha256=CbGzFWtYw5u38Z9-Hy3CbDXpKOOvPFAAkSqu2HZl8no,16862
+bartz/__init__.py,sha256=CcWbDdelQoObG4MW7_jit0Y9fQumU2_Yv3LYtBi1Ag8,1480
+bartz/_version.py,sha256=qBs4HqYsPn6yUHWEuCN4rGUC05gABbl800d8LBd8h9w,22
+bartz/debug.py,sha256=9ZH-JfwZVu5OPhHBEyXQHAU5H9KIu1vxLK7yNv4m4Ew,5314
+bartz/grove.py,sha256=x_6NK_l-hrXfy1PhssYNJkX41-w_WqjDziww0E7YRS8,8500
+bartz/jaxext.py,sha256=UiAHw6pQMz8rjtn_JChFAIdPgUF393S5m-Dduyndnbo,11409
+bartz/mcmcloop.py,sha256=lKDszvniNXka99X3e9RCrTgvEAZHA7ZbVXEgxUYvKMY,7634
+bartz/mcmcstep.py,sha256=diI9vHXHMvu_Lk_bSJ-a038OnEbXDpNEikVPhRcxEys,54987
+bartz/prepcovars.py,sha256=mMgfL-LGJ_8QpOL6iy7yfkL8A7FrT7Zfn5M3voyNwSQ,5818
+bartz-0.4.1.dist-info/LICENSE,sha256=heuIJZQK9IexJYC-fYHoLUrgj8HG8yS3G072EvKh-94,1073
+bartz-0.4.1.dist-info/METADATA,sha256=nOOpuP_b-eAJSR1sgyzJaGK4I1iuEWVFMJ9D3JMkY5M,2890
+bartz-0.4.1.dist-info/WHEEL,sha256=fGIA9gx4Qxk2KDKeNJCbOEwSrmLtjWCwzBz351GyrPQ,88
+bartz-0.4.1.dist-info/RECORD,,

{bartz-0.4.0.dist-info → bartz-0.4.1.dist-info}/WHEEL

@@ -1,4 +1,4 @@
 Wheel-Version: 1.0
-Generator: poetry-core 1.9.0
+Generator: poetry-core 2.1.2
 Root-Is-Purelib: true
 Tag: py3-none-any

bartz-0.4.0.dist-info/METADATA

@@ -1,77 +0,0 @@
-Metadata-Version: 2.1
-Name: bartz
-Version: 0.4.0
-Summary: A JAX implementation of BART
-Home-page: https://github.com/Gattocrucco/bartz
-License: MIT
-Author: Giacomo Petrillo
-Author-email: info@giacomopetrillo.com
-Requires-Python: >=3.10,<4.0
-Classifier: License :: OSI Approved :: MIT License
-Classifier: Programming Language :: Python :: 3
-Classifier: Programming Language :: Python :: 3.10
-Classifier: Programming Language :: Python :: 3.11
-Classifier: Programming Language :: Python :: 3.12
-Requires-Dist: jax (>=0.4.23,<0.5.0)
-Requires-Dist: jaxlib (>=0.4.23,<0.5.0)
-Requires-Dist: numpy (>=1.25.2,<2.0.0)
-Requires-Dist: scipy (>=1.11.4,<2.0.0)
-Project-URL: Bug Tracker, https://github.com/Gattocrucco/bartz/issues
-Project-URL: Repository, https://github.com/Gattocrucco/bartz
-Description-Content-Type: text/markdown
-
-[![PyPI](https://img.shields.io/pypi/v/bartz)](https://pypi.org/project/bartz/)
-
-# BART vectoriZed
-
-A branchless vectorized implementation of Bayesian Additive Regression Trees (BART) in JAX.
-
-BART is a nonparametric Bayesian regression technique. Given predictors $X$ and responses $y$, BART finds a function to predict $y$ given $X$. The result of the inference is a sample of possible functions, representing the uncertainty over the determination of the function.
-
-This Python module provides an implementation of BART that runs on GPU, to process large datasets faster. It is also good on CPU. Most other implementations of BART are for R, and run on CPU only.
-
-On CPU, bartz runs at the speed of dbarts (the fastest implementation I know of), but using half the memory. On GPU, the speed premium depends on sample size; with 50000 datapoints and 5000 trees, on an Nvidia Tesla V100 GPU it's 12 times faster than an Apple M1 CPU, and this factor is linearly proportional to the number of datapoints.
-
-## Links
-
-- [Documentation (latest release)](https://gattocrucco.github.io/bartz/docs)
-- [Documentation (development version)](https://gattocrucco.github.io/bartz/docs-dev)
-- [Repository](https://github.com/Gattocrucco/bartz)
-- [Code coverage](https://gattocrucco.github.io/bartz/coverage)
-
-## Other BART packages
-
-- [stochtree](https://github.com/StochasticTree) C++ library with R and Python bindings taylored to researchers who want to make their own BART variants
-- [bnptools](https://github.com/rsparapa/bnptools) Feature-rich R packages for BART and some variants
-- [dbarts](https://github.com/vdorie/dbarts) Fast R package
-- [bartMachine](https://github.com/kapelner/bartMachine) Fast R package, supports missing predictors imputation
-- [SoftBART](https://github.com/theodds/SoftBART) R package with a smooth version of BART
-- [bcf](https://github.com/jaredsmurray/bcf) R package for a version of BART for causal inference
-- [flexBART](https://github.com/skdeshpande91/flexBART) Fast R package, supports categorical predictors
-- [flexBCF](https://github.com/skdeshpande91/flexBCF) R package, version of bcf optimized for large datasets
-- [XBART](https://github.com/JingyuHe/XBART) R/Python package, XBART is a faster variant of BART
-- [BART](https://github.com/JingyuHe/BART) R package, BART warm-started with XBART
-- [XBCF](https://github.com/socket778/XBCF)
-- [BayesTree](https://cran.r-project.org/package=BayesTree) R package, original BART implementation
-- [bartCause](https://github.com/vdorie/bartCause) R package, pre-made BART-based workflows for causal inference
-- [stan4bart](https://github.com/vdorie/stan4bart)
-- [VCBART](https://github.com/skdeshpande91/VCBART)
-- [monbart](https://github.com/jaredsmurray/monbart)
-- [mBART](https://github.com/remcc/mBART_shlib)
-- [SequentialBART](https://github.com/mjdaniels/SequentialBART)
-- [sparseBART](https://github.com/cspanbauer/sparseBART)
-- [pymc-bart](https://github.com/pymc-devs/pymc-bart)
-- [semibart](https://github.com/zeldow/semibart)
-- [CSP-BART](https://github.com/ebprado/CSP-BART)
-- [AMBARTI](https://github.com/ebprado/AMBARTI)
-- [MOTR-BART](https://github.com/ebprado/MOTR-BART)
-- [bcfbma](https://github.com/EoghanONeill/bcfbma)
-- [bartBMAnew](https://github.com/EoghanONeill/bartBMAnew)
-- [BART-BMA](https://github.com/BelindaHernandez/BART-BMA) (superseded by bartBMAnew)
-- [gpbart](https://github.com/MateusMaiaDS/gpbart)
-- [GPBART](https://github.com/nchenderson/GPBART)
-- [bartpy](https://github.com/JakeColtman/bartpy)
-- [BayesTreePrior](https://github.com/AlexiaJM/BayesTreePrior)
-- [BayesTree.jl](https://github.com/mathcg/BayesTree.jl)
-- [longbet](https://github.com/google/longbet)
-

bartz-0.4.0.dist-info/RECORD

@@ -1,13 +0,0 @@
-bartz/BART.py,sha256=CbGzFWtYw5u38Z9-Hy3CbDXpKOOvPFAAkSqu2HZl8no,16862
-bartz/__init__.py,sha256=E96vsP0bZ8brejpZmEmRoXuMsUdinO_B_SKUUl1rLsg,1448
-bartz/_version.py,sha256=2eiWQI55fd-roDdkt4Hvl9WzrTJ4xQo33VzFud6D03U,22
-bartz/debug.py,sha256=9ZH-JfwZVu5OPhHBEyXQHAU5H9KIu1vxLK7yNv4m4Ew,5314
-bartz/grove.py,sha256=x_6NK_l-hrXfy1PhssYNJkX41-w_WqjDziww0E7YRS8,8500
-bartz/jaxext.py,sha256=RcVWTCGS8lXF7GBsNbKrpuA4MTcokItq0CpWm3s7CGk,12033
-bartz/mcmcloop.py,sha256=lKDszvniNXka99X3e9RCrTgvEAZHA7ZbVXEgxUYvKMY,7634
-bartz/mcmcstep.py,sha256=diI9vHXHMvu_Lk_bSJ-a038OnEbXDpNEikVPhRcxEys,54987
-bartz/prepcovars.py,sha256=mMgfL-LGJ_8QpOL6iy7yfkL8A7FrT7Zfn5M3voyNwSQ,5818
-bartz-0.4.0.dist-info/LICENSE,sha256=heuIJZQK9IexJYC-fYHoLUrgj8HG8yS3G072EvKh-94,1073
-bartz-0.4.0.dist-info/METADATA,sha256=K86CVXT6ayPnc2hjhreYGMEeYWfYJIZdDkKuBB0-FYA,4500
-bartz-0.4.0.dist-info/WHEEL,sha256=sP946D7jFCHeNz5Iq4fL4Lu-PrWrFsgfLXbbkciIZwg,88
-bartz-0.4.0.dist-info/RECORD,,