Ero sivun ”Radonin terveysvaikutukset” versioiden välillä

Opasnet Suomista
Siirry navigaatioon Siirry hakuun
Ei muokkausyhteenvetoa
Rivi 124: Rivi 124:
radon <- radon[, 1:3]
radon <- radon[, 1:3]
colnames(radon)[3] <- "r"
colnames(radon)[3] <- "r"
head(radon)
#head(radon)
radon <- radon[substr(radon$Havainto, 1, 7) != "kaikki.", ]
radon <- radon[substr(radon$Havainto, 1, 7) != "kaikki.", ]
radon <- radon[radon$Maakunta != "Yhteensä", ]
radon <- radon[radon$Maakunta != "Yhteensä", ]
Rivi 212: Rivi 212:
out <- merge(out, tausta.sairastuvuus)
out <- merge(out, tausta.sairastuvuus)
out <- merge(out, tupakoivat)
out <- merge(out, tupakoivat)
# for(i in 1:nrow(out)){out$Tupakka[i] <- rbinom(1, 1, out$Osuus[i])} JÄRKYTTÄVÄN HIDAS
out$Asukkaita <- out$Asukkaita/2 * ifelse(out$Tupakka == "Kyllä", out$Tupakoivien.osuus, 1-out$Tupakoivien.osuus)
#head(out)
 
#out[is.na(out$Radon), ]
out$Sairastuvuus <- out$Sairastuvuus * (1 - out$Tupakoivien.osuus * (rr.tupakka - 1) / (out$Tupakoivien.osuus * (rr.tupakka - 1) +1))/100000
#out
 
out$Taustariski <- out$Sairastuvuus * (1 - out$Tupakoivien.osuus * (rr.tupakka - 1) / (out$Tupakoivien.osuus * (rr.tupakka - 1) +1))
out$Kokonaisvaikutus <- exp(log(rr.tupakka)*ifelse(out$Tupakka == "Kyllä", 1,0) + log(rr.radon)*out$Radon)
#head(out$Sairastuvuus)
out$Radonvaikutus  <- out$Kokonaisvaikutus - exp(log(rr.tupakka)*ifelse(out$Tupakka == "Kyllä", 1,0))
#head((1 - out$Tupakoivien.osuus * (rr.tupakka - 1) / (out$Tupakoivien.osuus * (rr.tupakka - 1) +1)))
 
#head(out$Tupakoivien.osuus)
out$Radonvaikutus <- out$Radonvaikutus * out$Sairastuvuus * out$Asukkaita
#head(rr.tupakka - 1)
out$Kokonaisvaikutus <- out$Kokonaisvaikutus * out$Sairastuvuus * out$Asukkaita
 
print(xtable(out[1:20, ]), type = 'html')
 
print("Yksikköriski")
print(xtable(as.data.frame(as.table(tapply(out$Radonvaikutus/out$Asukkaita, list(out$Maakunta, out$Tupakka, out$Talo), sum)))), type= 'html')


print("Radonin aiheuttama lisäriski")
print(xtable(as.data.frame(as.table(((tapply(out$Radonvaikutus, list(out$Talo, out$Tupakka), sum)))))), type = 'html')


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"
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)))))
 
out$Radonvaikutus <- ifelse(out$Radonvaikutus < 0, 0, out$Radonvaikutus)
out$Radonvaikutus <- ifelse(out$Radonvaikutus < 0, 0, out$Radonvaikutus)
#head(out)
#ggplot(out, aes(x = Radon, weight = Radonvaikutus)) + geom_density(alpha = 0.2) # , fill = Tupakka
#barplot(out$Radon)
temp <- out[, c("Radon", "Radonvaikutus", "Tupakka")]


temp$Radon <- round(as.numeric(as.character(temp$Radon)), -2)
out$Radonsata <- round(as.numeric(as.character(out$Radon)), -2)
head(temp)
 
#temp <- as.data.frame(as.table(((tapply(temp$Radonvaikutus, list(temp$Radon, temp$Tupakka), sum)))))
ggplot(out, aes(x = Radon, weight = Radonvaikutus, fill = Tupakka)) + geom_density(alpha = 1)
#head(temp)
ggplot(out, aes(x = Radonsata, weight = Radonvaikutus, fill = Tupakka)) + geom_bar(position = "dodge")
#temp$Var1
ggplot(temp, aes(x = Radon, weight = Radonvaikutus, fill = Tupakka)) + geom_density(alpha = 10)
</rcode>
</rcode>


Versio 7. joulukuuta 2011 kello 20.43




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>

Noudettu kohteesta ”https://fi.opasnet.org/index.php?title=Radonin_terveysvaikutukset&oldid=10277”