Analysis of the Influence of Socio-Demographic Variables and Some Nutrition and Lifestyle Habits on Beverage Consumption in the Spanish Population

Beverages and drinks play a significant role in maintaining the integral health of individuals. The aim of this study is to discover the pattern of beverage consumption in different groups of the Spanish population and to investigate its relationship with other nutritional variables and habits. To achieve the objectives, an observational, descriptive and cross-sectional study was conducted. For data collection, a questionnaire was designed and validated that explored different beverage and food consumption variables as well as socio-demographic and lifestyle variables. The instrument was disseminated, among the Spanish young adult population, through snowball sampling using social networks, collecting a sample of 17,541 valid surveys. Bivariate comparative analyses and correlation analyses were performed, and finally, the principal component analysis (PCA) method was used in order to study the relationships between variables related to drinking and health. The main results show significant differences in the pattern of beverage consumption between the socio-demographic variables of sex, age and educational level, as well as between different areas of Spain, while the PCA model shows the relationship between the consumption of sugar-sweetened beverages with the Healthy Nutrition Index of the population and sport practice. Based on the results of the study, the following conclusions were reached: the beverage consumption pattern of the Spanish population is affected by socio-demographic variables. Healthier drinking habits affect the nutrition and health of the population.


Introduction
Beverages and drinks are essential elements in maintaining the health and well-being of individuals [1].Staying adequately hydrated is essential for maintaining bodily functions such as digestion, circulation, temperature regulation, and the transport of nutrients and oxygen [2,3], and the type of beverages consumed can have a significant impact on overall health [4][5][6].Water is the most consumed drink and the primary source of hydration, but other beverages also contribute to fluid intake.
Body hydration is not the only health contribution that beverages make, many beverages, such as fruit juices and smoothies, can provide essential vitamins, minerals, and antioxidants [7,8] that support the body's immune system and overall vitality.Drinks like milk provide calcium and vitamin D [9], which are crucial for maintaining strong and healthy bones, and certain beverages, like herbal teas, can aid in digestion and soothe gastrointestinal discomfort [10,11].Certain beverages, like green tea and red wine (in Foods 2023, 12, 4310 2 of 19 moderation), contain antioxidants that may support heart health [12,13].Finally, beverages that contain carbohydrates, such as sports drinks, or drinks containing excitants, such as coffee, can provide a quick source of energy during physical activities or during the working day [14,15].
However, it is important to be mindful of the type and quantity of beverages consumed.Some beverages can be high in added sugars, artificial flavours, and empty calories, which may contribute to health issues such as obesity, diabetes, and dental problems [4,16,17].Choosing low-calorie, nutrient-dense beverages can help with weight management and control calorie intake [1,18].It is advisable to prioritise the consumption of water relative to that of other beverages to keep the body hydrated, or choose beverages that provide nutritional benefits without excessive sugar or unhealthy additives [19].
The types of drinks people consume are often influenced by cultural and socioeconomic factors.These factors play a significant role in shaping dietary preferences and habits, including beverage choices.Another important factor that can influence the choice of one beverage over another is the availability of ingredients and the climatic conditions of the place where you live [20].Or cultural norms can influence whether certain beverages are considered acceptable or appropriate to consume in specific settings.
If we analyse the influence of the economic or cultural level, it is observed that people with higher incomes can afford to choose from a wider range of beverages and may be more likely to afford premium or specialty drinks [21], or that people with a higher level of education and culture may be better informed about healthier drink options [22].
The impact that beverage consumption can have on a person's health is significant, which is why it is important to understand the beverage consumption patterns of a population and to investigate how these relate to socio-cultural factors as well as to other nutritional habits and lifestyles.Understanding these influences is important for public health efforts, as it helps in developing strategies to promote healthier beverage choices that are culturally sensitive and accessible to people from different socioeconomic backgrounds.
Despite the growing concern about population nutrition, there are still few studies that focus specifically on the impact of beverages on nutrition and health, and that analyse the contribution of different types of beverages.
While there are numerous studies exploring beverage consumption, most of them focus only on measuring the health impact of consuming a specific beverage.There are, for example, studies conducted on energy drinks in Denmark [23] or the United States [24,25], studies focusing on sugary drinks such as those conducted in the UK [26], Peru [27], India [28] or Sweden [29], or studies focused on exploring the consumption of carbonated beverages [30].In contrast, there are few studies that focus on exploring the demographic determinants of beverage consumption in general.In this respect, the study conducted on the French population in 2016, which includes a very large sample, is worth mentioning [31].
With regard to Spain, studies have been carried out analysing beverage consumption in specific regions of the territory, such as Catalonia [32], the Balearic Islands [33] or in the Community of Madrid [34].Furthermore, there are other studies that focus only on the consumption of specific beverages such as energy drinks [35] or alcohol [36][37][38][39] that have been studied in depth.Spanish studies analysing the impact of different types of beverages on health and their socio-demographic determinants are scarce and the sample used is not very large.The ANIBES study [40] has a sample of 2285 participants, while the study conducted by Cíntia Ferreira-Pêgo et al. [41] has a sample of 1262 adults.
This study, on the beverage consumption of the Spanish population, aims to contribute to the evidence in this field.The large number of participants is undoubtedly one of its strengths and may allow extrapolable results to be obtained.
The objectives of this study are as follows: 1.
To explore the associations between beverage consumption and sociodemographic factors.Specifically, the relationship between age, sex, educational level, purchasing power, place and type of residence of the person and beverage consumption was analysed.

2.
To analyse the relationship between the frequency of consumption of a given beverage and some nutritional, lifestyle and health variables.
Specifically, the variables explored were as follows: the Healthy Nutrition Index, body mass index, self-perceived level of health, sedentary lifestyle, physical activity, and time and quality of sleep.

Study and Sampling
A prospective cross-sectional study was carried out on the Spanish population aged between 18 and 45 years living in Spain.We decided to limit the sample to those ages to minimise the possible bias caused by the greater difficulty in accessing the internet and social networks that most of the older population has.Subjects were excluded if, at the time of answering the survey, they had any pathology or limitation that could affect their diet, such as hospitalisation or confinement.

Ethical Approval
The study was carried out with the approval of the ethics committee of the Universidad Catolica de Valencia (UCV/2019-2020/152).All data collected came with the informed consent of the participants in the study, and the ethical recommendations contained in the Declaration of Helsinki [42] were followed at all times.

Instrument
The instrument used for data collection was a self-developed questionnaire consisting of several parts: the first part collected information on anthropometric data such as weight, height and socio-demographic variables, the second part collected information on different food groups and their respective frequencies of consumption, the third part collected information on drinking habits, and the fourth and last part collected data on sports and lifestyle habits related to health.Rigorous methodological approaches were followed for instrument development and psychometric testing.
The questionnaire was validated through a pilot group of 53 people, with similar characteristics to the study population, and through a nominal group of 7 health experts, namely a nutritionist, two family doctors, two psychologists, a social educator and an expert from the world of communication.

Data Collection
The questionnaire was disseminated using non-probabilistic snowball sampling.The support used was the Google Form program and the channels used were mostly telematic; firstly, the Instagram account @elretonutricional was used, but so were the researchers' personal social networks and mailings to various associations.
On the other hand, we also tried to disseminate the questionnaire physically, in order to facilitate access to those people who are less fond of networks and technology, using posters hung in different shops and businesses where many different people go every day.The data collection period was from August 2020 to November 2021.

Variables
The socio-demographic variables used were as follows: sex (which was analysed in binary form, distinguishing between male and female), age (which was categorised between young people who were (18-30) years old and adults who were (31)(32)(33)(34)(35)(36)(37)(38)(39)(40)(41)(42)(43)(44)(45) years old), level of studies (which was categorised into basic studies (no studies, only primary or secondary studies, vocational training or baccalaureate) and higher studies (degree, master's and PhD)), type of residence (which was categorised into those living in the family home and those living outside), shared residence (where a distinction was made between those living alone and those living with others) and place of residence (where a distinction was made between the different regions of Spain).
To analyse nutrition, we decided to use the IASE index (Índice de Alimentación Saludable para la población Española) [43] by adapting a reduced version of it.The adapted index reports a maximum score of 73 points, in accordance with the degree of compliance with the recommendations proposed by the Spanish Society of Community Nutrition (SENC) [44].According to the IASE score obtained, the nutrition habits of the population can be classified into 3 categories: 58.4 < IASE < 73: "Healthy"; 36.5 < IASE < 58.4: "Needs changes"; IASE < 36.5:"Unhealthy".Table 1 shows the categorisation of variables used for the IASE.The beverage consumption variables were as follows: water, sugary soft drinks, juice, coffee and energy drinks, drinks at mealtimes, alcohol and binge drinking.The health and lifestyle variables were a sedentary lifestyle, sleeping hours and sleep quality.and were categorised with a score from 1 to 4 according to the indicated criteria in Table 2, except for body mass index (BMI) and minutes of sport, which was used as a numerical.

Statistical Analysis
To begin with, the data were cleaned by eliminating erroneous data, as a consequence of incorrect entry by the respondent, or atypical data on both height and weight, trying to eliminate BMI values below 14 and above 40.Once the information had been cleaned, the corresponding descriptive and inferential statistical analysis was carried out using Excel 2013 for Windows and SPSS v.23 (SPSS Inc., Chicago, IL, USA).
The non-parametric Mann-Whitney U test for independent samples (to analyse dichotomous variables) and the Kruskal-Wallis test (for comparison between more than two groups) were used to compare means and study significance.The choice was motivated by the fact that the normality criterion had previously been tested using the Kolmogorov-Smirnov test with Lilliefors correction and the results had indicated that the data did not follow a normal distribution.Both tests were applied with a significance level of 0.05.
Performing an analysis for each pair of variables studied is not functional, so we decided to use the principal component analysis (PCA) method to explore the relationships between the ordinal numerical variables, since this method allows the dimensionality of a model to be reduced.Two principal components were chosen that explain a percentage of variability of about 24%, as can be seen in Figure 1.
The first dimension (or principal component) is mainly explained by sugary drinks, IASE, meal beverages, and sleep quality, as can be seen in Figure 2. Figure 3 shows how the variables alcohol and getting drunk are the main contributors to the second dimension.This method was carried out with RStudio 4.3.0.
between the ordinal numerical variables, since this method allows the dimensionality of a model to be reduced.Two principal components were chosen that explain a percentage of variability of about 24%, as can be seen in Figure 1.The first dimension (or principal component) is mainly explained by sugary drinks, IASE, meal beverages, and sleep quality, as can be seen in Figure 2. Figure 3 shows how the variables alcohol and getting drunk are the main contributors to the second dimension.This method was carried out with RStudio 4.3.0.The first dimension (or principal component) is mainly explained by sugary drinks, IASE, meal beverages, and sleep quality, as can be seen in Figure 2. Figure 3 shows how the variables alcohol and getting drunk are the main contributors to the second dimension.This method was carried out with RStudio 4.3.0.

Results
Table 3 shows the socio-demographic characteristics of the respondents.After data collection and subsequent filtering, we finally obtained 17,541 valid surveys representative of people resident in the whole of Spain, including the islands.The behaviour of the sample was studied in relation to the frequency of consumption of water, sugary soft drinks, juices, coffee and energy drinks, and alcohol, and the frequency of binge drinking, distinguishing this via the socio-demographic variables of sex,

Results
Table 3 shows the socio-demographic characteristics of the respondents.After data collection and subsequent filtering, we finally obtained 17,541 valid surveys representative of people resident in the whole of Spain, including the islands.The behaviour of the sample was studied in relation to the frequency of consumption of water, sugary soft drinks, juices, coffee and energy drinks, and alcohol, and the frequency of binge drinking, distinguishing this via the socio-demographic variables of sex, age, educational level, type of residence and shared residence (Table 4).Significant differences in the frequency of consumption of the same measures were also analysed with respect to the region of residence in Spain (Table 5).As all p-values for each Kruskal-Wallis between the drinking variables and the autonomous communities of Spain (regions of Spain) are significant, it could be interesting to see the differences in a map (Figure 4) particularly with respect to sugar-sweetened beverages, since, together with alcohol, they are the ones that have the greatest impact on health.
As all p-values for each Kruskal-Wallis between the drinking variables and the autonomous communities of Spain (regions of Spain) are significant, it could be interesting to see the differences in a map (Figure 4) particularly with respect to sugar-sweetened beverages, since, together with alcohol, they are the ones that have the greatest impact on health.
It can be observed how the consumption of sugar-sweetened beverages is higher in the northern part of Spain than in the centre-south, particularly in the autonomous communities of Cantabria and Asturias, followed by Galicia and the Balearic Islands.Before applying a PCA model, it is essential to study the correlation between the variables implied in the model.A Pearson correlation graph can be seen in Figure 5.It can be observed that there are no high correlations between the variables, except for get drunk and alcohol (0.295), as is logical, sleeping hours and sleeping quality (0.366), as is logical too, and more interestingly between sugary drinks and the IASE (0.365).
People who do not drink alcohol or drink alcohol infrequently (high punctuation in alcohol) tend not to get drunk (high punctuation in get drunk).Furthermore, survey respondents who do not usually drink sugary drinks also have a good IASE.Finally, those who have good sleep quality sleep more hours.However, it is worth noting that all these correlations are not high.It can be observed how the consumption of sugar-sweetened beverages is higher in the northern part of Spain than in the centre-south, particularly in the autonomous communities of Cantabria and Asturias, followed by Galicia and the Balearic Islands.
Before applying a PCA model, it is essential to study the correlation between the variables implied in the model.A Pearson correlation graph can be seen in Figure 5.It can be observed that there are no high correlations between the variables, except for get drunk and alcohol (0.295), as is logical, sleeping hours and sleeping quality (0.366), as is logical too, and more interestingly between sugary drinks and the IASE (0.365).
People who do not drink alcohol or drink alcohol infrequently (high punctuation in alcohol) tend not to get drunk (high punctuation in get drunk).Furthermore, survey respondents who do not usually drink sugary drinks also have a good IASE.Finally, those who have good sleep quality sleep more hours.However, it is worth noting that all these correlations are not high.
Figure 6 shows the PCA model and helps to explore the relationships between all variables related to drinks and those related with health.
In the PCA, via the colour and length of the arrows, following the colour scale that marks the level of contribution of the variables, it can be seen how the variables sedentary lifestyle, and frequency of consumption of juice, coffee and energy drinks have the least influence on the model; in fact, they are coloured in blue and their lengths are very short.The variables educational level, water consumption, sport and BMI, which are coloured green, also have a limited contribution, while the variables related to sleep, binge drinking and IASE have a medium contribution, with their arrow being coloured in yellow and of considerable length.Finally, the consumption of alcohol and that of sugar-sweetened beverages are the variables with the greatest influence on health, with their arrows being the longest and coloured orange-red in the PCA model.Looking now at the orientation of the arrows, alcohol is more related to get drunk and coffee and energy drinks.However, this last variable contributes less to the model.Additionally, sleeping hours and sleep quality are highly positively correlated and negatively correlated with BMI.Therefore, the more time spent and the greater quality of sleep, the smaller the punctuation in BMI.Finally, as already observed in Pearson's correlation model, IASE and sugary drinks are correlated.Then, those who drink less sugary drinks have a higher IASE.

Discussion
One of the main findings of this study is that, for the young adult Spanish population, the frequency of beverage consumption shows a clear dependence on socio-demographic variables, particularly with regard to gender, age and educational attainment (see Table 4 and Appendix A).This is in line with previous studies in the literature that focus on specific beverages in different regions of the world, such as the one conducted by Jozaa Z AlTamimi [29], the other one conducted by Manu Raj Mathur [45], both focusing on sugar-sweetened beverages, several that focus on alcohol consumption [46][47][48], or those conducted on the consumption of energy drinks [23,24] or caffeinated beverages [49].This dependence of beverage consumption on socio-demographic factors is of fundamental importance in public health.Knowing that gender, age, level of education or place of residence determines people's choice or prioritisation of one drink or another can help to design and disseminate more targeted and therefore more effective public health actions.Training or awareness-raising campaigns can be targeted at a particular gender or age group, or personalised according to the place of residence and therefore the climatic conditions or customs of that territory.We can also use a simpler and plainer message for those with basic education and more scientific and specific language for those with higher education.Although it is certainly necessary to carry out more studies focused on specific socio-demographic variables, which will allow a further refinement of the results, so that the public health measures adopted can be effective, the results obtained can be a starting point.
On the other hand, the data allow us to observe that living alone as opposed to living with other people, or living in a family home as opposed to living away from home, although there are significant differences between categories in beverage consumption, tare not as marked as in the case of other socio-demographic variables (see Figures A4 and A5 of Appendix A).This is important given the steady increase in recent years in Western society in people living alone [50][51][52].While it has been observed that living outside or away from the family home influences some nutritional patterns and health habits [53][54][55][56], the results seem to indicate that it has little influence on drinking.Perhaps it could be interesting to study, and therefore make a future line of work, whether this behaviour, observed for the young adult population, also holds for the older population, which is generally a population more at risk of suffering the detrimental effects of loneliness and living alone [57][58][59][60].Our research group, without going any further, plans to carry out a new round of dissemination of the questionnaire and a new data collection round, this time using the inclusion criterion of being over 45 years old.In this way, we will be able to obtain results from adults and elderly people that can be compared with those obtained from the younger age group.
Another result that may attract attention is the consumption of sugar-sweetened beverages according to the regions of Spain; although the differences in absolute values are not very marked, they are statistically significant.Figure 4 visually shows how it is higher in general in the northern part of Spain than in the centre-south, particularly in the autonomous communities of Cantabria and Asturias, followed by Galicia and the Balearic Islands.This could be explained by the fact that these regions have different cultural and gastronomic aspects compared to those in other parts of Spain, and this may also have a significant influence on the type of beverages consumed in that area.The relationship between culture, food and beverages is intrinsic and complex, and may be a reflection of history, climate, availability of ingredients and local traditions [61,62].It could be interesting in the future to carry out studies on food, nutrition and beverage consumption habits focused on specific regions of the Spanish territory in order to study in more detail the influence of gastronomy and customs and how they vary at the territorial level.Such studies, in a specific area, could perhaps use mixed techniques, both quantitative and qualitative.They could combine the dissemination of a questionnaire to the population, with experimental data collection in restaurants and establishments to analyse the type of diets and drinks offered and the consumption pattern of customers.They could also use the data collected in national surveys to analyse the food purchasing pattern of the population.
Another striking result of the study is the relationship between the practice of physical activity and the consumption of sugar-sweetened beverages, which can be seen visually in the PCA graph in Figure 6.In this figure, it can be seen how the lines representing the contribution of the two variables are almost parallel and point in the same direction.Increased physical activity increases the need for hydration due to the loss of fluids and electrolytes through sweat [63][64][65].The most common beverages for hydration during exercise are typically water, sports drinks designed to replenish electrolytes and, in some cases, carbohydrate drinks to provide additional energy.The choice of beverage depends on the duration and intensity of exercise, as well as personal preference.Sport has become an extremely popular activity in Western society and around the world in recent years.This popularity has led to a significant increase in the movement of money and a growing media interest in sport.All of these factors have led to increasingly aggressive marketing and advertising campaigns by brands seeking to capitalise on this trend.Not all brand marketing and advertising messages provide accurate information about the dietary and nutritional properties of a product and its healthfulness [66][67][68].While sugary drinks can be very appealing in taste, and in some cases effective thirst quenchers, they are not always the healthiest option for replenishing fluids after sporting activities [69,70].This opens up another important field of public health, which especially affects the young and young adult population, which is the one that generally practices sports on a more consistent basis.It is a priority to design effective training campaigns on the consumption of beverages before, during and after sports practice, to teach the population to practice healthy hydration and to choose products that are not only commercially attractive but also nutritionally adequate.However, in order to achieve this effectively, it will be necessary to rely on successive studies specifically designed to analyse the consumption of sugar-sweetened beverages in young people, sportsmen and sportswomen.

Strengths and Limitations
One of the weaknesses of the study is represented by the sampling technique used.Snowball sampling was used, which offers the advantage of reaching a large number of responses, but may result in a biased sample because the sample collected depends on participants recommending other potential participants.
Analysing our sample, we detected some self-selection bias in that the sample collected had a large number of people concerned about nutrition and health.This is mainly due to the fact that when the survey was distributed through social networks, several influencers and professionals working in these fields supported the dissemination of the survey.It can be assumed that the community of followers of these accounts are, for the most part, people interested in nutrition and health who may have higher-than-average knowledge of the subject and who generally seek to adopt healthier habits.
We tried to minimise self-selection bias by also disseminating the survey physically, outside of social networks.To this end, we contacted different establishments throughout Spain with a heterogeneous public (pharmacies, tobacconists, associations, etc.) and asked them to display a poster in their establishments explaining the objectives of the research and containing a QR code that referred to the survey.In this way, people visiting these establishments could access the survey if they wished to do so.
We should also highlight a gender bias, with 82.26% of those who responded to the survey being female.This result is not surprising because in general there is a greater concern among women about nutrition and health issues, which predisposes them to read and be more interested in these topics and, consequently, to be more inclined to participate in studies in this field.
To minimise, as much as possible, this sex bias, an effort was made to recruit male representatives for the study.Although it was not possible to match female participation, a broadly significant sample of 3111 males was achieved.
Finally, the diffusion of the sample through social networks may induce an age bias in the sample, presenting a higher response rate among the young population.Young people in general are more familiar with and attracted to the use of social networks, tend to devote more time to their use, and were therefore more likely to answer the survey.
After reflecting on this bias and in order to minimise it, it was finally decided to focus the study on the population between 18 and 45 years of age, ages in which the use of networks is widespread.
The greatest potential of this study undoubtedly lies in the large sample size, which provides great statistical power and reduces the study's margin of error.Having a large population allows a representation of different demographic groups, which increases the likelihood that the conclusions are representative of the general population.
Another noteworthy point of the present study is the geographic scope achieved, with a representative sample from every region of Spain, including the islands.Geographic diversity allows a more complete understanding of the subject under investigation and increases the external validity of the study, reducing the possibility of geographic bias.
Finally, the potential of the questionnaire obtained for data collection cannot be overlooked; the ease and speed with which it can be completed makes it suitable for wide dissemination through different channels, thus making it possible to achieve very large and significant samples in different groups.
A next step could be to replicate the research by focusing it on an older target population to explore the behaviour, in terms of beverage consumption and nutrition and lifestyle habits, of adults and the elderly.

Conclusions
The main results show significant differences in the pattern of beverage consumption between the socio-demographic variables of sex, age and educational level, as well as between different areas of Spain.Furthermore, it can be concluded that living alone or with others and living in the family home or away from home also affect drinking habits, although to a lesser extent than other socio-demographic variables do.
The PCA model shows the inverse relationship between the consumption of sugarsweetened beverages with the Healthy Nutrition Index of the population and the direct relationship between the consumption of sugar-sweetened beverages and the practice of physical activity.In addition, hours of sleep and sleep quality have been found to correlate with body mass index, leading to the conclusion that those who spend more time sleeping and who have a higher quality of sleep in general have a lower BMI.
It can be concluded that the beverage consumption pattern of the Spanish population is affected by socio-demographic variables.Healthier drinking habits affect the nutrition and health of the population.Figure A2 shows the distribution of subjects according to sex, with men in green and women in red.Significant differences can be seen visually between the two ellipses that show the mean of the groups of points.The ellipse for men is displaced in the quadrant Figure A2 shows the distribution of subjects according to sex, with men in green and women in red.Significant differences can be seen visually between the two ellipses that show the mean of the groups of points.The ellipse for men is displaced in the quadrant of the PCA where the variables alcohol consumption, binge drinking and consumption of sugary soft drinks are collected, indicating that men consume these beverages more frequently than women do.
Figure A3 shows the distribution of subjects according to the level of education, where green represents the subjects with basic education and red represents the respondents with higher education.Significant differences can be seen visually between the two ellipses that show the mean of the groups of points.The ellipse of those with basic studies is displaced in the quadrant of the PCA where the variables consumption of coffee and energy drinks, BMI, alcohol consumption and binge drinking are collected, indicating that people with basic studies in general have worse eating habits and a higher BMI compared to their counterparts with higher studies.
Figure A4 shows the distribution of subjects according to whether they live alone or accompanied, where green represents the subjects who live alone and red represents the respondents who live accompanied.Figure A5 shows the distribution of subjects according to whether they live in the family home or outside it, where green represents the subjects who live in the family home and red represents those who live outside it.For these last two variables, the differences are minimal and in both cases the ellipses of one group and the other almost overlap.
of the PCA where the variables alcohol consumption, binge drinking and consumption of sugary soft drinks are collected, indicating that men consume these beverages more frequently than women do. Figure A3 shows the distribution of subjects according to the level of education, where green represents the subjects with basic education and red represents the respondents with higher education.Significant differences can be seen visually between the two ellipses that show the mean of the groups of points.The ellipse of those with basic studies is displaced in the quadrant of the PCA where the variables consumption of coffee and energy drinks, BMI, alcohol consumption and binge drinking are collected, indicating that people with basic studies in general have worse eating habits and a higher BMI compared to their counterparts with higher studies.Figure A3 shows the distribution of subjects according to the level of education, where green represents the subjects with basic education and red represents the respondents with higher education.Significant differences can be seen visually between the two ellipses that show the mean of the groups of points.The ellipse of those with basic studies is displaced in the quadrant of the PCA where the variables consumption of coffee and energy drinks, BMI, alcohol consumption and binge drinking are collected, indicating that people with basic studies in general have worse eating habits and a higher BMI compared to their counterparts with higher studies.accompanied, where green represents the subjects who live alone and red represents the respondents who live accompanied.Figure A5 shows the distribution of subjects according to whether they live in the family home or outside it, where green represents the subjects who live in the family home and red represents those who live outside it.For these last two variables, the differences are minimal and in both cases the ellipses of one group and the other almost overlap.accompanied, where green represents the subjects who live alone and red represents the respondents who live accompanied.Figure A5 shows the distribution of subjects according to whether they live in the family home or outside it, where green represents the subjects who live in the family home and red represents those who live outside it.For these last two variables, the differences are minimal and in both cases the ellipses of one group and the other almost overlap.

Figure 1 .
Figure 1.Scree plot of the PCA Model.NOTE: The red line indicates the average percentage of explained variance (if all dimensions are explained as the same).

Figure 2 .
Figure 2. Contribution of variables to the first dimension.NOTE: The red line indicates the average percentage of contribution if all dimensions contribute the same.

Figure 1 .
Figure 1.Scree plot of the PCA Model.NOTE: The red line indicates the average percentage of explained variance (if all dimensions are explained as the same).

Figure 1 .
Figure 1.Scree plot of the PCA Model.NOTE: The red line indicates the average percentage of explained variance (if all dimensions are explained as the same).

Figure 2 .
Figure 2. Contribution of variables to the first dimension.NOTE: The red line indicates the average percentage of contribution if all dimensions contribute the same.

Figure 2 .
Figure 2. Contribution of variables to the first dimension.NOTE: The red line indicates the average percentage of contribution if all dimensions contribute the same.

Figure 3 .
Figure 3. Contribution of variables to the second dimension.NOTE: The red line indicates the average percentage of contribution if all dimensions contribute the same.

Figure 3 .
Figure 3. Contribution of variables to the second dimension.NOTE: The red line indicates the average percentage of contribution if all dimensions contribute the same.

Figure 4 .
Figure 4. Map of autonomous communities of Spain coloured depending on the punctuation of sugar drinks.

Figure 4 .
Figure 4. Map of autonomous communities of Spain coloured depending on the punctuation of sugar drinks.

FoodsFoodsFigure 5 .
Figure 5. Pearson correlation plot between "the drinking variables" and health variablesFigure6shows the PCA model and helps to explore the relationships be variables related to drinks and those related with health.

Figure 6 .
Figure 6.Variables plot of the PCA model.

Figure 5 .
Figure 5. Pearson correlation plot between "the drinking variables" and health variables.

Figure 5 .
Figure 5. Pearson correlation plot between "the drinking variables" and health variables.

Figure 6
Figure6shows the PCA model and helps to explore the relationships between all variables related to drinks and those related with health.

Figure 6 .
Figure 6.Variables plot of the PCA model.Figure 6. Variables plot of the PCA model.

Figure 6 .
Figure 6.Variables plot of the PCA model.Figure 6. Variables plot of the PCA model.

Figure A1 .
Figure A1.Variable plot of the age distribution.

Figure A1 .
Figure A1.Variable plot of the age distribution.

Figure A2 .
Figure A2.Variable plot of the sex distribution.

Figure A3 .
Figure A3.Variable plot of the study-level distribution.

Figure A2 .
Figure A2.Variable plot of the sex distribution.

Figure A2 .
Figure A2.Variable plot of the sex distribution.

Figure A3 .
Figure A3.Variable plot of the study-level distribution.

Figure A3 .
Figure A3.Variable plot of the study-level distribution.

Figure A4 .
Figure A4.Variable plot distribution of the shared residence.

Figure A5 .
Figure A5.Variable plot distribution of living at home.

Figure A4 .
Figure A4.Variable plot distribution of the shared residence.

Figure A4 .
Figure A4.Variable plot distribution of the shared residence.

Figure A5 .
Figure A5.Variable plot distribution of living at home.

Figure A5 .
Figure A5.Variable plot distribution of living at home.

Table 1 .
Conversion table applied to convert qualitative questions into quantitative questions of the nutritional index.

Table 2 .
Categorisation of beverage consumption variables and the health and lifestyle variables.

Table 3 .
Data collected and their socio-demographic characteristics.

Table 3 .
Data collected and their socio-demographic characteristics.

Table 4 .
Comparison of drinking habits differentiated by sociodemographic variables.

Table 5 .
Comparison of drinking habits differentiated by autonomous communities of Spain.