Mercury concentrations in fish in Finland: Difference between revisions

Revision as of 14:49, 5 October 2017

[show] This page is a knowledge crystal of subtype variable. The page identifier is Op_en4004
Moderator:Jouni (see all)

Upload data {{#opasnet_base_link:Op_en4004}}

Question

What are the current total mercury concentrations in different fish species in Finland?

Answer

Simulated data from a previous study.

Rationale

Data

File:Mercury vs omega-3.ANA

The table below is data from EU-kalat 1 project (2003), and the values are about Baltic herring. (Other fish species were measured but data are not here).[1]

Mercury in Baltic herring(mg /kg f.w.)
Obs	Age	Size	Concentration
1	1	1	0.0005
2	2	1	0.0005
3	3	1	0.0005
4	2	1	0.005
5	2	4	0.005
6	5	2	0.006
7	2	1	0.007
8	1	1	0.008
9	2	1	0.008
10	3	2	0.008
11	3	2	0.008
12	2	2	0.009
13	3	2	0.009
14	3	3	0.01
15	2	2	0.011
16	3	2	0.011
17	1	1	0.012
18	2	1	0.014
19	4	2	0.015
20	4	3	0.015
21	1	1	0.018
22	3	2	0.018
23	2	1	0.02
24	2	2	0.021
25	2	2	0.022
26	4	2	0.022
27	6	4	0.022
28	8	4	0.024
29	6	3	0.025
30	5	2	0.026
31	4	3	0.027
32	5	4	0.027
33	8	4	0.027
34	6	4	0.028
35	5	3	0.03
36	2	1	0.031
37	7	5	0.031
38	3	2	0.033
39	10	4	0.033
40	6	3	0.038
41	10	5	0.038
42	12	5	0.039
43	6	4	0.04
44	5	5	0.043
45	7	5	0.043
46	7	4	0.048
47	6	3	0.049
48	6	3	0.049
49	3	2	0.051
50	9	5	0.054
51	9	4	0.066
52	7	5	0.069
53	3	5	0.072
54	11	4	0.083
55	11	5	0.107

Calculations

+ Show code - Hide code

library(OpasnetUtils)

concentration <- Ovariable("concentration", ddata = "Op_en4004.mercury_in_baltic_herring")

objects.store(concentration)

cat("Ovariable concentration stored.\n")

+ Show code - Hide code

# This is code Op_en4004/MeHgmean on page [[Mercury concentrations in fish in Finland]]

library(OpasnetUtils)

mehg <- opbase.data(
  "Op_en4004",
  subset = "Kerty database"
#  include = list(Näytelaji = c("lohi", "silakka")) # makes an error due to an unknown reason
)

mehg2 <- mehg[mehg$Näytelaji %in% c("lohi","silakka"),c(8,9,12,15,16,17,20,21,22)] 
mehg2$Conc <- as.numeric(as.character(mehg2$Tulos))
mehg2$Year <- as.numeric(as.character(substr(mehg2$Alkupäivä, 7,8)))
mehg2 <- mehg2[!(is.na(mehg2$Conc) | (mehg2$Year > 13 & mehg2$Näytelaji == "silakka")) , ] 
mehg2[mehg2$Näytelaji == "lohi" , "Helcom-alue"] <- ""

oprint(aggregate(mehg2["Conc"], by = mehg2[c("Helcom-alue", "Näytelaji")], FUN=length))
oprint(aggregate( 
  mehg2["Conc"], 
  by = mehg2[c("Näytelaji", "Year")],
  FUN=length
))

MeHgmean <- aggregate(mehg2["Conc"], by = mehg2[c("Helcom-alue", "Näytelaji")], FUN=mean)
MeHgsd <- aggregate(mehg2["Conc"], by = mehg2[c("Helcom-alue", "Näytelaji")], FUN=sd)

objects.store(MeHgmean, MeHgsd)
cat("Data.frames MeHgmean and MeHgsd stored.\n")

References