I would like to use the following distribution as a part of my model import pymc as pm from pymc.math import log, exp def _podium_logp(value, mu, sigma): dist = pm.Gamma.dist(mu=mu, sigma=sigma) density1 = exp(pm.logcdf(dist, value + 3)) - exp(pm.logcdf(dist, value - 2)) density2 = ex…

Implementing rounding (by manual integration) more efficiently

ricardoV94 February 19, 2024, 10:33am 5

I mean pytensor.tensor.vectorize which works just like numpy.

Topic		Replies	Views
Dynamic shaping, "round" function, JAX, and a "few" more questions v5 modeling , jax , pytensor	28	1641	September 29, 2023
Custom discrete distribution for zero truncated Poisson Questions	10	2608	December 19, 2017
pm.DensityDist problem under PMYC 4 v5	17	2478	June 14, 2022
DensityDist from scipy.stats distribution Questions	5	2287	May 20, 2018
Why would scipy/numpy wrapped in @as_op cause much faster sampling than using pytensor operations? Questions	10	1038	April 29, 2023