cancerRates: Download cancer incidence rates from the International Agency...

View source: R/cancerRates.R

cancerRatesR 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

"M" or "F" for male or female rates only, c("M","F") (the default) for both sexes.

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)
}


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