Difference between revisions of "Assessment of building policies' effect on dampness and asthma in Europe"

From Testiwiki
Jump to: navigation, search
m (R code for detailed analysis)
m (R code for detailed analysis)
Line 38: Line 38:
 
final <- apply(array, c(6,4,8,7,3), gmean)
 
final <- apply(array, c(6,4,8,7,3), gmean)
 
final <- apply(final, c(1,2,3,4), sum, na.rm=TRUE)
 
final <- apply(final, c(1,2,3,4), sum, na.rm=TRUE)
final <- final[,,"0",]-final[,,"1",]
+
final <- final[,,"1",]-final[,,"0",]
 
finaldf <- as.data.frame(as.table(final))
 
finaldf <- as.data.frame(as.table(final))
 
barplot(finaldf[,"Freq"])
 
barplot(finaldf[,"Freq"])

Revision as of 13:17, 16 December 2010



Scope

Purpose

What are the effects of different building policies on dampness and asthma prevalence in Europe?

Boundaries etc.

Boundaries, scenarios, intended users, and participants are the same as in the Mega case study.

Definition

Decision variables
Other variables
Indicators

R code for detailed analysis

BP <- op_baseGetData("opasnet_base", "Op_en4068")
asthma <- op_baseGetData("opasnet_base", "Op_en4723")
BParray <- DataframeToArray(BP)
asthmaxBP <- IntArray(asthma, BParray, "Successful")
array <- DataframeToArray(asthmaxBP)
gmean <- function(x) exp(mean(log(x[is.na(x)==FALSE])))
final <- apply(array, c(6,4,8,7,3), gmean)
final <- apply(final, c(1,2,3,4), sum, na.rm=TRUE)
final <- final[,,"1",]-final[,,"0",]
finaldf <- as.data.frame(as.table(final))
barplot(finaldf[,"Freq"])
#final <- data.frame(Policy=rep(levels(factor(asthmaxBP[,"policy"])), each=2), Successful=rep(c(1,0),4), 
#	Mean=c(mean(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[1]&asthmaxBP[,"Successful"]==1,"Result"]), 
#	mean(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[1]&asthmaxBP[,"Successful"]==0,"Result"]),
#	mean(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[2]&asthmaxBP[,"Successful"]==1,"Result"]),
#	mean(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[2]&asthmaxBP[,"Successful"]==0,"Result"]),
#	mean(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[3]&asthmaxBP[,"Successful"]==1,"Result"]), 
#	mean(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[3]&asthmaxBP[,"Successful"]==0,"Result"]),
#	mean(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[4]&asthmaxBP[,"Successful"]==1,"Result"]), 
#	mean(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[4]&asthmaxBP[,"Successful"]==0,"Result"])),
#	GMean=c(exp(mean(log(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[1]&asthmaxBP[,"Successful"]==1,"Result"]))), 
#	exp(mean(log(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[1]&asthmaxBP[,"Successful"]==0,"Result"]))),
#	exp(mean(log(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[2]&asthmaxBP[,"Successful"]==1,"Result"]))),
#	exp(mean(log(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[2]&asthmaxBP[,"Successful"]==0,"Result"]))),
#	exp(mean(log(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[3]&asthmaxBP[,"Successful"]==1,"Result"]))), 
#	exp(mean(log(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[3]&asthmaxBP[,"Successful"]==0,"Result"]))),
#	exp(mean(log(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[4]&asthmaxBP[,"Successful"]==1,"Result"]))), 
#	exp(mean(log(asthmaxBP[asthmaxBP[,"policy"]==levels(factor(asthmaxBP[,"policy"]))[4]&asthmaxBP[,"Successful"]==0,"Result"])))))
#barplot(final[,4])

Result

Asthma incidence increase due to building dampness in Europe: Show results

Extra asthma cases due to building dampness
Policy Successful Mean GMean
Behavioural changes 1 124627.8 42716.72
Behavioural changes 0 123267.9 42221.99
EPBD 2008 1 122072.8 42204.77
EPBD 2008 0 125699.5 42716.64
More efficient urban structures 1 122790.4 42540.45
More efficient urban structures 0 125024.9 42401.92
More renewables 1 123290.7 42495.84
More renewables 0 124626.3 42440.56

See also

Keywords

Dampness, indoor air, asthma, Europe

References