Ero sivun ”Järvisedimenttien metallipitoisuudet” versioiden välillä
Siirry navigaatioon
Siirry hakuun
(datanpäivityskoodi kirjoitettu ja dataa siistitty) |
(koodit päivitetty vastaamaan uutta, parempaa datarakennetta) |
||
Rivi 18: | Rivi 18: | ||
name:säde|description:Määritä säde (km; järvet säteen määräämän ympyränmuotoisen alueen sisäpuolella huomioidaan)|default:100 | name:säde|description:Määritä säde (km; järvet säteen määräämän ympyränmuotoisen alueen sisäpuolella huomioidaan)|default:100 | ||
"> | "> | ||
library(OpasnetUtils) | library(OpasnetUtils) | ||
library(xtable) | library(xtable) | ||
library(ggplot2) | library(ggplot2) | ||
data <- opbase.data("Op_en5462") | data <- opbase.data("Op_en5462") | ||
data$Easting <- as.numeric(as.character(data$Easting)) | |||
data$Northing <- as.numeric(as.character(data$Northing)) | |||
data <- | |||
out <- data | out <- data | ||
out$Distance <- sqrt((out$Easting - easting)^2 + (out$Northing - northing)^2) | |||
out$ | out$Distance <- out$Distance <= säde*1000 | ||
colnames( | out2 <- as.data.frame(as.table(tapply(out$Result, out[c("Metal", "Distance")], mean))) | ||
colnames(out2) <- gsub("Freq", "Pitoisuus", colnames(out2)) | |||
cat("Metallipitoisuuden keskiarvo ", säde, " km säteellä pisteestä ", easting, "E, ", northing, " N.\n") | cat("Metallipitoisuuden keskiarvo ", säde, " km säteellä pisteestä ", easting, "E, ", northing, " N.\n") | ||
print(xtable( | print(xtable(out2[out2$Distance == TRUE, ]), type = 'html') | ||
cat("Pitoisuudet metallille", metalli, ".\n") | cat("Pitoisuudet metallille", metalli, ".\n") | ||
ggplot(out[out$Metal == metalli, ], aes(x = Easting, y = Northing, size = Result)) +geom_point() # Concentrations on a "map". | |||
</rcode> | </rcode> | ||
Rivi 104: | Rivi 53: | ||
library(xtable) | library(xtable) | ||
data <- tidy(opbase.data("Op_en5462"), objname = "conc | data <- tidy(opbase.data("Op_en5462"), objname = "conc") | ||
data <- data[data$Metal %in% metalli & data$Lake.name == Järvi, ] | |||
data <- data[data$ | |||
conc <- new("ovariable", | conc <- new("ovariable", | ||
name = "conc", | name = "conc", | ||
Rivi 130: | Rivi 76: | ||
==Perustelut== | ==Perustelut== | ||
Datan rakennetta (ei sisältöä) parannettu tällä koodilla: [http://fi.opasnet.org/fi_wiki/index.php?title=J%C3%A4rvisedimenttien_metallipitoisuudet&oldid=16253#Update_data_structure]. | |||
== Katso myös == | == Katso myös == |
Versio 19. joulukuuta 2012 kello 09.07
Edistymisluokitus |
---|
Opasnetissa lukuisat sivut ovat työn alla eri vaiheissa. Niiden tietosisältöön pitää siis suhtautua harkiten. Tämän sivun sisällön edistyminen on arvioitu:
|
Moderaattori:Marjo (katso kaikki)
Sivun edistymistä ei ole arvioitu. Arvostuksen määrää ei ole arvioitu (ks. peer review). |
Lisää dataa
|
Kysymys
Miten lasketaan suomalaisen järven sedimentin metallipitoisuus?
Vastaus
1) Lasketaan järven sedimentin metallipitoisuus perustuen ympäröivän, ympyrän säteen määräämän alueen järvien keskipitoisuuteen.
2) Lasketaan järven sedimentin metallipitoisuus läheisen/sopivan suuren järven pitoisuuden perusteella.
--#: Koodia voisi siistiä seuraavasti:
- joka järvi omalle rivilleen, niin on helpompi lukea (sieltä puuttuu muutamia lainausmerkkejä, mikä kaataa koodin jos valitsee sellaisen järven.
- Koodista voisi tehdä ovariablen nimeltä sediment.concentration (ei conc), joka ottaa syötteeksi tiedot Järvi, metalli, coeff.pH ja coeff.landuse (pitäisi olla johdonmukaisia kielen kanssa).
- Data kaivetaan kannasta formulan sisällä.
- Jos syötteet coeff.pH ja coeff.landuse voidaan johtaa jostain, sen johtamisen voisi tehdä käyttämällä syötteenä pH ja landuse ja laskemalla kertoimet.
--Jouni Tuomisto 18. joulukuuta 2012 kello 17.09 (EET)
Perustelut
Datan rakennetta (ei sisältöä) parannettu tällä koodilla: [5].
Katso myös
Perustuu GTK:n dataan. GTK.
Viitteet
Lähdeluettelo tulostuu tänne automaattisesti kun viittauksissa käytetään <ref></ref> -tagia. Lähteet voi myös lisätä halutessaan manuaalisesti.