For this data exploration I used a data set from the Tidy Tuesday project (https://github.com/rfordatascience/tidytuesday/blob/master/data/2020/2020-02-18/readme.md), February 18, 2020. I also reviewed Andrew Couch’s EDA of this dataset, available on YouTube (https://www.youtube.com/watch?v=VKCPYet9qLM). The variables of interest for this analysis are consumption of food by food category (kg/person/year) and CO2 emission by food category (kg CO2/person/year).
Around the world people eat more plant-based products (2:1) than meat, when measuring consumption by kilograms per person per year. Plant-based products emit less co2 emissions than most meat-based products, especially beef, lamb, and goat. For some countries, most of their consumption is plant-based (i.e., Liberia’s consumption is 82% plant based) while others, mostly animal-based (i.e., Finland, 86%). When we look at foods in a lacto-ova vegetarian diet (eggs, milk & cheese + plants), we find that India’s consumption, for example, is 96% plant, egg, and dairy-based. If we go a step further and create a group representing an ova vegetarian diet (eggs + plants), Liberia’s consumption of these products comes in the highest at 84%. Since eggs are an efficient food as well in terms of co2 emissions, comparing an ova-vegetarian group to lacto-ova vegetarian and vegan may be meaningful.
There are many limitations to this analysis. For one, food consumption (“consumption”) is measured as kilograms of food per person per year. There could be some significant bias here, in that vegetables may weigh less that animal products. For example, the nutritional value of a kilogram of milk or cheese is very different from that of a kilogram of spinach. Second, analysis is done in a global accurate and does not take into account cultural differences across countries that account for population’s food choices. For example, pork consumption is nil or very low in countries with large Muslim populations.
Future analysis: Rather than looking at consumption measured in kilograms per person, a more meaningful analysis might be to consider consumption as calories consumed per person per year by food product. It would also be interesting to see what the relationship is between a country’s GDP and the proportion of their diet that is plant-based. The intention behind this data may have been to shine light on the fact that animal products are less energy efficient than plant-based products. However, this propensity towards a plant-based diet may have more to do with a country’s income than climate awareness.
| Variable | Class | Description |
|---|---|---|
| country | character | Country Name |
| food_category | character | Food Category |
| consumption | double | Consumption (kg/person/year) |
| co2_emission | double | Co2 Emission (kg CO2/person/year) |
I conducted a cursory review to check for NA values and determine number of food categories per country. Each of the 130 countries has a row for each of the 11 food categories. There were no NA values. I left all 0 values in the data set. ‘0’ values are present for 31 countries and the food categories to have a 0 value are soybeans (22), pork (8), and Lamb & Goat (1). Logically this makes sense to me, as soybeans are not a common food staple in all countries and pork is not consumed in some countries, especially those with large Muslim populations.
Here is a list of the 11 food categories.
## [1] "Pork" "Poultry"
## [3] "Beef" "Lamb & Goat"
## [5] "Fish" "Eggs"
## [7] "Milk - inc. cheese" "Wheat and Wheat Products"
## [9] "Rice" "Soybeans"
## [11] "Nuts inc. Peanut Butter"
Analysis: t-test
I conducted a t-test comparing means between co2 emission related to animal products versus non-animal (vegan) products. I created new categories to assign food products to a “Meat” or “Vegan” group. In the “Meat” group is Pork, Poultry, Beef, Lamb & Goat, Fish, Milk inc. cheese, and Eggs. The remaining food products make up the “Vegan” group.
These t-tests suggest that there is not a significant difference in consumption levels between meat and vegan groups but a very significant difference between the meat group and the vegan group in terms of co2 emissions. Overall, more vegan foods are consumed world-wide than meat foods by about 2:1 (13,769.92 meat, 26427.96 vegan).
## # A tibble: 2 × 3
## food_type mean_consumption mean_co2
## <chr> <dbl> <dbl>
## 1 Meat 29.0 108.
## 2 Vegan 26.5 14.7
##
## Two Sample t-test
##
## data: consumption by food_type
## t = 0.93514, df = 1428, p-value = 0.3499
## alternative hypothesis: true difference in means between group Meat and group Vegan is not equal to 0
## 95 percent confidence interval:
## -2.811294 7.933492
## sample estimates:
## mean in group Meat mean in group Vegan
## 29.04171 26.48062
##
## Two Sample t-test
##
## data: co2_emission by food_type
## t = 11.734, df = 1428, p-value < 0.00000000000000022
## alternative hypothesis: true difference in means between group Meat and group Vegan is not equal to 0
## 95 percent confidence interval:
## 78.0623 109.4009
## sample estimates:
## mean in group Meat mean in group Vegan
## 108.4682 14.7366
## # A tibble: 2 × 2
## food_type total_consumption
## <chr> <dbl>
## 1 Meat 26428.
## 2 Vegan 13770.
| Country | Percent of consumption plant-based, by kg per person |
|---|---|
| Liberia | 0.82 |
| Bangladesh | 0.81 |
| Guinea | 0.74 |
| Cambodia | 0.74 |
| Indonesia | 0.74 |
| Country | Percent of consumption meat-based, by kg per person |
|---|---|
| Finland | 0.86 |
| Netherlands | 0.85 |
| Iceland | 0.83 |
| Sweden | 0.83 |
| Mexico | 0.82 |
| Country | Percent of consumption of a lacto-ova vegetarian diet, by kg per person |
|---|---|
| India | 0.96 |
| Pakistan | 0.95 |
| Algeria | 0.93 |
| Nepal | 0.93 |
| Ethiopia | 0.93 |