Ero sivun ”Radonin terveysvaikutukset” versioiden välillä

Opasnet Suomista
Siirry navigaatioon Siirry hakuun
(→‎Kaava: toteutettu jakaumat fraktiileista)
(→‎Kaava: melkein jo toimii)
Rivi 11: Rivi 11:
=== Kaava ===
=== Kaava ===


<rcode variables="name:n|default:100" graphics="1">
<rcode variables="name:n|default:10" graphics="1">
library(sorvi)
library(sorvi)
library(OpasnetBaseUtils)
library(OpasnetBaseUtils)
library(xtable)
library(xtable)
kunnat <- data.frame(kunta = c("Akaa", "Alajärvi", "Alavieska", "Alavus", "Asikkala",
library(ggplot2)
kunnat <- data.frame(Kunta = c("Akaa", "Alajärvi", "Alavieska", "Alavus", "Asikkala",
"Askola", "Aura", "Brändö", "Eckerö", "Enonkoski", "Enontekiö", "Espoo", "Eura",
"Askola", "Aura", "Brändö", "Eckerö", "Enonkoski", "Enontekiö", "Espoo", "Eura",
"Eurajoki", "Evijärvi", "Finström", "Forssa", "Föglö", "Geta", "Haapajärvi", "Haapavesi",
"Eurajoki", "Evijärvi", "Finström", "Forssa", "Föglö", "Geta", "Haapajärvi", "Haapavesi",
Rivi 63: Rivi 64:
"Äänekoski"),
"Äänekoski"),


maakunta = c("Pirkanmaa", "Etelä-Pohjanmaa", "Pohjois-Pohjanmaa", "Etelä-Pohjanmaa",
Maakunta = c("Pirkanmaa", "Etelä-Pohjanmaa", "Pohjois-Pohjanmaa", "Etelä-Pohjanmaa",
"Päijät-Häme", "Uusimaa", "Varsinais-Suomi", "Ahvenanmaa", "Ahvenanmaa", "Etelä-Savo",
"Päijät-Häme", "Uusimaa", "Varsinais-Suomi", "Ahvenanmaa", "Ahvenanmaa", "Etelä-Savo",
"Lappi", "Uusimaa", "Satakunta", "Satakunta", "Etelä-Pohjanmaa", "Ahvenanmaa",
"Lappi", "Uusimaa", "Satakunta", "Satakunta", "Etelä-Pohjanmaa", "Ahvenanmaa",
Rivi 196: Rivi 197:
out <- rbind(out, radon2)
out <- rbind(out, radon2)
}
}
radon2
#radon2
head(out)
#head(out)
 
talot <- op_baseGetData("opasnet_base", "Op_fi2761")[, -c(1, 2, 8)]
talot <- talot[talot$Talotyyppi != "Kaikki talotyypit" & talot$Kunta !="KOKO MAA - HELA LANDET" & talot$Asuntokunta != "Kaikki asuntokunnat", ]
talot$Talotyyppi <- ifelse(talot$Talotyyppi=="Asuinkerrostalo", "Kerrostalo", "Pientalo")
#head(talot)
väki <- data.frame(Asuntokunta = c("1 henk.", "2 henk.", "3 henk.", "4 henk.", "5 henk.", "6 henk.", "7+ henk."), n = 1:7)
talot <- merge(talot, väki)
#head(talot)
talot$Result <- talot$Result * talot$n
talot <- as.data.frame(as.table(tapply(talot$Result, list(talot$Kunta, talot$Talotyyppi), sum)))
colnames(talot) <- c("Kunta", "Talo", "Asukkaita")
#head(talot)
talot <- merge(talot, kunnat)
#head(talot)
out <- merge(out, talot)
#head(out)
 
rr.radon <- 1.0016
rr.tupakka <- 20
tausta.altistus <- 5
tausta.sairastuvuus <- data.frame( Maakunta = rep(c("Uusimaa", "Itä-Uusimaa", "Varsinais-Suomi", "Satakunta", "Kanta-Häme",
"Pirkanmaa", "Päijät-Häme", "Kymenlaakso", "Etelä-Karjala", "Etelä-Savo", "Itä-Savo", "Pohjois-Karjala", "Pohjois-Savo",
"Keski-Suomi", "Pohjanmaa", "Vaasa", "Keski-Pohjanmaa", "Pohjois-Pohjanmaa", "Kainuu", "Länsi-Pohja", "Lappi", "Ahvenanmaa", "Koko maa"), 2),
Sukupuoli = (rep( c("Miehet", "Naiset"), each = 23)), Sairastuvuus = c(29.5, 29.5, 26.3, 30.2, 30.3, 26.5, 26.9, 26.6, 29.0, 26.1, 28.7, 29.0, 26.7, 24.7, 26.1, 26.1, 29.8, 30.9, 28.3, 33.7, 34.3, 27.9, 28.2, 11.8, 11.8, 9.2, 7.7, 9.0, 7.7, 7.7, 9.5, 7.2, 7.4, 6.6, 5.8, 5.4, 6.0, 5.9, 8.7, 9.3, 8.7, 7.9, 10.1, 11.3, 13.9, 8.8))
tupakoivat <- data.frame(Sukupuoli = rep(c("Miehet", "Naiset"), 2), Tupakoivien.osuus = rep(c(0.22, 0.16), 2), Tupakka = rep(c("Kyllä", "Ei"), each = 2))
out <- merge(out, tausta.sairastuvuus)
out <- merge(out, tupakoivat)
# for(i in 1:nrow(out)){out$Tupakka[i] <- rbinom(1, 1, out$Osuus[i])} JÄRKYTTÄVÄN HIDAS
#head(out)
 
out$Taustariski <- out$Sairastuvuus * (1 - out$Tupakoivien.osuus * (rr.tupakka - 1) / (out$Tupakoivien.osuus * (rr.tupakka - 1) +1))
#head(out$Sairastuvuus)
#head((1 - out$Tupakoivien.osuus * (rr.tupakka - 1) / (out$Tupakoivien.osuus * (rr.tupakka - 1) +1)))
#head(out$Tupakoivien.osuus)
#head(rr.tupakka - 1)
 
out$Kokonaisvaikutus <- out$Taustariski/100000 * out$Asukkaita/2 * ifelse(out$Tupakka == "Kyllä", out$Tupakoivien.osuus, 1-out$Tupakoivien.osuus) # insidenssikerroin ja sukupuolen huomioonotto
out$Radonvaikutus  <- out$Kokonaisvaikutus * (exp(log(rr.tupakka)*ifelse(out$Tupakka == "Kyllä", 1,0) + log(rr.radon)*out$Radon) - exp(log(rr.tupakka)*ifelse(out$Tupakka == "Kyllä", 1,0)))
out$Kokonaisvaikutus <- out$Kokonaisvaikutus * exp(log(rr.tupakka)*ifelse(out$Tupakka == "Kyllä", 1,0) + log(rr.radon)*out$Radon)
#head(out)
print(xtable(out[1:60, ]), type = 'html')
out2 <- as.data.frame(as.table(((tapply(out$Kokonaisvaikutus, list(out$Maakunta, out$Talo, out$Tupakka), sum)))))
print(xtable(out2), type = 'html')
out2 <- as.data.frame(as.table(((tapply(out$Radonvaikutus, list(out$Maakunta, out$Talo, out$Tupakka), sum)))))
print(xtable(out2), type = 'html')
out2 <- as.data.frame(as.table(((tapply(out$Radonvaikutus, list(out$Talo, out$Tupakka), sum)))))
print(xtable(out2), type = 'html')
ggplot(out, aes(x = Maakunta, weight = Kokonaisvaikutus, fill = Tupakka)) + geom_bar(position = "stack") # position = "dodge"
temp <- as.data.frame(as.table(((tapply(out$Radonvaikutus, list(out$obs, out$Radon, out$Radonvaikutus), sum)))))
ggplot(out, aes(x = Radonvaikutus)) + geom_density(alpha = 0.2) # , fill = Scenario # TÄSSÄ ON NYT YMMÄRRYKSEN PUUTE
 
</rcode>
</rcode>
* Ikävakioitu (maailman standardiväestöön) syöpäkuolleisuus 100 000 henkilövuotta kohden vuosina 2005-2009 primaaripaikoittain ja sairaanhoitopiireittäin. Keuhkosyöpä miehet [http://stats.cancerregistry.fi/stats/fin/vfin0013m0.html] naiset [http://stats.cancerregistry.fi/stats/fin/vfin0014m0.html]
* Päivittäin tupakoivien osuus Suomessa: Tupakoivat [http://www.stat.fi/til/tup/2009/tup_2009_2010-12-22_tie_001.html]
* Tupakan annosvaste: UK: päivittäin tupakoivia 22 % Miehet, 20 % Naiset. Tämä aiheuttaa 88 % ja 84 % keuhkosyövistä vastaavasti. [http://info.cancerresearchuk.org/cancerstats/causes/lifestyle/tobacco/] Tämän perusteella voidaan laskea riskisuhde RR = AF / (EF(RR-1)+1), missä AF on attributable fraction eli altisteen aiheuttama osuus koko tautikuormasta ja EF on altistuneiden osuus koko väestöstä. Tämän perusteella päivittäisen tupakoinnin riskisuhteeksi keuhkosyövälle saadaan 20 - 30. (Laskennassa käytetään 20:tä.)


==Katso myös==
==Katso myös==

Versio 7. joulukuuta 2011 kello 15.32




Kysymys

Mitkä ovat radonin terveysvaikutukset Suomessa?

Vastaus

Perustelut

Kaava

n:

+ Näytä koodi

  • Ikävakioitu (maailman standardiväestöön) syöpäkuolleisuus 100 000 henkilövuotta kohden vuosina 2005-2009 primaaripaikoittain ja sairaanhoitopiireittäin. Keuhkosyöpä miehet [1] naiset [2]
  • Päivittäin tupakoivien osuus Suomessa: Tupakoivat [3]
  • Tupakan annosvaste: UK: päivittäin tupakoivia 22 % Miehet, 20 % Naiset. Tämä aiheuttaa 88 % ja 84 % keuhkosyövistä vastaavasti. [4] Tämän perusteella voidaan laskea riskisuhde RR = AF / (EF(RR-1)+1), missä AF on attributable fraction eli altisteen aiheuttama osuus koko tautikuormasta ja EF on altistuneiden osuus koko väestöstä. Tämän perusteella päivittäisen tupakoinnin riskisuhteeksi keuhkosyövälle saadaan 20 - 30. (Laskennassa käytetään 20:tä.)

Katso myös

Viitteet


Aiheeseen liittyviä tiedostoja

<mfanonymousfilelist></mfanonymousfilelist>