In a previous post (), I discussed the frequently reported lowest-mortality body mass index (BMI), which is about 26. The empirical results reviewed in that post suggest that fat-free mass plays an important role in that context. Keep in mind that this "BMI=26 phenomenon" is often reported in studies of populations from developed countries, which are likely to be relatively sedentary. This is important for the point made in this post.
A lowest-mortality BMI of 26 is somehow at odds with the fact that many healthy and/or long-living populations have much lower BMIs. You can clearly see this in the distribution of BMIs among males in Kitava and Sweden shown in the graph below, from a study by Lindeberg and colleagues (). This distribution is shifted in such a way that would suggest a much lower BMI of lowest-mortality among the Kitavans, assuming a U-curve shape similar to that observed in studies of populations from developed countries ().
Another relevant example comes from the China Study II (see, e.g., ), which is based on data from 8000 adults. The `verage BMI in the China Study II dataset, with data from the 1980s, is approximately 21; for an average weight that is about 116 lbs. That BMI is relatively uniform across Chinese counties, including those with the lowest mortality rates. No county has an average BMI that is 26; not even close. This also supports the idea that Chinese people were, at least during that period, relatively thin.
Now take a look at the graph below, also based on the China Study II dataset, from a previous post (), relating total daily calorie intake with longevity. I should note that the relationship between total daily calorie intake and longevity depicted in this graph is not really statistically significant. Still, the highest longevity seems to be in the second tercile of total daily calorie intake.
Again, the average weight in the dataset is about 116 lbs. A conservative estimate of the number of calories needed to maintain this weight without any physical activity would be about 1740. Add about 700 calories to that, for a reasonable and healthy level of physical activity, and you get 2440 calories needed daily for weight maintenance. That is right in the middle of the second tercile, the one with the highest longevity.
What does this have to do with the lowest-mortality BMI of 26 from studies of samples from developed countries? Populations in these countries are likely to be relatively sedentary, at least on average, in which case a low BMI will be associated with a low total calorie intake. And a low total calorie intake will lead to a low intake of nutrients needed by the body to fight disease.
And don’t think you can fix this problem by consuming lots of vitamin and mineral pills. When I refer here to a higher or lower nutrient intake, I am not talking only about micronutrients, but also about macronutrients (fatty and amino acids) in amounts that are needed by your body. Moreover, important micronutrients, such as fat-soluble vitamins, cannot be properly absorbed without certain macronutrients, such as fat.
Industrial nutrient isolation for supplementation use has not been a very successful long-term strategy for health optimization (). On the other hand, this type of supplementation has indeed been found to have had modest-to-significant success in short-term interventions aimed at correcting acute health problems caused by severe nutritional deficiencies ().
So the "BMI=26 phenomenon" may be a reflection not of a direct effect of high muscularity on health, but of an indirect effect mediated by a high intake of needed nutrients among sedentary folks. This may be so even though the lowest mortality is for the combination of that BMI with a relatively small waist (), which suggests some level of muscularity, but not necessarily serious bodybuilder-level muscularity. High muscularity, of the serious bodybuilder type, is not very common; at least not enough to significantly sway results based on the analysis of large samples.
The combination of a BMI=26 with a relatively small waist is indicative of more muscle and less body fat. Having more muscle and less body fat has an advantage that is rarely discussed. It allows for a higher total calorie intake, and thus a higher nutrient intake, without an unhealthy increase in body fat. Muscle mass increases one's caloric requirement for weight maintenance, more so than body fat. Body fat also increases that caloric requirement, but it also acts like an organ, secreting a number of hormones into the bloodstream, and becoming pro-inflammatory in an unhealthy way above a certain level.
Clearly having a low body fat percentage is associated with lower incidence of degenerative diseases, but it will likely lead to a lower intake of nutrients relative to one’s needs unless other factors are present, e.g., being fairly muscular or physically active. Chronic low nutrient intake tends to get people closer to the afterlife like nothing else ().
In this sense, having a BMI=26 and being relatively sedentary (without being skinny-fat) has an effect that is similar to that of having a BMI=21 and being fairly physically active. Both would lead to consumption of more calories for weight maintenance, and thus more nutrients, as long as nutritious foods are eaten.
Showing posts with label longevity. Show all posts
Showing posts with label longevity. Show all posts
Monday, July 2, 2012
Monday, June 18, 2012
The lowest-mortality BMI: What is its relationship with fat-free mass?
Do overweight folks live longer? It is not uncommon to see graphs like the one below, from the Med Journal Watch blog (), suggesting that, at least as far as body mass index (BMI) is concerned (), overweight folks (25 < BMI < 30) seem to live longer. The graph shows BMI measured at a certain age, and risk of death within a certain time period (e.g., 20 years) following the measurement. The lowest-mortality BMI is about 26, which is in the overweight area of the BMI chart.
Note that mortality risk, relative to the mortality risk of people in the same age group, increases less steeply in response to weight variations as one becomes older. An older person increases the risk of dying to a lesser extent by weighing more or less than does a younger person. This seems to be particularly true for gains in weight.
The table below is from a widely cited 2002 article by Allison and colleagues (), where they describe a study of 10,169 males aged 25-75. Almost all of the participants, ninety-eight percent, were followed up for many years after measurement; a total of 3,722 deaths were recorded.
Take a look at the two numbers circled in red. The one on the left is the lowest-mortality BMI not adjusting for fat mass or fat-free mass: a reasonably high 27.4. The one on the right is the lowest-mortality BMI adjusting for fat mass and fat-free mass: a much lower 21.6.
I know this may sound confusing, but due to possible statistical distortions this does not mean that you should try to bring your BMI to 21.6 if you want to reduce your risk of dying. What this means is that fat mass and fat-free mass matter. Moreover, all of the participants in this study were men. The authors concluded that: “…marked leanness (as opposed to thinness) has beneficial effects.”
Then we have an interesting 2003 article by Bigaard and colleagues () reporting on a study of 27,178 men and 29,875 women born in Denmark, 50 to 64 years of age. The table below summarizes deaths in this study, grouping them by BMI and waist circumference.
These are raw numbers; no complex statistics here. Circled in green is the area with samples that appear to be large enough to avoid “funny” results. Circled in red are the lowest-mortality percentages; I left out the 0.8 percentage because it is based on a very small sample.
As you can see, they refer to men and women with BMIs in the 25-29.9 range (overweight), but with waist circumferences in the lower-middle range: 90-96 cm for men and 74-82 cm for women; or approximately 35-38 inches for men and 29-32 inches for women.
Women with BMIs in the 18.5-24.9 range (normal) and the same or lower waists also died in small numbers. Underweight men and women had the highest mortality percdntages.
A relatively small waist (not a wasp waist), together with a normal or high BMI, is an indication of more fat-free mass, which is retained together with some body fat. It is also an indication of less visceral body fat accumulation.
Note that mortality risk, relative to the mortality risk of people in the same age group, increases less steeply in response to weight variations as one becomes older. An older person increases the risk of dying to a lesser extent by weighing more or less than does a younger person. This seems to be particularly true for gains in weight.
The table below is from a widely cited 2002 article by Allison and colleagues (), where they describe a study of 10,169 males aged 25-75. Almost all of the participants, ninety-eight percent, were followed up for many years after measurement; a total of 3,722 deaths were recorded.
Take a look at the two numbers circled in red. The one on the left is the lowest-mortality BMI not adjusting for fat mass or fat-free mass: a reasonably high 27.4. The one on the right is the lowest-mortality BMI adjusting for fat mass and fat-free mass: a much lower 21.6.
I know this may sound confusing, but due to possible statistical distortions this does not mean that you should try to bring your BMI to 21.6 if you want to reduce your risk of dying. What this means is that fat mass and fat-free mass matter. Moreover, all of the participants in this study were men. The authors concluded that: “…marked leanness (as opposed to thinness) has beneficial effects.”
Then we have an interesting 2003 article by Bigaard and colleagues () reporting on a study of 27,178 men and 29,875 women born in Denmark, 50 to 64 years of age. The table below summarizes deaths in this study, grouping them by BMI and waist circumference.
These are raw numbers; no complex statistics here. Circled in green is the area with samples that appear to be large enough to avoid “funny” results. Circled in red are the lowest-mortality percentages; I left out the 0.8 percentage because it is based on a very small sample.
As you can see, they refer to men and women with BMIs in the 25-29.9 range (overweight), but with waist circumferences in the lower-middle range: 90-96 cm for men and 74-82 cm for women; or approximately 35-38 inches for men and 29-32 inches for women.
Women with BMIs in the 18.5-24.9 range (normal) and the same or lower waists also died in small numbers. Underweight men and women had the highest mortality percdntages.
A relatively small waist (not a wasp waist), together with a normal or high BMI, is an indication of more fat-free mass, which is retained together with some body fat. It is also an indication of less visceral body fat accumulation.
Monday, March 19, 2012
The 2012 red meat-mortality study (Arch Intern Med): The data suggests that red meat is protective
I am not a big fan of using arguments such as “food questionnaires are unreliable” and “observational studies are worthless” to completely dismiss a study. There are many reasons for this. One of them is that, when people misreport certain diet and lifestyle patterns, but do that consistently (i.e., everybody underreports food intake), the biasing effect on coefficients of association is minor. Measurement errors may remain for this or other reasons, but regression methods (linear and nonlinear) assume the existence of such errors, and are designed to yield robust coefficients in their presence. Besides, for me to use these types of arguments would be hypocritical, since I myself have done several analyses on the China Study data (), and built what I think are valid arguments based on those analyses.
My approach is: Let us look at the data, any data, carefully, using appropriate analysis tools, and see what it tells us; maybe we will find evidence of measurement errors distorting the results and leading to mistaken conclusions, or maybe not. With this in mind, let us take a look at the top part of Table 3 of the most recent (published online in March 2012) study looking at the relationship between red meat consumption and mortality, authored by Pan et al. (Frank B. Hu is the senior author) and published in the prestigious Archives of Internal Medicine (). This is a prominent journal, with an average of over 270 citations per article according to Google Scholar. The study has received much media attention recently.
Take a look at the area highlighted in red, focusing on data from the Health Professionals sample. That is the multivariate-adjusted cardiovascular mortality rate, listed as a normalized percentage, in the highest quintile (Q5) of red meat consumption from the Health Professionals sample. The non-adjusted percentages are 1.4 percent mortality in Q5 and 1.13 in Q1 (from Table 1 of the same article); so the multivariate adjustment-normalization changed the values of the percentages somewhat, but not much. The highlighted 1.35 number suggests that for each group of 100 people who consumed a lot of red meat (Q5), when compared with a group of 100 people who consumed little red meat (Q1), there were on average 0.35 more deaths over the same period of time (more than 20 years).
The heavy red meat eaters in Q5 consumed 972.37 percent more red meat than those in Q1. This is calculated with data from Table 1 of the same article, as: (2.36-0.22)/0.22. In Q5, the 2.36 number refers to the number of servings of red meat per day, with each serving being approximately 84 g. So the heavy red meat eaters ate approximately 198 g per day (a bit less than 0.5 lb), while the light red meat eaters ate about 18 g per day. In other words, the heavy red meat eaters ate 9.7237 times more, or 972.37 percent more, red meat.
So, just to be clear, even though the folks in Q5 consumed 972.37 percent more red meat than the folks in Q1, in each matched group of 100 you would not find a single additional death over the same time period. If you looked at matched groups of 1,000 individuals, you would find 3 more deaths among the heavy red meat eaters. The same general pattern, of a minute difference, repeats itself throughout Table 3. As you can see, all of the reported mortality ratios are 1-point-something. In fact, this same pattern repeats itself in all mortality tables (all-cause, cardiovascular, cancer). This is all based on a multivariate analysis that according to the authors controlled for a large number of variables, including baseline history of diabetes.
Interestingly, looking at data from the same sample (Health Professionals), the incidence of diabetes is 75 percent higher in Q5 than in Q1. The same is true for the second sample (Nurses Health), where the Q5-Q1 difference in incidence of diabetes is even greater - 81 percent. This caught my eye, being diabetes such a prototypical “disease of affluence”. So I entered the whole data reported in the article into HCE () and WarpPLS (), and conducted some analyses. The graphs below are from HCE. The data includes both samples – Health Professionals and Nurses Health.
HCE calculates bivariate correlations, and so does WarpPLS. But WarpPLS stores numbers with a higher level of precision, so I used WarpPLS for calculating coefficients of association, including correlations. I also double-checked the numbers with other software, just in case (e.g., SPSS and MATLAB). Here are the correlations calculated by WarpPLS, which refer to the graphs above: 0.030 for red meat intake and mortality; 0.607 for diabetes and mortality; and 0.910 for food intake and diabetes. Yes, you read it right, the correlation between red meat intake and mortality is a very low and non-significant 0.030 in this dataset. Not a big surprise when you look at the related HCE graph, with the line going up and down almost at random. Note that I included the quintiles data from both the Health Professionals and Nurses Health samples in one dataset.
Those folks in Q5 had a much higher incidence of diabetes, and yet the increase in mortality for them was significantly lower, in percentage terms. A key difference between Q5 and Q1 being what? The Q5 folks ate a lot more red meat. This looks suspiciously suggestive of a finding that I came across before, based on an analysis of the China Study II data (). The finding was that animal food consumption (and red meat is an animal food) was protective, actually reducing the negative effect of wheat flour consumption on mortality. That analysis actually suggested that wheat flour consumption may not be so bad if you eat 221 g or more of animal food daily.
So, I built the model below in WarpPLS, where red meat intake (RedMeat) is hypothesized to moderate the relationship between diabetes incidence (Diabetes) and mortality (Mort). Below I am also including the graphs for the direct and moderating effects; the data is standardized, which reduces estimation error, particularly in moderating effects estimation. I used a standard linear algorithm for the calculation of the path coefficients (betas next to the arrows) and jackknifing for the calculation of the P values (confidence = 1 – P value). Jackknifing is a resampling technique that does not require multivariate normality and that tends to work well with small samples; as is the case with nonparametric techniques in general.
The direct effect of diabetes on mortality is positive (0.68) and almost statistically significant at the P < 0.05 level (confidence of 94 percent), which is noteworthy because the sample size here is so small – only 10 data points, 5 quintiles from the Health Professionals sample and 5 from the Nurses Health sample. The moderating effect is negative (-0.11), but not statistically significant (confidence of 61 percent). In the moderating effect graphs (shown side-by-side), this negative moderation is indicated by a slightly less steep inclination of the regression line for the graph on the right, which refers to high red meat intake. A less steep inclination means a less strong relationship between diabetes and mortality – among the folks who ate the most red meat.
Not too surprisingly, at least to me, the results above suggest that red meat per se may well be protective. Although we should consider a least two other possibilities. One is that red meat intake is a marker for consumption of some other things, possibly present in animal foods, that are protective - e.g., choline and vitamin K2. The other possibility is that red meat is protective in part by displacing other less healthy foods. Perhaps what we are seeing here is a combination of these.
Whatever the reason may be, red meat consumption seems to actually lessen the effect of diabetes on mortality in this sample. That is, according to this data, the more red meat is consumed, the fewer people die from diabetes. The protective effect might have been stronger if the participants had eaten more red meat, or more animal foods containing the protective factors; recall that the threshold for protection in the China Study II data was consumption of 221 g or more of animal food daily (). Having said that, it is also important to note that, if you eat excess calories to the point of becoming obese, from red meat or any other sources, your risk of developing diabetes will go up – as the earlier HCE graph relating food intake and diabetes implies.
Please keep in mind that this post is the result of a quick analysis of secondary data reported in a journal article, and its conclusions may be wrong, even though I did my best not to make any mistake (e.g., mistyping data from the article). The authors likely spent months, if not more, in their study; and have the support of one of the premier research universities in the world. Still, this post raises serious questions. I say this respectfully, as the authors did seem to try their best to control for all possible confounders.
I should also say that the moderating effect I uncovered is admittedly a fairly weak effect on this small sample and not statistically significant. But its magnitude is apparently greater than the reported effects of red meat on mortality, which are not only minute but may well be statistical artifacts. The Cox proportional hazards analysis employed in the study, which is commonly used in epidemiology, is nothing more than a sophisticated ANCOVA; it is a semi-parametric version of a special case of the broader analysis method automated by WarpPLS.
Finally, I could not control for confounders because, given the small sample, inclusion of confounders (e.g., smoking) leads to massive collinearity. WarpPLS calculates collinearity estimates automatically, and is particularly thorough at doing that (calculating them at multiple levels), so there is no way to ignore them. Collinearity can severely distort results, as pointed out in a YouTube video on WarpPLS (). Collinearity can even lead to changes in the signs of coefficients of association, in the context of multivariate analyses - e.g., a positive association appears to be negative. The authors have the original data – a much, much larger sample - which makes it much easier to deal with collinearity.
Moderating effects analyses () – we need more of that in epidemiological research eh?
My approach is: Let us look at the data, any data, carefully, using appropriate analysis tools, and see what it tells us; maybe we will find evidence of measurement errors distorting the results and leading to mistaken conclusions, or maybe not. With this in mind, let us take a look at the top part of Table 3 of the most recent (published online in March 2012) study looking at the relationship between red meat consumption and mortality, authored by Pan et al. (Frank B. Hu is the senior author) and published in the prestigious Archives of Internal Medicine (). This is a prominent journal, with an average of over 270 citations per article according to Google Scholar. The study has received much media attention recently.
Take a look at the area highlighted in red, focusing on data from the Health Professionals sample. That is the multivariate-adjusted cardiovascular mortality rate, listed as a normalized percentage, in the highest quintile (Q5) of red meat consumption from the Health Professionals sample. The non-adjusted percentages are 1.4 percent mortality in Q5 and 1.13 in Q1 (from Table 1 of the same article); so the multivariate adjustment-normalization changed the values of the percentages somewhat, but not much. The highlighted 1.35 number suggests that for each group of 100 people who consumed a lot of red meat (Q5), when compared with a group of 100 people who consumed little red meat (Q1), there were on average 0.35 more deaths over the same period of time (more than 20 years).
The heavy red meat eaters in Q5 consumed 972.37 percent more red meat than those in Q1. This is calculated with data from Table 1 of the same article, as: (2.36-0.22)/0.22. In Q5, the 2.36 number refers to the number of servings of red meat per day, with each serving being approximately 84 g. So the heavy red meat eaters ate approximately 198 g per day (a bit less than 0.5 lb), while the light red meat eaters ate about 18 g per day. In other words, the heavy red meat eaters ate 9.7237 times more, or 972.37 percent more, red meat.
So, just to be clear, even though the folks in Q5 consumed 972.37 percent more red meat than the folks in Q1, in each matched group of 100 you would not find a single additional death over the same time period. If you looked at matched groups of 1,000 individuals, you would find 3 more deaths among the heavy red meat eaters. The same general pattern, of a minute difference, repeats itself throughout Table 3. As you can see, all of the reported mortality ratios are 1-point-something. In fact, this same pattern repeats itself in all mortality tables (all-cause, cardiovascular, cancer). This is all based on a multivariate analysis that according to the authors controlled for a large number of variables, including baseline history of diabetes.
Interestingly, looking at data from the same sample (Health Professionals), the incidence of diabetes is 75 percent higher in Q5 than in Q1. The same is true for the second sample (Nurses Health), where the Q5-Q1 difference in incidence of diabetes is even greater - 81 percent. This caught my eye, being diabetes such a prototypical “disease of affluence”. So I entered the whole data reported in the article into HCE () and WarpPLS (), and conducted some analyses. The graphs below are from HCE. The data includes both samples – Health Professionals and Nurses Health.
HCE calculates bivariate correlations, and so does WarpPLS. But WarpPLS stores numbers with a higher level of precision, so I used WarpPLS for calculating coefficients of association, including correlations. I also double-checked the numbers with other software, just in case (e.g., SPSS and MATLAB). Here are the correlations calculated by WarpPLS, which refer to the graphs above: 0.030 for red meat intake and mortality; 0.607 for diabetes and mortality; and 0.910 for food intake and diabetes. Yes, you read it right, the correlation between red meat intake and mortality is a very low and non-significant 0.030 in this dataset. Not a big surprise when you look at the related HCE graph, with the line going up and down almost at random. Note that I included the quintiles data from both the Health Professionals and Nurses Health samples in one dataset.
Those folks in Q5 had a much higher incidence of diabetes, and yet the increase in mortality for them was significantly lower, in percentage terms. A key difference between Q5 and Q1 being what? The Q5 folks ate a lot more red meat. This looks suspiciously suggestive of a finding that I came across before, based on an analysis of the China Study II data (). The finding was that animal food consumption (and red meat is an animal food) was protective, actually reducing the negative effect of wheat flour consumption on mortality. That analysis actually suggested that wheat flour consumption may not be so bad if you eat 221 g or more of animal food daily.
So, I built the model below in WarpPLS, where red meat intake (RedMeat) is hypothesized to moderate the relationship between diabetes incidence (Diabetes) and mortality (Mort). Below I am also including the graphs for the direct and moderating effects; the data is standardized, which reduces estimation error, particularly in moderating effects estimation. I used a standard linear algorithm for the calculation of the path coefficients (betas next to the arrows) and jackknifing for the calculation of the P values (confidence = 1 – P value). Jackknifing is a resampling technique that does not require multivariate normality and that tends to work well with small samples; as is the case with nonparametric techniques in general.
The direct effect of diabetes on mortality is positive (0.68) and almost statistically significant at the P < 0.05 level (confidence of 94 percent), which is noteworthy because the sample size here is so small – only 10 data points, 5 quintiles from the Health Professionals sample and 5 from the Nurses Health sample. The moderating effect is negative (-0.11), but not statistically significant (confidence of 61 percent). In the moderating effect graphs (shown side-by-side), this negative moderation is indicated by a slightly less steep inclination of the regression line for the graph on the right, which refers to high red meat intake. A less steep inclination means a less strong relationship between diabetes and mortality – among the folks who ate the most red meat.
Not too surprisingly, at least to me, the results above suggest that red meat per se may well be protective. Although we should consider a least two other possibilities. One is that red meat intake is a marker for consumption of some other things, possibly present in animal foods, that are protective - e.g., choline and vitamin K2. The other possibility is that red meat is protective in part by displacing other less healthy foods. Perhaps what we are seeing here is a combination of these.
Whatever the reason may be, red meat consumption seems to actually lessen the effect of diabetes on mortality in this sample. That is, according to this data, the more red meat is consumed, the fewer people die from diabetes. The protective effect might have been stronger if the participants had eaten more red meat, or more animal foods containing the protective factors; recall that the threshold for protection in the China Study II data was consumption of 221 g or more of animal food daily (). Having said that, it is also important to note that, if you eat excess calories to the point of becoming obese, from red meat or any other sources, your risk of developing diabetes will go up – as the earlier HCE graph relating food intake and diabetes implies.
Please keep in mind that this post is the result of a quick analysis of secondary data reported in a journal article, and its conclusions may be wrong, even though I did my best not to make any mistake (e.g., mistyping data from the article). The authors likely spent months, if not more, in their study; and have the support of one of the premier research universities in the world. Still, this post raises serious questions. I say this respectfully, as the authors did seem to try their best to control for all possible confounders.
I should also say that the moderating effect I uncovered is admittedly a fairly weak effect on this small sample and not statistically significant. But its magnitude is apparently greater than the reported effects of red meat on mortality, which are not only minute but may well be statistical artifacts. The Cox proportional hazards analysis employed in the study, which is commonly used in epidemiology, is nothing more than a sophisticated ANCOVA; it is a semi-parametric version of a special case of the broader analysis method automated by WarpPLS.
Finally, I could not control for confounders because, given the small sample, inclusion of confounders (e.g., smoking) leads to massive collinearity. WarpPLS calculates collinearity estimates automatically, and is particularly thorough at doing that (calculating them at multiple levels), so there is no way to ignore them. Collinearity can severely distort results, as pointed out in a YouTube video on WarpPLS (). Collinearity can even lead to changes in the signs of coefficients of association, in the context of multivariate analyses - e.g., a positive association appears to be negative. The authors have the original data – a much, much larger sample - which makes it much easier to deal with collinearity.
Moderating effects analyses () – we need more of that in epidemiological research eh?
Monday, February 20, 2012
The “pork paradox”? National pork consumption and obesity
In my previous post () I discussed some country data linking pork consumption and health, analyzed with WarpPLS (). One of the datasets used, the most complete, contained data from Nationmaster.com () for the following countries: Australia, Brazil, Canada, China, Denmark, France, Germany, Hong Kong, Hungary, Japan, Mexico, Poland, Russia, Singapore, Spain, Sweden, United Kingdom, and United States. That previous post also addressed a study by Bridges (), based on country-level data, suggesting that pork consumption may cause liver disease.
In this post we continue that analysis, but with a much more complex model containing the following country variables: wealth (PPP-adjusted GNP/person), pork consumption (lbs/person/year), alcohol consumption (liters/person/year), obesity (% of population), and life expectancy (years). The model and results, generated by WarpPLS, are shown on the figure below. (See notes at the end of this post.) These results are only for direct effects.
WarpPLS also calculates total effects, which are the effects of each variable on any other variable to which it is linked directly and/or indirectly. Two variables may be linked indirectly, through various paths, even if they are not linked directly (i.e., have an arrow directly connecting them). Another set of outputs generated by the software are effect sizes, which are calculated as Cohen’s f-squared coefficients. The figure below shows the total effects table. The values underlined in red are for total effects that are both statistically significant and also above the effect size threshold recommended by Cohen to be considered relevant (f-squared > 0.02).
As I predicted in my previous post, wealth is positively associated with pork consumption. So is alcohol consumption, and more strongly than wealth; which is consistent with a study by Jeanneret and colleagues showing a strong association between alcohol consumption and protein rich diets (). The inclusion of wealth in the model, compared with the model without wealth in the previous post, renders the direct and total effects of alcohol and pork consumption on life expectancy statistically indistinguishable from zero. (This often happens when a confounder is added to a model.)
Pork consumption is negatively associated with obesity, which is interesting. So is alcohol consumption, but much less strongly than pork consumption. This does not mean that if you eat 20 doughnuts every day, together with 1 lb of pork, you are not going to become obese. What this does suggest is that maybe countries where pork is consumed more heavily are somewhat more resistant to obesity. Here it should be noted that pork is very popular in Asian countries, which are becoming increasingly wealthy, but without the widespread obesity that we see in the USA.
But it is not the inclusion of Asian countries in the dataset that paints such a positive picture for pork consumption vis-à-vis obesity, and even weakens the association between wealth and obesity so much as to make it statistically non-significant. Denmark is a wealthy country that has very low levels of obesity. And it happens to have the highest level of pork consumption in the whole dataset: 142.6 lbs/person/year. So we are not talking about an “Asian paradox” here.
More like a “pork paradox”.
Finally, as far as life expectancy is concerned, the key factors seem to be wealth and obesity. Wealth has a major positive effect on life expectancy, while obesity has a much weaker negative effect. Well, access to sanitation, medical services, and other amenities of civilization, still trumps obesity in terms of prolonging life; however miserable life may turn out to be. The competing effects of these two variables (i.e., wealth and obesity) were taken into consideration, or controlled for, in the calculation of total effects and effect sizes.
The fact that pork consumption is negatively associated with obesity goes somewhat against the idea that pork is inherently unhealthy; even though pork certainly can cause disease if not properly prepared and/or cooked, which is true for many other plant and animal foods. The possible connection with liver problems, alluded to in the previous post, is particularly suspicious in light of these results. Liver diseases often impair that organ’s ability to make glycogen based on carbohydrates and protein; that is, liver diseases frequently lead to liver insulin resistance. And obesity frequently follows from liver insulin resistance.
Given that pork consumption appears to be negatively associated with obesity, it would be surprising if it was causing widespread liver disease, unless its relationship with liver dhsease was found to be nonlinear. (Alcohol consumption seems to be nonlinearly associated with liver disease.) Still, most studies that suggest the existence of a causal link between pork consumption and liver disease, like Bridges’s (), hint at a linear and dose-dependent relationship.
Notes
- Country-level data is inherently problematic, particularly when simple models are used (e.g., a model with only two variables). There are just too many possible confounders that may lead to the appearance of causal associations.
- More complex models ameliorate the above situation somewhat, but bump into another problem associated with country-level data – small sample sizes. We used data from 18 countries in this analysis, which is more than in the Bridges study. Still, the effective sample size here (N=18) is awfully small.
- There were some missing values in this dataset, which were handled by WarpPLS employing the most widely used approach in these cases – i.e., by replacing the missing values with the mean of each column. The percentages of missing values per variable (i.e., column) were: alcohol consumption: 27.78%; life expectancy: 5.56%; and obesity: 33.33%.
In this post we continue that analysis, but with a much more complex model containing the following country variables: wealth (PPP-adjusted GNP/person), pork consumption (lbs/person/year), alcohol consumption (liters/person/year), obesity (% of population), and life expectancy (years). The model and results, generated by WarpPLS, are shown on the figure below. (See notes at the end of this post.) These results are only for direct effects.
WarpPLS also calculates total effects, which are the effects of each variable on any other variable to which it is linked directly and/or indirectly. Two variables may be linked indirectly, through various paths, even if they are not linked directly (i.e., have an arrow directly connecting them). Another set of outputs generated by the software are effect sizes, which are calculated as Cohen’s f-squared coefficients. The figure below shows the total effects table. The values underlined in red are for total effects that are both statistically significant and also above the effect size threshold recommended by Cohen to be considered relevant (f-squared > 0.02).
As I predicted in my previous post, wealth is positively associated with pork consumption. So is alcohol consumption, and more strongly than wealth; which is consistent with a study by Jeanneret and colleagues showing a strong association between alcohol consumption and protein rich diets (). The inclusion of wealth in the model, compared with the model without wealth in the previous post, renders the direct and total effects of alcohol and pork consumption on life expectancy statistically indistinguishable from zero. (This often happens when a confounder is added to a model.)
Pork consumption is negatively associated with obesity, which is interesting. So is alcohol consumption, but much less strongly than pork consumption. This does not mean that if you eat 20 doughnuts every day, together with 1 lb of pork, you are not going to become obese. What this does suggest is that maybe countries where pork is consumed more heavily are somewhat more resistant to obesity. Here it should be noted that pork is very popular in Asian countries, which are becoming increasingly wealthy, but without the widespread obesity that we see in the USA.
But it is not the inclusion of Asian countries in the dataset that paints such a positive picture for pork consumption vis-à-vis obesity, and even weakens the association between wealth and obesity so much as to make it statistically non-significant. Denmark is a wealthy country that has very low levels of obesity. And it happens to have the highest level of pork consumption in the whole dataset: 142.6 lbs/person/year. So we are not talking about an “Asian paradox” here.
More like a “pork paradox”.
Finally, as far as life expectancy is concerned, the key factors seem to be wealth and obesity. Wealth has a major positive effect on life expectancy, while obesity has a much weaker negative effect. Well, access to sanitation, medical services, and other amenities of civilization, still trumps obesity in terms of prolonging life; however miserable life may turn out to be. The competing effects of these two variables (i.e., wealth and obesity) were taken into consideration, or controlled for, in the calculation of total effects and effect sizes.
The fact that pork consumption is negatively associated with obesity goes somewhat against the idea that pork is inherently unhealthy; even though pork certainly can cause disease if not properly prepared and/or cooked, which is true for many other plant and animal foods. The possible connection with liver problems, alluded to in the previous post, is particularly suspicious in light of these results. Liver diseases often impair that organ’s ability to make glycogen based on carbohydrates and protein; that is, liver diseases frequently lead to liver insulin resistance. And obesity frequently follows from liver insulin resistance.
Given that pork consumption appears to be negatively associated with obesity, it would be surprising if it was causing widespread liver disease, unless its relationship with liver dhsease was found to be nonlinear. (Alcohol consumption seems to be nonlinearly associated with liver disease.) Still, most studies that suggest the existence of a causal link between pork consumption and liver disease, like Bridges’s (), hint at a linear and dose-dependent relationship.
Notes
- Country-level data is inherently problematic, particularly when simple models are used (e.g., a model with only two variables). There are just too many possible confounders that may lead to the appearance of causal associations.
- More complex models ameliorate the above situation somewhat, but bump into another problem associated with country-level data – small sample sizes. We used data from 18 countries in this analysis, which is more than in the Bridges study. Still, the effective sample size here (N=18) is awfully small.
- There were some missing values in this dataset, which were handled by WarpPLS employing the most widely used approach in these cases – i.e., by replacing the missing values with the mean of each column. The percentages of missing values per variable (i.e., column) were: alcohol consumption: 27.78%; life expectancy: 5.56%; and obesity: 33.33%.
Monday, February 13, 2012
Does pork consumption cause cirrhosis? Perhaps, if people become obese from eating pork
The idea that pork consumption may cause cirrhosis has been around for a while. A fairly widely cited 1985 study by Nanji and French () provides one of the strongest indictments of pork: “In countries with low alcohol consumption, no correlation was obtained between alcohol consumption and cirrhosis. However, a significant correlation was obtained between cirrhosis and pork.”
Recently Paul Jaminet wrote a blog post on the possible link between pork consumption and cirrhosis (). Paul should be commended for bringing this topic to the fore, as the implications are far-reaching and very serious. One of the key studies mentioned in Paul’s post is a 2009 article by Bridges (), from which the graphs below were taken.
The graphs above show a correlation between cirrhosis and alcohol consumption of 0.71, and a correlation between cirrhosis and pork consumption of 0.83. That is, the correlation between cirrhosis and pork consumption is the stronger of the two! Combining this with the Nanji and French study, we have evidence that: (a) in countries with low alcohol consumption we can find a significant correlation between cirrhosis and pork consumption; and (b) in countries where both alcohol and pork are consumed, pork consumption has the strongest correlation with cirrhosis.
Do we need anything else to ban pork from our diets? Yes, we do, as there is more to this story.
Clearly alcohol and pork consumption are correlated as well, as we can see from the graphs above. That is, countries where alcohol is consumed more heavily also tend to have higher levels of pork consumption. If alcohol and pork consumption are correlated, then a multivariate analysis of their effects should be conducted, as one of the hypothesized effects (of alcohol or pork) on cirrhosis may even disappear after controlling for the other effect.
I created a dataset, as best as I could, based on the graphs from the Bridges article. (I could not get the data online.) I then entered it into WarpPLS (). I wanted to run a moderating effect analysis, which is a form of nonlinear multivariate analysis. This is important, because the association between alcohol consumption and disease in general is well known to be nonlinear.
In fact, the relationship between alcohol consumption and disease is often used as a classic example of hormesis (), and its characteristic J-curve shape. Since correlation is a measure of linear association, the lower correlation between alcohol consumption and cirrhosis, when compared with pork consumption, may be just a “mirage of linearity”. In multivariate analyses, this mirage of linearity may lead to what are known as type I and II errors, at the same time ().
I should note that the Bridges study did something akin to a moderating effect analysis; through an analysis of the interaction between alcohol and pork consumption. However, in that analysis the values of the variables that were multiplied to create a “dummy” interaction variable were on their original scales, which can be a major source of bias. A more advisable way to conduct an interaction effect analysis is to first make the variables dimensionless, by standardizing them, and then creating a dummy interaction variable as a product of the two variables. That is what WarpPLS does for moderating effects’ estimation.
One more detour, leading to an important implication, and then we will get to the results. In a 1988 article, Jeanneret and colleagues show evidence of a strong and possibly causal association between alcohol consumption and protein-rich diets (). One possible implication of this is that in countries where pork is a dietary staple, like Denmark and Germany, alcohol consumption should be strongly and causally associated with pork consumption. (I guess Anthony Bordain would agree with this eh?)
Below are the results of a multivariate analysis on a model that incorporates the above implication, by including a link between alcohol and pork consumption. The model also explores the role of pork consumption as a moderator of the relationship between alcohol and cirrhosis, as well as the direct effect of pork consumption on cirrhosis. Finally, the total effects of alcohol and pork consumption on cirrhosis are also investigated; they are shown on the left.
The total effects are both statistically significant, with the total effect of alcohol consumption being 94 percent stronger than the total effect of pork consumption on cirrhosis. Looking at the model, alcohol consumption is strongly associated with pork consumption (which is consistent with Jeanneret and colleagues’s study). Alcohol consumption is also strongly associated with cirrhosis, through a direct effect; much more so than pork. Finally, pork consumption seems to strengthen the relationship between alcohol consumption and cirrhosis (the moderating effect).
As we can see the relationship between pork consumption and cirrhosis is still there, in moderating and direct effects, even though it seems to be a lot weaker than that between alcohol consumption and cirrhosis. Why does pork seem to influence cirrhosis at all in this dataset?
Well, there is another factor that is strongly associated with cirrhosis, and that is obesity (). In fact, obesity is associated with just about any major disease, including various types of cancer ().
And in countries where pork is a dietary staple, isn’t it reasonable to assume that pork consumption will play a role in obesity? Often folks who consume a lot of addictive industrial foods (e.g., bread, candy, regular sodas) also eat plenty of foods with saturated fat; and the latter end up showing up in disease statistics, misleadingly supporting the lipid hypothesis. The phenomenon involving pork and cirrhosis may well be similar.
But you may find the above results and argument not convincing enough. Maybe you want to see some evidence that pork is actually good for one’s health. The results above suggest that it may not be bad at all, if you buy into the obesity angle, but not that it can be good.
So I downloaded the most recent data from Nationmaster.com () on the following variables: pork consumption, alcohol consumption, and life expectancy. The list of countries was a bit larger than and different from that in the Bridges study; the following countries were included: Australia, Brazil, Canada, China, Denmark, France, Germany, Hong Kong, Hungary, Japan, Mexico, Poland, Russia, Singapore, Spain, Sweden, United Kingdom, and United States. Below are the results of a simple multivariate analysis with WarpPLS.
As with the Bridges dataset, there is a strong multivariate association between alcohol and pork consumption (0.43). The multivariate association between alcohol consumption and life expectancy is negative (-0.14). The multivariate association between pork consumption and life expectancy is positive (0.36). Neither association is statistically significant, although the association involving pork consumption gets close to significance with a P=0.11 (a confidence level of 89 percent; calculated through jackknifing, a nonparametric technique). The graphs show the plots for the associations and the best-fitting lines; the blue dashed arrows indicate the multivariate associations to which the graphs refer. So, in this second dataset from Nationmaster.com, the more pork is consumed in a country, the longer is the life expectancy in that country.
In other words, for each 1 standard deviation variation in pork consumption, there is a 0.36 standard deviation variation in life expectancy, after we control for alcohol consumption. The standard deviation for pork consumption is 36.281 lbs/person/year, or 45.087 g/person/day; for life expectancy, it is 4.677 years. Working the numbers a bit more, the results above suggest that each extra gram of pork consumed per person per day is associated with approximately 13 additional days of overall life expectancy in a country! This is calculated as: 4.677/45.087*0.36*365 = 13.630.
Does this prove that eating pork will make you live longer? No single study will “prove” something like that. Pork consumption is also likely a marker for wealth in a country; and wealth is strongly and positively associated with life expectancy at the country level. Moreover, when you aggregate dietary and disease incidence data by country, often the statistical effects are caused by those people in the dietary extremes (e.g., alcohol abuse, not moderate consumption). Finally, if people avoid death from certain diseases, they will die in higher quantities from other diseases, which may bias statistical results toward what may look like a higher incidence of those other diseases.
What the results summarized in this post do suggest is that pork consumption may not be a problem at all, unless you become obese from eating it. How do you get obese from eating pork? Eating it together with industrial foods that are addictive would probably help.
Recently Paul Jaminet wrote a blog post on the possible link between pork consumption and cirrhosis (). Paul should be commended for bringing this topic to the fore, as the implications are far-reaching and very serious. One of the key studies mentioned in Paul’s post is a 2009 article by Bridges (), from which the graphs below were taken.
The graphs above show a correlation between cirrhosis and alcohol consumption of 0.71, and a correlation between cirrhosis and pork consumption of 0.83. That is, the correlation between cirrhosis and pork consumption is the stronger of the two! Combining this with the Nanji and French study, we have evidence that: (a) in countries with low alcohol consumption we can find a significant correlation between cirrhosis and pork consumption; and (b) in countries where both alcohol and pork are consumed, pork consumption has the strongest correlation with cirrhosis.
Do we need anything else to ban pork from our diets? Yes, we do, as there is more to this story.
Clearly alcohol and pork consumption are correlated as well, as we can see from the graphs above. That is, countries where alcohol is consumed more heavily also tend to have higher levels of pork consumption. If alcohol and pork consumption are correlated, then a multivariate analysis of their effects should be conducted, as one of the hypothesized effects (of alcohol or pork) on cirrhosis may even disappear after controlling for the other effect.
I created a dataset, as best as I could, based on the graphs from the Bridges article. (I could not get the data online.) I then entered it into WarpPLS (). I wanted to run a moderating effect analysis, which is a form of nonlinear multivariate analysis. This is important, because the association between alcohol consumption and disease in general is well known to be nonlinear.
In fact, the relationship between alcohol consumption and disease is often used as a classic example of hormesis (), and its characteristic J-curve shape. Since correlation is a measure of linear association, the lower correlation between alcohol consumption and cirrhosis, when compared with pork consumption, may be just a “mirage of linearity”. In multivariate analyses, this mirage of linearity may lead to what are known as type I and II errors, at the same time ().
I should note that the Bridges study did something akin to a moderating effect analysis; through an analysis of the interaction between alcohol and pork consumption. However, in that analysis the values of the variables that were multiplied to create a “dummy” interaction variable were on their original scales, which can be a major source of bias. A more advisable way to conduct an interaction effect analysis is to first make the variables dimensionless, by standardizing them, and then creating a dummy interaction variable as a product of the two variables. That is what WarpPLS does for moderating effects’ estimation.
One more detour, leading to an important implication, and then we will get to the results. In a 1988 article, Jeanneret and colleagues show evidence of a strong and possibly causal association between alcohol consumption and protein-rich diets (). One possible implication of this is that in countries where pork is a dietary staple, like Denmark and Germany, alcohol consumption should be strongly and causally associated with pork consumption. (I guess Anthony Bordain would agree with this eh?)
Below are the results of a multivariate analysis on a model that incorporates the above implication, by including a link between alcohol and pork consumption. The model also explores the role of pork consumption as a moderator of the relationship between alcohol and cirrhosis, as well as the direct effect of pork consumption on cirrhosis. Finally, the total effects of alcohol and pork consumption on cirrhosis are also investigated; they are shown on the left.
The total effects are both statistically significant, with the total effect of alcohol consumption being 94 percent stronger than the total effect of pork consumption on cirrhosis. Looking at the model, alcohol consumption is strongly associated with pork consumption (which is consistent with Jeanneret and colleagues’s study). Alcohol consumption is also strongly associated with cirrhosis, through a direct effect; much more so than pork. Finally, pork consumption seems to strengthen the relationship between alcohol consumption and cirrhosis (the moderating effect).
As we can see the relationship between pork consumption and cirrhosis is still there, in moderating and direct effects, even though it seems to be a lot weaker than that between alcohol consumption and cirrhosis. Why does pork seem to influence cirrhosis at all in this dataset?
Well, there is another factor that is strongly associated with cirrhosis, and that is obesity (). In fact, obesity is associated with just about any major disease, including various types of cancer ().
And in countries where pork is a dietary staple, isn’t it reasonable to assume that pork consumption will play a role in obesity? Often folks who consume a lot of addictive industrial foods (e.g., bread, candy, regular sodas) also eat plenty of foods with saturated fat; and the latter end up showing up in disease statistics, misleadingly supporting the lipid hypothesis. The phenomenon involving pork and cirrhosis may well be similar.
But you may find the above results and argument not convincing enough. Maybe you want to see some evidence that pork is actually good for one’s health. The results above suggest that it may not be bad at all, if you buy into the obesity angle, but not that it can be good.
So I downloaded the most recent data from Nationmaster.com () on the following variables: pork consumption, alcohol consumption, and life expectancy. The list of countries was a bit larger than and different from that in the Bridges study; the following countries were included: Australia, Brazil, Canada, China, Denmark, France, Germany, Hong Kong, Hungary, Japan, Mexico, Poland, Russia, Singapore, Spain, Sweden, United Kingdom, and United States. Below are the results of a simple multivariate analysis with WarpPLS.
As with the Bridges dataset, there is a strong multivariate association between alcohol and pork consumption (0.43). The multivariate association between alcohol consumption and life expectancy is negative (-0.14). The multivariate association between pork consumption and life expectancy is positive (0.36). Neither association is statistically significant, although the association involving pork consumption gets close to significance with a P=0.11 (a confidence level of 89 percent; calculated through jackknifing, a nonparametric technique). The graphs show the plots for the associations and the best-fitting lines; the blue dashed arrows indicate the multivariate associations to which the graphs refer. So, in this second dataset from Nationmaster.com, the more pork is consumed in a country, the longer is the life expectancy in that country.
In other words, for each 1 standard deviation variation in pork consumption, there is a 0.36 standard deviation variation in life expectancy, after we control for alcohol consumption. The standard deviation for pork consumption is 36.281 lbs/person/year, or 45.087 g/person/day; for life expectancy, it is 4.677 years. Working the numbers a bit more, the results above suggest that each extra gram of pork consumed per person per day is associated with approximately 13 additional days of overall life expectancy in a country! This is calculated as: 4.677/45.087*0.36*365 = 13.630.
Does this prove that eating pork will make you live longer? No single study will “prove” something like that. Pork consumption is also likely a marker for wealth in a country; and wealth is strongly and positively associated with life expectancy at the country level. Moreover, when you aggregate dietary and disease incidence data by country, often the statistical effects are caused by those people in the dietary extremes (e.g., alcohol abuse, not moderate consumption). Finally, if people avoid death from certain diseases, they will die in higher quantities from other diseases, which may bias statistical results toward what may look like a higher incidence of those other diseases.
What the results summarized in this post do suggest is that pork consumption may not be a problem at all, unless you become obese from eating it. How do you get obese from eating pork? Eating it together with industrial foods that are addictive would probably help.
Subscribe to:
Posts (Atom)












