Empirical Bayes Estimates of Relative Risk

Share:

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

1
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 mean 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

See Also

scotland, mapvariable

Examples

1
2
3
4
5
6
7
8
9
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)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.