Het ggd tool (url=ggd-rekentool-luchtkwaliteit-en-gezondheid - Academische Werkplaats Gezonde Leefomgeving ) is een spreadsheet model waarin je kan zien wat voor invloed gemiddelde fijnstof en stikstofdioxide jaarconcentraties op de gezondheid van een nederlandse populatie hebben.
Op bovenstaande website staan zowel de spreadsheet als enige toelichtende artikelen.
Je kunt het tool als een black box gebruiken, en zo zien wat een langdurige blootstelling aan en gemiddelde concentratie betekent.
Maar als je de “hidden tabs” van de spreadsheet zichtbaar maakt, dan kun je het Rekenmodel bekijken en de formules reversed engineeren.Het helpt daarbij als je de artikelen op die site leest.
In dit eerste artikeltje proberen we de onderliggende methode kwantitatief te doorgronden.
We nemen als voorbeeld het effect dat fijnstof heeft op het jaarlijks aantal beroertes in NL voor 35+
In het model is het aantal beroertes alleen afhankelijk van pm25
RelativeRisk (RR) geeft de sterkte van dat effect.
De eerste stap is dat we met een python programmatje hetzelfde resultaat hopen te halen als de spreadsheet (heel nederland, pm25=10 ug/m3, 35+ beroertes = 7409)
Dat is gelukt met deze code :
De output is:
Voor andere jaarlijkse gezondheidsproblemen, zoals vermeld in het ggdtool, moeten we andere constantes nemen.Maar het principe werkt !
In volgende arikelen uit deze reeks gaan we in op dirty details, onzekerheden , leeftijdsopbouw afhankelijkheden en evt vragen op het forum nav dit stukje