getRates {diseasemapping}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)[1], breaks = NULL) 

Arguments

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
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.

Author(s)

Patrick Brown

Examples

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

[Package diseasemapping version 0.4 Index]