Question
What are metal concentrations in lake sediments? The focus is on Finland.
Answer
Show results
+ Show code- Hide code
library(OpasnetBaseUtils)
library(xtable)
library(ggplot2)
data <- (op_baseGetData("opasnet_base", "Op_en5462"))
data <- tidy(data)
#out <- data[data$Observation %in% metalli, ]
#colnames(out) <- c("DateTime", "TagName", "Value")
#out$Value <- as.numeric(out$Value)
#out$DateTime <- as.numeric(out$DateTime)
#cat("Metalli: ", metalli, "\n")
data$Result <- as.numeric(data$Result)
data <- reshape(data, timevar = "Observation", idvar = "ID", direction = "wide")
colnames(data) <- gsub("Result.", "", colnames(data))
data$distance <- sqrt((as.numeric(data$Easting) - easting)^2 + (as.numeric(data$Northing) - northing)^2) <= säde*1000
objects.put(data)
#data$Co <- as.numeric(data$Co)
#data$Cr <- as.numeric(data$Cr)
#data$Cu <- as.numeric(data$Cu)
#data$Fe <- as.numeric(data$Fe)
#data$Mn <- as.numeric(data$Mn)
#data$Ni <- as.numeric(data$Ni)
#data$Pb <- as.numeric(data$Pb)
#data$S <- as.numeric(data$S)
#data$V <- as.numeric(data$V)
#data$Zn <- as.numeric(data$Zn)
#as.data.frame(as.table(tapply(data[, metalli], data[c("distance")], mean)))
plot(data[data$distance, metalli])
qplot(Easting, Northing, data = data, size = Co)
library(reshape)
data <- melt(data, measure.vars = c("Co", "Cr", "Cu", "Fe", "Mn", "Ni", "Pb", "S", "V", "Zn"))
#print(xtable(data), type = 'html')
out <- as.data.frame(as.table(tapply(data$value, data[c("variable", "distance")], mean)))
cat("Metallipitoisuuden keskiarvo ", säde, " km säteellä pisteestä ", easting, "E, ", northing, " N.\n")
print(xtable(out[out$distance == TRUE, ]), type = 'html')
qplot(Easting, Northing, data = data[data$Pollutant = metalli, size = value)
| |
Rationale
Based on data from GTK, Finland.
Dependencies
Formula
See also
Keywords
References
Related files
<mfanonymousfilelist></mfanonymousfilelist>
- Metal concentrations in lake sediments. Opasnet . [1]. Accessed 02 Jun 2024.