# getRates: Calculate the estimated coefficients of age and sex group... In diseasemapping: Modelling Spatial Variation in Disease Risk for Areal Data

 getRates R Documentation

## Calculate the estimated coefficients of age and sex group from the glm model

### Description

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.

### Usage

``````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), fit.numeric=NULL, breaks = NULL)
``````

### Arguments

 `casedata` A data frame of case data, with columns corresponding to variables in `formula`. Assumed to be one row per case, unless a column called `y` or `cases` or `count` is included, in which case this column gives the number of cases per row. `popdata` population data set `formula` the glm model you want to fit. ie. `~age*sex` `family` the distribution to fit the model `minimumAge` the lower boundary of the age, default is 0 `maximumAge` the higher boundary of the age, default is 100 `S` vector of sexes to include in the analysis. Defaults to both "M" and "F" `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

### Details

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.

### Value

A summary of the glm model contains set of estimated coefficients for different age and sex groups.

Patrick Brown

### Examples

``````
data('casedata')
data('popdata')
popdata = terra::unwrap(popdata)
therates = getRates(casedata, popdata, ~sex*age,
breaks=c(seq(0, 80, by=10), Inf))
therates

``````

diseasemapping documentation built on Sept. 22, 2023, 1:07 a.m.