Difference between revisions of "Benefit-risk assessment of Baltic herring and salmon intake"
(→Calculations) |
(→Calculations: bayes model and graphs separated. Intermediate results stored as arrays) |
||
Line 108: | Line 108: | ||
* Model run 1.3.2017 [http://en.opasnet.org/en-opwiki/index.php?title=Special:RTools&id=3Xu19vkWK1lyWVg3] | * Model run 1.3.2017 [http://en.opasnet.org/en-opwiki/index.php?title=Special:RTools&id=3Xu19vkWK1lyWVg3] | ||
− | <rcode name="bayes" label=" | + | <rcode name="bayes" label="Sample Bayes model (for developers only)"> |
− | # This is code Op7748/ on page [[Benefit-risk assessment of Baltic herring and salmon intake]] | + | # This is code Op7748/bayes on page [[Benefit-risk assessment of Baltic herring and salmon intake]] |
library(OpasnetUtils) | library(OpasnetUtils) | ||
library(reshape2) | library(reshape2) | ||
− | |||
library(rjags) | library(rjags) | ||
Line 138: | Line 137: | ||
# Total amount eaten of each fish per modelled individual is finally calculated. | # Total amount eaten of each fish per modelled individual is finally calculated. | ||
− | objects.latest("Op_en7749", "preprocess") | + | objects.latest("Op_en7749", "preprocess") # [[Goherr: Fish consumption study]]: survey, surcol |
# The following three rows should be in the preprocessing code. | # The following three rows should be in the preprocessing code. | ||
survey$Row <- 1:nrow(survey) | survey$Row <- 1:nrow(survey) | ||
− | survey$Weighting <- as.double( | + | survey$Weighting <- as.double(as.character(survey$Weighting)) |
survey$Ages <- ifelse(as.numeric(as.character(survey$Age)) < 46, "18-45",">45") | survey$Ages <- ifelse(as.numeric(as.character(survey$Age)) < 46, "18-45",">45") | ||
# Interesting fish eating questions | # Interesting fish eating questions | ||
surv <- survey[c(1,3,4,16,29,30,31,46:49,95:98,158)] | surv <- survey[c(1,3,4,16,29,30,31,46:49,95:98,158)] | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
agel <- unique(survey$Ages) | agel <- unique(survey$Ages) | ||
Line 234: | Line 188: | ||
for(i in 1:F) { # F = fish species | for(i in 1:F) { # F = fish species | ||
for(j in 1:C) { # C = congener | for(j in 1:C) { # C = congener | ||
− | mu[i,j] ~ dunif(-3,3)#0.001 | + | mu[i,j] ~ dunif(-3,3) # Why does this not work with dnorm(0, 0.001)? |
tau[i,j] <- pow(sigma[i,j], -2) | tau[i,j] <- pow(sigma[i,j], -2) | ||
sigma[i,j] ~ dunif(0, 10) | sigma[i,j] ~ dunif(0, 10) | ||
Line 272: | Line 226: | ||
Q = ncol(ans), | Q = ncol(ans), | ||
ans = ans, | ans = ans, | ||
− | qlen = qlen, | + | qlen = qlen-1, # Subtract 1 because starts from 0. |
age = match(survey$Ages, agel), | age = match(survey$Ages, agel), | ||
country = match(survey$Country, countryl), | country = match(survey$Country, countryl), | ||
Line 292: | Line 246: | ||
#temp <- adply(samps$mu, c(1,2,3,4)) | #temp <- adply(samps$mu, c(1,2,3,4)) | ||
− | + | pcd.pred <- array( | |
− | + | samps$pcd.pred, | |
− | |||
− | |||
dim = c(length(conl), length(fisl), 1000, 4), | dim = c(length(conl), length(fisl), 1000, 4), | ||
dimnames = list( | dimnames = list( | ||
Line 303: | Line 255: | ||
Seed = c("S1","S2","S3","S4") | Seed = c("S1","S2","S3","S4") | ||
) | ) | ||
− | )) | + | ) |
+ | |||
+ | mu.pred <- array( | ||
+ | samps$mu, | ||
+ | dim = c(length(conl), length(fisl), 1000, 4), | ||
+ | dimnames = list( | ||
+ | Fish = fisl, | ||
+ | Congener = conl, | ||
+ | Iter = 1:1000, | ||
+ | Seed = c("S1","S2","S3","S4") | ||
+ | ) | ||
+ | ) | ||
+ | |||
+ | ans.pred <- array( | ||
+ | samps$ans.pred, | ||
+ | dim = c( | ||
+ | length(agel), | ||
+ | length(countryl), | ||
+ | length(genderl), | ||
+ | ncol(ans), | ||
+ | 1000, # iterations | ||
+ | 4 # Seeds | ||
+ | ), | ||
+ | dimnames = list( | ||
+ | Age = agel, | ||
+ | Country = countryl, | ||
+ | Gender = genderl, | ||
+ | Question = colnames(ans), | ||
+ | Iter = 1:1000, | ||
+ | Seed = c("S1","S2","S3","S4") | ||
+ | ) | ||
+ | ) | ||
+ | |||
+ | objects.store(pcd.pred, mu.pred, ans.pred) | ||
+ | cat("Arrays pcd.pred, mu.pred, ans.pred stored.\n") | ||
+ | </rcode> | ||
+ | |||
+ | === Plot concentrations and survey === | ||
+ | |||
+ | * Model run 1.3.2017 [http://en.opasnet.org/en-opwiki/index.php?title=Special:RTools&id=ewXNdifWMuu4emJf] | ||
+ | |||
+ | <rcode label="Plot concentrations and answers" graphics=1> | ||
+ | #This is code Op_en7748/ on page [[Benefit-risk assessment of Baltic herring and salmon intake]] | ||
+ | |||
+ | library(OpasnetUtils) | ||
+ | library(ggplot2) | ||
+ | library(reshape2) | ||
+ | |||
+ | objects.latest("Op_en7748", code_name = "bayes") #: pcd.pred, ans.pred, mu.pred | ||
+ | objects.latest("Op_en3104", code_name = "preprocess") # [[EU-kalat]]: eu | ||
+ | |||
+ | pl <- melt(pcd.pred) | ||
+ | mul <- melt(mu.pred) | ||
+ | ql <- melt(ans.pred) | ||
+ | |||
+ | ############## PCDDF concentrations plotted | ||
ggplot(eu@output, aes(x = euResult, colour = Fish))+geom_density(size = 1) + | ggplot(eu@output, aes(x = euResult, colour = Fish))+geom_density(size = 1) + | ||
Line 336: | Line 343: | ||
facet_wrap(~ Fish, scales = "free_y") + coord_cartesian(xlim = c(1E-3,1E2))+ | facet_wrap(~ Fish, scales = "free_y") + coord_cartesian(xlim = c(1E-3,1E2))+ | ||
labs(title = "Selected fish species") | labs(title = "Selected fish species") | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
########## Questionnaire plotted | ########## Questionnaire plotted | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
for(i in colnames(ans)) { | for(i in colnames(ans)) { | ||
− | print(ggplot(ql[ql$Question == i,], aes(x = value, | + | print(ggplot(ql[ql$Question == i,], aes(x = value, fill = Gender))+ |
− | + | geom_bar(position = "dodge") + | |
facet_grid(Country ~ Age)+ | facet_grid(Country ~ Age)+ | ||
labs(title = paste(i, "from bayes survey model"))) | labs(title = paste(i, "from bayes survey model"))) | ||
} | } | ||
− | # | + | ##################### Mu plotted |
− | # | ||
− | |||
− | # | ||
− | # | ||
− | # | ||
− | # | ||
ggplot(mul, aes(x = value, colour = Congener))+stat_ecdf(size = 1) + | ggplot(mul, aes(x = value, colour = Congener))+stat_ecdf(size = 1) + | ||
Line 398: | Line 365: | ||
labs(title = "Mixing of model. Runs 1..1000, no thinning") | labs(title = "Mixing of model. Runs 1..1000, no thinning") | ||
+ | #tef <- Ovariable("tef", ddata = "Op_en4017", subset = "TEF values") | ||
+ | #tef <- EvalOutput(tef) | ||
+ | |||
+ | #levels(eu$Congener) <- gsub("HCDD", "HxCDD", levels(eu$Congener)) | ||
+ | #levels(eu$Congener) <- gsub("HCDF", "HxCDF", levels(eu$Congener)) | ||
+ | #levels(eu$Congener) <- gsub("CoPCB", "PCB", levels(eu$Congener)) | ||
+ | #euteq <- eu * tef | ||
</rcode> | </rcode> | ||
Revision as of 12:46, 1 March 2017
This page is a assessment.
The page identifier is Op_en7748 |
---|
Moderator:Arja (see all) |
This page is a stub. You may improve it into a full page, and then a rating bar will appear here. |
Upload data
|
Scope
This assessment is part of the WP5 work in Goherr project. Purpose is to evaluate health benefits and risks caused of eating Baltic herring and salmon in four Baltic sea countries (Denmark, Estonia, Finland and Sweden). This assessment is currently on-going.
Question
What are the current population level health benefits and risks of eating Baltic herring and salmon in Finland, Estonia, Denmark and Sweden? How would the health effects change in the future, if consumption of Baltic herring and salmon changes due to actions caused by different management scenarios of Baltic sea fish stocks?
Intended use and users
Results of this assessment are used to inform policy makers about the health impacts of fish. Further, this assessment will be combined with the results of the other Goherr WPs to produce estimates of future health impacts of Baltic fish related to different policy options. Especially, results of this assessment will be used as input in the decision support model built in Goherr WP6.
Participants
- National institute for health and welfare (THL)
- Goherr project group
Boundaries
- Four baltic sea countries (Denmark, Estonia, Finland, Sweden)
- Current situation (fish use year 2016, pollutant levels in fish year 2010?)
- Estimation for future (year 2020?)
Decisions and scenarios
Management scenarios developed in Goherr WP3 frames the following boundaries to the use and consumption of Baltic herring and salmon as human food. Effect of these scenarios to the dioxin levels and the human food use will be evalauted quantitatively and feed into the health benefit-risk model to assess the health effect changes.
- Scenario 1: “Transformation to sustainability”
- Hazardous substances, including dioxins, are gradually flushed out and the dioxin levels in Baltic herring are below or close to the maximum allowable level.
- Fish stocks are allowed to recover to levels, which makes maximum sustainable yield possible and increases the total catches of wild caught fish. The catches of salmon by commercial fisheries has stabilized at low level, while the share of recreational catch increases slightly.
- The use of the Baltic herring catch for food increases. A regional proactive management plan for the use of catch has increased the capacity of the fishing fleets to fish herring for food and through product development and joint marketing, have increased consumer demand for Baltic herring.
- Scenario 2: “Business-as-usual”
- The commercial catches of salmon continue to decrease. The demand for top predatory species, such as salmon and cod remains high, while the demand for herring decreased further as a result of demographic changes.
- Most of the herring catch are used for fish meal and oil production in the region.
- The use of Baltic herring from the southern parts of the Baltic Sea where the dioxin contents are not likely to exceed the maximum allowable level, are prioritised for human consumption. In the absence of the demand in many of the Baltic Sea countries, majority of the herring intended for direct human consumption are exported to Russia.
- Scenario 3: “Inequality”
- The nutrient and dioxins levels continue to decrease slowly.
- The commercial catches of salmon have decreased further as the general attitudes favour recreational fishing, which has also resulted in decreased demand.
- The herring catches have increased slightly, but the availability of herring suitable for human consumption remains low due to both, dioxin levels that remain above the maximum allowable limit in the northern Baltic Sea and the poor capacity to fish for food.
- The use of the catch varies between countries. In Estonia, for example, where the whole catch has been traditionally used for human consumption, there is no significant change in this respect, but in Finland, Sweden and Denmark, herring fishing is predominantly feed directed.
- Scenario 4: “Transformation to protectionism”
- The level of hazardous substances also increases as emission sources are not adequately addressed.
- Commercial salmon fisheries disappears almost completely from the Baltic Sea, although restocking keeps small scale fisheries going.
- Many of the Baltic herring stocks are also fished above the maximum sustainable yield and total catches are declining.
- Owing to the growing dioxin levels detected in herring, majority of the catch is used for aquaculture.
Timing
- Model development during 2016 and 2017
- First set of results in March 2017, draft publication in March 2018
Answer
This section will be updated as soon as preliminary results are available
Results
Conclusions
Rationale
Stakeholders
- Policy makers
- Food safety authorities
- Fisheries management
- Researchers
- Food safety
- Health
- NGO's
- WWF
- Active consumers
- Marine Stewardship Council
- Baltic sea fishers and producers?
Dependencies
Content of all variable pages below are preliminary and will be updated when this assessment proceeds.
- Toxic equivalency factor (TEF)
- Consumption of fish
- Pollutant and fatty acid levels
- Persistent organic pollutant levels in Baltic herring
- Persistent organic pollutant levels in Baltic salmon
- Nutrient levels in fish --# : Includes three tables, content modelled with Analytica years ago. How to update these? --Arja (talk) 07:20, 20 October 2016 (UTC)
- Risk functions of health effects
- Background health data (incidences and DALY)
- Disability weights of health effects
- Length of disease
Analyses
Indices
- Country (Denmark, Estonia, Finland, Sweden)
- Year (current, future)
- Sex
- Age (categories?)
- Fish species (Baltic herring, Baltic salmon)
- Health end-point (ICD-10 code?, name? something else? --# : This need to be matched for disability weight, duration of disease, background incidence and disease burden data, dose-responses --Arja (talk) 09:39, 22 September 2016 (UTC)
- Compound (pollutant, fatty acid) --# : This needs to be matched with dose-responses --Arja (talk) 09:39, 22 September 2016 (UTC)
Calculations
This section will have the actual health benefit-risk model (schematically described in the above figure) written with R. The code will utilise all variables listed in the above Dependencies section. Model results are presented as tables and figures when those are available.
- Model run 28.2.2017 [1]
- Model run 28.2.2017 with corrected survey model [2]
- Model run 28.2.2017 with Mu estimates [3]
- Model run 1.3.2017 [4]
Plot concentrations and survey
- Model run 1.3.2017 [5]
References
Keywords
See also
- Risk and Benefit Assessment of Herring and Salmonid Fish from the Baltic Sea Area
- Swedish Market Basket 2010
- Riksmaten 2010
- Danskernes kostvaner 2011-2013