Ero sivun ”Hiekkalaatikko” versioiden välillä

Opasnet Suomista
Siirry navigaatioon Siirry hakuun
Ei muokkausyhteenvetoa
Ei muokkausyhteenvetoa
Rivi 1: Rivi 1:
<p>Päivitetty
Päivitetty
</p><p><span class="texhtml">123 + &#945; &#8722; 555 &#8722; 666</span>
 
</p><p><br />
<math>123+\alpha -555-666</math>
</p>
 
<h3> Terveysvaikutukset </h3>
[[op_en:Sandbox]]
<p><b>HUOM! Tällä hetkellä malli laskee ainoastaan murskausprosessin pölypäästöt ja niiden leviämisen 10-15 km:n säteellä olevaan väestöön. Mallia päivitetään jatkuvasti.</b>
 
</p><p>&lt;img src="/fi_wiki/images/c/cc/Harjoituskuva.png" _fck_mw_filename="Harjoituskuva.png" _fck_mw_width="200" _fck_mw_type="frame" alt="" class="fck_mw_frame fck_mw_right" /&gt;
 
</p><p>&lt;img src="/fi_wiki/images/0/06/Koala.jpg" _fck_mw_filename="Koala.jpg" _fck_mw_width="300" _fck_mw_height="600" _fck_mw_type="frame" alt="" class="fck_mw_frame fck_mw_right" /&gt;
=== Terveysvaikutukset ===
</p><p><span class="fck_mw_special">fckLRlibrary(OpasnetBaseUtils)fckLRlibrary(ggplot2)fckLRfckLRearth.radius &lt;- 6372.8 # quadratic mean or root mean square approximation of the average great-circle fckLR # circumference derives a radius of about 6372.8 km (Wikipedia)fckLRcentral.angle &lt;- function(s.la, s.lo, f.la, f.lo) 2 * asin((sin((s.la - f.la) / 2)^2 + cos(s.la) * cos(f.la) * sin((s.lo - f.lo) / 2)^2)^0.5)fckLRfckLRdtheta.y &lt;- 1/earth.radius*180/pi # central angle increase per 1 kilometer north from a given point assuming no displacement on x axisfckLRdtheta.x &lt;- 2*asin(sin(1/(2*earth.radius))/cos(LA/180*pi))*180/pi # central angle increase per 1 kilometer east from a given point fckLR# - assuming no displacement on y axisfckLRfckLR# Populaatio datafckLRfckLRpop.locs &lt;- op_baseGetLocs(&quot;heande_base&quot;, &quot;Heande3182&quot;, apply.utf8 = FALSE)fckLRhead(pop.locs)fckLRpop.slice.la &lt;- pop.locs[pop.locs$ind == &quot;Latitude&quot;, &quot;loc_id&quot;][pop.locs[pop.locs$ind == &quot;Latitude&quot;, &quot;loc&quot;] &lt; LA + 10.5 * dtheta.y &amp; fckLR pop.locs[pop.locs$ind == &quot;Latitude&quot;, &quot;loc&quot;] &gt; LA - 10.5 * dtheta.y]fckLRpop.slice.lo.inverse &lt;- pop.locs[pop.locs$ind == &quot;Longitude&quot;, &quot;loc_id&quot;][pop.locs[pop.locs$ind == &quot;Longitude&quot;, &quot;loc&quot;] &gt; LO + 10.5 * dtheta.x | fckLR pop.locs[pop.locs$ind == &quot;Longitude&quot;, &quot;loc&quot;] &lt; LO - 10.5 * dtheta.x]fckLRfckLRpop &lt;- op_baseGetData(&quot;heande_base&quot;, &quot;Heande3182&quot;, include = pop.slice.la, exclude = pop.slice.lo.inverse)fckLRfckLRhead(pop)fckLRfckLRpop$Longitude &lt;- as.numeric(as.character(pop$Longitude))fckLRpop$Latitude &lt;- as.numeric(as.character(pop$Latitude))fckLRfckLRpop$LObin &lt;- cut(pop$Longitude, breaks = LO + ((-11:10) + 0.5 ) * dtheta.x)fckLRpop$LAbin &lt;- cut(pop$Latitude, breaks = LA + ((-11:10) + 0.5 ) * dtheta.y)fckLRfckLR# Pitoisuus datafckLRfckLRpitoisuus &lt;- function(n, paasto, L.matrix) { #, X_coord, Y_coord) {fckLR ID.list &lt;- tapply(1:nrow(L.matrix), L.matrix[,c(&quot;Kaupunki&quot;, &quot;Vuosi&quot;, &quot;Tyyppi&quot;)], list)fckLR ID.list.samples &lt;- sample(ID.list, n, replace = TRUE)fckLR ID.vec &lt;- unlist(ID.list.samples)fckLRprint(ID.vec)fckLRprint(ID.list)fckLR #c.matrix &lt;- pitoisuus(L.matrix, Paasto) #, X_coord, Y_coord)fckLR l.matrix &lt;- L.matrix[ID.vec,]fckLR l.matrix$obs &lt;- rep(1:n, each = length(ID.vec)/n)fckLRprint(head(l.matrix)) fckLR c.matrix &lt;- merge(l.matrix, paasto)fckLR c.matrix &lt;- model.frame(I(Paasto * k) ~., data = c.matrix)fckLR colnames(c.matrix)[1] &lt;- &quot;Pitoisuus&quot;fckLRprint(head(c.matrix))fckLR return(c.matrix)fckLR}fckLRfckLRPILTTI.matrix &lt;- op_baseGetData(&quot;heande_base&quot;, &quot;Heande3181&quot;)[,-c(1,2,9)] # unit: ugm^-3/Mga^-1fckLRfckLRPILTTI.matrix$dy &lt;- as.numeric(as.character(PILTTI.matrix$dy))fckLRPILTTI.matrix$dx &lt;- as.numeric(as.character(PILTTI.matrix$dx))fckLRfckLRcolnames(PILTTI.matrix)[colnames(PILTTI.matrix)==&quot;Result&quot;] &lt;- &quot;k&quot;fckLRfckLRhead(PILTTI.matrix)fckLRfckLRPaasto &lt;- data.frame(Paasto = murskaus.maara * paasto.kerroin / 1e6) # unit: Mga^-1fckLRfckLRhead(Paasto)fckLRfckLRC.matrix &lt;- pitoisuus(N, Paasto, PILTTI.matrix) # unit: ugm^-3fckLRNfckLRhead(C.matrix)fckLRfckLRC.matrix$LObin &lt;- cut(C.matrix$dx / 1000 * dtheta.x + LO, breaks = LO + ((-11:10) + 0.5 ) * dtheta.x)fckLRC.matrix$LAbin &lt;- cut(C.matrix$dy / 1000 * dtheta.y + LA, breaks = LA + ((-11:10) + 0.5 ) * dtheta.y)fckLRfckLR# YhdistysfckLRfckLRhead(C.matrix)fckLRfckLRpop.paasto &lt;- merge(C.matrix, pop[,-c(1,2,6)])fckLRfckLRpop.paasto.korjaus &lt;- data.frame(Pitoisuus = pop.paasto$Pitoisuus, Vaesto = pop.paasto$Result / N)fckLRfckLRhead(pop.paasto.korjaus)fckLRfckLRplot1 &lt;- ggplot(pop.paasto.korjaus, aes(x = Pitoisuus, weight = Vaesto)) + geom_histogram(binwidth = fckLR(max(pop.paasto.korjaus$Pitoisuus[pop.paasto.korjaus$Vaesto!=0]) - fckLRmin(pop.paasto.korjaus$Pitoisuus[pop.paasto.korjaus$Vaesto!=0]))/100) + fckLRxlim(min(pop.paasto.korjaus$Pitoisuus[pop.paasto.korjaus$Vaesto!=0]) - 0.000001, fckLRmax(pop.paasto.korjaus$Pitoisuus[pop.paasto.korjaus$Vaesto!=0]) + 0.000001)fckLRfckLRplot1 + geom_vline(xintercept = 40, colour = &quot;red&quot;) # OhjearvofckLRfckLRexposure.pop &lt;- tapply(pop.paasto$Pitoisuus * pop.paasto$Result, pop.paasto[,c(&quot;obs&quot;)], sum)fckLRfckLRbg.mort &lt;- 45182/5203826 # same values as used in PILTTIfckLRfckLRerf &lt;- 0.0097 # J. T. Tuomisto, A. Wilson, et al. Uncertainty in mortality response to airborne fine particulate matter... 2008fckLR# unit: 1/ugm^-3fckLRfckLRmort.out &lt;- erf * bg.mort * exposure.popfckLRfckLRqplot(mort.out, geom=&quot;density&quot;)fckLRfckLRcat(&quot;Odotusarvo kuolemille vuodessa:&quot;, mean(mort.out), &quot;\n&quot;)fckLRfckLRcat(&quot;Ohjearvon 40 ugm^-3 mukaisen altistusrajan ylitti&quot;, sum(pop.paasto.korjaus$Vaesto[pop.paasto.korjaus$Pitoisuus&gt;40]), &quot;asukasta.\n&quot;)fckLR</span>
 
</p>
'''HUOM! Tällä hetkellä malli laskee ainoastaan murskausprosessin pölypäästöt ja niiden leviämisen 10-15 km:n säteellä olevaan väestöön. Mallia päivitetään jatkuvasti.
<ul><li>Oletusarvot Siilinjärven Yaralle, yhdelle tonnille murskattua malmia vuodessa olettaen että murskaustapa on kuivaseulonta ilman avulla (suurin päästökerroin). Katso: &lt;a _fcknotitle="true" href=":op fi:Metallimalmin murskausprosessin kokonais- (TSP) ja hengittyvien hiukkasten (PM10) päästökertoimet"&gt;:op fi:Metallimalmin murskausprosessin kokonais- (TSP) ja hengittyvien hiukkasten (PM10) päästökertoimet&lt;/a&gt;.
 
</li><li>Havainnollisia lähtöarvojen muutoksia:
[[Image:Harjoituskuva.png|600px|Hassu kuva]]
<ul><li> Luikonlahden koordinaatit (pienempi väestö kaivoksen ympärillä): latitude=62.936836, longitude=28.70749  
 
</li><li> Hienomurskaus kostealle malmille: päästökerroin 9.1 (oletus: kuiva malmi 72.6)
[[Image:Harjoituskuva.png|200px|Hassu kuva]]
</li></ul>
 
</li></ul>
[[Image:Harjoituskuva.png|600x200px|Hassumpi kuva]]
<h2> R-testi </h2>
 
<p><span class="fck_mw_special">fckLR######## Haetaan R-koodi generic sivulta Projektinhallinta. Sisältää funktiot dropall ja PTable.fckLR######## Ladataan tarvittavat paketitfckLRlibrary(OpasnetBaseUtils)fckLRlibrary(ggplot2)fckLRlibrary(xtable)fckLRfckLRprint(&quot;Haetaan tarvittava data Opasnet-kannasta&quot;)fckLRsaanto.siemenet &lt;- op_baseGetData(&quot;opasnet_base&quot;, &quot;Op_fi2633&quot;)[,-c(1,2,7)] # Jatropan siementen saanto viljelystäfckLRsaanto.öljy    &lt;- op_baseGetData(&quot;opasnet_base&quot;, &quot;Op_fi2634&quot;)[,-c(1,2,5)] # Öljyn saanto jatropan siemenistäfckLRsaanto.diesel  &lt;- op_baseGetData(&quot;opasnet_base&quot;, &quot;Op_fi2632&quot;)[,-c(1,2,5)] # Biodieselin saanto jatropaöljystäfckLRviljelyala      &lt;- op_baseGetData(&quot;opasnet_base&quot;, &quot;Op_fi2642&quot;)[,-c(1,2)] # Jatropan viljelyalueetfckLRpäästö.ilmasto  &lt;- op_baseGetData(&quot;opasnet_base&quot;, &quot;Op_fi2547&quot;)[,-c(1,2)] # Jatropan viljelyn ilmastovaikutuksetfckLRpäästö.sosiaali &lt;- op_baseGetData(&quot;opasnet_base&quot;, &quot;Op_fi2552&quot;)[,-c(1,2)] # Jatropan viljelyn sosiaaliset vaikutuksetfckLRpäästö.ekosyst  &lt;- op_baseGetData(&quot;opasnet_base&quot;, &quot;Op_fi2548&quot;)[,-c(1,2)] # Jatropan viljelyn ekosysteemivaikutuksetfckLRP              &lt;- op_baseGetData(&quot;opasnet_base&quot;, &quot;Op_fi2539&quot;)[,-c(1,2,7)] # Jatropan käyttö bioenergian lähteenäfckLRfckLRprint(&quot;Ajetaan malli&quot;)fckLRfckLR### Muutetaan sarakkeiden nimet sopiviksi yhdistämistä varten.fckLRcolnames(saanto.siemenet)[4] &lt;- &quot;siemenet&quot;fckLRcolnames(saanto.öljy)[2] &lt;- &quot;öljy&quot;fckLRcolnames(saanto.diesel)[2] &lt;- &quot;diesel&quot;fckLRfckLR## Yhdistetään tiedot toisiinsa yhdeksi data.frameksi. Lasketaan saanto.fckLRsaanto &lt;- merge(saanto.siemenet, saanto.öljy)fckLRsaanto &lt;- merge(saanto, saanto.diesel)fckLRsaanto[,9] &lt;- saanto$siemenet * saanto$öljy * saanto$diesel * alafckLRcolnames(saanto)[9] &lt;- &quot;saanto (kg/a)&quot;fckLRfckLR## Muutetaan tulos jakaumaksi Monte Carlolla.fckLRP &lt;- PTable(P, n)fckLRsaanto &lt;- merge(P, saanto)fckLRfckLR## Lasketaan tulostaulu.fckLRif(length(divisions)&gt;1) divisions &lt;- as.list(saanto[, divisions]) else divisions &lt;- saanto[, divisions]fckLRout1 &lt;- as.data.frame(as.table(tapply(saanto[, 10], divisions, mean))) fckLRout1 &lt;- dropall(out1[!is.na(out1$Freq), ])fckLRfckLRprint(xtable(out1), type = 'html')fckLRfckLR## Lasketaan tuloskuvaaja.fckLRout2 &lt;- as.data.frame(as.table(tapply(saanto[, 10], list(saanto[, divisions2], saanto$obs), mean))) fckLRout2 &lt;- dropall(out2[!is.na(out2$Freq), ])fckLRggplot(out2, aes(x = Freq, weight = 1, fill = Var1)) +geom_density() fckLR## Jostain syystä vain osa kuvista piirtyy oikein, riippuen mitä parametreja valitaan. En ymmärrä syytä.fckLR</span>
[[Tiedosto:koala.jpg|300x600px|Hassuin kuva]]
</p><p><br />
 
<p><span class="fck_mw_template"><span class="fck_mw_template">{{bluebox|fckLR==Komia loota==fckLR* sininenfckLRfckLRon komia...fckLR}}</span></span>
<rcode variables = "LA|Leveyspiiri|63.110577|LO|Pituuspiiri|27.735929|murskaus.maara|Murskattavan malmin määrä (Mg/a)|10000000|paasto.kerroin|Murskaukseen käytettävän prosessin päästökerroin (g/Mg)|72.6|N|Iteraatioiden määrä|100" graphics = "1">
</p>
library(OpasnetBaseUtils)
</p><p><span class="fck_mw_special">fckLR2001|Finland|1001fckLR2002|Finland|judanssfckLR2001|Sweden|900fckLR2002|Sweden|800fckLR</span>
library(ggplot2)
</p><p><span class="fck_mw_template"><span class="fck_mw_template">{{greenbox|fckLRJEPJEJPKEPKEfckLRfckLRJOOfckLR}}</span></span>
 
</p><p>fwae
earth.radius <- 6372.8 # quadratic mean or root mean square approximation of the average great-circle  
</p><p>JEPULISTA
# circumference derives a radius of about 6372.8 km (Wikipedia)
<p><span class="fck_mw_special">fckLRcategory = THL:n julkaisuja vuonna 2010fckLRmode = orderedfckLRinclude = {julkaisu}fckLRordermethod=titlefckLRformat  = ,\n&lt;big&gt;<b><a href="%PAGE%">%TITLE%</a></b>&lt;/big&gt;,,fckLRnotcategory = glossaryfckLRorder = ascendingfckLR</span>
central.angle <- function(s.la, s.lo, f.la, f.lo) 2 * asin((sin((s.la - f.la) / 2)^2 + cos(s.la) * cos(f.la) * sin((s.lo - f.lo) / 2)^2)^0.5)
</p>
 
</p><p><br />
dtheta.y <- 1/earth.radius*180/pi # central angle increase per 1 kilometer north from a given point assuming no displacement on x axis
<p><span class="fck_mw_template"><span class="fck_mw_template">{{budjettirajapintafckLR| työpaketit        = 'TP1';TP1;'TP2';TP2;'TP3';TP3;'TP4';TP4;'TP5';TP5fckLR| vuodet            = '2012';2012;'2013';2013;'2014;2014fckLR| sivu              = Op_fi2682fckLR| rahoittaja        = STMfckLR| projekti          = TekaisufckLR| työpaketti        = TP1fckLR| työpaketti.add    = TP1fckLR| rahoittaja.add    = STMfckLR| työpaketit_default = 'TP1';'TP2';'TP3';'TP5'fckLR| rahoittaja_default = STMfckLR}}</span></span>
dtheta.x <- 2*asin(sin(1/(2*earth.radius))/cos(LA/180*pi))*180/pi # central angle increase per 1 kilometer east from a given point  
</p>
# - assuming no displacement on y axis
</p>&lt;a _fcknotitle="true" href="Category:Testi"&gt;Testi&lt;/a&gt;
 
<p>&lt;a href="op en:Sandbox"&gt;op en:Sandbox&lt;/a&gt;
# Populaatio data
</p>
 
pop.locs <- op_baseGetLocs("heande_base", "Heande3182", apply.utf8 = FALSE)
head(pop.locs)
pop.slice.la <- pop.locs[pop.locs$ind == "Latitude", "loc_id"][pop.locs[pop.locs$ind == "Latitude", "loc"] < LA + 10.5 * dtheta.y &  
pop.locs[pop.locs$ind == "Latitude", "loc"] > LA - 10.5 * dtheta.y]
pop.slice.lo.inverse <- pop.locs[pop.locs$ind == "Longitude", "loc_id"][pop.locs[pop.locs$ind == "Longitude", "loc"] > LO + 10.5 * dtheta.x |  
pop.locs[pop.locs$ind == "Longitude", "loc"] < LO - 10.5 * dtheta.x]
 
pop <- op_baseGetData("heande_base", "Heande3182", include = pop.slice.la, exclude = pop.slice.lo.inverse)
 
head(pop)
 
pop$Longitude <- as.numeric(as.character(pop$Longitude))
pop$Latitude <- as.numeric(as.character(pop$Latitude))
 
pop$LObin <- cut(pop$Longitude, breaks = LO + ((-11:10) + 0.5 ) * dtheta.x)
pop$LAbin <- cut(pop$Latitude, breaks = LA + ((-11:10) + 0.5 ) * dtheta.y)
 
# Pitoisuus data
 
pitoisuus <- function(n, paasto, L.matrix) { #, X_coord, Y_coord) {
ID.list <- tapply(1:nrow(L.matrix), L.matrix[,c("Kaupunki", "Vuosi", "Tyyppi")], list)
ID.list.samples <- sample(ID.list, n, replace = TRUE)
ID.vec <- unlist(ID.list.samples)
print(ID.vec)
print(ID.list)
#c.matrix <- pitoisuus(L.matrix, Paasto) #, X_coord, Y_coord)
l.matrix <- L.matrix[ID.vec,]
l.matrix$obs <- rep(1:n, each = length(ID.vec)/n)
print(head(l.matrix))
c.matrix <- merge(l.matrix, paasto)
c.matrix <- model.frame(I(Paasto * k) ~., data = c.matrix)
colnames(c.matrix)[1] <- "Pitoisuus"
print(head(c.matrix))
return(c.matrix)
}
 
PILTTI.matrix <- op_baseGetData("heande_base", "Heande3181")[,-c(1,2,9)] # unit: ugm^-3/Mga^-1
 
PILTTI.matrix$dy <- as.numeric(as.character(PILTTI.matrix$dy))
PILTTI.matrix$dx <- as.numeric(as.character(PILTTI.matrix$dx))
 
colnames(PILTTI.matrix)[colnames(PILTTI.matrix)=="Result"] <- "k"
 
head(PILTTI.matrix)
 
Paasto <- data.frame(Paasto = murskaus.maara * paasto.kerroin / 1e6) # unit: Mga^-1
 
head(Paasto)
 
C.matrix <- pitoisuus(N, Paasto, PILTTI.matrix) # unit: ugm^-3
N
head(C.matrix)
 
C.matrix$LObin <- cut(C.matrix$dx / 1000 * dtheta.x + LO, breaks = LO + ((-11:10) + 0.5 ) * dtheta.x)
C.matrix$LAbin <- cut(C.matrix$dy / 1000 * dtheta.y + LA, breaks = LA + ((-11:10) + 0.5 ) * dtheta.y)
 
# Yhdistys
 
head(C.matrix)
 
pop.paasto <- merge(C.matrix, pop[,-c(1,2,6)])
 
pop.paasto.korjaus <- data.frame(Pitoisuus = pop.paasto$Pitoisuus, Vaesto = pop.paasto$Result / N)
 
head(pop.paasto.korjaus)
 
plot1 <- ggplot(pop.paasto.korjaus, aes(x = Pitoisuus, weight = Vaesto)) + geom_histogram(binwidth =  
(max(pop.paasto.korjaus$Pitoisuus[pop.paasto.korjaus$Vaesto!=0]) -  
min(pop.paasto.korjaus$Pitoisuus[pop.paasto.korjaus$Vaesto!=0]))/100) +  
xlim(min(pop.paasto.korjaus$Pitoisuus[pop.paasto.korjaus$Vaesto!=0]) - 0.000001,  
max(pop.paasto.korjaus$Pitoisuus[pop.paasto.korjaus$Vaesto!=0]) + 0.000001)
 
plot1 + geom_vline(xintercept = 40, colour = "red") # Ohjearvo
 
exposure.pop <- tapply(pop.paasto$Pitoisuus * pop.paasto$Result, pop.paasto[,c("obs")], sum)
 
bg.mort <- 45182/5203826 # same values as used in PILTTI
 
erf <- 0.0097 # J. T. Tuomisto, A. Wilson, et al. Uncertainty in mortality response to airborne fine particulate matter... 2008
# unit: 1/ugm^-3
 
mort.out <- erf * bg.mort * exposure.pop
 
qplot(mort.out, geom="density")
 
cat("Odotusarvo kuolemille vuodessa:", mean(mort.out), "\n")
 
cat("Ohjearvon 40 ugm^-3 mukaisen altistusrajan ylitti", sum(pop.paasto.korjaus$Vaesto[pop.paasto.korjaus$Pitoisuus>40]), "asukasta.\n")
</rcode>
 
*Oletusarvot Siilinjärven Yaralle, yhdelle tonnille murskattua malmia vuodessa olettaen että murskaustapa on kuivaseulonta ilman avulla (suurin päästökerroin). Katso: [[:op_fi:Metallimalmin murskausprosessin kokonais- (TSP) ja hengittyvien hiukkasten (PM10) päästökertoimet]].
*Havainnollisia lähtöarvojen muutoksia:
** Luikonlahden koordinaatit (pienempi väestö kaivoksen ympärillä): latitude=62.936836, longitude=28.70749  
** Hienomurskaus kostealle malmille: päästökerroin 9.1 (oletus: kuiva malmi 72.6)
 
== R-testi ==
 
<rcode
graphics="1"  
include="page:Funktioita_R-toolsiin|name:generic"
variables="name:ala|default:900000|description:Jatropan viljelyala (ha)|
name:n|default:10|
name:divisions|description:Mitkä tekijät halua eritellä tuloksessa?|type:checkbox|options:'Katalyytin määrä';Katalyytin määrä;'Ikä';Ikä;'Kastelu';Kastelu;'Käytetty puristin';Käytetty puristin|default:'Kastelu'|
name:divisions2|description:Minkä yhden tekijän halua eritellä kuvaajassa?|type:selection|options:'Katalyytin määrä';Katalyytin määrä;'Ikä';Ikä;'Kastelu';Kastelu;'Käytetty puristin';Käytetty puristin|default:'Kastelu'
">
######## Haetaan R-koodi generic sivulta Projektinhallinta. Sisältää funktiot dropall ja PTable.
######## Ladataan tarvittavat paketit
library(OpasnetBaseUtils)
library(ggplot2)
library(xtable)
 
print("Haetaan tarvittava data Opasnet-kannasta")
saanto.siemenet <- op_baseGetData("opasnet_base", "Op_fi2633")[,-c(1,2,7)] # Jatropan siementen saanto viljelystä
saanto.öljy    <- op_baseGetData("opasnet_base", "Op_fi2634")[,-c(1,2,5)] # Öljyn saanto jatropan siemenistä
saanto.diesel  <- op_baseGetData("opasnet_base", "Op_fi2632")[,-c(1,2,5)] # Biodieselin saanto jatropaöljystä
viljelyala      <- op_baseGetData("opasnet_base", "Op_fi2642")[,-c(1,2)] # Jatropan viljelyalueet
päästö.ilmasto  <- op_baseGetData("opasnet_base", "Op_fi2547")[,-c(1,2)] # Jatropan viljelyn ilmastovaikutukset
päästö.sosiaali <- op_baseGetData("opasnet_base", "Op_fi2552")[,-c(1,2)] # Jatropan viljelyn sosiaaliset vaikutukset
päästö.ekosyst  <- op_baseGetData("opasnet_base", "Op_fi2548")[,-c(1,2)] # Jatropan viljelyn ekosysteemivaikutukset
P              <- op_baseGetData("opasnet_base", "Op_fi2539")[,-c(1,2,7)] # Jatropan käyttö bioenergian lähteenä
 
print("Ajetaan malli")
 
### Muutetaan sarakkeiden nimet sopiviksi yhdistämistä varten.
colnames(saanto.siemenet)[4] <- "siemenet"
colnames(saanto.öljy)[2] <- "öljy"
colnames(saanto.diesel)[2] <- "diesel"
 
## Yhdistetään tiedot toisiinsa yhdeksi data.frameksi. Lasketaan saanto.
saanto <- merge(saanto.siemenet, saanto.öljy)
saanto <- merge(saanto, saanto.diesel)
saanto[,9] <- saanto$siemenet * saanto$öljy * saanto$diesel * ala
colnames(saanto)[9] <- "saanto (kg/a)"
 
## Muutetaan tulos jakaumaksi Monte Carlolla.
P <- PTable(P, n)
saanto <- merge(P, saanto)
 
## Lasketaan tulostaulu.
if(length(divisions)>1) divisions <- as.list(saanto[, divisions]) else divisions <- saanto[, divisions]
out1 <- as.data.frame(as.table(tapply(saanto[, 10], divisions, mean)))  
out1 <- dropall(out1[!is.na(out1$Freq), ])
 
print(xtable(out1), type = 'html')
 
## Lasketaan tuloskuvaaja.
out2 <- as.data.frame(as.table(tapply(saanto[, 10], list(saanto[, divisions2], saanto$obs), mean)))
out2 <- dropall(out2[!is.na(out2$Freq), ])
ggplot(out2, aes(x = Freq, weight = 1, fill = Var1)) +geom_density()
## Jostain syystä vain osa kuvista piirtyy oikein, riippuen mitä parametreja valitaan. En ymmärrä syytä.
</rcode>
 
 
{{bluebox|
==Komia loota==
* sininen
 
on komia...
}}
 
<t2b index="Year,Country" unit="#" obs="Mortality">
2001|Finland|1001
2002|Finland|judanss
2001|Sweden|900
2002|Sweden|800
</t2b>
 
{{greenbox|
JEPJEJPKEPKE
 
JOO
}}
 
fwae
 
JEPULISTA
<dpl>
category = THL:n julkaisuja vuonna 2010
mode = ordered
include = {julkaisu}
ordermethod=title
format  = ,\n<big>'''[[%PAGE%|%TITLE%]]'''</big>,,
notcategory = glossary
order = ascending
</dpl>
[[Luokka:testi]]
 
 
{{budjettirajapinta
| työpaketit        = 'TP1';TP1;'TP2';TP2;'TP3';TP3;'TP4';TP4;'TP5';TP5
| vuodet            = '2012';2012;'2013';2013;'2014;2014
| sivu              = Op_fi2682
| rahoittaja        = STM
| projekti          = Tekaisu
| työpaketti        = TP1
| työpaketti.add    = TP1
| rahoittaja.add    = STM
| työpaketit_default = 'TP1';'TP2';'TP3';'TP5'
| rahoittaja_default = STM
}}

Versio 12. tammikuuta 2012 kello 11.53

Päivitetty

<math>123+\alpha -555-666</math>


Terveysvaikutukset

HUOM! Tällä hetkellä malli laskee ainoastaan murskausprosessin pölypäästöt ja niiden leviämisen 10-15 km:n säteellä olevaan väestöön. Mallia päivitetään jatkuvasti.

Hassu kuva

Hassu kuva

Hassumpi kuva

Hassuin kuva

Leveyspiiri (LA):

Pituuspiiri (LO):

Murskattavan malmin määrä (Mg/a) (murskaus.maara):

Murskaukseen käytettävän prosessin päästökerroin (g/Mg) (paasto.kerroin):

Iteraatioiden määrä (N):

+ Näytä koodi

R-testi

Jatropan viljelyala (ha):

n:

Mitkä tekijät halua eritellä tuloksessa?:
Katalyytin määrä
Ikä
Kastelu
Käytetty puristin

Minkä yhden tekijän halua eritellä kuvaajassa?:

+ Näytä koodi


Komia loota

  • sininen

on komia...

Ero sivun ”Hiekkalaatikko” versioiden välillä(#)
ObsYearCountryMortality
12001Finland1001
22002Finlandjudanss
32001Sweden900
42002Sweden800

JEPJEJPKEPKE

JOO

fwae

JEPULISTA

Apua typeryyksien keksijöille!


Jouni T. Tuomisto: Apua typeryyksien keksijöille!. Ajatelma toimintatapojen muuttamisesta. Opasnet 2010. Viitattu 09.01.2025.

Argumenta


Sami Majaniemi, Mikko Pohjola, Jouni T. Tuomisto: Argumenta. Opasnet 2010. Viitattu 09.01.2025.

Arseenialtistuminen Suomessa


Hannu Komulainen: Arseeni altistuminen Suomi. Opasnet 2010. Viitattu 09.01.2025.

Arviointi THL:n budjetista


Jouni T. Tuomisto: Arviointi THL:n budjetista. Opasnet 2010. Viitattu 09.01.2025.

Arviointi kirjailijoiden tulojen turvaamisesta


Jouni T. Tuomisto: Arviointi kirjailijoiden tulojen turvaamisesta. Opasnet 2010. Viitattu 09.01.2025.

Arviointi pohjavesilaitoksen mikrobiologisista riskeistä


Päivi Meriläinen: Arviointi pohjavesilaitoksen mikrobiologisista riskeistä. Opasnet 2010. Viitattu 09.01.2025.


Päivi Meriläinen: Arviointi pohjavesilaitoksen mikrobiologisista riskeistä. Opasnet 2010. Viitattu 09.01.2025.

Arvostuspeli


Jouni T. Tuomisto: Arvostuspeli. Opasnet 2010. Viitattu 09.01.2025.

Arvostusteoria. Ehdotus markkinatalouden aiheuttaman riippuvuuden hoitoon


Jouni T. Tuomisto: Arvostusteoria. Ehdotus markkinatalouden aiheuttaman riippuvuuden hoitoon. Luento. Opasnet 2010. Viitattu 09.01.2025.

Avoin arviointi erityistilanteessa


Jouni T. Tuomisto: Avoin arviointi erityistilanteessa. Ohjeistus. Opasnet 2010. Viitattu 09.01.2025.

Energiantuotannon kasvihuonekaasupäästöt Suomessa


Matleena Tuomisto, Noora Koponen, Teemu Rintala: Energiantuotannon kasvihuonekaasupäästöt Suomessa. Opasnet 2010. Viitattu 09.01.2025.

Fluoridi altistuminen


Hannu Komulainen, Olli Leino: Fluoridi altistuminen. Muuttuja. Opasnet 2010. Viitattu 09.01.2025.

Formaldehydi altistus Suomi


Eero Priha, Otto Hänninen, Olli Leino: Formaldehydi altistus Suomi. Opasnet 2010. Viitattu 09.01.2025.

Heande-projektinhallinta


Jouni T. Tuomisto: Heande-projektinhallinta. Projektinhallintajärjestelmän esittely. Opasnet 2010. Viitattu 09.01.2025.

Hyöty-riskiarviointi ydinvoimasta Suomessa


Teemu Rintala, Sami Majaniemi, Minna-Mari Hämäläinen, Noora Koponen, Matleena Tuomisto, Aino Paakkinen, Jouni T. Tuomisto: Hyöty-riskiarviointi ydinvoimasta Suomessa. Opasnet 2010. Viitattu 09.01.2025.

Ilmastokäräjät


Jouni T. Tuomisto: Ilmastokäräjät. Tutkimussuunnitelma. Opasnet 2010. Viitattu 09.01.2025.

Ilmastonmuutos


Matleena Tuomisto, Jouni T. Tuomisto: Ilmastonmuutos. Opasnet 2010. Viitattu 09.01.2025.

Ilmastopolitiikkaa voisi ohjata tieteen menetelmin


Jouni T. Tuomisto, Mikko V. Pohjola: Ilmastopolitiikkaa voisi ohjata tieteen menetelmin. Mielipidekirjoitus. Opasnet 2010. Viitattu 09.01.2025.

Iskuryhmä ilmastonmuutosta vastaan


Jouni T. Tuomisto: Iskuryhmä ilmastonmuutosta vastaan. Pamfletti. Opasnet 2010. Viitattu 09.01.2025.

Islannin tulivuorenpurkaus 2010


Jouni T. Tuomisto, Pauliina Ahtoniemi: Islannin tulivuorenpurkaus 2010. Opasnet 2010. Viitattu 09.01.2025.

Kalaruoan terveysvaikutukset


Olli Leino: Kalaruoan terveysvaikutukset. Opasnet 2010. Viitattu 09.01.2025.

Kaupungin verkkosivujen määrittelyt


Jouni T. Tuomisto: Kaupungin verkkosivujen määrittelyt. Opasnet 2010. Viitattu 09.01.2025.

Kaupunkien lämpösaarekkeiden lieventämisen terveysvaikutukset EuroopassaMalline:Julkaisu.default Kaupunkirakenne - kansanterveysMalline:Julkaisu.default Kloorauksen sivutuotteet altistus Suomi


Päivi Meriläinen, Olli Leino: Kloorauksen sivutuotteet altistus Suomi. Opasnet 2010. Viitattu 09.01.2025.

Komiteoita, tiedostoja ja muita entisajan keksintöjä


Jouni T. Tuomisto: Komiteoita, tiedostoja ja muita entisajan keksintöjä. Opasnet 2010. Viitattu 09.01.2025.

Kosteusvauriot altistus


Ulla Haverinen-Shaughnessy, Olli Leino: Kosteusvauriot altistus. Opasnet 2010. Viitattu 09.01.2025.

Kuntakone


Sami Majaniemi: Kuntakone. Opasnet 2010. Viitattu 09.01.2025.

Lyijy ympäristö altistus


Otto Hänninen, Olli Leino: Lyijy ympäristö altistus. Opasnet 2010. Viitattu 09.01.2025.

Matkojenhallintajärjestelmä


Jouni T. Tuomisto: Matkojenhallintajärjestelmä. Ehdotus parannuksesta. Opasnet 2010. Viitattu 09.01.2025.

Metyylielohopea altistus Suomi


Olli Leino, Juha Villman: Metyylielohopea altistus Suomi. Opasnet 2010. Viitattu 09.01.2025.

Narkolepsia Suomessa


Noora Koponen, Minna-Mari Hämäläinen: Narkolepsia Suomessa. Muuttuja. Opasnet 2010. Viitattu 09.01.2025.

Nettiajan kansalaisyhteiskunta


Sami Majaniemi: Nettiajan kansalaisyhteiskunta. Opasnet 2010. Viitattu 09.01.2025.

Opasnet-uskomusjärjestelmä


Jouni T. Tuomisto: Opasnet-uskomusjärjestelmä. Opasnet 2010. Viitattu 09.01.2025.

Opasnetin kehitystarpeita


Jouni T. Tuomisto: Opasnetin kehitystarpeita. Yhteistyömahdollisuuksien kartoitus. Opasnet 2010. Viitattu 09.01.2025.

Opasnetin käyttökoulutus


Jouni T. Tuomisto, Mikko Pohjola, Juha Villman: Opasnetin käyttökoulutus. Koulutustilaisuus. Opasnet 2010. Viitattu 09.01.2025.

Opastopian toimintaympäristö


Sami Majaniemi: Opastopian toimintaympäristö. Opasnet 2010. Viitattu 09.01.2025.

Postilaki


Jouni T. Tuomisto: Postilaki. Ehdotus parannuksesta. Opasnet 2010. Viitattu 09.01.2025.

Päättäjän tietomonopoli on murrettava


Jouni T. Tuomisto: Päättäjän tietomonopoli on murrettava. Pamfletti. Opasnet 2010. Viitattu 09.01.2025.

Radioaktiivisen säteilyn terveysvaikutukset Suomessa


Teemu Rintala, Noora Koponen, Minna-Mari Hämäläinen, Jouni T. Tuomisto: Radioaktiivisen säteilyn terveysvaikutukset Suomessa. Opasnet 2010. Viitattu 09.01.2025.

Seurakuntavaalit 2010


Jouni T. Tuomisto: Seurakuntavaalit 2010. Tutkimussuunnitelma. Opasnet 2010. Viitattu 09.01.2025.

Sikainfluenssa


Noora Koponen, Jouni T. Tuomisto: Sikainfluenssa. Luonnos arvioinniksi. Opasnet 2010. Viitattu 09.01.2025.

Sikainfluenssarokotteen terveyshaitat


Noora Koponen, Jouni T. Tuomisto: Sikainfluenssarokotteen terveyshaitat. Opasnet 2010. Viitattu 09.01.2025.

Suomalainen energiapolitiikka


Teemu Rintala, Noora Koponen, Matleena Tuomisto, Sami Majaniemi, Aino Paakkinen, Jouni T. Tuomisto: Suomalainen energiapolitiikka. Opasnet 2010. Viitattu 09.01.2025.

Suomalaiset ydinvoimalat


Teemu Rintala, Minna-Mari Hämäläinen, Noora Koponen, Jouni T. Tuomisto: Suomalaiset ydinvoimalat. Muuttuja. Opasnet 2010. Viitattu 09.01.2025.

Suomen energiantuotannon terveysvaikutusten arviointi


Teemu Rintala, Sami Majaniemi: Suomen energiantuotannon terveysvaikutusten arviointi. Opasnet 2010. Viitattu 09.01.2025.

Syövät työperäiset lukumäärä


Antti Karjalainen, Olli Leino: Syövät työperäiset lukumäärä. Opasnet 2010. Viitattu 09.01.2025.

Sähköinen sairauskertomus


Jouni T. Tuomisto: Sähköinen sairauskertomus. Opasnet 2010. Viitattu 09.01.2025.

THL:n kannanottoja ydinvoimasta


Matti Jantunen, Jouko Tuomisto, Raimo Salonen: THL:n kannanottoja ydinvoimasta. Opasnet 2010. Viitattu 09.01.2025.

Terveyden ja hyvinvoinnin ajatushautomo


Jouni T. Tuomisto, Sami Majaniemi: Terveyden ja hyvinvoinnin ajatushautomo. Opasnet 2010. Viitattu 09.01.2025.

Terveydenhuollon tietojärjestelmä


Jouni T. Tuomisto: Terveydenhuollon tietojärjestelmä. Alustava arviointi. Opasnet 2012. Viitattu 09.01.2025.

Tiede 2.0


Jouni T. Tuomisto: Tiede 2.0. Opasnet-muuttuja. Opasnet 2010. Viitattu 09.01.2025.

Tieteellis-poliittinen vallankumous


Jouni T. Tuomisto: Tieteellis-poliittinen vallankumous. Ajatelma. Opasnet 2010. Viitattu 09.01.2025.

Tieteen ja politiikan vuorovaikutusstrategiaMalline:Julkaisu.default Tietoteknisiä ratkaisuja Opasnetissä


Jouni T. Tuomisto, Juha Villman: Tietoteknisiä ratkaisuja Opasnetissä. Opasnet 2010. Viitattu 09.01.2025.

Tietotyö


Jouni T. Tuomisto: Tietotyö. Opasnet 2010. Viitattu 09.01.2025.

Tutkimuslaitoksen tietopolitiikka


Jouni T. Tuomisto, Sami Majaniemi: Tutkimuslaitoksen tietopolitiikka. Opasnet 2010. Viitattu 09.01.2025.

Uutisraivaaja 2011


Sami Majaniemi: Uutisraivaaja 2011. Suunnitelma. Opasnet 2010. Viitattu 09.01.2025.

Valtion palkkausjärjestelmä


Jouni T. Tuomisto: Valtion palkkausjärjestelmä. Arviointi. Opasnet 2010. Viitattu 09.01.2025.

Verovähennyskortti


Jouni T. Tuomisto: Verovähennyskortti. Verojärjestelmän parannusehdotus. Opasnet 2010. Viitattu 09.01.2025.

Wikipuolue


Jouni T. Tuomisto: Wikipuolue. Opasnet 2010. Viitattu 09.01.2025.

YVA


Mikko Pohjola, Erkki Kuusisto, Jouni T. Tuomisto: YVA. Kommentteja YVA-lain arvioinnista. Opasnet 2010. Viitattu 09.01.2025.

Ydinvoimaan liittyvien riskien kokeminen Suomessa


Noora Koponen, Minna-Mari Hämäläinen, Matleena Tuomisto, Teemu Rintala, Aino Paakkinen, Sami Majaniemi, Jouni T. Tuomisto: Ydinvoimaan liittyvien riskien kokeminen Suomessa. Opasnet 2010. Viitattu 09.01.2025.

Ympäristömelun vaikutuslaskentatiedot Suomi


Erkki Kuusisto: Ympäristömelun vaikutuslaskentatiedot Suomi. Opasnet 2010. Viitattu 09.01.2025.

Ympäristön tupakansavu altistus Suomi


Otto Hänninen, Erkki Kuusisto, Olli Leino: Ympäristön tupakansavu altistus Suomi. Opasnet 2010. Viitattu 09.01.2025.

Ytrips OyMalline:Julkaisu.default


Tämä koodi tuottaa vain tapahtumalistauksen ilman yhteenvetoja.

Voit rajata tarkastelusi yhteen työpakettiin tai vuoteen.

Tarkasteltava työpaketti:

Tarkasteltava vuosi:

+ Näytä koodi


Tämä tuottaa halutut yhteenvedot mutta ei tapahtumalistausta.

Voit valita näytettävät yksityiskohdat ja myös rajata tarkastelusi yhteen työpakettiin tai vuoteen.

Mitkä jaottelut haluat näyttää yhteenvedossa?:
Projekti
Työpaketti
Vuosi
Laitos
Rahoittaja
Kustannuslaji

Tarkasteltavat työpaketit:
TP1
TP2
TP3
TP4
TP5

Tarkasteltava vuosi:

Voit testata budjettia tällä tilapäisellä lisätapahtumalla.

Työpaketti:

Rahoittaja:

Kustannuslaji:

Budjettiin lisättävä summa:

Kuvaus:

+ Näytä koodi


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