EBpostthresh: Produce the probabilities of exceeding a threshold given a...

Description Usage Arguments Value Author(s) See Also Examples

View source: R/EBpostthresh.R

Description

This function produces the posterior probabilities of exceeding a threshold given a gamma distributions with parameters (alpha+Y, (alpha+E*mu)/mu) where mu = exp(x beta). This model arises from Y being Poisson with mean theta times E where theta is the relative risk and E are the expected numbers. The prior on theta is gamma with parameters alpha and beta. The parameters alpha and beta may be estimated using empirical Bayes.

Usage

1
EBpostthresh(Y, E, alpha, beta, Xrow = NULL, rrthresh)

Arguments

Y

observed disease counts

E

expected disease counts

alpha
beta
Xrow
rrthresh

Value

Posterior probabilities of exceedence are returned.

Author(s)

Jon Wakefield

See Also

eBayes

Examples

1
2
3
4
5
6
7
8
data(scotland)
Y <- scotland$data$cases
E <- scotland$data$expected
ebresults <- eBayes(Y,E)
# Find probabilities of exceedence of 3
thresh3 <- EBpostthresh(Y, E, alpha=ebresults$alpha, beta=ebresults$beta,
rrthresh=3)
mapvariable(thresh3, scotland$spatial.polygon)

SpatialEpi documentation built on May 30, 2017, 5:23 a.m.