cancerRates {diseasemapping} | R Documentation |
Download cancer incidence rates from the International Agency for Research on Cancer (IARC)
Description
Rates by age and sex group are retreived from http://ci5.iarc.fr/CI5plus/ci5plus.htm
Usage
cancerRates(area = "canada", year=2000, sex=c("M", "F"), site="Lung")
Arguments
area |
Region to retrieve rates from, |
year |
year or sequence of years to retrieve data from, within the period 1978 to 2002 |
site |
a vector of cancer sites, see details |
sex |
|
Details
area
must be one of Canada, Norway,
Latvia,
Lithuania,
Iceland,
Finland,
Estonia,
Denmark,
"Czech Republic",
"Costa Rica",
USA,
Iowa,
"New Mexico"
or the Canadian provinces of
Newfoundland, Prince Edward Island,
Nova Scotia,
Ontario,
Manitoba,
Saskatchewan,
Alberta, and
British Columbia. Alternately an integer specifying a registry code from http://ci5.iarc.fr.
site
must be one or more of
All Sites,
Oral Cavity and Pharynx,
Oesophagus.
Stomach,
Colon,
Rectum and Anus,
Liver,
Gallbladder,
Pancreas,
Larynx,
Lung,
Bone,
Melanoma of skin,
Prostate (Males only),
Testis (Males only),
Breast (Females only),
Cervix uteri (Females only),
Corpus uteri (Females only),
Ovary and other uterine adnexa (Females only),
Kidney,
Bladder,
Eye,
Brain and Central Nervous System,
Thyroid,
Non-Hodgkin Lymphoma,
Hodgkin Lymphoma,
Multiple myeloma,
Leukaemia.
Value
vector of cancer rates, by age and sex group
Examples
# won't work if offline or if the iarc web site is down
qcLungF=try(cancerRates(area="canada",
year=2001:2002, site="lung", sex="F"), silent=TRUE)
if(length(grep("try-error", class(qcLungF)))) {
qcLungF = structure(c(0, 5e-06, 0, 0, 5e-06, 1e-05, 0, 3.4e-05, 9.6e-05,
0.000211, 0.000559, 0.001289, 0.002003, 0.002508, 0.002728, 0.003189,
0.002792, 0.001905), .Names = c("F_0", "F_5", "F_10", "F_15",
"F_20", "F_25", "F_30", "F_35", "F_40", "F_45", "F_50", "F_55",
"F_60", "F_65", "F_70", "F_75", "F_80", "F_85"), site = "Lung",
area = "Canada", year = "2001-2002")
}
qcLungF
data('popdata')
popdata = terra::unwrap(popdata)
qcLungExp = getSMR(popdata, qcLungF)
names(qcLungExp)
if(require('mapmisc', quietly=TRUE)) {
mycol = colourScale(qcLungExp$expected,
breaks=12, dec=0, style='quantile')
plot(popdata[1:400,])
plot(qcLungExp, col=mycol$plot, border='#00000040',add=TRUE)
legendBreaks('topright', mycol)
}