Lung cancer cases due to radon in Europe: Difference between revisions

Revision as of 15:02, 26 January 2011

This page is a knowledge crystal of subtype variable. The page identifier is Op_en4715
Moderator:Teemu R (see all)
This page is a stub. You may improve it into a full page.
Upload data {{#opasnet_base_link:Op_en4715}}

Scope

Mortality due to indoor radon concentrations.

Spatial: Europe
Temporal: years 2010-2050

Definition

Lung cancer cases calculated from radon concentration using <math>cases = relative risk * background rate of cases * concentration * population</math>. RR taken from HEIMTSA and INTARESE (Darby 2004; Darby 2005).
- Background rate of new lung cancers is assumed constant for the whole Europe: 58.2 cases per 100,000 population (Globocan 2008, crude rate).
- Linear no threshold model assumed.

Data

Dependencies

Radon concentrations in European residences
Population of Europe
Lung cancer mortality in Europe (from WHO mortality data)
ERF of radon exposure on lung cancer (heande:ERFs of several pollutants)

Unit

cases per year

Formula

This code features R functions described on pages Opasnet Base Connection for R and Operating intelligently with multidimensional arrays in R.

conc <- op_baseGetData("opasnet_base", "Op_en4713")
pop <- op_baseGetData("opasnet_base", "Op_en4691", include = 1367, exclude = c(1435, 1436))
erf <- op_baseGetData("opasnet_base", "Erac2987", include = 48847)
countries <- c("Austria", "Belgium", "Bulgaria", "Switzerland", "Cyprus", "Czech Republic", "Germany", "Denmark", "Estonia", "Spain", 
	"Finland", "France", "Greece", "Hungary", "Ireland", "Iceland", "Italy", "Lithuania", "Luxembourg", "Latvia", "Malta", "Netherlands", 
	"Norway", "Poland", "Portugal", "Romania", "Sweden", "Slowenia", "Slovakia", "United Kingdom")
levels(pop[,"CountryID"]) <- countries #Country IDs converted to names for compatibility
colnames(pop)[4] <- "Country"
colnames(pop)[7] <- "Population"
concarray <- DataframeToArray(conc)
popxconc <- IntArray(pop, concarray, "Concentration")#[1:1000,,,]
erfarray <- DataframeToArray(erf[,c("obs","Result")])
popxconcxerf <- IntArray(popxconc, erfarray, "ERF")
lungmortality <- data.frame(popxconcxerf[,c("obs","Country","policy","Year","Age","Sex")], Result = popxconcxerf[,"ERF"] * 58.2 * 
popxconcxerf[,"Concentration"] / 100 * popxconcxerf[, "Population"] / 100000)

Result

**Lung cancer cases in Europe due to indoor radon in residences (mean and 95% confidence interval).**
	Year
Policy	2010	2020	2030	2050
BAU	42691 (22042-95542)	51138 (24733-111979)	57618 (25785-137181)	62461 (27023-143632)
All	NA	53385 (25402-130453)	68338 (28863-171718)	81873 (31372-239696)
Biomass	NA	NA	NA	81849 (32497-218247)
Insulation	NA	NA	NA	80534 (32445-236962)
Renovation	NA	NA	NA	94075 (35396-288323)

**Lung cancer cases attributable to indoor radon in residences in Europe, year 2010.**
Country of observation	Mean	SD
Austria	1105	1146
Belgium	857	796
Bulgaria	303	382
Switzerland	1636	2835
Cyprus	10	14
Germany	4754	4305
Denmark	382	408
Estonia	222	239
Spain	8245	15249
Finland	805	792
France	8204	12176
Greece	857	977
Hungary	1583	1985
Ireland	533	580
Italy	5035	5117
Lithuania	223	233
Luxembourg	68	64
Latvia	229	289
Malta	43	50
Netherlands	504	338
Norway	566	667
Poland	2175	1864
Portugal	1183	1208
Romania	1233	1357
Sweden	1353	1475
Slovakia	586	661
Total	42691

Keywords

References

Related files

@@ Line 11: / Line 11: @@
 *Lung cancer cases calculated from radon concentration using <nowiki><math>cases = relative risk * background rate of cases * concentration * population</math></nowiki>. RR taken from HEIMTSA and INTARESE (Darby 2004; Darby 2005).
 ** Background rate of new lung cancers is assumed constant for the whole Europe: 58.2 cases per 100,000 population (Globocan 2008, crude rate).
-** Mathematical formula is a simplification. Linear no threshold model assumed.
+** Linear no threshold model assumed.
 === Data ===
@@ Line 31: / Line 31: @@
   <nowiki>
+conc <- op_baseGetData("opasnet_base", "Op_en4713")
 pop <- op_baseGetData("opasnet_base", "Op_en4691", include = 1367, exclude = c(1435, 1436))
+erf <- op_baseGetData("opasnet_base", "Erac2987", include = 48847)
 countries <- c("Austria", "Belgium", "Bulgaria", "Switzerland", "Cyprus", "Czech Republic", "Germany", "Denmark", "Estonia", "Spain",
 	"Finland", "France", "Greece", "Hungary", "Ireland", "Iceland", "Italy", "Lithuania", "Luxembourg", "Latvia", "Malta", "Netherlands",
 	"Norway", "Poland", "Portugal", "Romania", "Sweden", "Slowenia", "Slovakia", "United Kingdom")
-conc <- op_baseGetData("opasnet_base", "Op_en4713")
+levels(pop[,"CountryID"]) <- countries #Country IDs converted to names for compatibility
-levels(pop[,"CountryID"]) <- countries
 colnames(pop)[4] <- "Country"
-colnames(pop)[8] <- "Population"
+colnames(pop)[7] <- "Population"
-pop <- pop[,c(3,4,5,6,7,8)]
-levels(pop[,"Age"])[1] <- "All Ages" #Fixed to match mortality format
 concarray <- DataframeToArray(conc)
-popxconc <- IntArray(pop, concarray[1:1000,,,], "Concentration")
+popxconc <- IntArray(pop, concarray, "Concentration")#[1:1000,,,]
-k <- rnorm(nrow(popxconc), 0.16, (0.31-0.05)/3.92)*58.2 #RR * background rate
+erfarray <- DataframeToArray(erf[,c("obs","Result")])
-lungmortality <- data.frame(popxconc[,c("obs","Country","policy","Year","Age","Rate","Sex")], Result = k * popxconc[,
+popxconcxerf <- IntArray(popxconc, erfarray, "ERF")
-"Concentration"] / 100 * popxconc[, "Population"] / 100000)</nowiki>
+lungmortality <- data.frame(popxconcxerf[,c("obs","Country","policy","Year","Age","Sex")], Result = popxconcxerf[,"ERF"] * 58.2 *
+popxconcxerf[,"Concentration"] / 100 * popxconcxerf[, "Population"] / 100000)</nowiki>
 == Result ==

Lung cancer cases due to radon in Europe: Difference between revisions

Revision as of 15:02, 26 January 2011

Contents

Scope

Definition

Data

Dependencies

Unit

Formula

Result

See also

Keywords

References

Related files

Navigation menu