I tried to sample RVs from a customized likelihood whose parameters are known. In other word, I want to do something like p0 = 0.3 #true bernoulli parameter obs = stats.bernoulli.rvs(p0, size=1000) #obs comes from bernoulli distribution but I need to modify two things: (1) use a customized likeli…

Can not correctly sample high dimensional RVs from customized likelihood

aakhmetz October 23, 2022, 1:59am 3

Just curious, if you simply write

c = pm.ZeroInflatedPoisson('c',psi,mu,observed=x)

would something be changed?

1 Like

Topic		Replies	Views
How to set up a custom likelihood function for two variables Questions	8	13306	March 9, 2018
Sample posterior predictive error Questions	3	1209	October 8, 2019
How to model product of multiple random variables Questions	8	4252	June 28, 2019
Likelihood given by product of two distributions? Questions	6	2215	August 4, 2018
Can pymc3 get samples from a (not posterior) complex distribution using MCMC? Questions	14	1622	July 2, 2019