A function that determines the posterior expectations
DS.posterior.reduce(DS.GF.obj, exposure)
Object resulting from running DS.prior function on a data set.
In the case of the Poisson family with exposure, represents the exposure values for the count data.
Returns
Mukhopadhyay, S. and Fletcher, D., 2018. "Generalized Empirical Bayes via Frequentist Goodness of Fit," Nature Scientific Reports, 8(1), p.9983, https://www.nature.com/articles/s41598-018-28130-5 .
# NOT RUN {
data(rat)
rat.start <- gMLE.bb(rat$y, rat$n)$estimate
rat.ds <- DS.prior(rat, max.m = 4, rat.start, family = "Binomial")
DS.posterior.reduce(rat.ds)
# }
Run the code above in your browser using DataLab