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

From Opasnet
Jump to navigation Jump to search
(→‎An alternative code for Finland: more fancy stuff. Population and disease from summary tables; redundant columns shown as summary)
mNo edit summary
 
(3 intermediate revisions by 2 users not shown)
Line 1: Line 1:
{{variable|moderator=Teemu R|stub=Yes}}
{{variable|moderator=Teemu R|stub=Yes}}
[[Category:Code under inspection]]


== Scope ==
== Scope ==
Line 63: Line 64:
library(xtable)
library(xtable)


conc <- op_baseGetData("opasnet_base", "Op_en1619") # Summary table Category:Exposures
conc <- summary.bring("Op_en1619")[show, ] # Summary table Category:Exposures
erf <- op_baseGetData("opasnet_base", "Op_en1914") # Summary table Category:Exposure-response functions
erf <- summary.bring("Op_en1914") # Summary table Category:Exposure-response functions
pop <- op_baseGetData("opasnet_base", "Op_en2154") # Summary table Category:Population characteristics
pop <- summary.bring("Op_en2154") # Summary table Category:Population characteristics
bg <- op_baseGetData("opasnet_base", "Op_en1620") # Summary table Category:Health effects
bg <- summary.bring("Op_en1620") # Summary table Category:Health effects


cat("Adjust data for the model.\n")
cat("Adjust data for the model.\n")


conc <- summary.bring(conc)[show, ]
print(xtable(conc), type = 'html')
print(xtable(conc), type = 'html')
erf <- summary.bring(erf)
print(xtable(erf), type = 'html')
print(xtable(erf), type = 'html')
pop <- summary.bring(pop)
print(xtable(pop), type = 'html')
print(xtable(pop), type = 'html')
bg <- summary.bring(bg)
print(xtable(bg), type = 'html')
print(xtable(bg), type = 'html')


Line 85: Line 79:
out <- merge(out, pop)
out <- merge(out, pop)
out <- merge(out, bg)
out <- merge(out, bg)
#print(xtable(out), type = 'html')


out$ERF <- as.numeric(substr(out$ERF, 1, 6))
out$ERF <- as.numeric(substr(out$ERF, 1, 6))
Line 92: Line 85:


cat("Run the actual model.\n")
cat("Run the actual model.\n")
out$Out <- (out$ERF - 1) / out$ERF * out$Exposure * out$Response * out$Amount / 100000
#print(xtable(out), type = 'html')


out <- dropall(out)
out <- dropall(out)
print(xtable(t(out[1, ])), type = 'html')
tst <- rep(TRUE, ncol(out)) # tst: are all values in a column the same?
for(i in 1:ncol(out)){tst[i] <- nlevels(as.factor(out[, i])) == 1}
print(xtable(t(out[1, tst])), type = 'html')


tst <- 1:ncol(out)
out$Out <- (out$ERF - 1) / out$ERF * out$Exposure * out$Response * out$Amount / 100000
for(i in 1:ncol(out)){tst[i] <- nlevels(as.factor(out[, i])) == 1}
for(i in 1:ncol(out)){tst[i] <- nlevels(as.factor(out[, i])) == 1}
print(xtable(out[, !tst]), type = 'html')
print(xtable(out[, !tst]), type = 'html')
</rcode>
</rcode>
{{attack|# |For some reason Population in the data.frame pop is named "x", after with the matching goes wrong (not a critical error. However, the reason for this must be found out.|--[[User:Jouni|Jouni]] 15:58, 3 January 2012 (EET)}}
{{attack|# |Response in data.frame bg is text (or factor) instead of numeric. Why?|--[[User:Jouni|Jouni]] 15:58, 3 January 2012 (EET)}}


== Result ==
== Result ==

Latest revision as of 10:58, 26 August 2013

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

Which exposure parameters do you want to see?:
Actual population average in Finland
Guidance value for new apartments
Guidance value for old apartments

+ 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=30647"