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

From Opasnet
Jump to navigation Jump to search
Line 31: Line 31:
*This code features [[R]] functions described on pages [[Opasnet Base Connection for R]] and [[Operating intelligently with multidimensional arrays in R]].
*This code features [[R]] functions described on pages [[Opasnet Base Connection for R]] and [[Operating intelligently with multidimensional arrays in R]].


<nowiki>
<rcode>
library(OpasnetBaseUtils)
conc <- op_baseGetData("opasnet_base", "Op_en4713")
conc <- op_baseGetData("opasnet_base", "Op_en4713")
pop <- op_baseGetData("opasnet_base", "Op_en4691", include = 1367, exclude = c(1435, 1436))
pop <- op_baseGetData("opasnet_base", "Op_en4691", include = 1367, exclude = c(1435, 1436))
Line 46: Line 47:
popxconcxerf <- IntArray(popxconc, erfarray, "ERF")
popxconcxerf <- IntArray(popxconc, erfarray, "ERF")
lungmortality <- data.frame(popxconcxerf[,c("obs","Country","policy","Year","Age","Sex")], Result = popxconcxerf[,"ERF"] * 58.2 *  
lungmortality <- data.frame(popxconcxerf[,c("obs","Country","policy","Year","Age","Sex")], Result = popxconcxerf[,"ERF"] * 58.2 *  
popxconcxerf[,"Concentration"] / 100 * popxconcxerf[, "Population"] / 100000)</nowiki>
popxconcxerf[,"Concentration"] / 100 * popxconcxerf[, "Population"] / 100000)
</rcode>
 
====An alternative code for Finland====
 
This code uses the summary table approach for obtaining data. The data sources are the same except [[ERF for long-term indoor exposure to radon and lung cancer]] (summary table [[:Category:Exposure-response functions]]) and [[Exposure to indoor radon in Finland]] (summary table [[:Category:Exposures]]).
 
<rcode include="page:OpasnetBaseUtils|name:generic">
 
############### Bring parts of summary table
#  data is the summary table.
summary.bring <- function(data){
pages <- levels(data$Page)
 
out <- op_baseGetData("opasnet_base", pages[1])
if(nlevels(pages) > 1){
for(i in 2:nlevels(pages)){
out[[i]] <- op_baseGetData("opasnet_base", pages[i])
}}
return(out)
}
##########
 
cat("Get data from Opasnet Base.\n")
library(OpasnetBaseUtils)
conc <- op_baseGetData("opasnet_base", "Op_en1619") # Summary table Category:Exposures
#pop <- op_baseGetData("opasnet_base", "Op_en4691", include = 1367, exclude = c(1435, 1436)) #
erf <- op_baseGetData("opasnet_base", "Op_en1914") # Summary table Category:Exposure-response functions
 
exposures <- summary.bring(conc)
exposures
 
erfs <- summary.bring(erf)
erfs
 
 
cat("Adjust data for the model.\n")
 
#conc <- conc[conc$Page == "Op_en5450", ] # Slice exposures to radon
conc
temp <- reshape(conc, idvar = c("Page", "Index"), timevar = "index", conc[conc$Observation == "Location", "Result.Text"]
temp
names(temp) <- conc[conc$Observation == "Location", "Index"]
as.data.frame(temp)
t(as.data.frame(temp))
 
 
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"
 
cat("Run the actual model.\n")
 
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)
</rcode>


== Result ==
== Result ==

Revision as of 12:39, 2 January 2012


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.
    • In a given iteration, radon concentrations are assumed equal in all residences in a country.

Data

Dependencies

Unit

cases per year

Formula

+ Show code

An alternative code for Finland

This code uses the summary table approach for obtaining data. The data sources are the same except ERF for long-term indoor exposure to radon and lung cancer (summary table Category:Exposure-response functions) and Exposure to indoor radon in Finland (summary table Category:Exposures).

+ Show code

Result

{{#opasnet_base_link:Op_en4715}}


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

See also

Keywords

References


Related files

<mfanonymousfilelist></mfanonymousfilelist>

Retrieved from "https://en.opasnet.org/index.php?title=Lung_cancer_cases_due_to_radon_in_Europe&oldid=23215"