Goherr: Fish consumption study

This page is a knowledge crystal of subtype variable. The page identifier is Op_en7749
Moderator:Arja (see all)
This page is a stub. You may improve it into a full page.
Upload data {{#opasnet_base_link:Op_en7749}}

Question

How Baltic herring and salmon are used as human food in Baltic sea countries? Which determinants affect on people’s eating habits of these fish species?

Answer

Survey data will be analysed during winter 2016-2017 and results will be updated here.

• Model run 13.3.2017

+ Show code- Hide code library(OpasnetUtils) library(ggplot2) library(reshape2) objects.latest("Op_en7749", "preprocess") for(i in c(5:6, 16, 29:30, 46:49, 85:86, 95:98, 135)) { temp <- survey[!is.na(survey[[i]]),] p <- ggplot(temp, aes(x = temp[[i]])) + geom_bar() + theme_gray(base_size = 18) + theme(axis.text.x=element_text(angle=90,hjust=1,vjust=0.4)) + labs(title = colnames(temp[i])) + xlab("") + facet_wrap(~ Country) print(p) } temp <- melt(survey, measure.vars = 145:153, variable.name = "Changing.factor", value.name = "Impact") levs <- c("-5 strongly disagree", "-4", "-3", "-2", "-1", "0 Neutral", "1", "2", "3", "4", "5 strongly agree", "I don't know") temp$Impact <- factor(temp$Impact, levels = levs, labels = levs, ordered = TRUE) ggplot(temp, aes(x = Impact, fill = Country)) + geom_histogram(position = "dodge") + facet_wrap(~ Changing.factor)+ theme_gray(base_size = 24) + theme(axis.text.x=element_text(angle=90,hjust=1,vjust=0.4)) surveytemp <- subset(survey, survey[[31]] == "Yes") temp <- melt(surveytemp, measure.vars = 38:43, variable.name = "Variable", value.name = "Value") ggplot(temp, aes(x = Value, fill = Country)) + geom_histogram(position = "dodge") + facet_wrap(~ Variable)+ theme_gray(base_size = 24) + labs(title = "Baltic salmon sources") temp <- melt(surveytemp, measure.vars = 50:59, variable.name = "Variable", value.name = "Value") ggplot(temp, aes(x = Value, fill = Country)) + geom_histogram(position = "dodge") + facet_wrap(~ Variable)+ theme_gray(base_size = 24) + labs(title = "Reasons to eat Baltic salmon") temp <- melt(surveytemp, measure.vars = 73:82, variable.name = "Variable", value.name = "Value") ggplot(temp, aes(x = Value, fill = Country)) + geom_histogram(position = "dodge") + facet_wrap(~ Variable)+ theme_gray(base_size = 24) + theme(axis.text.x=element_text(angle=90,hjust=1,vjust=0.4)) + labs(title = "Baltic salmon actions") surveytemp <- subset(survey, survey[[31]] == "No") temp <- melt(surveytemp, measure.vars = 62:70, variable.name = "Variable", value.name = "Value") ggplot(temp, aes(x = Value, fill = Country)) + geom_histogram(position = "dodge") + facet_wrap(~ Variable)+ theme_gray(base_size = 24) + labs(title = "Reasons not to eat Baltic salmon") surveytemp <- subset(survey, survey[[86]] == "Yes") temp <- melt(surveytemp, measure.vars = 87:92, variable.name = "Variable", value.name = "Value") ggplot(temp, aes(x = Value, fill = Country)) + geom_histogram(position = "dodge") + facet_wrap(~ Variable)+ theme_gray(base_size = 24) + labs(title = "Baltic herring sources") temp <- melt(surveytemp, measure.vars = 99:108, variable.name = "Variable", value.name = "Value") ggplot(temp, aes(x = Value, fill = Country)) + geom_histogram(position = "dodge") + facet_wrap(~ Variable)+ theme_gray(base_size = 24) + labs(title = "Reasons to eat Baltic herring") temp <- melt(surveytemp, measure.vars = 123:132, variable.name = "Variable", value.name = "Value") ggplot(temp, aes(x = Value, fill = Country)) + geom_histogram(position = "dodge") + facet_wrap(~ Variable)+ theme_gray(base_size = 24) + theme(axis.text.x=element_text(angle=90,hjust=1,vjust=0.4)) + labs(title = "Baltic herring actions") surveytemp <- subset(survey, survey[[86]] == "No") temp <- melt(surveytemp, measure.vars = 111:120, variable.name = "Variable", value.name = "Value") ggplot(temp, aes(x = Value, fill = Country)) + geom_histogram(position = "dodge") + facet_wrap(~ Variable)+ theme_gray(base_size = 24) + labs(title = "Reasons not to eat Baltic herring")

Rationale

Survey of eating habits of Baltic herring and salmon in Denmark, Estonia, Finland and Sweden has been done in September 2016 by Taloustutkimus oy. Content of the questionnaire can be accessed in Google drive. The actual data will be uploaded to Opasnet base on Octobere 2016.

The R-code to analyse the survey data will be provided on this page later on.

Data

Original datafile File:Goherr fish consumption.csv

Preprocessing

This code is used to preprocess the original questionnaire data from the above .csv file and to store the data as a usable variable to Opasnet base.

+ Show code- Hide code library(OpasnetUtils) ############## Generic functions and objects are defined first. ### webropol.convert converts a csv file from Webropol into a useful data.frame. webropol.convert <- function( data, # Data.frame created from a Webropol csv file. The first row should contain headings. rowfact, # Row number where the factor levels start (in practice, last row + 3) textmark = "Other open" # The text that is shown in the heading if there is an open sub-question. ) { out <- dropall(data[2:(rowfact - 3) , ]) subquestion <- t(data[1 , ]) subquestion <- gsub("\xa0", " ", subquestion) subquestion <- gsub("\xb4", " ", subquestion) subquestion <- gsub("\n", " ", subquestion) # subquestion <- gsub("\\(", " ", subquestion) # subquestion <- gsub("\\)", " ", subquestion) textfield <- regexpr(textmark, subquestion) != -1 subquestion <- strsplit(subquestion, ":") # Divide the heading into a main question and a subquestion. subqtest <- 0 # The previous question name. for(i in 1:ncol(out)) { #print(i) if(subquestion[[i]][1] != subqtest) { # If part of previous question, use previous fact. fact <- as.character(data[rowfact:nrow(data) , i]) # Create factor levels from the end of Webropol file. fact <- fact[fact != ""] # Remove empty rows fact <- gsub("\xa0", " ", fact) fact <- gsub("\xb4", " ", fact) fact <- gsub("\n", " ", fact) fact <- strsplit(fact, " = ") # Separate value (level) and interpretation (label) } if(length(fact) != 0 & !textfield[i]) { # Do this only if the column is not a text type column. out[[i]] <- factor( out[[i]], levels = unlist(lapply(fact, function(x) x[1])), labels = unlist(lapply(fact, function(x) x[2])), ordered = TRUE ) } subqtest <- subquestion[[i]][1] } return(out) } # merge.questions takes a multiple checkbox question and merges that into a single factor. # First levels in levs have priority over others, if several levels apply to a row. merge.questions <- function( dat, # data.frame with questionnaire data cols, # list of vectors of column names or numbers to be merged into one level in the factor levs, # vector (with the same length as cols) of levels of factors into which questions are merged. name # text string for the name of the new factor column in the data. ) { for(i in length(cols):1) { temp <- FALSE for(j in rev(cols[[i]])) { temp <- temp | !is.na(dat[[j]]) } dat[[name]][temp] <- levs[i] } dat[[name]] <- factor(dat[[name]], levels = levs, ordered = TRUE) return(dat) } ############# Data preprocessing # Get the data either from Opasnet or your own hard drive. #Survey original file: N:/Ymal/Projects/Goherr/WP5/Goherr_fish_consumption.csv survey <- opasnet.csv("5/57/Goherr_fish_consumption.csv", sep = ";", fill = TRUE, quote = "\"") surcol <- t(survey[1,]) survey <- webropol.convert(survey, 2121, textmark = ":Other open") # Data file is converted to data.frame using levels at row 1269. # Take the relevant columnames from the table on the page. colframe <- opbase.data("Op_en7749", subset = "Questions in the Goherr questionnaire") colnames(survey) <- colframe$Result survey$Row <- 1:nrow(survey) survey$Weighting <- as.double(levels(survey$Weighting))[survey$Weighting] survey$Ages <- ifelse(as.numeric(as.character(survey$Age)) < 46, "18-45",">45") freqlist <- c( "less than once a year", "A few times a year", "1 - 3 times per month", "once a week", "2 - 4 times per week", "5 or more times per week" ) amlist <- c( "1/6 plate or below (50 grams)", "1/3 plate (100 grams)", "1/2 plate (150 grams)", "2/3 plate (200 grams)", "5/6 plate (250 grams)", "full plate (300 grams)", "overly full plate (500 grams)" # "Not able to estimate" ) sidel <- c( "1/6 plate or below (50 grams)", "1/4 plate (100 grams)", "1/2 plate (150 grams)", "2/3 plate (200 grams)", "5/6 plate (250 grams)" # "Not able to estimate" ) ansl <- list( `How often eat fish` = freqlist, `How often Baltic salmon` = c("Never", freqlist), `How much Baltic salmon` = amlist, `How often side Baltic salmon` = c("Never", freqlist), `How much side Baltic salmon` = sidel, `How often Baltic herring` = c("Never", freqlist), `How much Baltic herring` = amlist, `How often side Baltic herring` = c("Never", freqlist), `How much side Baltic herring` = sidel ) for (i in names(ansl)) { survey[[i]] <- factor(survey[[i]], levels = ansl[[i]]) } objects.store(survey, surcol) cat("objects survey, surcol were stored.\n")

Bayes model

• Model run 3.3.2017. All variables assumed independent. [1]
• Model run 3.3.2017. p has more dimensions. [2]

+ Show code- Hide code # This is code Op_en7749/bayes on page [[Goherr: Fish consumption study]] library(OpasnetUtils) library(reshape2) library(rjags) library(car) # Fish intake in humans # Data from data.frame survey from page [[Goherr: Fish consumption study]] # Start with salmon questions 46:49 (amounts eaten) # Some preprocessing should be moved away from this code. # Model assumes that questionnaire data follows binomial distribution with parameter p and max value levels. # p depends on country (4 groups), sex (2), age (2) that are known individually. # We should test whether sex makes a difference or whether it can be omitted. # Predicted p for each question, determined by the group, are produced from the bayes model. # p is used to sample answers, which are then combined with estimates of amounts related to each answer. # Total amount eaten of each fish per modelled individual is finally calculated. objects.latest("Op_en7749", "preprocess") # [[Goherr: Fish consumption study]]: survey, surcol agel <- as.character(unique(survey$Ages)) countryl <- sort(as.character(unique(survey$Country))) genderl <- sort(as.character(unique(survey$Gender))) fisl <- c("Salmon", "Herring") # Interesting fish eating questions surv <- survey[c(1,3,158,16,29,30,31,46:49,86,95:98)] colnames(surv) #[1] "Country" "Gender" #[3] "Ages" "Fish eating" #[5] "How often eat fish" "Salmon eating" #[7] "Baltic salmon" "How often Baltic salmon" #[9] "How much Baltic salmon" "How often side Baltic salmon" #[11] "How much side Baltic salmon" "Eat Baltic herring" #[13] "How often Baltic herring" "How much Baltic herring" #[15] "How often side Baltic herring" "How much side Baltic herring" surv$Ages <- match(surv$Ages, agel) # Not a factor, coerce to integer surv <- as.data.frame(lapply(surv, FUN = function(x) as.integer(x))) # Coerce to integers surv[is.na(surv[[12]]) | surv[[12]] == 3 , 12] <- 1 # Eat Baltic herring: I don't know --> No surv[4:16] <- lapply(surv[4:16], FUN = function(x) x-1) # Scale to start from zero # Row numbers for respondents that have eaten fish, Baltic salmon, and Baltic herring eatfish <- (1:nrow(surv))[surv[[4]] %in% 1] eatsalm <- (1:nrow(surv))[!(surv[[7]] %in% 0 | is.na(surv[[7]]))] eatherr <- (1:nrow(surv))[surv[[12]] %in% 1 & !is.na(rowSums(surv[13:16]))] # Oletetaan, että covarianssimatriisi on vakio kaikille maille ja sukupuolille yms # mutta keskiarvo on spesifi näille ja kysymykselle. qlen <- c(4,2,2,2,6,2,2,7,7,7,5,2,7,7,7,5) # Number of options in each question agel countryl genderl fisl ############ Version with both binomial and multivariate normal mod <- textConnection(" model{ for(i in 1:S) { for(j in 1:4) { surv[i,j] ~ dbin(phi(z[j]), qlen[j]) } } z[1:4] ~ dmnorm(mu[], Omega[,]) mu[1:4] ~ dmnorm(mu0[1:4], S2[1:4,1:4]) Omega[1:4,1:4] ~ dwish(S3[1:4,1:4],4) } ") jags <- jags.model( mod, data = list( surv = surv[eatherr,13:16], qlen = qlen[13:16], S = length(eatherr), mu0 = rep(0.5,4), S2 = diag(4)/100000, S3 = diag(4)/10000 ), n.chains = 4, n.adapt = 100 ) update(jags, 100) samps.j <- jags.samples(jags, c('mu','Omega', 'z'), 1000) j <- array( samps.j$mu, dim = c(4,1000,4), dimnames = list(Question = 1:4, Iter = 1:1000, Seed = 1:4) ) scatterplotMatrix(t(j[,,1])) samps.c <- coda.samples(jags, c('mu','Omega','z'), 1000) plot(samps.c) ############################## ############ Version with multivariate normal mod <- textConnection(" model{ for(i in 1:S) { surv[i,1:4] ~ dmnorm(mu[], Omega[,]) } mu[1:4] ~ dmnorm(mu0[1:4], S2[1:4,1:4]) Omega[1:4,1:4] ~ dwish(S3[1:4,1:4],4) } ") jags <- jags.model( mod, data = list( surv = surv[eatherr,13:16], S = length(eatherr), mu0 = rep(0.5,4), S2 = diag(4)/100000, S3 = diag(4)/10000 ), n.chains = 4, n.adapt = 100 ) update(jags, 100) samps.j <- jags.samples(jags, c('mu','Omega', 'z'), 1000) j <- array( samps.j$mu, dim = c(4,1000,4), dimnames = list(Question = 1:4, Iter = 1:1000, Seed = 1:4) ) scatterplotMatrix(t(j[,,1])) samps.c <- coda.samples(jags, c('mu','Omega','z'), 1000) plot(samps.c) ############################## ########### MODEL WITH BINOMIAL DISTRIBUTION ASSUMPTION # Interesting fish eating questions surv <- survey[c(1,3,4,16,29,30,31,46:49,86,95:98,158)] colnames(surv) consume <- data.matrix(survey[c( "Baltic salmon", "Eat Baltic herring" )]) consume <- 1 - (consume!=2 | is.na(consume)) # Yes=1, other=0 agel <- as.character(unique(survey$Ages)) countryl <- sort(as.character(unique(survey$Country))) genderl <- sort(as.character(unique(survey$Gender))) fisl <- c("Salmon", "Herring") ans.sal <- survey[ grep("Yes", survey[["Baltic salmon"]]) , c( "How often Baltic salmon", "How much Baltic salmon", "How often side Baltic salmon", "How much side Baltic salmon", "Ages", "Country", "Gender" )] ans.her <- survey[ grep("Yes", survey[["Eat Baltic herring"]]) , c( "How often Baltic herring", "How much Baltic herring", "How often side Baltic herring", "How much side Baltic herring", "Ages", "Country", "Gender" )] quesl <- c( "How often", "How much", "How often side", "How much side", "Ages", "Country", "Gender" ) colnames(ans.sal) <- quesl colnames(ans.her) <- quesl ans <- rbind(ans.sal, ans.her) head(ans) qlen <- c(7,7,7,5) # Number of options in each question questionl <- colnames(ans.sal) questionl.her <- colnames(ans.her) agel countryl genderl fisl mod <- textConnection(" model{ # Data1: Full questionnaire for(f in 1:2) { # Salmon or herring for(i in 1:U) { # Full questionnaire consume[i,f] ~ dbern(p[f,ageU[i],countryU[i],genderU[i]]) } } # Data2: Those responses that reported either Baltic herring or salmon for(q in 1:Q) { # Q = # of fish-specific questions for(i in 1:S) { # S = # respondents who eat Baltic salmon or herring ans[i,q] ~ dbin(p2[fis[i],age[i],country[i],gender[i],q], qlen[q]) } } # Priors (multidimensional) for(i in 1:A) { # Age groups for(j in 1:C) { # Countries for(k in 1:G) { #Genders for(f in 1:2) { # Fish species p[f,i,j,k] ~ dbeta(1,1) for(q in 1:Q) { # Q = # of fish-specific questions p2[f,i,j,k,q] ~ dbeta(1,1) ans.pred[f,i,j,k,q] ~ dbin(p2[f,i,j,k,q], qlen[q]) } } } } } } ") jags <- jags.model( mod, data = list( U = nrow(surv), Q = 4, S = nrow(ans), consume = consume, ans = data.matrix(ans[1:4]) - 1, # Zero is the first option qlen = qlen-1, # Subtract 1 because starts from 0. age = match(ans$Ages, agel), country = match(ans$Country, countryl), gender = match(ans$Gender, genderl), ageU = match(survey$Ages, agel), countryU = match(survey$Country, countryl), genderU = match(survey$Gender, genderl), fis = c(rep(1, nrow(ans.sal)), rep(2, nrow(ans.her))), A = length(agel), C = length(countryl), G = length(genderl) ), n.chains = 4, n.adapt = 100 ) update(jags, 100) samps <- jags.samples(jags, c('ans.pred','p','p2'), 1000) #samps.coda <- coda.samples(jags, c('mu', 'pcd.pred', 'ans.pred'), 1000) ans.pred <- array( samps$ans.pred, dim = c(length(fisl), length(agel), length(countryl), length(genderl), length(quesl), 1000, 4), dimnames = list( Fish = fisl, Age = agel, Country = countryl, Gender = genderl, Question = quesl, Iter = 1:1000, Seed = c("S1","S2","S3","S4") ) ) p.pred <- array( samps$p, dim = c(length(fisl), length(agel), length(countryl), length(genderl), 1000, 4), dimnames = list( Fish = fisl, Age = agel, Country = countryl, Gender = genderl, Iter = 1:1000, Seed = c("S1","S2","S3","S4") ) ) p2.pred <- array( samps$p2, dim = c(length(fisl), length(agel), length(countryl), length(genderl), length(quesl), 1000, 4), dimnames = list( Fish = fisl, Age = agel, Country = countryl, Gender = genderl, Question = quesl, Iter = 1:1000, Seed = c("S1","S2","S3","S4") ) ) ansm <- melt(ans.pred) pm <- melt(p.pred) p2m <- melt(p2.pred) scatterplotMatrix(t(ans.pred[1,1,,1,1,,1])) scatterplotMatrix(t(p2.pred[1,1,,1,1,,1])) scatterplotMatrix(t(p2.pred[,1,3,1,1,,1])) scatterplotMatrix(t(p2.pred[1,,3,1,1,,1])) scatterplotMatrix(t(p2.pred[1,1,3,,1,,1])) scatterplotMatrix(t(p2.pred[1,1,3,1,,,1])) scatterplotMatrix(t(p.pred[1,1,,1,,1])) scatterplotMatrix(t(p.pred[,1,3,1,,1])) scatterplotMatrix(t(p.pred[1,,3,1,,1])) scatterplotMatrix(t(p.pred[1,1,3,,,1])) objects.store(ans.pred, p.pred, p2.pred) cat("Arrays ans.pred (answer), p.pred (probability to eat), p2.pred (p of answer) stored.\n")

Calculations

This code calculates how much (g/day) Baltic herring and salmon are eaten based on an Bayesian model build up based on the questionnaire data.

+ Show code- Hide code # This is code Op_en7749/calculate_amount on page [[Goherr: Fish consumption study]] library(OpasnetUtils) library(ggplot2) library(reshape2) library(rjags) objects.latest("Op_en7749", code_name = "bayes") #: ans.pred, p.pred, p2.pred ql <- melt(ans.pred) ##Assumptions for amount and frequencies assumptions <- opbase.data("Op_en7749.assumptions_for_calculations") salmon_dishtimesassump <- Ovariable("salmon_dishtimesassump", data = assumptions[assumptions$Variable == "freq" , c("value", "Result")]) salmon_dishamountassump <- Ovariable("salmon_dishamountassump", data = assumptions[assumptions$Variable == "amdish" , c("value", "Result")]) salmon_ingridient <- Ovariable("salmon_ingridient", data = assumptions[assumptions$Variable == "ingridient", ]["Result"]) salmon_sidetimesassump <- Ovariable("salmon_wholetimesassump", data = assumptions[assumptions$Variable == "freq" , c("value", "Result")]) salmon_sideamountassump <- Ovariable("salmon_wholeamountassump", data = assumptions[assumptions$Variable == "amside" , c("value", "Result")]) ##Build up the variables to calculate amount of Baltic salmon eating salmon_dishtimes <- Ovariable("salmon_dishtimes", data = data.frame( subset(ql, ql$Fish == "Salmon" & ql$Question == "How often" & ql$Seed == "S1"), Result = 1 )) salmon_dishtimes@data$Question <- NULL salmon_dishamount <- Ovariable("salmon_dishamount", data = data.frame( subset(ql, ql$Fish == "Salmon" & ql$Question == "How much" & ql$Seed == "S1"), Result = 1 )) #salmon with dish amount per one serving salmon_dishamount@data$Question <- NULL salmon_sidetimes <- Ovariable("salmon_sidetimes", data = data.frame( subset(ql, ql$Fish == "Salmon" & ql$Question == "How often side" & ql$Seed == "S1"), Result = 1 )) #salmon as such times / year salmon_sidetimes@data$Question <- NULL salmon_sideamount <- Ovariable("salmon_sideamount", data = data.frame( subset(ql, ql$Fish == "Salmon" & ql$Question == "How much side" & ql$Seed == "S1"), Result = 1 )) #salmon as whole amount per one serving salmon_sideamount@data$Question <- NULL amount <- Ovariable("amount", dependencies = data.frame(Name = c( "salmon_dishamount", "salmon_dishamountassump", "salmon_ingridient", "salmon_dishtimes", "salmon_dishtimesassump", "salmon_sideamount", "salmon_sideamountassump", "salmon_sidetimes", "salmon_sidetimesassump" )), formula = function(...) { out1 <- (salmon_dishamount * salmon_dishamountassump * salmon_ingridient) out1@output$value <- NULL out1 <- ((salmon_dishtimes * salmon_dishtimesassump) * out1) / 365 out1@output$value <- NULL out2 <- (salmon_sideamount * salmon_sideamountassump) out2@output$value <- NULL out2 <- ((salmon_sidetimes * salmon_sidetimesassump) * out2) / 365 out2@output$value <- NULL out <- out1 + out2 # Per vuosi -> per d #result(out)[result(out) == 0] <- 0.01 # Ei jätetä nollia saantiin return(out) } ) amountSalmon <- EvalOutput(amount, forceEval = TRUE, N = 1000) temp <- amountSalmon@output temp <- aggregate(temp$amountResult, by=list(temp$Country, temp$Age, temp$Gender), FUN = mean) colnames(temp) <- c("Country","Age","Gender","Salmon amount") oprint(temp) ##Build up the variables to calculate amount of Baltic herring eating herring_dishtimes <- Ovariable("herring_dishtimes", data = data.frame( subset(ql, ql$Fish == "Herring" & ql$Question == "How often" & ql$Seed == "S1"), Result = 1 )) herring_dishtimes@data$Question <- NULL herring_dishamount <- Ovariable("herring_dishamount", data = data.frame( subset(ql, ql$Fish == "Herring" & ql$Question == "How much" & ql$Seed == "S1"), Result = 1 )) #herring with dish amount per one serving herring_dishamount@data$Question <- NULL herring_sidetimes <- Ovariable("herring_sidetimes", data = data.frame( subset(ql, ql$Fish == "Herring" & ql$Question == "How often side" & ql$Seed == "S1"), Result = 1 )) #herring as such times / year herring_sidetimes@data$Question <- NULL herring_sideamount <- Ovariable("herring_sideamount", data = data.frame( subset(ql, ql$Fish == "Herring" & ql$Question == "How much side" & ql$Seed == "S1"), Result = 1 )) #herring as whole amount per one serving herring_sideamount@data$Question <- NULL herring_dishtimesassump <- Ovariable("herring_dishtimesassump", data = assumptions[assumptions$Variable == "freq" , c("value", "Result")]) herring_dishamountassump <- Ovariable("herring_dishamountassump", data = assumptions[assumptions$Variable == "amdish" , c("value", "Result")]) herring_ingridient <- Ovariable("herring_ingridient", data = assumptions[assumptions$Variable == "ingridient", ]["Result"]) herring_sidetimesassump <- Ovariable("herring_wholetimesassump", data = assumptions[assumptions$Variable == "freq" , c("value", "Result")]) herring_sideamountassump <- Ovariable("herring_wholeamountassump", data = assumptions[assumptions$Variable == "amside" , c("value", "Result")]) amount <- Ovariable("amount", dependencies = data.frame(Name = c( "herring_dishamount", "herring_dishamountassump", "herring_ingridient", "herring_dishtimes", "herring_dishtimesassump", "herring_sideamount", "herring_sideamountassump", "herring_sidetimes", "herring_sidetimesassump" )), formula = function(...) { out1 <- (herring_dishamount * herring_dishamountassump * herring_ingridient) out1@output$value <- NULL out1 <- ((herring_dishtimes * herring_dishtimesassump) * out1) / 365 out1@output$value <- NULL out2 <- (herring_sideamount * herring_sideamountassump) out2@output$value <- NULL out2 <- ((herring_sidetimes * herring_sidetimesassump) * out2) / 365 out2@output$value <- NULL out <- out1 + out2 # Per vuosi -> per d #result(out)[result(out) == 0] <- 0.01 # Ei jätetä nollia saantiin return(out) } ) amountHerring <- EvalOutput(amount, forceEval = TRUE, N = 1000) temp <- amountHerring@output temp <- aggregate(temp$amountResult, by=list(temp$Country, temp$Age, temp$Gender), FUN = mean) colnames(temp) <- c("Country","Age","Gender","Herring amount") oprint(temp) objects.store(amountHerring, amountSalmon) cat("objects amountHerring, amountSalmon were stored.\n")

Assumptions

The following assumptions are used:

Assumptions for calculations(-)

Obs	Variable	value	Explanation	Result
1	freq	6	times per year	260 - 364
2	freq	5	times per year	104 - 208
3	freq	4	times per year	52
4	freq	3	times per year	12 - 36
5	freq	2	times per year	2 - 5
6	freq	1	times per year	0.5 - 0.9
7	freq	0	times per year	0
8	amdish	0	grams / serving	20 - 50
9	amdish	1	grams / serving	70 - 100
10	amdish	2	grams / serving	120 - 150
11	amdish	3	grams / serving	170 - 200
12	amdish	4	grams / serving	220 - 250
13	amdish	5	grams / serving	270 - 300
14	amdish	6	grams / serving	450 - 500
15	ingridient		fraction	0.1 - 0.3
16	amside	0	grams / serving	20 - 50
17	amside	1	grams / serving	70 - 100
18	amside	2	grams / serving	120 - 150
19	amside	3	grams / serving	170 - 200
20	amside	4	grams / serving	220 - 250

Questionnaire

Show details

Questions in the Goherr questionnaire(-) Obs Index Full question Question 1 V1 Country Country 2 V2 1. What is the postal code of your residence? Postal code 3 V3 2. What is your gender? Gender 4 V4 2. How old are you? Age 5 V5 3. What is the level of your highest education? Education 6 V6 4. How would you rate the purchasing power of your household? Purchasing power 7 V7 5. Which of the following food supplements (as a separate product or as a part of some multicompound) do you eat often at least part of the year?:Fatty acids (Omega 3, DHA, EPA etc.) Omega3 8 V8 5. Which of the following food supplements (as a separate product or as a part of some multicompound) do you eat often at least part of the year?:D vitamin Vitamin D 9 V9 5. Which of the following food supplements (as a separate product or as a part of some multicompound) do you eat often at least part of the year?:Other vitamins Other vitamins 10 V10 5. Which of the following food supplements (as a separate product or as a part of some multicompound) do you eat often at least part of the year?:None of the above None above 11 V11 5. Which of the following food supplements (as a separate product or as a part of some multicompound) do you eat often at least part of the year?:Don't Know Don't know 12 V12 6. Does someone in your household fish?:No No fishing 13 V13 6. Does someone in your household fish?:Yes, for recreational purposes Recreational 14 V14 6. Does someone in your household fish?:Yes, for household use Household use 15 V15 6. Does someone in your household fish?:Yes, as a professional fisherman Professional 16 V16 7. Do you eat fish? Fish eating 17 V17 8. Why don't you eat fish?:Allergy Allergy 18 V18 8. Why don't you eat fish?:I don't like the taste Don't like taste 19 V19 8. Why don't you eat fish?:I don't like fish bones Don't like bones 20 V20 8. Why don't you eat fish?:Difficult to cook Cooking difficult 21 V21 8. Why don't you eat fish?Health risks Risky to health 22 V22 8. Why don't you eat fish?:I'm not used to Not used to 23 V23 8. Why don't you eat fish?:I'm concerned about the sustainability of fish stocks Sustainability of stock 24 V24 8. Why don't you eat fish?:It is ethically wrong to eat animals Ethically wrong 25 V25 8. Why don't you eat fish?:I'm a vegetarian Vegetarian 26 V26 8. Why don't you eat fish?:Don't Know Don't know 27 V27 8. Why don't you eat fish?:Other Other 28 V28 8. Why don't you eat fish?:Other open Open 29 V29 9. How often do you eat fish? How often eat fish 30 V30 10. Do you eat salmon or trout at least sometimes? Salmon eating 31 V31 11. Which of the following salmonids do you eat at least sometimes?:Wild Baltic salmon Baltic salmon 32 V32 11. Which of the following salmonids do you eat at least sometimes?:Wild Atlantic salmon Atlantic salmon 33 V33 11. Which of the following salmonids do you eat at least sometimes?:Norwegian cultured salmon Norwegian salmon 34 V34 11. Which of the following salmonids do you eat at least sometimes?:Rainbow trout Rainbow trout 35 V35 11. Which of the following salmonids do you eat at least sometimes?:I never know which salmon I eat Unknown salmon 36 V36 11. Which of the following salmonids do you eat at least sometimes?:Some other salmon, what? Other salmon 37 V37 11. Which of the following salmonids do you eat at least sometimes?:Other open Open 38 V38 12. Where do you usually get the wild Baltic salmon that you eat?:Someone in our hosehold catches it Own catch 39 V39 12. Where do you usually get the wild Baltic salmon that you eat?:Buy from a grocery store Grocery store 40 V40 12. Where do you usually get the wild Baltic salmon that you eat?:Buy from a market place/fishmonger or similar Market place 41 V41 12. Where do you usually get the wild Baltic salmon that you eat?:Eat in a restaurant/cafe/canteen or similar Restaurant 42 V42 12. Where do you usually get the wild Baltic salmon that you eat?:Directly from a fisher From fisher 43 V43 12. Where do you usually get the wild Baltic salmon that you eat?:From a relative, friend, or similar From friend 44 V44 12. Where do you usually get the wild Baltic salmon that you eat?:Other Other 45 V45 12. Where do you usually get the wild Baltic salmon that you eat?:Other open Open 46 V46 13. How often do you on average eat wild Baltic salmon as an ingredient in a dish (i.e. as part of a soup, gravy, sushi, casserole etc.)? How often Baltic salmon 47 V47 14. How much on average do you eat of a dish with wild Baltic salmon in it? Think of a medium sized lunch plate and the amount of the whole dish, not just the salmon. How much Baltic salmon 48 V48 15. How often do you on average eat wild Baltic salmon as a side dish (i.e as a separate product such as fish balls, stake, fillets, graved, canned etc.)? How often side Baltic salmon 49 V49 16. What is the average amount of wild Baltic salmon as a side dish? Think of a medium sized lunch plate. How much side Baltic salmon 50 V50 17. Select up to three most important reasons for you to eat wild Baltic salmon?:It tastes good Tastes good 51 V51 17. Select up to three most important reasons for you to eat wild Baltic salmon?:I caught it I caught it 52 V52 17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's easy to cook Easy to cook 53 V53 17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's quick to cook Quick to cook 54 V54 17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's readily available Readily available 55 V55 17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's healthy Healthy 56 V56 17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's inexpensive Inexpensive 57 V57 17. Select up to three most important reasons for you to eat wild Baltic salmon?:My family likes it Family likes it 58 V58 17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's environmentally/climate friendly choice Environmentally friendly 59 V59 17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's a traditional dish Traditional dish 60 V60 17. Select up to three most important reasons for you to eat wild Baltic salmon?:Some other reasons, what? Other 61 V61 17. Select up to three most important reasons for you to eat wild Baltic salmon?:Other open Open 62 V62 18. You answered in the previous question that you don't eat wild Baltic salmon. Select up to three most important reasons for why not eating it.:I don't like the taste Don't like taste 63 V63 18. You answered in the previous question that you don't eat wild Baltic salmon. Select up to three most important reasons for why not eating it.:I am not used to eating it / not offered at home Not used to 64 V64 18. You answered in the previous question that you don't eat wild Baltic salmon. Select up to three most important reasons for why not eating it.:I don't know how to cook it Don't know how to cook 65 V65 18. You answered in the previous question that you don't eat wild Baltic salmon. Select up to three most important reasons for why not eating it.:It is not easily available Not available 66 V66 18. You answered in the previous question that you don't eat wild Baltic salmon. Select up to three most important reasons for why not eating it.:I'm worried about the possible health risks caused by the harmful chemicals in it Health risks 67 V67 18. You answered in the previous question that you don't eat wild Baltic salmon. Select up to three most important reasons for why not eating it.:I doubt the quality/freshness of the products Quality issues 68 V68 18. You answered in the previous question that you don't eat wild Baltic salmon. Select up to three most important reasons for why not eating it.:I'm worried about the sustainability of the stocks Sustainability of stock 69 V69 18. You answered in the previous question that you don't eat wild Baltic salmon. Select up to three most important reasons for why not eating it.:It doesn't belong to traditional dishes of my home area Not traditional 70 V70 18. You answered in the previous question that you don't eat wild Baltic salmon. Select up to three most important reasons for why not eating it.:It is too expensive Expensive 71 V71 18. You answered in the previous question that you don't eat wild Baltic salmon. Select up to three most important reasons for why not eating it.:Some other reason, what? Other 72 V72 18. You answered in the previous question that you don't eat wild Baltic salmon. Select up to three most important reasons for why not eating it.:Other open Open 73 V73 19. How would the following factors influence on your consumption of wild Baltic salmon compared to your current use? Select one for each option:Lower price Lower price 74 V74 19. How would the following factors influence on your consumption of wild Baltic salmon compared to your current use? Select one for each option:Higher price Higher price 75 V75 19. How would the following factors influence on your consumption of wild Baltic salmon compared to your current use? Select one for each option:Better availability Better availability 76 V76 19. How would the following factors influence on your consumption of wild Baltic salmon compared to your current use? Select one for each option:Availability of wider variety of ready meals / refined products Ready meals 77 V77 19. How would the following factors influence on your consumption of wild Baltic salmon compared to your current use? Select one for each option:Cooking suggestions and recipes provided within the store or in the package. Recipes 78 V78 19. How would the following factors influence on your consumption of wild Baltic salmon compared to your current use? Select one for each option:Better information (e.g. package markings or in the store) on the catch date, area, fisher and/or the refiner Package markings 79 V79 19. How would the following factors influence on your consumption of wild Baltic salmon compared to your current use? Select one for each option:Improvement of the stocks (eco-labelling) Eco-labelling 80 V80 19. How would the following factors influence on your consumption of wild Baltic salmon compared to your current use? Select one for each option:Lower level of harmful chemicals in the fish Less chemicals 81 V81 19. How would the following factors influence on your consumption of wild Baltic salmon compared to your current use? Select one for each option:National food safety authorities publish a recommendation to eat it Recommendation 82 V82 19. How would the following factors influence on your consumption of wild Baltic salmon compared to your current use? Select one for each option:National food safety authorities publish a recommendation to limit eating it Limit recommendation 83 V83 19. How would the following factors influence on your consumption of wild Baltic salmon compared to your current use? Select one for each option:Some other reason, what? Other 84 V84 19. How would the following factors influence on your consumption of wild Baltic salmon compared to your current use? Select one for each option:Other open Open 85 V85 20. Do you eat any herring (Baltic or other) at least sometimes? Eat herring 86 V86 21. Do you eat Baltic herring? This means herring caught from Baltic sea, not from Atlantic ocean Eat Baltic herring 87 V87 22. Where do you usually get the wild Baltic herring that you eat?:Someone in our hosehold catches it Own catch 88 V88 22. Where do you usually get the wild Baltic herring that you eat?:Buy from a grocery store Grocery store 89 V89 22. Where do you usually get the wild Baltic herring that you eat?:Buy from a market place/fishmonger or similar Market place 90 V90 22. Where do you usually get the wild Baltic herring that you eat?:Eat in a restaurant/cafe/canteen or similar Restaurant 91 V91 22. Where do you usually get the wild Baltic herring that you eat?:Directly from a fisher From fisher 92 V92 22. Where do you usually get the wild Baltic herring that you eat?:From a relative, friend, or similar From friend 93 V93 22. Where do you usually get the wild Baltic herring that you eat?:Other Other 94 V94 22. Where do you usually get the wild Baltic herring that you eat?:Other open Open 95 V95 23. How often do you on average eat Baltic herring as ingredient of a dish? (i.e as part of a dish such as in a casserole or some other food including Baltic herring)? How often Baltic herring 96 V96 24. How much on average do you eat of a dish with Baltic herring in it? Think of a medium sized lunch plate and the amount of the whole dish, not just the amount of Baltic herring How much Baltic herring 97 V97 25. How often do you on average eat whole Baltic herrings or fillets as a side dish (fried, smoked, rolls, pickled, salted or otherwise prepared)? How often side Baltic herring 98 V98 26. How much on average of Baltic herring as a side dish (think of a regular sized lunch plate)? How much side Baltic herring 99 V99 27. Select up to three most important reasons for you to eat Baltic herring?:It tastes good Tastes good 100 V100 27. Select up to three most important reasons for you to eat Baltic herring?:I caught it I caught it 101 V101 27. Select up to three most important reasons for you to eat Baltic herring?:It's easy to cook Easy to cook 102 V102 27. Select up to three most important reasons for you to eat Baltic herring?:It's quick to cook Quick to cook 103 V103 27. Select up to three most important reasons for you to eat Baltic herring?:It's readily available Readily available 104 V104 27. Select up to three most important reasons for you to eat Baltic herring?:It's healthy Healthy 105 V105 27. Select up to three most important reasons for you to eat Baltic herring?:It's inexpensive Inexpensive 106 V106 27. Select up to three most important reasons for you to eat Baltic herring?:My family likes it Family likes it 107 V107 27. Select up to three most important reasons for you to eat Baltic herring?:It's environmentally/climate friendly choice Environmentally friendly 108 V108 27. Select up to three most important reasons for you to eat Baltic herring?:It's a traditional dish Traditional dish 109 V109 27. Select up to three most important reasons for you to eat Baltic herring?:Some other reasons, what? Other 110 V110 27. Select up to three most important reasons for you to eat Baltic herring?:Other open Open 111 V111 28. You answered in the previous question that you don't eat Baltic herring. Select up to three most important reasons for why not eating it.:I don't like the taste and/or smell of it Don't like taste 112 V112 28. You answered in the previous question that you don't eat Baltic herring. Select up to three most important reasons for why not eating it.:I am not used to eating it / not offered at home Not used to 113 V113 28. You answered in the previous question that you don't eat Baltic herring. Select up to three most important reasons for why not eating it.:I don't know how to cook it Don't know how to cook 114 V114 28. You answered in the previous question that you don't eat Baltic herring. Select up to three most important reasons for why not eating it.:It is difficult to cook Difficult to cook 115 V115 28. You answered in the previous question that you don't eat Baltic herring. Select up to three most important reasons for why not eating it.:It is not easily available Not available 116 V116 28. You answered in the previous question that you don't eat Baltic herring. Select up to three most important reasons for why not eating it.:I'm worried about the possible risks caused by the harmful chemicals in it Health risks 117 V117 28. You answered in the previous question that you don't eat Baltic herring. Select up to three most important reasons for why not eating it.:It is more suitable for animal feed than human food Better for feed 118 V118 28. You answered in the previous question that you don't eat Baltic herring. Select up to three most important reasons for why not eating it.:I doubt the quality/freshness of the products Quality issues 119 V119 28. You answered in the previous question that you don't eat Baltic herring. Select up to three most important reasons for why not eating it.:I'm worried about the sustainability of the stocks Sustainability of stock 120 V120 28. You answered in the previous question that you don't eat Baltic herring. Select up to three most important reasons for why not eating it.:It doesn't belong to traditional dishes of my home area Not traditional 121 V121 28. You answered in the previous question that you don't eat Baltic herring. Select up to three most important reasons for why not eating it.:Some other reason, what? Other 122 V122 28. You answered in the previous question that you don't eat Baltic herring. Select up to three most important reasons for why not eating it.:Other open Open 123 V123 29. How would the following factors influence your consumption of Baltic herring compared to your current use? Select one option for each row.:Lower price Lower price 124 V124 29. How would the following factors influence your consumption of Baltic herring compared to your current use? Select one option for each row.:Higher price Higher price 125 V125 29. How would the following factors influence your consumption of Baltic herring compared to your current use? Select one option for each row.:Better availability Better availability 126 V126 29. How would the following factors influence your consumption of Baltic herring compared to your current use? Select one option for each row.:Availability of wider variety of ready meals / processed products Ready meals 127 V127 29. How would the following factors influence your consumption of Baltic herring compared to your current use? Select one option for each row.:Cooking suggestions and recipes provided in the store or in the package. Recipes 128 V128 29. How would the following factors influence your consumption of Baltic herring compared to your current use? Select one option for each row.:Better information (e.g. package markings or in the store) on the catch date, area, fisher and/or processing company. Package markings 129 V129 29. How would the following factors influence your consumption of Baltic herring compared to your current use? Select one option for each row.:Improved stocks (eco-labelling) Eco-labelling 130 V130 29. How would the following factors influence your consumption of Baltic herring compared to your current use? Select one option for each row.:Lower level of harmful chemicals in the fish Less chemicals 131 V131 29. How would the following factors influence your consumption of Baltic herring compared to your current use? Select one option for each row.:National food safety authorities publish a recommendation to eat it Recommendation 132 V132 29. How would the following factors influence your consumption of Baltic herring compared to your current use? Select one option for each row.:National food safety authorities publish a recommendation to limit eating it Limit recommendation 133 V133 29. How would the following factors influence your consumption of Baltic herring compared to your current use? Select one option for each row.:Some other reason, what? Other 134 V134 29. How would the following factors influence your consumption of Baltic herring compared to your current use? Select one option for each row.:Other open Open 135 V135 30. Are you familiar with the possible eating recommendations of Baltic herring and/or wild Baltic salmon in your country? Recommendation awareness 136 V136 31. How strongly do you agree or disagree with the following statements on Baltic salmon? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It contributes to the economic wellbeing of the society Economic wellbeing 137 V137 31. How strongly do you agree or disagree with the following statements on Baltic salmon? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It is a natural resource for humans to utilize Natural resource 138 V138 31. How strongly do you agree or disagree with the following statements on Baltic salmon? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It belongs to all Baltic Sea citizens equally Belongs to Baltic Sea 139 V139 31. How strongly do you agree or disagree with the following statements on Baltic salmon? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It involves long traditions related to fishing and fish eating Long traditions 140 V140 31. How strongly do you agree or disagree with the following statements on Baltic salmon? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It inspires to search experiences from nature, cultural sites, fish markets or from different dishes Inspiring experiences 141 V141 31. How strongly do you agree or disagree with the following statements on Baltic salmon? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It is among the most renowned fish species of the Baltic Sea Renowned species 142 V142 31. How strongly do you agree or disagree with the following statements on Baltic salmon? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It is part of the Baltic Sea ecosystem Part of ecosystem 143 V143 31. How strongly do you agree or disagree with the following statements on Baltic salmon? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:I have an emotional bond to it Emotional bond 144 V144 31. How strongly do you agree or disagree with the following statements on Baltic salmon? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It helps me to achieve my personal aims Personal aims 145 V145 32. How strongly do you agree or disagree with the following statements on Baltic herring? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It contributes to the economic wellbeing of the society Economic wellbeing 146 V146 32. How strongly do you agree or disagree with the following statements on Baltic herring? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It is a natural resource for humans to utilize Natural resource 147 V147 32. How strongly do you agree or disagree with the following statements on Baltic herring? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It belongs to all Baltic Sea citizens equally Belongs to Baltic Sea 148 V148 32. How strongly do you agree or disagree with the following statements on Baltic herring? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It involves long traditions related to fishing and fish eating Long traditions 149 V149 32. How strongly do you agree or disagree with the following statements on Baltic herring? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It inspires to search experiences from nature, cultural sites, fish markets or from different dishes Inspiring experiences 150 V150 32. How strongly do you agree or disagree with the following statements on Baltic herring? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It is among the most renowned fish species of the Baltic Sea Renowned species 151 V151 32. How strongly do you agree or disagree with the following statements on Baltic herring? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It is part of the Baltic Sea ecosystem Part of ecosystem 152 V152 32. How strongly do you agree or disagree with the following statements on Baltic herring? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:I have an emotional bond to it Emotional bond 153 V153 32. How strongly do you agree or disagree with the following statements on Baltic herring? Please indicate your position by selecting a number from -5 (strongly disagree) to 5 (strongly agree).:It helps me to achieve my personal aims Personal aims 154 V154 Weighting Weighting 155 V155 How old are you? Age category FI Age category FI 156 V156 2. How old are you? Age category SWE, DK, EST Age category other

Dependencies

The survey data will be used as input in the benefit-risk assessment of Baltic herring and salmon intake, which is part of the WP5 work in Goherr-project.

Related files

<mfanonymousfilelist></mfanonymousfilelist>

Goherr: Fish consumption study. Opasnet . [3]. Accessed 15 May 2025.