Lung cancer cases due to radon in Europe: Difference between revisions
Jump to navigation
Jump to search
(→An alternative code for Finland: code cleaned, functions moved to OpasnetBaseUtils) |
(→An alternative code for Finland: more fancy stuff. Population and disease from summary tables; redundant columns shown as summary) |
||
| Line 65: | Line 65: | ||
conc <- op_baseGetData("opasnet_base", "Op_en1619") # Summary table Category:Exposures | conc <- op_baseGetData("opasnet_base", "Op_en1619") # Summary table Category:Exposures | ||
erf <- op_baseGetData("opasnet_base", "Op_en1914") # Summary table Category:Exposure-response functions | erf <- op_baseGetData("opasnet_base", "Op_en1914") # Summary table Category:Exposure-response functions | ||
pop <- op_baseGetData("opasnet_base", "Op_en2154") # Summary table Category:Population characteristics | |||
bg <- op_baseGetData("opasnet_base", "Op_en1620") # Summary table Category:Health effects | |||
cat("Adjust data for the model.\n") | cat("Adjust data for the model.\n") | ||
| Line 73: | Line 75: | ||
erf <- summary.bring(erf) | erf <- summary.bring(erf) | ||
print(xtable(erf), type = 'html') | print(xtable(erf), type = 'html') | ||
pop <- summary.bring(pop) | |||
print(xtable(pop), type = 'html') | |||
bg <- summary.bring(bg) | |||
print(xtable(bg), type = 'html') | |||
out <- merge(conc, erf) | out <- merge(conc, erf) | ||
print(xtable(out), type = 'html') | out <- merge(out, pop) | ||
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)) | ||
out$Exposure <- as.numeric(substr(out$Exposure, 1, 3)) | out$Exposure <- as.numeric(substr(out$Exposure, 1, 3)) | ||
out$Response <- as.numeric(out$Response) | |||
cat("Run the actual model.\n") | cat("Run the actual model.\n") | ||
out$ | out$Out <- (out$ERF - 1) / out$ERF * out$Exposure * out$Response * out$Amount / 100000 | ||
out | #print(xtable(out), type = 'html') | ||
out | out <- dropall(out) | ||
print(xtable(out), type = 'html') | print(xtable(t(out[1, ])), type = 'html') | ||
tst <- 1:ncol(out) | |||
for(i in 1:ncol(out)){tst[i] <- nlevels(as.factor(out[, i])) == 1} | |||
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 == | ||
Revision as of 13:58, 3 January 2012
| Moderator:Teemu R (see all) |
| This page is a stub. You may improve it into a full page. |
| Upload data
|
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
- 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.
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).
⇤--#: . 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. --Jouni 15:58, 3 January 2012 (EET) (type: truth; paradigms: science: attack)
⇤--#: . Response in data.frame bg is text (or factor) instead of numeric. Why? --Jouni 15:58, 3 January 2012 (EET) (type: truth; paradigms: science: attack)
Result
{{#opasnet_base_link:Op_en4715}}
| 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) |
| 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>