Difference between revisions of "Asthma prevalence due to building dampness in Europe"

From Opasnet
Jump to: navigation, search
(Short summary added)
m (Formula)
 
(15 intermediate revisions by 2 users not shown)
Line 1: Line 1:
{{variable|moderator=Teemu R|stub=Yes}}
+
[[Category:Indoor air]]
 +
[[Category:Mega case study]]
 +
[[Category:Health impact]]
 +
{{variable|moderator=Teemu R}}
  
 
== Scope ==
 
== Scope ==
Line 7: Line 10:
 
== Definition ==
 
== Definition ==
  
Initial Prevalence data from GINA<ref name=GINA>[http://www.ginasthma.com/download.asp?intId=29 Global Initiative for Asthma report: Global Burden of Asthma]</ref>.  
+
This variable is calculated entirely from upstream variables listed under [[#Dependencies]]. Mathematical method described under [[#Formula]].
 
 
Formula calculates number of asthma cases which can be attributed to dampness in Europe in years 2010-2050 by using the following:
 
#<nowiki><math>extra cases = (prevalence_total - prevalence_nondamp) population</math></nowiki>
 
#<nowiki><math>prevalence_total = prevalence_nondamp(%damp * OR + 1 - %damp)</math></nowiki>
 
#<nowiki><math>prevalence_nondamp = \frac{prevalence_0}{%damp * OR + 1 - %damp}</math></nowiki>
 
 
 
Population data[http://en.opasnet.org/w/Population_of_Europe] included population growth scenarios so they too are included.  
 
  
 
=== Data ===
 
=== Data ===
Line 22: Line 18:
 
*[[:heande:Moisture damage]]
 
*[[:heande:Moisture damage]]
 
*[[Population of Europe by Country|Population of Europe]]
 
*[[Population of Europe by Country|Population of Europe]]
*Prevelance of Clinical Asthma<ref name=GINA></ref>
+
*[[Asthma prevalence]]
 
*[[ERF of indoor dampness on respiratory health effects]]
 
*[[ERF of indoor dampness on respiratory health effects]]
  
Line 31: Line 27:
 
=== Formula ===
 
=== Formula ===
  
#<nowiki><math>extra cases = (prevalence_total - prevalence_nondamp) population</math></nowiki>
+
Formula calculates number of asthma cases which can be attributed to dampness in Europe in years 2010-2050 by using the following:
#<nowiki><math>prevalence_total = prevalence_nondamp(%damp * OR + 1 - %damp)</math></nowiki>
+
#<math>extra cases = (prevalencetotal - prevalencenondamp) * population</math>
#<nowiki><math>prevalence_nondamp = \frac{prevalence_0}{%damp * OR + 1 - %damp}</math></nowiki>
+
#<math>prevalencetotal = prevalencenondamp * (%damp * OR + 1 - %damp)</math>
 +
#<math>prevalence_0 = %damp * OR * prevalencenondamp + (1 - %damp) * prevalencenondamp</math>
 +
::<math>prevalencenondamp = \frac{prevalence_0}{%damp * OR + 1 - %damp}</math>
  
*ERF approximated as that for current asthma (1.56).  
+
*Prevalences for nondamp and damp homes are assumed constant, for a given iteration in a given country.  
  
 
  <nowiki>
 
  <nowiki>
Line 43: Line 41:
 
"Finland", "France", "Greece", "Hungary", "Ireland", "Iceland", "Italy", "Lithuania", "Luxembourg", "Latvia", "Malta", "Netherlands",  
 
"Finland", "France", "Greece", "Hungary", "Ireland", "Iceland", "Italy", "Lithuania", "Luxembourg", "Latvia", "Malta", "Netherlands",  
 
"Norway", "Poland", "Portugal", "Romania", "Sweden", "Slowenia", "Slovakia", "United Kingdom")
 
"Norway", "Poland", "Portugal", "Romania", "Sweden", "Slowenia", "Slovakia", "United Kingdom")
levels(pop[,"CountryID"]) <- countries
+
levels(pop[,"CountryID"]) <- countries #IDs converted to actual names, for compatibility with other data
colnames(pop)[4] <- "Country"
+
colnames(pop)[c(4,7)] <- c("Country","Population")
colnames(pop)[8] <- "Population"
+
asthma <- op_baseGetData("opasnet_base", "Op_en4789")
#asthma <- read.csv("C:/Documents and Settings/tris/My Documents/Asthma prevalence.csv", sep = ";")
+
erf <- op_baseGetData("opasnet_base", "Op_en4716")
asthma <- data.frame(Country=c('Scotland','Jersey','Guernsey','Wales','Isle of Man','England','New Zealand','Australia','Republic of Ireland',
 
'Canada','Peru','Trinidad & Tobago','Costa Rica','Brazil','United States of America','Fiji','Paraguay','Uruguay','Israel','Barbados','Panama',
 
'Kuwait','Ukraine','Ecuador','South Africa','Czech Republic','Finland','Malta','Ivory Coast','Colombia','Turkey','Lebanon','Kenya','Germany',
 
'France','Norway','Japan','Sweden','Thailand','Hong Kong','Philippines','United Arab Emirates','Belgium','Austria','Spain','Saudi Arabia',
 
'Argentina','Iran','Estonia','Nigeria','Chile','Singapore','Malaysia','Portugal','Uzbekistan','FYR Macedonia','Italy','Oman','Pakistan',
 
'Tunisia','Cape Verde','Latvia','Poland','Algeria','South Korea','Bangladesh','Morocco','Occupied Territory of Palestine','Mexico','Ethiopia',
 
'Denmark','India','Taiwan','Cyprus','Switzerland','Russia','China','Greece','Georgia','Nepal','Romania','Albania','Indonesia','Macau'),
 
InitialPrevalence=c(18.4,17.6,17.5,16.8,16.7,15.3,15.1,14.7,14.6,14.1,13,12.6,11.9,11.4,10.9,10.5,9.7,9.5,9,8.9,8.8,8.5,8.3,8.2,8.1,8,8,8,7.8,7.4,
 
7.4,7.2,7,6.9,6.8,6.8,6.7,6.5,6.5,6.2,6.2,6.2,6,5.8,5.7,5.6,5.5,5.5,5.4,5.4,5.1,4.9,4.8,4.8,4.6,4.5,4.5,4.5,4.3,4.3,4.2,4.2,4.1,3.9,3.9,3.8,
 
3.8,3.6,3.3,3.1,3,3,2.6,2.4,2.3,2.2,2.1,1.9,1.8,1.5,1.5,1.3,1.1,0.7))
 
#asthma <- data.frame(Country=asthma[1:26,1], Casesper1000=(asthma[1:26,2]+asthma[1:26,3]))
 
erf <- 1.56
 
 
poparray <- DataframeToArray(pop, "Population")
 
poparray <- DataframeToArray(pop, "Population")
 
dampxpop <- IntArray(dampness, poparray, "Population")
 
dampxpop <- IntArray(dampness, poparray, "Population")
asthmaarray <- DataframeToArray(asthma, "InitialPrevalence")
+
asthmaarray <- DataframeToArray(asthma)
 
dampxpopxasthma <- IntArray(dampxpop, asthmaarray, "InitialPrevalence")
 
dampxpopxasthma <- IntArray(dampxpop, asthmaarray, "InitialPrevalence")
p_nd <- data.frame(dampxpopxasthma[dampxpopxasthma[,"Year"]=="2010", c(2,3,6,7,8)], p_nd=dampxpopxasthma[dampxpopxasthma[,"Year"]=="2010",  
+
erfarray <- DataframeToArray(erf)
"InitialPrevalence"] * 100 /(dampxpopxasthma[dampxpopxasthma[,"Year"]=="2010", "Result"] * erf + 100 - dampxpopxasthma[dampxpopxasthma[,"Year"]==
+
dampxpopxasthmaxerf <- IntArray(dampxpopxasthma, erfarray, "ERF")
"2010", "Result"]))
+
p_nd <- data.frame(dampxpopxasthmaxerf[dampxpopxasthmaxerf[,"Year"]=="2010", c("obs","Country","Age","Sex")],  
 +
p_nd=dampxpopxasthmaxerf[dampxpopxasthmaxerf[,"Year"]=="2010", "InitialPrevalence"] * 100 /(dampxpopxasthmaxerf[
 +
dampxpopxasthmaxerf[,"Year"]=="2010", "Result"] * dampxpopxasthmaxerf[dampxpopxasthmaxerf[,"Year"]=="2010",
 +
"ERF"] + 100 - dampxpopxasthmaxerf[dampxpopxasthmaxerf[,"Year"]=="2010", "Result"]))
 
p_ndarray <- DataframeToArray(p_nd, "p_nd")
 
p_ndarray <- DataframeToArray(p_nd, "p_nd")
dampxpopxasthmaxp_nd <- IntArray(dampxpopxasthma, p_ndarray, "p_nd")
+
dampxpopxasthmaxerfxp_nd <- IntArray(dampxpopxasthmaxerf, p_ndarray, "p_nd")
final <- data.frame(dampxpopxasthmaxk[,c(2,3,4,6,7,8)], Result=dampxpopxasthmaxp_nd[,"Population"] * dampxpopxasthmaxp_nd[,"p_nd"] *  
+
final <- data.frame(dampxpopxasthmaxerfxp_nd[,c("obs","Country","policy","Year","Age","Sex")],  
dampxpopxasthmaxp_nd[,"Result"] * (erf - 1) / 10000)</nowiki>
+
Result=dampxpopxasthmaxerfxp_nd[,"Population"] * dampxpopxasthmaxerfxp_nd[,"p_nd"] *  
 +
dampxpopxasthmaxerfxp_nd[,"Result"] * (dampxpopxasthmaxerfxp_nd[,"ERF"] - 1) / 10000)
 +
###Fancy alternative below, might be better, requires more testing
 +
#final <- model.frame(I(Result*Population*p_nd*(ERF-1)/10000)~ obs + Country + policy + Age + Sex + Outcome,
 +
#data = dampxpopxasthmaxerfxp_nd)</nowiki>
  
 
== Result ==
 
== Result ==
  
 
{{resultlink}}
 
{{resultlink}}
 +
 +
{|{{prettytable}}
 +
|+'''Asthma cases (prevalence) in Europe due to residential building dampness (mean and 95% confidence interval).'''
 +
! !!colspan="4"|Year
 +
|----
 +
!Policy!!2010!!2020!!2030!!2050
 +
|----
 +
|BAU || 1715846 (794208-2918407) || 2069089 (929518-3645690) || 2300513 (1007103-4193891) || 2417413 (1016202-4559645)
 +
|----
 +
|All || NA || 2071501 (940391-3650210) || 2634778 (1139578-4745158) || 3009693 (1251020-5519308)
 +
|----
 +
|Biomass || NA || NA || NA || 2998888 (1249803-5529395)
 +
|----
 +
|Insulation || NA || NA || NA || 3002498 (1239186-5524389)
 +
|----
 +
|Renovation || NA || NA || NA || 3416010 (1443227-6233562)
 +
|----
 +
|}
 +
 +
{| {{prettytable}}
 +
|+ '''Asthma cases (prevalence) attributable to residential building dampness in Europe in 2010.
 +
! Country of observation!! Mean!! SD
 +
|----
 +
|| Austria || 23661 || 22103
 +
|----
 +
|| Belgium || 46341 || 30025.2
 +
|----
 +
|| Cyprus || 2988 || 1251
 +
|----
 +
|| Czech Republic || 65025 || 39220
 +
|----
 +
|| Denmark|| 9051 || 7432
 +
|----
 +
|| Estonia|| 7828 || 3876
 +
|----
 +
|| Finland|| 10929 || 18613
 +
|----
 +
|| France|| 302344 || 201280
 +
|----
 +
|| Germany|| 375145 || 265165
 +
|----
 +
|| Greece|| 20343 || 10627
 +
|----
 +
|| Italy|| 276619 || 142459
 +
|----
 +
|| Latvia|| 11875 || 5228
 +
|----
 +
|| Poland|| 267934 || 105393
 +
|----
 +
|| Portugal|| 47961 || 25133
 +
|----
 +
|| Spain|| 227574 || 126839
 +
|----
 +
|| Sweden|| 20225 || 27113
 +
|----
 +
! Total!! 1715846 ||
 +
|----
 +
|}
  
 
==See also==
 
==See also==

Latest revision as of 12:23, 11 April 2011



Scope

Europe in 2010, 2020, 2030, 2050.

Definition

This variable is calculated entirely from upstream variables listed under #Dependencies. Mathematical method described under #Formula.

Data

Dependencies

Unit

#

Formula

Formula calculates number of asthma cases which can be attributed to dampness in Europe in years 2010-2050 by using the following:

  1. extra cases = (prevalencetotal - prevalencenondamp) * population
  2. prevalencetotal = prevalencenondamp * (%damp * OR + 1 - %damp)
  3. prevalence_0 = %damp * OR * prevalencenondamp + (1 - %damp) * prevalencenondamp
prevalencenondamp = \frac{prevalence_0}{%damp * OR + 1 - %damp}
  • Prevalences for nondamp and damp homes are assumed constant, for a given iteration in a given country.
dampness <- op_baseGetData("opasnet_base", "Erac2988")
pop <- op_baseGetData("opasnet_base", "Op_en4691", include = 1367, exclude = c(1435, 1436))
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 #IDs converted to actual names, for compatibility with other data
colnames(pop)[c(4,7)] <- c("Country","Population")
asthma <- op_baseGetData("opasnet_base", "Op_en4789")
erf <- op_baseGetData("opasnet_base", "Op_en4716")
poparray <- DataframeToArray(pop, "Population")
dampxpop <- IntArray(dampness, poparray, "Population")
asthmaarray <- DataframeToArray(asthma)
dampxpopxasthma <- IntArray(dampxpop, asthmaarray, "InitialPrevalence")
erfarray <- DataframeToArray(erf)
dampxpopxasthmaxerf <- IntArray(dampxpopxasthma, erfarray, "ERF")
p_nd <- data.frame(dampxpopxasthmaxerf[dampxpopxasthmaxerf[,"Year"]=="2010", c("obs","Country","Age","Sex")], 
	p_nd=dampxpopxasthmaxerf[dampxpopxasthmaxerf[,"Year"]=="2010", "InitialPrevalence"] * 100 /(dampxpopxasthmaxerf[
	dampxpopxasthmaxerf[,"Year"]=="2010", "Result"] * dampxpopxasthmaxerf[dampxpopxasthmaxerf[,"Year"]=="2010", 
	"ERF"] + 100 - dampxpopxasthmaxerf[dampxpopxasthmaxerf[,"Year"]=="2010", "Result"]))
p_ndarray <- DataframeToArray(p_nd, "p_nd")
dampxpopxasthmaxerfxp_nd <- IntArray(dampxpopxasthmaxerf, p_ndarray, "p_nd")
final <- data.frame(dampxpopxasthmaxerfxp_nd[,c("obs","Country","policy","Year","Age","Sex")], 
Result=dampxpopxasthmaxerfxp_nd[,"Population"] * dampxpopxasthmaxerfxp_nd[,"p_nd"] * 
dampxpopxasthmaxerfxp_nd[,"Result"] * (dampxpopxasthmaxerfxp_nd[,"ERF"] - 1) / 10000)
###Fancy alternative below, might be better, requires more testing
#final <- model.frame(I(Result*Population*p_nd*(ERF-1)/10000)~ obs + Country + policy + Age + Sex + Outcome, 
#data = dampxpopxasthmaxerfxp_nd)

Result

Show results


Asthma cases (prevalence) in Europe due to residential building dampness (mean and 95% confidence interval).
Year
Policy 2010 2020 2030 2050
BAU 1715846 (794208-2918407) 2069089 (929518-3645690) 2300513 (1007103-4193891) 2417413 (1016202-4559645)
All NA 2071501 (940391-3650210) 2634778 (1139578-4745158) 3009693 (1251020-5519308)
Biomass NA NA NA 2998888 (1249803-5529395)
Insulation NA NA NA 3002498 (1239186-5524389)
Renovation NA NA NA 3416010 (1443227-6233562)
Asthma cases (prevalence) attributable to residential building dampness in Europe in 2010.
Country of observation Mean SD
Austria 23661 22103
Belgium 46341 30025.2
Cyprus 2988 1251
Czech Republic 65025 39220
Denmark 9051 7432
Estonia 7828 3876
Finland 10929 18613
France 302344 201280
Germany 375145 265165
Greece 20343 10627
Italy 276619 142459
Latvia 11875 5228
Poland 267934 105393
Portugal 47961 25133
Spain 227574 126839
Sweden 20225 27113
Total 1715846

See also

Keywords

Asthma, indoor air, dampness, Europe

References


Related files

<mfanonymousfilelist></mfanonymousfilelist>