Computes the parameters and mixing probabilities for a Poisson sampling problem where the prior is a discrete mixture of gamma densities.
poisson.gamma.mix(probs,gammapar,data)
vector of probabilities of the gamma components of the prior
matrix where each row contains the shape and rate parameters for a gamma component of the prior
list with components y, vector of counts, and t, vector of time intervals
vector of probabilities of the gamma components of the posterior
matrix where each row contains the shape and rate parameters for a gamma component of the posterior
# NOT RUN { probs=c(.5, .5) gamma.par1=c(1,1) gamma.par2=c(10,2) gammapar=rbind(gamma.par1,gamma.par2) y=c(1,3,2,4,10); t=c(1,1,1,1,1) data=list(y=y,t=t) poisson.gamma.mix(probs,gammapar,data) # }
Run the code above in your browser using DataLab