python (12.9k questions)

javascript (9.2k questions)

reactjs (4.7k questions)

java (4.2k questions)

java (4.2k questions)

c# (3.5k questions)

c# (3.5k questions)

html (3.3k questions)

Questions - mcmc

MCMC for estimating negative binomial distribution

I want to estimate parameters of negative binomial distribution using MCMC Metropolis-Hastings algorithm. In other words, I have sample: set.seed(42) y <- rnbinom(20, size = 3, prob = 0.2) and I w...
test-img

John

r

optimization

random

mcmc

Votes: 0

Answers: 1

Latest Answer

Change dnorm in loglikelihood to dnbinom and fix the proposal for prob so it doesn't go outside (0,1): set.seed(42) y <- rnbinom(20, size = 3, prob = 0.2) prior_r <- function(r) { return(dpoi...
test-img

jblood94

Sampling a product of 4 independent functions where two of them are constant - with emcee

I need to sample from the following function: p = p1* p2* p3*p4 Where p1 and p2 are Gaussian functions with the same sigma but different means. But both p3 and p4 are equal to 1/(2*(pi)), both consta...
test-img

Matthew D.

python

sampling

mcmc

emcee

Votes: 0

Answers: 0

Python: MCMC fit to data with variable uncertainties (using lmfit)

Please bear with me as I try to explain what I'm trying to do. I'm trying fit an arctangent model to some data. I have two independent measurements in my dataset; one of these has unknown uncertaintie...
test-img

Ekbal Hussain

python

mcmc

lmfit

emcee

Votes: 0

Answers: 1

Latest Answer

First, please give a more complete example. As the lmfit documentation shows, you can provide uncertainties for the data. If you think you don't know the uncertainties, then try setting them to "...
test-img

M Newville

Posts

Questions

Blogs

Jobs