getRates {diseasemapping} | R Documentation |
The getRates function calculates the estimated coefficient of the age and sex group from the case and population data set. It fits a glm model with poisson distribution by default.
getRates(casedata, popdata, formula, family = poisson, minimumAge = 0, maximumAge = 100, S = c("M", "F"), years = NULL, year.range = NULL, case.years = grep("^year$", names(casedata), ignore.case = TRUE, value = TRUE)[1], fit.numeric=NULL, breaks = NULL)
casedata |
case data set name |
popdata |
population data set name |
formula |
the glm model you want to fit. ie. $~age*sex$~ |
family |
the distribution to fit the model |
minimumAge |
the lower boundry of the age, default is 0 |
maximumAge |
the higher boundry of the age, default is 100 |
S |
the group of sex you want to choose, default is both male and female. "M" denotes male; "F" denotes female |
years |
a vector of census years |
year.range |
study period: a vector of two elements, starting dates and ending dates |
case.years |
variable name in the case data which contains time |
fit.numeric |
the variables which needed to be changed from factor to numeric |
breaks |
the age breaks |
It fits a glm model with poisson or binomial distribution over case and population data sets. If there is no data set in some age and sex group, an NA will show there.
A summary of the glm model contains set of estimated coefficients for different age and sex groups.
Patrick Brown
data(casedata) data(popdata) therates = getRates(casedata, popdata, ~age*sex,breaks=c(seq(0, 80, by=10), Inf)) therates