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

From Testiwiki
Jump to: navigation, search
m (Unit)
(Better data and revised formula)
Line 18: Line 18:
 
*[[:heande:Moisture damage]]
 
*[[:heande:Moisture damage]]
 
*[[Population of Europe by Country|Population of Europe]]
 
*[[Population of Europe by Country|Population of Europe]]
*Asthma incidence from the [http://data.euro.who.int/hmdb/index.php WHO: European hospital morbidity database (HMDB)]
+
*Prevelance of Clinical Asthma<ref>[http://www.ginasthma.com/download.asp?intId=29 Global Initiative for Asthma report: Global Burden of Asthma]</ref>
 
*[[ERF of indoor dampness on respiratory health effects]]
 
*[[ERF of indoor dampness on respiratory health effects]]
  
Line 27: Line 27:
 
=== Formula ===
 
=== Formula ===
  
*Asthma prevalence is estimated by summing in-patients and day-cases from the hospital morbidity database.
+
#<nowiki><math>extra cases = \left ( \frac{OR - 1}{\frac{1}{%damp}-1} \right ) d</math></nowiki>
*Wheeze ERF used, as it approximately corresponds to the type of cases likely to be registered by the database.  
+
#<nowiki><math>d = \frac{cases}{\frac{OR}{\frac{1}{%damp}-1}+1}</math></nowiki>
 +
 
 +
*ERF approximated as that for current asthma (1.56).  
  
 
  <nowiki>
 
  <nowiki>
Line 39: Line 41:
 
colnames(pop)[4] <- "Country"
 
colnames(pop)[4] <- "Country"
 
colnames(pop)[8] <- "Population"
 
colnames(pop)[8] <- "Population"
asthma <- read.csv("C:/Documents and Settings/tris/My Documents/Asthma cases hmdb.csv", sep = ";") #File constructed from hmdb data
+
asthma <- read.csv("C:/Documents and Settings/tris/My Documents/Asthma prevalence.csv", sep = ";")
asthma <- data.frame(Country=asthma[1:26,1], Casesper1000=(asthma[1:26,2]+asthma[1:26,3]))
+
#asthma <- data.frame(Country=asthma[1:26,1], Casesper1000=(asthma[1:26,2]+asthma[1:26,3]))
erf <- 1.50
+
erf <- 1.56
 
poparray <- DataframeToArray(pop, "Population")
 
poparray <- DataframeToArray(pop, "Population")
 
dampxpop <- IntArray(dampness, poparray, "Population")
 
dampxpop <- IntArray(dampness, poparray, "Population")
asthmaarray <- DataframeToArray(asthma, "Casesper1000")
+
asthmaarray <- DataframeToArray(asthma, "Prevalence")
dampxpopxasthma <- IntArray(dampxpop, asthmaarray, "Casesper1000")
+
dampxpopxasthma <- IntArray(dampxpop, asthmaarray, "Prevalence")
final <- data.frame(dampxpopxasthma[,c(2,3,4,6,7,8)], Result=(dampxpopxasthma[,"Result"] * dampxpopxasthma[,"Population"] * (erf - 1) *
+
final <- data.frame(dampxpopxasthma[,c(2,3,4,6,7,8)], Result=(dampxpopxasthma[,"Prevalence"] / 100 * dampxpopxasthma[,"Population"] / (erf / (100 /
dampxpopxasthma[,"Casesper1000"] / 1000))</nowiki>
+
dampxpopxasthma[,"Result"] - 1) + 1) * (erf - 1) / (100 / dampxpopxasthma[,"Result"] - 1)))</nowiki>
  
 
== Result ==
 
== Result ==

Revision as of 14:47, 13 December 2010


Scope

Europe in 2010, 2020, 2030, 2050.

Definition

Data

Description of the data used for obtaining the value of the variable
(e.g. measurement data; mathematical method and its parameters). 
Please include references (preferably using the <ref> </ref> tags) and links to original data, as appropriate.

Dependencies

Unit

#

Formula

  1. <math>extra cases = \left ( \frac{OR - 1}{\frac{1}{%damp}-1} \right ) d</math>
  2. <math>d = \frac{cases}{\frac{OR}{\frac{1}{%damp}-1}+1}</math>
  • ERF approximated as that for current asthma (1.56).
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
colnames(pop)[4] <- "Country"
colnames(pop)[8] <- "Population"
asthma <- read.csv("C:/Documents and Settings/tris/My Documents/Asthma prevalence.csv", sep = ";")
#asthma <- data.frame(Country=asthma[1:26,1], Casesper1000=(asthma[1:26,2]+asthma[1:26,3]))
erf <- 1.56
poparray <- DataframeToArray(pop, "Population")
dampxpop <- IntArray(dampness, poparray, "Population")
asthmaarray <- DataframeToArray(asthma, "Prevalence")
dampxpopxasthma <- IntArray(dampxpop, asthmaarray, "Prevalence")
final <- data.frame(dampxpopxasthma[,c(2,3,4,6,7,8)], Result=(dampxpopxasthma[,"Prevalence"] / 100 * dampxpopxasthma[,"Population"] / (erf / (100 / 
dampxpopxasthma[,"Result"] - 1) + 1) * (erf - 1) / (100 / dampxpopxasthma[,"Result"] - 1)))

Result

Show results


See also

Keywords

Asthma, indoor air, dampness, Europe

References

Related files

<mfanonymousfilelist></mfanonymousfilelist>