eBayes: Empirical Bayes Estimates of Relative Risk

In SpatialEpi: Methods and Data for Spatial Epidemiology

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

Y	a length n vector of observed cases
E	a length n vector of expected number of cases
Xmat	n x p dimension matrix of covariates

Value

A list with 5 elements:

RR	the ecological relative risk posterior mean estimates
RRmed	the ecological relative risk posterior median estimates
beta	the MLE's of the regression coefficients
alpha	the MLE of negative binomial dispersion parameter
SMR	the 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

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)

SpatialEpi documentation built on March 7, 2023, 8 p.m.