Goherr: Fish consumption study

From Opasnet
Jump to: navigation, search


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?


Original questionnaire analysis results

  • 13.3.2017 #: . These should be presented somewhere --Arja (talk) 07:39, 26 April 2017 (UTC) (type: truth; paradigms: science: relevant comment)

Consumption amount estimates

Do you want to use directly the survey data rather than modelled data?:

+ Show code


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 can be found from the link below (see Data).



Original datafile File:Goherr fish consumption.csv.


The following assumptions are used:

Assumptions for calculations(-)
1freq1times /a0Never
2freq2times /a0.5 - 0.9less than once a year
3freq3times /a2 - 5A few times a year
4freq4times /a12 - 361 - 3 times per month
5freq5times /a52once a week
6freq6times /a104 - 2082 - 4 times per week
7freq7times /a260 - 3645 or more times per week
8amdish1g /serving20 - 701/6 plate or below (50 grams)
9amdish2g /serving70 - 1301/3 plate (100 grams)
10amdish3g /serving120 - 1801/2 plate (150 grams)
11amdish4g /serving170 - 2302/3 plate (200 grams)
12amdish5g /serving220 - 2805/6 plate (250 grams)
13amdish6g /serving270 - 400full plate (300 grams)
14amdish7g /serving400 - 550overly full plate (500 grams)
15ingredientfraction0.1 - 0.3Fraction of fish in the dish
16amside1g /serving20 - 701/6 plate or below (50 grams)
17amside2g /serving70 - 1301/4 plate (100 grams)
18amside3g /serving120 - 1801/2 plate (150 grams)
19amside4g /serving170 - 2302/3 plate (200 grams)
20amside5g /serving220 - 2805/6 plate (250 grams)
21change1fraction-1 - -0.8Decrease it to zero
22change2fraction-0.9 - -0.5Decrease it to less than half
23change3fraction-0.6 - -0.1Decrease it a bit
24change4fraction0No effect
25change5fraction0.1 - 0.6Increase it a bit
26change6fraction0.5 - 0.9Increase it over by half
27change7fraction0.8 - 1.3Increase it over to double
28change8fraction-0.3 - 0.3Don't know


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. The code stores a data.frame named survey.

  • Model run 11.7.2018 [2]

+ Show code


  • Sketches about modelling determinants of eating (spring 2018) [3]

Figures, tables and stat analyses for the first manuscript

  • Model run 15.6.2018 [4]
  • Model run 16.6.2018 [5]
  • Model run 20.6.2018 [6]
  • Model run 11.7.2018 [7]
  • Model run 9.8.2018 [8]
  • Model run 12.8.2018 [9]
  • Model run 12.9.2018 with country-specific regression analyses, open-ended answers and fig about eating any herring [10]
  • Model run 5.2.2019 [11]
  • Model run 6.2.2019 [12]
  • Model run 7.2.2019 with systematic analyses of reasons to eat and not to eat Baltic herring [13]
  • Model run 20.3.2019 with country-specific regression analyses [14]: printregr(countries=TRUE)
  • Model run 22.3.2019 with better colour patterns [15]

Methodological concerns:

+ Show code

Luke data about fish consumption in Finland [16][17]

Fish consumption as food in Finland(kg/a per person)
1domestic fishTotal6.
2domestic fishFarmed rainbow trout1.
3domestic fishBaltic herring0.
4domestic fishPike0.
5domestic fishPerch0.
6domestic fishVendace0.
7domestic fishEuropean whitefish0.
8domestic fishPike perch0.
9domestic fishOther domestic fish0.
10imported fishTotal6.
11imported fishFarmed rainbow trout0.
12imported fishFarmed salmon1.
13imported fishTuna (prepared and preserved)
14imported fishSaithe (frozen fillet)
15imported fishShrimps0.
16imported fishHerring and Baltic herring (preserved)
17imported fishOther imported fish2.

Descriptive statistics

Correlation matrix of all questions in the survey (answers converted to numbers).

Model must contain predictors such as country, gender, age etc. Maybe we should first study what determinants are important? Model must also contain determinants that would increase or decrease fish consumption. This should be conditional on the current consumption. How? Maybe we should look at principal coordinates analysis with all questions to see how they behave.

Also look at correlation table to see clusters.

Some obvious results:

  • If reports no fish eating, many subsequent answers are NA.
  • No vitamins correlates negatively with vitamin intake.
  • Unknown salmon correlates negatively with the types of salmon eaten.
  • Different age categories correlate with each other.

However, there are also meaningful negative correlations:

  • Country vs allergy
  • Country vs Norwegian salmon and Rainbow trout
  • Country vs not traditional.
  • Country vs recommendation awareness
  • Allergy vs economic wellbeing
  • Baltic salmon use (4 questions) vs Don't like taste and Not used to
  • All questions between Easy to cook ... Traditional dish

Meaningful positive correlations:

  • All questions between Baltic salmon ... Rainbow trout
  • How often Baltic salmon/herring/side salmon/side herring
  • How much Baltic salmon/herring/side salmon/side herring
  • Better availability ... Recommendation
  • All questions between Economic wellbeing...Personal aims
  • Omega3, Vitamin D, and Other vitamins

Model runs

+ Show code

Bayes model

  • Model run 3.3.2017. All variables assumed independent. [19]
  • Model run 3.3.2017. p has more dimensions. [20]
  • Model run 25.3.2017. Several model versions: strange binomial+multivarnormal, binomial, fractalised multivarnormal [21]
  • Model run 27.3.2017 [22]
  • Other models except multivariate normal were archived and removed from active code 29.3.2017.
  • Model run 29.3.2017 with raw data graphs [23]
  • Model run 29.3.2017 with salmon and herring ovariables stored [24]
  • Model run 13.4.2017 with first version of coordinate matrix and principal coordinate analysis [25]
  • Model run 20.4.2017 [26] code works but needs a safety check against outliers
  • Model run 21.4.2017 [27] some model results plotted
  • Model run 21.4.2017 [28] ovariables produced by the model stored.
  • Model run 18.5.2017 [29] small updates
  • 13.2.2018 old model run but with new Opasnet [30]

+ Show code

Initiate ovariables

jsp taken directly from data WITHOUT salmpling.

+ Show code

Amount estimated from a bayesian model.

  • Model run 24.5.2017 [31]

+ Show code

Amount estimates directly from data rather than from a bayesian model.

  • Initiation run 18.5.2017 [32]
  • Initiation run 24.2.2018: sampling from survey rather than each respondent once [33]

+ Show code

Initiate other ovariables

  • Code stores ovariables assump, often, much, oftenside, muchside, amount.
  • Model run 19.5.2017 [34]
  • Initiation run 24.5.2017 without jsp [35]
  • Model run 8.6.2017 [36]

+ Show code

Other code

This is code for analysing EFSA food intake data about fish for BONUS GOHERR manuscript Ronkainen L, Lehikoinen A, Haapasaari P, Tuomisto JT. 2019.

+ Show code


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.

See also



Related files

Goherr: Fish consumption study. Opasnet . [37]. Accessed 23 Mar 2019.