eBayes: Empirical Bayes Estimates of Relative Risk

Description

The computes empirical Bayes estimates of relative risk of study region with n areas, given observed and expected numbers of counts of disease and covariate information.

Usage

eBayes(Y, E, Xmat = NULL)

Arguments

a length n vector of observed cases

a length n vector of expected number of cases

Xmat

n x p dimension matrix of covariates

Value

A list with 5 elements:
RRthe ecological relative risk posterior mean estimates
RRmedthe ecological relative risk posterior mean estimates
betathe MLE's of the regression coefficients
alphathe MLE of negative binomial dispersion parameter
SMRthe standardized mortality/morbidity ratio Y/E

References

Clayton D. and Kaldor J. (1987) Empirical Bayes estimates of age-standardized relative risks for use in disease mapping. Biometrics, 43, 671--681

Examples

Run this code

data(scotland)
data <- scotland$data

x <- data$AFF
Xmat <- cbind(x,x^2)
results <- eBayes(data$cases,data$expected,Xmat)

scotland.map <- scotland$spatial.polygon
mapvariable(results$RR, scotland.map)