The article Total Meat Intake is Associated with Life Expectancy: A Cross-Sectional Data Analysis of 175 Contemporary Populations by You [*] examined the association between meat intake and life expectancy at a population level based on ecological data published by the United Nations agencies. Their results were
Worldwide, bivariate correlation analyses revealed that meat intake is positively correlated with life expectancies. This relationship remained significant when influences of caloric intake, urbanization, obesity, education and carbohydrate crops were statistically controlled. Stepwise linear regression selected meat intake, not carbohydrate crops, as one of the significant predictors of life expectancy
They also mentioned:
Thirdly, GDP PPP may be a comprehensive life expectancy contributor. For instance, populations with greater GDP PPP may have higher meat affordability, better medical service and better education level. Each factor may contribute to life expectancy in its unique way, but it is impossible to collect all these data and include them as the potential separate confounders in the data analyses to remove their competing effects on life expectancy.
I tried to study the influence of better medical service and better education level.
Method
Supplementary data sources have been procured and are duly cited in the sources section. The analysis was executed using Python in conjunction with statsmodels and Streamlit, and the resultant page is accessible at the following location: https://rcsmit.streamlit.app/?choice=26
At first some scatter plots were made. The original study contained only the first one.
A correlation heatmap with the additional datasources was made
A standard multiple regression is made. Like in the original study the values has been log normalized. After that the data was standardized by Z-score normalization to the selected columns.
First it is done with all the countries like in the study. The results aligned with those presented in Table 3 of the original study.
Subsequently, health index and education index were included in the analysis.
The coefficients obtained were as follows: meat consumption (0.1916), health index (0.3330), and education index (0.5085). These coefficients imply that health and education have a more significant influence on life expectancy than meat consumption.
In contrast to the original study, a weighted linear regression was performed in this analysis. The weight factor is determined by the population size.
The obtained coefficients were as follows: meat consumption (0.0574), health index (17.2756), and education index (0.150622). These coefficients indicate that health and education exert an even more pronounced influence on life expectancy compared to meat consumption.
Conclusion
In conclusion, this study indicates that healthcare quality and education play a more substantial role in influencing life expectancy when compared to meat consumption.
Discussion
However, it’s important to note that multiple linear regression necessitates the assumption that all variables are not correlated. In this case, meat consumption, health index, and education index are related to overall development.
Additionally, the final regression analysis was conducted with a reduced number of countries due to incomplete and outdated data. To enhance the robustness of the findings, it is recommended that this analysis be repeated with more comprehensive and up-to-date data.
Remarks
It is worth mentioning that upon reviewing the authors’ previous work, no other publications on a similar topic were identified. This suggests that this study may have been conducted independently and on a relatively ad hoc basis.
Links
For a critical examination of the methodology employed in this study, please refer to the provided link:
MEAT INTAKE & LONGEVITY, OR FOOD SUPPLY & BIRTH LIFE EXPECTANCY?
Data Sources
Meat consumption etc. : https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8881926/ (appendix 1)
- newborn life expectancy (e(0)), life expectancy at 5 years of life (e(5)) and intakes of meat, and carbohydrate crops, respectively. The established risk factors to life expectancy — caloric intake, urbanization, obesity and education levels — were included as the potential confounders.
Gapminder data set, values from 2018 : https://www.kaggle.com/datasets/albertovidalrod/gapminder-dataset?resource=download
- Country (country). Describes the country name
- Continent (continent). Describes the continent to which the country belongs
- Year (year). Describes the year to which the data belongs
- Life expectancy (life_exp). Describes the life expectancy for a given country in a given year
- Human Development Index (hdi_index). Describes the HDI index value for a given country in a given year
- — — The Human Development Index (HDI) is a statistical composite index of life expectancy, education (mean years of schooling completed and expected years of schooling upon entering the education system), and per capita income indicators, which is used to rank countries into four tiers of human development.
- CO2 emissions per person(co2_consump). Describes the CO2 emissions in tonnes per person for a given country in a given year
- Gross Domestic Product per capita (gdp). Describes the GDP per capita in dollars for a given country in a given year
- % Service workers (services). Describes the the % of service workers for a given country in a given year
Health Efficiency Index : https://web.archive.org/web/20200313135813/https://www.who.int/healthinfo/paper30.pdf
- Factors : 25% for health (DALE), 25% for health inequality, 12.5% for the level of responsiveness, 12.5% for the distribution of responsiveness, and 25% for fairness in financing.
Education: https://ourworldindata.org/human-development-index
- Mean years of schooling estimates the average number of years of total schooling adults aged 25 years and older have received. This data extends back to the year 1870 and is based on the combination of data from Lee and Lee (2016); Barro-Lee (2018); and the UN Development Programme. Fifteen is the projected maximum of this indicator for 2025.
- Expected years of schooling measures the number of years of schooling that a child of school entrance age can expect to receive if the current age-specific enrollment rates persist throughout the child’s life by country. Eighteen is equivalent to achieving a master’s degree in most countries.
- education_index = ((schooling_expected / 18) + (schooling_mean / 15) )/2
Country codes and info: https://countrycode.org/
renesmit