Difference between revisions of "Goherr: Fish consumption study"

From Testiwiki
Jump to: navigation, search
(Bayes model)
(technical edits)
Line 5: Line 5:
  
 
== Answer ==
 
== Answer ==
Survey data will be analysed during winter 2016-2017 and results will be updated here.
 
  
* [http://en.opasnet.org/en-opwiki/index.php?title=Special:RTools&id=QaMJZqUX0cPaTfOF Model run 13.3.2017]
 
  
<rcode label="Show results" embed=0 graphics=1>
 
 
 
library(OpasnetUtils)
 
library(ggplot2)
 
library(reshape2)
 
  
objects.latest("Op_en7749", "preprocess")
+
== Rationale ==
  
for(i in c(5:6, 16, 29:30, 46:49, 85:86, 95:98, 135)) {
+
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 [https://docs.google.com/document/d/1IHyQzzdleM7IVY3PaoVImjxKXVWla6ahL62YuXjKe24 Google drive]. The actual data can be found from the link below (see Data).
  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")
+
=== Data ===
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)+
+
==== Questionnaire ====
  theme_gray(base_size = 24) + theme(axis.text.x=element_text(angle=90,hjust=1,vjust=0.4))
 
  
surveytemp <- subset(survey, survey[[31]] == "Yes")
+
Original datafile [[File:Goherr_fish_consumption.csv]].
  
temp <- melt(surveytemp, measure.vars = 38:43, variable.name = "Variable", value.name = "Value")
+
{{hidden|
ggplot(temp, aes(x = Value, fill = Country)) + geom_histogram(position = "dodge") + facet_wrap(~ Variable)+
+
<t2b name="Questions in the Goherr questionnaire" index="Index,Full question" obs="Question" unit="-">
  theme_gray(base_size = 24) + labs(title = "Baltic salmon sources")
+
V1|Country|Country
 +
V2|1. What is the postal code of your residence?|Postal code
 +
V3|2. What is your gender?|Gender
 +
V4|2. How old are you?|Age
 +
V5|3. What is the level of your highest education?|Education
 +
V6|4. How would you rate the purchasing power of your household?|Purchasing power
 +
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
 +
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
 +
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
 +
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|No vitamins
 +
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|Unknown vitamins
 +
V12|6. Does someone in your household fish?:No|No fishing
 +
V13|6. Does someone in your household fish?:Yes, for recreational purposes|Recreational
 +
V14|6. Does someone in your household fish?:Yes, for household use|Household use
 +
V15|6. Does someone in your household fish?:Yes, as a professional fisherman|Professional fishing
 +
V16|7. Do you eat fish?|Eat fish
 +
V17|8. Why don't you eat fish?:Allergy|Allergy
 +
V18|8. Why don't you eat fish?:I don't like the taste|Bad taste
 +
V19|8. Why don't you eat fish?:I don't like fish bones|Bones
 +
V20|8. Why don't you eat fish?:Difficult to cook|Cooking difficult
 +
V21|8. Why don't you eat fish?Health risks|Health risks fish
 +
V22|8. Why don't you eat fish?:I'm not used to|Not used to
 +
V23|8. Why don't you eat fish?:I'm concerned about the sustainability of fish stocks|Sustainability fish
 +
V24|8. Why don't you eat fish?:It is ethically wrong to eat animals|Ethically wrong
 +
V25|8. Why don't you eat fish?:I'm a vegetarian|Vegetarian
 +
V26|8.  Why don't you eat fish?:Don't Know|Why not fish unknown
 +
V27|8. Why don't you eat fish?:Other|Why not fish other
 +
V28|8. Why don't you eat fish?:Other open|Why not fish open
 +
V29|9. How often do you eat fish?|How often fish
 +
V30|10. Do you eat salmon or trout at least sometimes?|Eat salmon
 +
V31|11. Which of the following salmonids do you eat at least sometimes?:Wild Baltic salmon|Baltic salmon
 +
V32|11. Which of the following salmonids do you eat at least sometimes?:Wild Atlantic salmon|Atlantic salmon
 +
V33|11. Which of the following salmonids do you eat at least sometimes?:Norwegian cultured salmon|Norwegian salmon
 +
V34|11. Which of the following salmonids do you eat at least sometimes?:Rainbow trout|Rainbow trout
 +
V35|11. Which of the following salmonids do you eat at least sometimes?:I never know which salmon I eat|Unknown salmon
 +
V36|11. Which of the following salmonids do you eat at least sometimes?:Some other salmon, what?|Other salmon
 +
V37|11. Which of the following salmonids do you eat at least sometimes?:Other open|Salmon open
 +
V38|12. Where do you usually get the wild Baltic salmon that you eat?:Someone in our hosehold catches it|Own catch BS
 +
V39|12. Where do you usually get the wild Baltic salmon that you eat?:Buy from a grocery store|Grocery store BS
 +
V40|12. Where do you usually get the wild Baltic salmon that you eat?:Buy from a market place/fishmonger or similar|Market place BS
 +
V41|12. Where do you usually get the wild Baltic salmon that you eat?:Eat in a restaurant/cafe/canteen or similar|Restaurant BS
 +
V42|12. Where do you usually get the wild Baltic salmon that you eat?:Directly from a fisher|Fisher BS
 +
V43|12. Where do you usually get the wild Baltic salmon that you eat?:From a relative, friend, or similar|Friend BS
 +
V44|12. Where do you usually get the wild Baltic salmon that you eat?:Other|Where BS other
 +
V45|12. Where do you usually get the wild Baltic salmon that you eat?:Other open|Where BS open
 +
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 BS
 +
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 BS
 +
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 BS
 +
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 BS
 +
V50|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It tastes good|Tastes good BS
 +
V51|17. Select up to three most important reasons for you to eat wild Baltic salmon?:I caught it|Self caught BS
 +
V52|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's easy to cook|Easy to cook BS
 +
V53|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's quick to cook|Quick to cook BS
 +
V54|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's readily available|Readily available BS
 +
V55|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's healthy|Healthy BS
 +
V56|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's inexpensive|Inexpensive BS
 +
V57|17. Select up to three most important reasons for you to eat wild Baltic salmon?:My family likes it|Family likes it BS
 +
V58|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's environmentally/climate friendly choice|Environmental BS
 +
V59|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's a traditional dish|Traditional BS
 +
V60|17. Select up to three most important reasons for you to eat wild Baltic salmon?:Some other reasons, what?|Why BS other
 +
V61|17. Select up to three most important reasons for you to eat wild Baltic salmon?:Other open|Why BS open
 +
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|Bad taste BS
 +
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 BS
 +
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|Cannot cook BS
 +
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 BS
 +
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 BS
 +
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 BS
 +
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 BS
 +
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 BS
 +
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 BS
 +
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?|Why not BS other
 +
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|Why not BS open
 +
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 BS
 +
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 BS
 +
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 BS
 +
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 BS
 +
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 BS
 +
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 BS
 +
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 BS
 +
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 BS
 +
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|Recommended BS
 +
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|Not recommended BS
 +
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?|Change BS other
 +
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|Change BS open
 +
V85|20. Do you eat any herring (Baltic or other) at least sometimes?|Eat herring
 +
V86|21. Do you eat Baltic herring? This means herring caught from Baltic sea, not from Atlantic ocean|Eat BH
 +
V87|22. Where do you usually get the wild Baltic herring that you eat?:Someone in our hosehold catches it|Own catch BH
 +
V88|22. Where do you usually get the wild Baltic herring that you eat?:Buy from a grocery store|Grocery store BH
 +
V89|22. Where do you usually get the wild Baltic herring that you eat?:Buy from a market place/fishmonger or similar|Market place BH
 +
V90|22. Where do you usually get the wild Baltic herring that you eat?:Eat in a restaurant/cafe/canteen or similar|Restaurant BH
 +
V91|22. Where do you usually get the wild Baltic herring that you eat?:Directly from a fisher|Fisher BH
 +
V92|22. Where do you usually get the wild Baltic herring that you eat?:From a relative, friend, or similar|Friend BH
 +
V93|22. Where do you usually get the wild Baltic herring that you eat?:Other|Where BH other
 +
V94|22. Where do you usually get the wild Baltic herring that you eat?:Other open|Where BH open
 +
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 BH
 +
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 BH
 +
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 BH
 +
V98|26. How much on average of Baltic herring as a side dish (think of a regular sized lunch plate)?|How much side BH
 +
V99|27. Select up to three most important reasons for you to eat Baltic herring?:It tastes good|Tastes good BH
 +
V100|27. Select up to three most important reasons for you to eat Baltic herring?:I caught it|Self caught BH
 +
V101|27. Select up to three most important reasons for you to eat Baltic herring?:It's easy to cook|Easy to cook BH
 +
V102|27. Select up to three most important reasons for you to eat Baltic herring?:It's quick to cook|Quick to cook BH
 +
V103|27. Select up to three most important reasons for you to eat Baltic herring?:It's readily available|Readily available BH
 +
V104|27. Select up to three most important reasons for you to eat Baltic herring?:It's healthy|Healthy BH
 +
V105|27. Select up to three most important reasons for you to eat Baltic herring?:It's inexpensive|Inexpensive BH
 +
V106|27. Select up to three most important reasons for you to eat Baltic herring?:My family likes it|Family likes it BH
 +
V107|27. Select up to three most important reasons for you to eat Baltic herring?:It's environmentally/climate friendly choice|Environmental BH
 +
V108|27. Select up to three most important reasons for you to eat Baltic herring?:It's a traditional dish|Traditional BH
 +
V109|27. Select up to three most important reasons for you to eat Baltic herring?:Some other reasons, what?|Why BH other
 +
V110|27. Select up to three most important reasons for you to eat Baltic herring?:Other open|Why BH open
 +
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|Bad taste BH
 +
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 BH
 +
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|Cannot cook BH
 +
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 BH
 +
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 BH
 +
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 BH
 +
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 BH
 +
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 BH
 +
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 BH
 +
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 BH
 +
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?|Why not BH other
 +
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|Why not BH open
 +
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 BH
 +
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 BH
 +
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 BH
 +
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 BH
 +
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 BH
 +
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 BH
 +
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 BH
 +
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 BH
 +
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|Recommended BH
 +
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|Not recommended BH
 +
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?|Change BH other
 +
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|Change BH open
 +
V135|30. Are you familiar with the possible eating recommendations of Baltic herring and/or wild Baltic salmon in your country?|Recommendation awareness
 +
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 BS
 +
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 BS
 +
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 all BS
 +
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 BS
 +
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 BS
 +
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 BS
 +
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 BS
 +
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 BS
 +
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
 +
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 BH
 +
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 BH
 +
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 all BH
 +
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 BH
 +
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 BH
 +
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 BH
 +
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 BH
 +
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 BH
 +
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 BH
 +
V154|Weighting|Weighting
 +
V155|How old are you? Age category FI|Age category FI
 +
V156|2. How old are you? Age category SWE, DK, EST|Age category other
 +
V157|Row number of respondent|Row
 +
V158|Age category "18-45" or ">45"|Ages
 +
</t2b>
 +
}}
  
temp <- melt(surveytemp, measure.vars = 50:59, variable.name = "Variable", value.name = "Value")
+
==== Assumptions ====
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")
+
The following assumptions are used:
ggplot(temp, aes(x = Value, fill = Country)) + geom_histogram(position = "dodge") + facet_wrap(~ Variable)+
+
<t2b name="Assumptions for calculations" index="Variable,Value,Unit" obs="Result" desc="Description" unit="-">
  theme_gray(base_size = 24) + theme(axis.text.x=element_text(angle=90,hjust=1,vjust=0.4)) + labs(title = "Baltic salmon actions")
+
freq|1|times /a|0|Never
 +
freq|2|times /a|0.5 - 0.9|less than once a year
 +
freq|3|times /a|2 - 5|A few times a year
 +
freq|4|times /a|12 - 36|1 - 3 times per month
 +
freq|5|times /a|52|once a week
 +
freq|6|times /a|104 - 208|2 - 4 times per week
 +
freq|7|times /a|260 - 364|5 or more times per week
 +
amdish|1|g /serving|20 - 70|1/6 plate or below (50 grams)
 +
amdish|2|g /serving|70 - 130|1/3 plate (100 grams)
 +
amdish|3|g /serving|120 - 180|1/2 plate (150 grams)
 +
amdish|4|g /serving|170 - 230|2/3 plate (200 grams)
 +
amdish|5|g /serving|220 - 280|5/6 plate (250 grams)
 +
amdish|6|g /serving|270 - 400|full plate (300 grams)
 +
amdish|7|g /serving|400 - 550|overly full plate (500 grams)
 +
ingredient||fraction|0.1 - 0.3|Fraction of fish in the dish
 +
amside|1|g /serving|20 - 70|1/6 plate or below (50 grams)
 +
amside|2|g /serving|70 - 130|1/4 plate (100 grams)
 +
amside|3|g /serving|120 - 180|1/2 plate (150 grams)
 +
amside|4|g /serving|170 - 230|2/3 plate (200 grams)
 +
amside|5|g /serving|220 - 280|5/6 plate (250 grams)
 +
</t2b>
  
surveytemp <- subset(survey, survey[[31]] == "No")
+
==== Preprocessing ====
 
 
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")
 
 
</rcode>
 
 
== 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 [https://docs.google.com/document/d/1IHyQzzdleM7IVY3PaoVImjxKXVWla6ahL62YuXjKe24 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. The code stores a data.frame named survey.
 
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 13.4.2017 [http://en.opasnet.org/en-opwiki/index.php?title=Special:RTools&id=1QsyOMpBHrvY0H2o]
 
* Model run 13.4.2017 [http://en.opasnet.org/en-opwiki/index.php?title=Special:RTools&id=1QsyOMpBHrvY0H2o]
 
* Model run 20.4.2017 [http://en.opasnet.org/en-opwiki/index.php?title=Special:RTools&id=JuLynAzoCuvu0U3g] (contains surv and helping vectors)
 
* Model run 20.4.2017 [http://en.opasnet.org/en-opwiki/index.php?title=Special:RTools&id=JuLynAzoCuvu0U3g] (contains surv and helping vectors)
 +
* Model run 21.4.2017 [http://en.opasnet.org/en-opwiki/index.php?title=Special:RTools&id=YdoQdrDj5uQnoogH] surv contains Eatfish, Eatherr and Eatsalm as columns.
  
 
<rcode name="preprocess2" label="Preprocess (only for developers)">
 
<rcode name="preprocess2" label="Preprocess (only for developers)">
Line 204: Line 331:
  
 
# Row numbers for respondents that have eaten fish, Baltic salmon, and Baltic herring
 
# Row numbers for respondents that have eaten fish, Baltic salmon, and Baltic herring
eatfish <- surv[[4]] %in% 1
+
surv$Eatfish <- surv[[4]] %in% 1
eatsalm <- !(surv[[7]] %in% 0 | is.na(rowSums(surv[7:11])))
+
surv$Eatsalm <- !(surv[[7]] %in% 0 | is.na(rowSums(surv[7:11])))
eatherr <- surv[[12]] %in% 1 & !is.na(rowSums(surv[13:16]))
+
surv$Eatherr <- surv[[12]] %in% 2 & !is.na(rowSums(surv[13:16]))
  
 
# Oletetaan, että covarianssimatriisi on vakio kaikille maille ja sukupuolille yms
 
# Oletetaan, että covarianssimatriisi on vakio kaikille maille ja sukupuolille yms
 
# mutta keskiarvo on spesifi näille ja kysymykselle.
 
# 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
+
#qlen <-  c(4,2,2,2,6,2,2,7,7,7,5,2,7,7,7,5) # Number of options in each question of surv
 
# qlen not needed when dbinom is not used.
 
# qlen not needed when dbinom is not used.
  
Line 219: Line 346:
 
fisl
 
fisl
  
objects.store(survey, surv, agel, countryl, genderl, fisl, eatfish, eatsalm, eatherr)
+
objects.store(survey, surv, agel, countryl, genderl, fisl)
cat("Data.frames survey and surv, and vectors agel, countryl, genderl, fisl, eatfish, eatsalm, eatherr were stored.\n")
+
cat("Data.frames survey and surv, and vectors agel, countryl, genderl and fisl were stored.\n")
 
</rcode>
 
</rcode>
  
Line 257: Line 384:
 
* Omega3, Vitamin D, and Other vitamins
 
* Omega3, Vitamin D, and Other vitamins
  
Study plan:
+
Model runs
* Determinants
+
* [http://en.opasnet.org/en-opwiki/index.php?title=Special:RTools&id=QaMJZqUX0cPaTfOF Model run 13.3.2017]
  
 
<rcode name="descrstat" label="Descriptive statistics" graphics=1>
 
<rcode name="descrstat" label="Descriptive statistics" graphics=1>
Line 269: Line 396:
 
library(vegan)
 
library(vegan)
  
objects.latest("Op_en7749", "preprocess2") # [[Goherr: Fish consumption study]]: survey, surcol
+
objects.latest("Op_en7749", "preprocess2") # [[Goherr: Fish consumption study]]: survey, surv
 +
 
 +
############################### From a previous code on Answer
 +
 
 +
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")
 +
 
 +
#####################################3
  
 
temp <- sapply(survey, as.numeric) # Can be done for surv to get a smaller matrix
 
temp <- sapply(survey, as.numeric) # Can be done for surv to get a smaller matrix
Line 299: Line 486:
  
 
## Eating frequencies of fish and Baltic salmon and herring with random noise, all
 
## Eating frequencies of fish and Baltic salmon and herring with random noise, all
pl <- surv[eatfish,c(5,8,10,13,15)]
+
pl <- surv[surv$Eatfish,c(5,8,10,13,15)]
 
scatterplotMatrix(
 
scatterplotMatrix(
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
Line 313: Line 500:
  
 
## Fish eating questions with some random noise, all
 
## Fish eating questions with some random noise, all
pl <- surv[eatfish,c(5,6,7,12)]
+
pl <- surv[surv$Eatfish,c(5,6,7,12)]
 
scatterplotMatrix(
 
scatterplotMatrix(
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
Line 319: Line 506:
  
 
## Fish eating questions with some random noise, FI
 
## Fish eating questions with some random noise, FI
pl <- surv[surv$Country == 1 & eatfish,c(5,6,7,12)]
+
pl <- surv[surv$Country == 1 & surv$Eatfish,c(5,6,7,12)]
 
scatterplotMatrix(
 
scatterplotMatrix(
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
Line 325: Line 512:
  
 
## Fish eating questions with some random noise, SWE
 
## Fish eating questions with some random noise, SWE
pl <- surv[surv$Country == 2 & eatfish,c(5,6,7,12)]
+
pl <- surv[surv$Country == 2 & surv$Eatfish,c(5,6,7,12)]
 
scatterplotMatrix(
 
scatterplotMatrix(
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
Line 331: Line 518:
  
 
## Fish eating questions with some random noise, Dk
 
## Fish eating questions with some random noise, Dk
pl <- surv[surv$Country == 3 & eatfish,c(5,6,7,12)]
+
pl <- surv[surv$Country == 3 & surv$Eatfish,c(5,6,7,12)]
 
scatterplotMatrix(
 
scatterplotMatrix(
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
Line 337: Line 524:
  
 
## Fish eating questions with some random noise, EST
 
## Fish eating questions with some random noise, EST
pl <- surv[surv$Country == 4 & eatfish,c(5,6,7,12)]
+
pl <- surv[surv$Country == 4 & surv$Eatfish,c(5,6,7,12)]
 
scatterplotMatrix(
 
scatterplotMatrix(
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
Line 343: Line 530:
  
 
## Baltic herring questions with some random noise
 
## Baltic herring questions with some random noise
pl <- surv[eatherr,13:16]
+
pl <- surv[surv$Eatherr,13:16]
 
scatterplotMatrix(
 
scatterplotMatrix(
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
Line 349: Line 536:
  
 
## Baltic salmon questions with some random noise
 
## Baltic salmon questions with some random noise
pl <- surv[eatsalm,8:11]
+
pl <- surv[surv$Eatsalm,8:11]
 
scatterplotMatrix(
 
scatterplotMatrix(
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
 
   data.matrix(pl) + runif(nrow(pl)*ncol(pl), -0.5, 0.5)
Line 416: Line 603:
 
* Model run 20.4.2017 [http://en.opasnet.org/en-opwiki/index.php?title=Special:RTools&id=QXRSZOliFNfJil39] code works but needs a safety check against outliers
 
* Model run 20.4.2017 [http://en.opasnet.org/en-opwiki/index.php?title=Special:RTools&id=QXRSZOliFNfJil39] code works but needs a safety check against outliers
  
<rcode name="bayes" label="Initiate Bayes model (for developers only)" graphics=1>
+
<rcode name="bayes" label="Initiate Bayes model (for developers only)" graphics=1 variables="
 +
name:store|description:Do you want to plot model results or store ovariables derived?|type:selection|options:
 +
FALSE;Plot results;
 +
TRUE;Store ovariables derived
 +
">
 
# This is code Op_en7749/bayes on page [[Goherr: Fish consumption study]]
 
# This is code Op_en7749/bayes on page [[Goherr: Fish consumption study]]
  
Line 432: Line 623:
  
 
objects.latest("Op_en7749", "preprocess2") # [[Goherr: Fish consumption study]]: survey, surv, ...
 
objects.latest("Op_en7749", "preprocess2") # [[Goherr: Fish consumption study]]: survey, surv, ...
 
# Move this to preprocess. Combine eatherr, eatsalm and eathfish to surv as columns
 
eatherr <- surv[[12]] %in% 2 & !is.na(rowSums(surv[13:16]))
 
  
 
cat("Version with multivariate normal.\n")
 
cat("Version with multivariate normal.\n")
  
 
mod <- textConnection("
 
mod <- textConnection("
                      model{
+
  model{
                      for(i in 1:S) {
+
    for(i in 1:S) {
                      survs[i,1:4] ~ dmnorm(mus[], Omegas[,])
+
      survs[i,1:4] ~ dmnorm(mus[], Omegas[,])
                      }
+
    }
                      for(j in 1:H) {
+
    for(j in 1:H) {
                      survh[j,1:4] ~ dmnorm(muh[], Omegah[,])
+
      survh[j,1:4] ~ dmnorm(muh[], Omegah[,])
                      }
+
    }
                      mus[1:4] ~ dmnorm(mu0[1:4], S2[1:4,1:4])
+
    mus[1:4] ~ dmnorm(mu0[1:4], S2[1:4,1:4])
                      Omegas[1:4,1:4] ~ dwish(S3[1:4,1:4],S)
+
    Omegas[1:4,1:4] ~ dwish(S3[1:4,1:4],S)
                      anss.pred ~ dmnorm(mus[], Omegas[,])
+
    anss.pred ~ dmnorm(mus[], Omegas[,])
                      muh[1:4] ~ dmnorm(mu0[1:4], S2[1:4,1:4])
+
    muh[1:4] ~ dmnorm(mu0[1:4], S2[1:4,1:4])
                      Omegah[1:4,1:4] ~ dwish(S3[1:4,1:4],H)
+
    Omegah[1:4,1:4] ~ dwish(S3[1:4,1:4],H)
                      ansh.pred ~ dmnorm(muh[], Omegah[,])
+
    ansh.pred ~ dmnorm(muh[], Omegah[,])
                      }
+
  }
                      ")
+
")
  
 
jags <- jags.model(
 
jags <- jags.model(
 
   mod,
 
   mod,
 
   data = list(
 
   data = list(
     survs = surv[eatsalm,c(8:11)],
+
     survs = surv[surv$Eatsalm,c(8:11)],
     S = sum(eatsalm),
+
     S = sum(surv$Eatsalm),
     survh = surv[eatherr,c(13:16)],
+
     survh = surv[surv$Eatherr,c(13:16)],
     H = sum(eatherr),
+
     H = sum(surv$Eatherr),
 
     mu0 = rep(2,4),
 
     mu0 = rep(2,4),
 
     S2 = diag(4)/100000,
 
     S2 = diag(4)/100000,
Line 500: Line 688:
 
)
 
)
  
######## fish.param contains expected values of the distribution parameters from the model
+
if(!store) {
 +
  # Mu for all questions about salmon
 +
  scatterplotMatrix(t(js[,,1,1]))
 +
  # All parameters for question 1 about salmon
 +
  scatterplotMatrix(t(js[1,,,1]))
 +
 
 +
  jsd <- melt(js)
 +
  ggplot(jsd, aes(x=value, colour=Question))+geom_density()+facet_grid(Parameter ~ Fish)
 +
  #ggplot(as.data.frame(js), aes(x = anss.pred, y = Sampled))+geom_point()+stat_ellipse()
  
fish.param <- list(
+
  coda.j <- coda.samples(
  mu = apply(js[,,1,], MARGIN = c(1,3), FUN = mean),
+
     jags,  
  Omega = lapply(
+
     c('mus', 'Omegas', 'anss.pred', 'mu0'),  
     1:2,
+
     1000
     FUN = function(x) {
 
      solve(apply(js[,,2:5,], MARGIN = c(1,3,4), FUN = mean)[,,x])
 
     } # solve matrix: precision->covariace
 
 
   )
 
   )
)
 
  
jsp <- Ovariable(
+
  plot(coda.j)
   "jsp",
+
 
   dependencies = data.frame(Name = "fish.param"),
+
} else {
  formula = function(...) {
+
 
     jsp <- lapply(
+
   ######## fish.param contains expected values of the distribution parameters from the model
 +
 
 +
   fish.param <- list(
 +
    mu = apply(js[,,1,], MARGIN = c(1,3), FUN = mean),
 +
     Omega = lapply(
 
       1:2,
 
       1:2,
 
       FUN = function(x) {
 
       FUN = function(x) {
         mvrnorm(openv$N, fish.param$mu[,x], fish.param$Omega[[x]])
+
         solve(apply(js[,,2:5,], MARGIN = c(1,3,4), FUN = mean)[,,x])
       }
+
       } # solve matrix: precision->covariace
 
     )
 
     )
      
+
  )
     jsp <- rbind(
+
 
      cbind(
+
  jsp <- Ovariable(
        Fish = "Salmon",
+
     "jsp",
        Iter = 1:nrow(jsp[[1]]),
+
     dependencies = data.frame(Name = "fish.param"),
        as.data.frame(jsp[[1]])
+
    formula = function(...) {
      ),
+
      jsp <- lapply(
      cbind(
+
        1:2,
        Fish = "Herring",
+
        FUN = function(x) {
        Iter = 1:nrow(jsp[[2]]),
+
          mvrnorm(openv$N, fish.param$mu[,x], fish.param$Omega[[x]])
        as.data.frame(jsp[[2]])
+
        }
 +
      )
 +
     
 +
      jsp <- rbind(
 +
        cbind(
 +
          Fish = "Salmon",
 +
          Iter = 1:nrow(jsp[[1]]),
 +
          as.data.frame(jsp[[1]])
 +
        ),
 +
        cbind(
 +
          Fish = "Herring",
 +
          Iter = 1:nrow(jsp[[2]]),
 +
          as.data.frame(jsp[[2]])
 +
        )
 
       )
 
       )
    )
+
      jsp <- melt(jsp, id.vars = c("Iter", "Fish"), variable.name = "Question", value.name = "Result")
    jsp <- melt(jsp, id.vars = c("Iter", "Fish"), variable.name = "Question", value.name = "Result")
+
      jsp <- Ovariable(output=jsp, marginal = colnames(jsp) %in% c("Iter", "Fish", "Question"))
    jsp <- Ovariable(output=jsp, marginal = colnames(jsp) %in% c("Iter", "Fish", "Question"))
+
      return(jsp)
    return(jsp)
+
    }
  }
+
  )
)
 
  
# Combine modelled survey answers with estimated amounts and frequencies by:
+
  # Combine modelled survey answers with estimated amounts and frequencies by:
# Rounding the modelled result and merging that with value in ovariable assump
+
  # Rounding the modelled result and merging that with value in ovariable assump
  
 
often <- Ovariable(
 
often <- Ovariable(
Line 628: Line 836:
 
)
 
)
  
jsp <- EvalOutput(jsp)
+
  objects.store(fish.param, jsp, assump, often, much, oftenside, muchside, amount)
often <- EvalOutput(often)
+
   cat("List fish.param and ovariables jsp, assump, often, much, oftenside, muchside, amount stored.\n")
much <- EvalOutput(much)
+
}
oftenside <- EvalOutput(oftenside)
 
muchside <- EvalOutput(muchside)
 
assump <- EvalOutput(assump)
 
 
 
amount <- EvalOutput(amount)
 
oprint(head(amount@output))
 
#scatterplotMatrix(jsp[[1]])
 
#scatterplotMatrix(t(js[,,6]))
 
 
 
#ggplot(melt(js), aes(x=value, colour=Var2))+geom_density()
 
#ggplot(as.data.frame(js), aes(x = anss.pred, y = Sampled))+geom_point()+stat_ellipse()
 
#coda.j <- coda.samples(
 
#  jags,
 
#  c('mus', 'Omegas', 'anss.pred', 'mu0'),
 
#  1000
 
#)
 
 
 
#plot(coda.j)
 
#jh <- array(
 
#  samps.j$ansh.pred,
 
#  dim = c(4,1000,4),
 
#  dimnames = list(Question = 1:4, Iter = 1:1000, Seed = 1:4)
 
#)
 
#scatterplotMatrix(t(jh[,,1]))
 
 
 
#objects.store(oftenh, muchh, ofsideh, musideh, oftens, muchs, ofsides, musides)
 
#cat("Ovariables oftenh, muchh, ofsideh, musideh, oftens, muchs, ofsides, musides stored.\n")</rcode>
 
 
 
==== 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.
 
 
 
<rcode label="Calculate amount" name=calculate_amount embed=1 graphics=1 store=1>
 
# 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")
 
 
</rcode>
 
</rcode>
  
==== Assumptions ====
+
=== Dependencies ===
 
 
The following assumptions are used:
 
<t2b name="Assumptions for calculations" index="Variable,Value,Unit" obs="Result" unit="-">
 
freq|6|times per year|260 - 364
 
freq|5|times per year|104 - 208
 
freq|4|times per year|52
 
freq|3|times per year|12 - 36
 
freq|2|times per year|2 - 5
 
freq|1|times per year|0.5 - 0.9
 
freq|0|times per year|0
 
amdish|0|grams / serving|20 - 50
 
amdish|1|grams / serving|70 - 100
 
amdish|2|grams / serving|120 - 150
 
amdish|3|grams / serving|170 - 200
 
amdish|4|grams / serving|220 - 250
 
amdish|5|grams / serving|270 - 300
 
amdish|6|grams / serving|450 - 500
 
ingredient||fraction|0.1 - 0.3
 
amside|0|grams / serving|20 - 50
 
amside|1|grams / serving|70 - 100
 
amside|2|grams / serving|120 - 150
 
amside|3|grams / serving|170 - 200
 
amside|4|grams / serving|220 - 250
 
</t2b>
 
 
 
==== Questionnaire ====
 
 
 
{{hidden|
 
<t2b name="Questions in the Goherr questionnaire" index="Index,Full question" obs="Question" unit="-">
 
V1|Country|Country
 
V2|1. What is the postal code of your residence?|Postal code
 
V3|2. What is your gender?|Gender
 
V4|2. How old are you?|Age
 
V5|3. What is the level of your highest education?|Education
 
V6|4. How would you rate the purchasing power of your household?|Purchasing power
 
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
 
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
 
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
 
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|No vitamins
 
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|Unknown vitamins
 
V12|6. Does someone in your household fish?:No|No fishing
 
V13|6. Does someone in your household fish?:Yes, for recreational purposes|Recreational
 
V14|6. Does someone in your household fish?:Yes, for household use|Household use
 
V15|6. Does someone in your household fish?:Yes, as a professional fisherman|Professional fishing
 
V16|7. Do you eat fish?|Eat fish
 
V17|8. Why don't you eat fish?:Allergy|Allergy
 
V18|8. Why don't you eat fish?:I don't like the taste|Bad taste
 
V19|8. Why don't you eat fish?:I don't like fish bones|Bones
 
V20|8. Why don't you eat fish?:Difficult to cook|Cooking difficult
 
V21|8. Why don't you eat fish?Health risks|Health risks fish
 
V22|8. Why don't you eat fish?:I'm not used to|Not used to
 
V23|8. Why don't you eat fish?:I'm concerned about the sustainability of fish stocks|Sustainability fish
 
V24|8. Why don't you eat fish?:It is ethically wrong to eat animals|Ethically wrong
 
V25|8. Why don't you eat fish?:I'm a vegetarian|Vegetarian
 
V26|8.  Why don't you eat fish?:Don't Know|Why not fish unknown
 
V27|8. Why don't you eat fish?:Other|Why not fish other
 
V28|8. Why don't you eat fish?:Other open|Why not fish open
 
V29|9. How often do you eat fish?|How often fish
 
V30|10. Do you eat salmon or trout at least sometimes?|Eat salmon
 
V31|11. Which of the following salmonids do you eat at least sometimes?:Wild Baltic salmon|Baltic salmon
 
V32|11. Which of the following salmonids do you eat at least sometimes?:Wild Atlantic salmon|Atlantic salmon
 
V33|11. Which of the following salmonids do you eat at least sometimes?:Norwegian cultured salmon|Norwegian salmon
 
V34|11. Which of the following salmonids do you eat at least sometimes?:Rainbow trout|Rainbow trout
 
V35|11. Which of the following salmonids do you eat at least sometimes?:I never know which salmon I eat|Unknown salmon
 
V36|11. Which of the following salmonids do you eat at least sometimes?:Some other salmon, what?|Other salmon
 
V37|11. Which of the following salmonids do you eat at least sometimes?:Other open|Salmon open
 
V38|12. Where do you usually get the wild Baltic salmon that you eat?:Someone in our hosehold catches it|Own catch BS
 
V39|12. Where do you usually get the wild Baltic salmon that you eat?:Buy from a grocery store|Grocery store BS
 
V40|12. Where do you usually get the wild Baltic salmon that you eat?:Buy from a market place/fishmonger or similar|Market place BS
 
V41|12. Where do you usually get the wild Baltic salmon that you eat?:Eat in a restaurant/cafe/canteen or similar|Restaurant BS
 
V42|12. Where do you usually get the wild Baltic salmon that you eat?:Directly from a fisher|Fisher BS
 
V43|12. Where do you usually get the wild Baltic salmon that you eat?:From a relative, friend, or similar|Friend BS
 
V44|12. Where do you usually get the wild Baltic salmon that you eat?:Other|Where BS other
 
V45|12. Where do you usually get the wild Baltic salmon that you eat?:Other open|Where BS open
 
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 BS
 
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 BS
 
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 BS
 
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 BS
 
V50|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It tastes good|Tastes good BS
 
V51|17. Select up to three most important reasons for you to eat wild Baltic salmon?:I caught it|Self caught BS
 
V52|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's easy to cook|Easy to cook BS
 
V53|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's quick to cook|Quick to cook BS
 
V54|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's readily available|Readily available BS
 
V55|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's healthy|Healthy BS
 
V56|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's inexpensive|Inexpensive BS
 
V57|17. Select up to three most important reasons for you to eat wild Baltic salmon?:My family likes it|Family likes it BS
 
V58|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's environmentally/climate friendly choice|Environmental BS
 
V59|17. Select up to three most important reasons for you to eat wild Baltic salmon?:It's a traditional dish|Traditional BS
 
V60|17. Select up to three most important reasons for you to eat wild Baltic salmon?:Some other reasons, what?|Why BS other
 
V61|17. Select up to three most important reasons for you to eat wild Baltic salmon?:Other open|Why BS open
 
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|Bad taste BS
 
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 BS
 
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|Cannot cook BS
 
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 BS
 
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 BS
 
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 BS
 
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 BS
 
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 BS
 
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 BS
 
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?|Why not BS other
 
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|Why not BS open
 
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 BS
 
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 BS
 
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 BS
 
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 BS
 
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 BS
 
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 BS
 
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 BS
 
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 BS
 
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|Recommended BS
 
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|Not recommended BS
 
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?|Change BS other
 
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|Change BS open
 
V85|20. Do you eat any herring (Baltic or other) at least sometimes?|Eat herring
 
V86|21. Do you eat Baltic herring? This means herring caught from Baltic sea, not from Atlantic ocean|Eat BH
 
V87|22. Where do you usually get the wild Baltic herring that you eat?:Someone in our hosehold catches it|Own catch BH
 
V88|22. Where do you usually get the wild Baltic herring that you eat?:Buy from a grocery store|Grocery store BH
 
V89|22. Where do you usually get the wild Baltic herring that you eat?:Buy from a market place/fishmonger or similar|Market place BH
 
V90|22. Where do you usually get the wild Baltic herring that you eat?:Eat in a restaurant/cafe/canteen or similar|Restaurant BH
 
V91|22. Where do you usually get the wild Baltic herring that you eat?:Directly from a fisher|Fisher BH
 
V92|22. Where do you usually get the wild Baltic herring that you eat?:From a relative, friend, or similar|Friend BH
 
V93|22. Where do you usually get the wild Baltic herring that you eat?:Other|Where BH other
 
V94|22. Where do you usually get the wild Baltic herring that you eat?:Other open|Where BH open
 
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 BH
 
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 BH
 
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 BH
 
V98|26. How much on average of Baltic herring as a side dish (think of a regular sized lunch plate)?|How much side BH
 
V99|27. Select up to three most important reasons for you to eat Baltic herring?:It tastes good|Tastes good BH
 
V100|27. Select up to three most important reasons for you to eat Baltic herring?:I caught it|Self caught BH
 
V101|27. Select up to three most important reasons for you to eat Baltic herring?:It's easy to cook|Easy to cook BH
 
V102|27. Select up to three most important reasons for you to eat Baltic herring?:It's quick to cook|Quick to cook BH
 
V103|27. Select up to three most important reasons for you to eat Baltic herring?:It's readily available|Readily available BH
 
V104|27. Select up to three most important reasons for you to eat Baltic herring?:It's healthy|Healthy BH
 
V105|27. Select up to three most important reasons for you to eat Baltic herring?:It's inexpensive|Inexpensive BH
 
V106|27. Select up to three most important reasons for you to eat Baltic herring?:My family likes it|Family likes it BH
 
V107|27. Select up to three most important reasons for you to eat Baltic herring?:It's environmentally/climate friendly choice|Environmental BH
 
V108|27. Select up to three most important reasons for you to eat Baltic herring?:It's a traditional dish|Traditional BH
 
V109|27. Select up to three most important reasons for you to eat Baltic herring?:Some other reasons, what?|Why BH other
 
V110|27. Select up to three most important reasons for you to eat Baltic herring?:Other open|Why BH open
 
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|Bad taste BH
 
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 BH
 
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|Cannot cook BH
 
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 BH
 
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 BH
 
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 BH
 
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 BH
 
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 BH
 
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 BH
 
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 BH
 
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?|Why not BH other
 
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|Why not BH open
 
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 BH
 
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 BH
 
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 BH
 
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 BH
 
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 BH
 
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 BH
 
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 BH
 
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 BH
 
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|Recommended BH
 
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|Not recommended BH
 
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?|Change BH other
 
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|Change BH open
 
V135|30. Are you familiar with the possible eating recommendations of Baltic herring and/or wild Baltic salmon in your country?|Recommendation awareness
 
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 BS
 
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 BS
 
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 all BS
 
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 BS
 
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 BS
 
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 BS
 
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 BS
 
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 BS
 
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
 
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 BH
 
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 BH
 
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 all BH
 
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 BH
 
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 BH
 
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 BH
 
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 BH
 
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 BH
 
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 BH
 
V154|Weighting|Weighting
 
V155|How old are you? Age category FI|Age category FI
 
V156|2. How old are you? Age category SWE, DK, EST|Age category other
 
V157|Row number of respondent|Row
 
V158|Age category "18-45" or ">45"|Ages
 
</t2b>
 
}}
 
  
=== Dependencies ===
 
 
The survey data will be used as input in the [[Benefit-risk_assessment_of_Baltic_herring_and_salmon_intake | benefit-risk assessment of Baltic herring and salmon intake]], which is part of the WP5 work in Goherr-project.
 
The survey data will be used as input in the [[Benefit-risk_assessment_of_Baltic_herring_and_salmon_intake | benefit-risk assessment of Baltic herring and salmon intake]], which is part of the WP5 work in Goherr-project.
 
=== Formula ===
 
  
 
==See also==
 
==See also==

Revision as of 10:34, 21 April 2017


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

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

Data

Questionnaire

Original datafile File:Goherr fish consumption.csv.



Assumptions

The following assumptions are used:

Assumptions for calculations(-)
ObsVariableValueUnitResultDescription
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)

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

  • Model run 13.4.2017 [1]
  • Model run 20.4.2017 [2] (contains surv and helping vectors)
  • Model run 21.4.2017 [3] surv contains Eatfish, Eatherr and Eatsalm as columns.

+ Show code

Analyses

Error creating thumbnail: Unable to save thumbnail to destination
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. [4]
  • Model run 3.3.2017. p has more dimensions. [5]
  • Model run 25.3.2017. Several model versions: strange binomial+multivarnormal, binomial, fractalised multivarnormal [6]
  • Model run 27.3.2017 [7]
  • Other models except multivariate normal were archived and removed from active code 29.3.2017.
  • Model run 29.3.2017 with raw data graphs [8]
  • Model run 29.3.2017 with salmon and herring ovariables stored [9]
  • Model run 13.4.2017 with first version of coordinate matrix and principal coordinate analysis [10]
  • Model run 20.4.2017 [11] code works but needs a safety check against outliers

Do you want to plot model results or store ovariables derived?:

+ Show code

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.

See also

Keywords

References


Related files

<mfanonymousfilelist></mfanonymousfilelist>

Goherr: Fish consumption study. Opasnet . [12]. Accessed 17 May 2024.