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

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)

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.