# NOT RUN {
set.seed(1234)
DataHist <- rpois(n = 100, lambda = 49)
set.seed(12345)
DataCur <- rpois(n = 100, lambda = 49)
PoissonNPP_MCMC(Data.Cur = DataCur, Data.Hist = DataHist,
CompStat = list(n0 = 20, mean0 = 10, n1 = 30, mean1 = 11),
prior = list(lambda.shape = 1/2, lambda.scale = 100,
delta.alpha = 1, delta.beta = 1),
MCMCmethod = 'RW', rw.logit.delta = 1,
ind.delta.alpha= 1, ind.delta.beta= 1,nsample = 10000,
control.mcmc = list(delta.ini = NULL, burnin = 2000, thin = 1))
# }
Run the code above in your browser using DataLab