Abstract

Background. Obesity and eating disorders are increasing in occidental countries and can undermine physical and psychological health. Therefore, preventing the insurgency of unhealthy eating habits in childhood is fundamental. Parents can play an important role in assisting pediatricians, psychiatrists, and clinical psychologists in the diagnosis of eating disorders because they have an active role in observing and assessing the quality of their children’s eating habits. Methods. In our study, we collected data from a sample of children (n = 125) and their parents (n = 161) without symptoms related to eating disorders. Parents assessed the eating habits, behavior problems, and mental health of their children and parental stress. In addition, we measured body mass index, anxiety, and lifestyle in children. Data were analyzed with bivariate correlation and MIMIC models. Results. Both mothers’ and fathers’ assessments of children’s eating habits are reliable. Unhealthy eating habits are connected with children’s behavioral problems and parental stress. We did not find significant differences in feeding styles and ways of assessing the quality of eating habits between mothers and fathers. Our study showed greater sensitivity of mothers toward the physical fitness of their children rather than fathers. Conclusions. Mothers and fathers both proved to be good observers of their children’s eating behavior, and they could cooperate with medical and psychological operators in preventing the risk of obesity.

1. Introduction

According to the National Health and Nutrition Examination Survey [1], during the period 2017–2018, the prevalence of severe obesity among USA adults was 9.2%. Women had a higher prevalence of severe obesity (11.5%) than men (6.9%). The prevalence was highest among adults aged 40–59 (11.5%), followed by adults aged 20–39 (9.1%) and adults aged 60 and over (5.8%). Moreover, in Europe, the problem of obesity is evident. To date, almost 60% of the population of the OECD countries is overweight, and 25% have severe obesity [2]. Average rates of adult obesity in OECD countries have increased from 21% in 2010 to 24% in 2016, so 50 million people are now obese [2]. In 2017, in Italy, obesity affects a total of 10.5% of the population over the age of 18, 14.5% among those over 65. This condition is found more frequently among residents of the South (11.8%), reaching the highest share in Sicily (12.6%); [3].

Correct eating habits prevent people from developing eating disorders, which, in turn, can undermine physical health. In fact, obesity is associated with a high risk of heart attack, arterial hypertension, hyperglycemia, high cholesterol level, and diabetic hyperinsulinemia [4]. In children, obesity can lead to chronic inflammation, hypothyroidism, and hypogonadism [4]. Obesity and eating disorders, however, affect also mental health and psychological well-being because a bad diet is associated with a low quality of life and a low level of mental health [59].

Therefore, it is important that people develop good eating habits, especially during their childhood and adolescence. Many studies showed the important impact that parents have on children’s eating habits. In particular, parents’ well-being and mental health resulted to have an impact on children’s eating behavior [10, 11]. Some studies showed that if children are left alone in selecting their food, they tended to select food with a high level of sugar [12]. Without the intervention of parents, particularly mothers, children tend to choose food more on the basis of its hedonistic characteristics than on the basis of its healthful effect on the body [12, 13]. The parental contribution to helping children find a healthy diet is indisputable [1416]. Parents have the responsibility for providing children with an array of healthful food, while children have the responsibility for how much she/he eats [17].

Recent research showed that public programs for interventions against childhood obesity through schools cannot be efficient without considering other sources of support for a healthy lifestyle, particularly family [18]. Parents play an important role in supporting and helping their children to keep healthy dietary habits. Some studies showed that low parental competencies and lack of control in children’s nutrition have negative effects on children’s dietary habits [19]. Other studies showed that parents’ underestimation of children’s obesity, parents’ pressure to eat, and low parent’s level of monitoring have negative effects on the selection of healthy food for children [20]. Parents’ attitude toward children’s dietary habits has an important effect on children’s nutrition quality [21].

Even if research showed that parents can play an important role in preventing the appearance of unhealthy eating habits in their children and in supporting healthy eating habits [15], it remains necessary to know their ability in assessing the quality of children’s eating habits, whether these habits have some relations with children’s mental health, and whether these habits can generate stress for parents.

Because other studies showed that children’s health conditions can affect parents’ attitude toward their eating habits [20, 21], we also took into consideration the psychological state of the children, particularly anxiety, and their physical conditions, particularly the body mass index (BMI), to see if there is a relationship between children’s physical and psychological characteristics and the assessment of their eating habits by parents. Children with unhealthy eating habits should have an abnormal BMI (BMI >30 for obesity and <18 for underweight) and a lower level of psychological well-being indicated by excessive anxiety [10].

We have formulated different hypotheses. Our first hypothesis (H1) is that children’s eating habits are affected by their physical and psychological conditions. In particular, this hypothesis is divided into two subhypotheses:H1a: Children’s eating habits are affected by children’s BMI (the higher the BMI, the more unhealthy are children’s eating habits).H1b: Children’s eating habits are affected by their anxiety (the higher the anxiety, the more unhealthy are children’s eating habits).

Our secondary hypothesis (H2) is that children’s unhealthy eating habits are connected to children’s mental health because children with unhealthy habits should show psychological or behavioral disorders more intensely [57] and connected to parental emotive reactions because children with unhealthy habits tend to increase the level of parental stress [10]. This hypothesis is divided into two subhypotheses:H2a: Children’s mental health is negatively affected by their unhealthy eating habits.H2b: Parental stress is increased by children’s unhealthy eating habits.

In addition, considering the fact that there are some studies showing differences between attitudes and reactions of mothers and fathers toward children’s eating habits [14, 16, 22], we decided to analyze whether there are some differences between mothers’ and fathers’ assessment of their children’s behavior in relation to eating habits.

2. Methods

2.1. Participants

For our study, we selected a sample of parents and their children from a population attending a local scholastic complex in the central part of Italy (Abruzzo) dedicated to primary education. The scholastic complex consisted of two school buildings, both of which consisted of primary classes. We have directed our attention to primary school children because not only this period is fundamental for the development of healthy eating habits [15, 17, 23] but also because children of this age can answer structured questionnaires. A total of 125 children (50.5% females) and 161 parents (57.1% females) participated in the research. Children’s age ranged from 8 to 11 years (mean age = 9.50; SD = 0.604); parents’ age ranged from 28 to 57 years (mean age = 43.10; SD = 5.322). The majority of parents were employed as dependent or autonomous workers (45%) or as employees in administrative or commercial jobs (21%) and had a high school degree (51%—see Table 1 for these and other demographic data).

2.2. Materials for Child Assessment

The materials used for the surveys on children consisted of tools for anthropometric and psychological assessment. The anthropometric assessment was done with a meter for measuring body stature and a scale for measuring body weight. These measures were then used to assess the BMI of children in relation to their age. Psychological assessment was done by asking children to compile a short questionnaire, in which we asked them if they played sports, if they liked eating, if they were judged fat by their classmates, and if they were satisfied with the quality of their relationship with their parents by asking how much time the mother and father spent with them. Finally, in the questionnaire, we included two scales taken from the Screen for Child Anxiety Related Emotional Disorders (SCARED) [24]. These subscales, each comprising five items, measured the children’s general anxiety level (Gen Anx) and separation anxiety (Sep Anx). The scores on the anxiety scales were then converted into values from 0 (no anxiety) and 2 (a lot of anxiety). Moreover, a global score of anxiety (Tot. Anx) was estimated for children.

2.3. Materials for Parent Assessment

In the questionnaire for parents, in addition to gender, age, and other personal data (educational qualification, level of income, type of profession, and marital status), some questions were asked about their lifestyle and diet and the way in which they fed their children. Regarding lifestyle, parents were asked if they played sports, if they followed a diet, and if they were overweight. Regarding the way in which they feed their children, they were asked if they encouraged their children to finish the meal, if they tended to cook the dishes preferred by their children, and if they let their children decide the food for lunch or dinner. After these questions, parents had to fill in three standardized questionnaires to assess the eating habits of their children, their level of stress, and the presence of behavioral and psychological disorders in their children.

2.4. The Children’s Eating Behavior Questionnaire (CEBQ) [25]

The CEBQ is a parent-report questionnaire designed to assess eating styles related to obesity risk. The questionnaire is composed of 35 items divided into 8 subscales, each of which assesses a particular eating style: food responsiveness (FR), emotional overeating (EOE), enjoyment of food (EF), desire to drink (DD), satiety responsiveness (SR), slowness in eating (SE), emotional undereating (EUE), and food fussiness (FF). FR evaluates if the child has a positive reaction to food (e.g., “My child’s always asking for food”); EF evaluates how much the child likes to eat (e.g., “My child enjoys eating”); EOE evaluates the child’s tendency to overeat when agitated or worried (e.g., “My child eats more when anxious”); DD evaluates the child’s tendency to drink (e.g., “My child is always asking for a drink”); SR evaluates the level of satiety reached by the child during meals (e.g., “My child gets full up easily”); SE evaluates how fast the child is in finishing eating (e.g. “My child eats slowly”); EUE evaluates the child’s tendency to eat when he is agitated, sad, or angry (e.g., “My child eats less when s/he is angry”); and FF evaluates the child’s tendency to be picky when eating (e.g., “My child is difficult to please with meals”). Responses are on a Likert scale from 0 (never true) to 2 (certainly true).

2.5. The Parenting Stress Index-Short Form (PSI-SF) [26]

The PSI-SF is composed of 36 items, with Likert scales ranging from 1 (strongly disagree) to 5 (strongly agree). The test is divided into 3 subscales, each made up of 12 items, which measure different aspects of parental stress: parental stress (PD), parent-child dysfunctional interaction (PCDI), and the difficult child (DC) subscale. The PD subscale measures the parent’s overall stress level; the PCDI subscale measures the quality of the parent-child interaction, in the sense that the higher the score, the lower the quality of the interaction; and the DC subscale measures the child’s level of problematic behavior. There is also a global score on the PSI-SF scale.

2.6. The Child Behavior Check List-6/18 (CBCL 6/18) Parent Form [27]

The CBCL 6/18 is a questionnaire composed of 113 items that explore various areas: activity, interest, attention, fear, play, interaction with peers and adults, anxiety, somatic conditions and problems, mood, aggression, and affective responsiveness and response to changes. CBCL is divided into 9 subscales that are: introversion (INT), psychosomatic disorders (SOM), anxiety and depression (ANX), socialization (SOC), ideation (IDE), attention (ATT), deviant behavior (DEV), aggression and hyperactivity (AGG), and other problems (OTHER). Scores are on a Likert scale from 0 (not true) to 2 (completely true). The sum of the INT, SOM, and ANX subscales generates the “internalizing” CBCL score (INT-CBCL), while the sum of the DEV and AGG subscales generates the “externalizing” CBCL score (EXT-CBCL). The sum of INT-CBCL, EXT-CBCL, and OTHER subscales gives the global CBCL score. For each subscale of the CBCL and for the internalizing and externalizing scales, there are cutoff values to determine if the child is normal, borderline, or clinical [28].

2.7. The Marlowe–Crowne Scale (MC) [29]

The MC is a control scale for assessing social desirability. If there is a positive correlation between the MC and other test scores, then subjects tend to overestimate their trait or characteristics; if there is a negative correlation, then subjects tend to underestimate their traits or characteristics. We used a short form of the MC [30], with nine items on a Likert scale from 1 (completely false) to 5 (completely true).

2.8. Procedure

Children and parents were recruited through a series of meetings between the research staff and families organized in collaboration with the teachers of the scholastic complex. In these meetings, the purposes and procedures of the research were explained both to parents and teachers of the 4th and 5th grade classrooms. In these meetings, we asked parents who wanted to participate to communicate their availability to the teachers who had the task to write on a list the names of the parents who accepted to participate. Participation was voluntary and not paid. After having collected the subscriptions from the research volunteers, some members of the research staff (A.S. and F.T.) went regularly to the school on specific days of the week, according to teachers’ consensus and availability, to measure the BMI of the children whose parents accepted to collaborate and to administer them a short questionnaire. When children ended the compilation, they were given another questionnaire to be filled out by their parents at home. Parents were also given instructions for returning the questionnaire within a certain period of time (two weeks). Parents’ compiled questionnaires were returned by the children to their teachers. Subsequently, the researchers moved on to the classes involved in the study to collect all the questionnaires. The total time required for administering and collecting all the questionnaires was approximately 3–4 months. The data collection occurred in the last months of 2019, shortly before the lockdown imposed by the COVID-19 pandemic. Parents signed informed consent for their participation. In the consent sheet, they were given information about respect for privacy, according to the Italian and European laws (Italian law no. 196/2003 and EU GDPR 679/2016, respectively). The authorization (ref: 185046, 10/04/2019) for the execution of the research was given to us by the head teacher manager of the Comprehensive Scholastic Institute “E. Fermi” of Alba Adriatica (TE) who, according to Italian law, is primarily responsible for guaranteeing the ethical rights of children attending school in concomitance with the authorization of the experimental protocol from our university department (27/07/2019).

2.9. Statistical Analyses

We performed different statistical analyses. Descriptives and frequencies for demographic characteristics of the children and parents sample, of children’s BMI and eating and life habits, and of parents’ lifestyle and children’s feeding styles. We reported descriptives and consistency measures (Cronbach’s alpha) for children’s and parents’ assessments on psychological tests, and bivariate correlations were estimated between mothers’ or fathers’ psychological measures. A multigroup confirmatory factor analysis (MG-CFA) with a multiple indicator multiple cause model (MIMIC) was computed for testing measurement invariance between mother and father. We tested the measurement invariance using two CFA models, one in which the criterion is a factor representing children’s mental state, measured by INT-CBCL, EXT-CBCL, and the OTHER, which are the components of the global score of the CBCL scale, and the other in which the criterion is a factor representing the parental stress measured by PD, PCDI, and DC, which are the components of the global score of the PSI scale. The predictors for the MIMIC models were the CEBQ subscales, which showed significant correlations with the total score of CBCL and PSI, both for mothers and fathers.

3. Results

Table 2 shows the frequencies and percentages of the values relative to the anthropometric data of the children, their lifestyle and food attitude, and the quality of their relationship with their parents. It also reports the frequency and percentages of lifestyle (sports, diet, and overweight) and feeding style of parents.

From the data in Table 2, it emerged that 78% of children played sports, 58% declared that eating is a pleasant experience, and 82% did not feel particularly overweight compared to their peers. On the basis of these declarations, the lifestyle of the 125 children can be considered healthy. This is also evidenced by the fact that only 13% of children had BMI values reporting obesity (only 5% of children were underweight). Regarding the quality of their relationship with their parents, most children reported that parents occasionally or often play with them (66% of mothers and 74% of fathers) and that parents were satisfied with them (69% of mothers and 68% of fathers were often satisfied with their children). Therefore, the relationship between parents and children showed a satisfactory quality. In relation to parents’ lifestyle, 65% of them said they practiced sports at least occasionally, 75% avoided overeating, 52% followed a diet, and 93% had a normal weight or were moderately overweight. Therefore, parents also showed a substantially healthy lifestyle. Regarding parents’ feeding habits, 27% of them often tried to convince their children to finish lunch or dinner, 82% cooked the dishes preferred by their children, and 62% declared that they let their children decide the food for lunch or dinner. Our data showed that parents’ feeding habits were sufficiently healthy. There were no substantial differences between mothers and fathers in their lifestyle and feeding habits.

Table 3 shows the descriptive quantitative data of children (BMI score, Gen Anx, and Sep Anx) and of CEBQ, PSI-SF, and CBCL subscales, and total scores divided for mother and father. Cronbach’s alpha and bivariate correlations with MC scores are reported to analyze the reliability of data.

Table 3 reports that the skewness and kurtosis values of the subscales of the various questionnaires, both for children and for parents, are included in the range of values −2 and +2. This means that data distributions of subjective responses are sufficiently normal [31]. Regarding reliability, most of Cronbach’s alpha values are good (>0.79) or sufficient (>0.60). Only for some scales, the reliability values are lower than 0.60, especially in some subscales of the CBCL [32]. It is necessary to say that the children in our sample do not present clinical, mental, or behavioral disorders (less than 7% of children show severe symptoms). Since the CBCL is a questionnaire designed to detect the presence of clinical symptoms, its scores have greater consistency in the presence of severe clinical subjects [33, 34]. Regarding the correlation with the MC scale, both the PSI-SF and CBCL subscales show significant correlations, although all of them were below 0.50, indicating a weak correlation [35]. Thus, there is a weak effect of social desirability on parents’ assessments of children’s mental health and their stress. Such correlations are lower in fathers.

Table 4 shows the bivariate correlations between CEBQ, PSI-SF, and CBCL subscales for mothers (upper triangle) and fathers (lower triangle). In addition, the correlations between pairs of parents (number of pairs = 61) of various quantitative measures are also reported to analyze whether or not there is an agreement in parents’ assessment of their children. The table also reports the correlations between BMI and children’s responses to the anxiety test (Gen Anx, Sep Anx, and Tot Anx).

Table 4 shows the existence of significant intracorrelations between the PSI-SF and CBCL subscales. There is a strong convergence between the subscales that measure children’s mental health and parental stress. There is also a strong intercorrelation between the CBCL subscale and the PSI-SF subscale. Other studies also showed that parents suffered from high levels of stress when their children had behavioral or psychological disorders [34]. With regard to the subscales of the CEBQ, there are no high intracorrelations between the subscales. The intercorrelations with the subscales of other questionnaires are significant for the CBCL and the PSI-SF subscales. Considering only the global score of the CBCL and PSI-SF, only the FR, EOE, DD, EUE, and FF subscales show significant correlations. Therefore, among the different eating habits measured by the CEBQ, the positive or negative reaction to food (FR or FF), the tendency to eat in relation to the emotional state (EOE or EUE), and the tendency to feel thirsty (DD) are the ones that are highly associated with children’s mental health and parental stress.

Regarding the correlations between parental stress, children’s mental health, and children's anxiety or BMI, no significant correlations emerged. It should be noted that only for mothers there are some significant correlations between some subscales of the CEBQ (FR, EOE, EF, SR, and SE) and the subscales of the PSI-SF with the BMI of the children.

Table 4 also shows the correlations between mothers’ and fathers’ assessments regarding the CEBQ, PSI-SF, and CBCL subscales. Almost all PSI-SF and CBCL subscales showed significant, although not very high, correlations (the highest correlation is 0.60 for INT-CBCL) both for mothers and fathers. The CEBQ subscales with significant correlations between mothers and fathers are DD, SE, EUE, and FF. Therefore, there is no strong convergence between mothers’ and fathers’ assessments of their children’s eating habits.

There are strong intracorrelations between children’s assessments of anxiety (Gen Anx, Sep Anx, and Tot Anx), but none of these are significantly correlated with BMI.

Bivariate correlations showed that FR, EOE, DD, EUE, and FF are the eating habits that have a significant correlation with children’s mental health and parental stress. All these subscales are indicators of unhealthy eating habits. Therefore, these subscales are used as predictors of the global level of children’s mental health and parental stress. Global mental health is measured by the indices of internalizing (INT-CBCL) and externalized (EXT-CBCL) symptoms and of other symptoms (OTHER) that make up the global CBCL score, while the global parental stress is measured by the three subscales of the PSI-SF (PD, PCDI, and DC). Figure 1 shows the MIMIC models tested in our analyses.

The measurement invariance of the MIMIC models was tested using the chi-square difference between the model with free parameters and fixed parameters. The first model (Mfree) is the comparison model with free parameters for both groups. The second model tested loading invariance between groups (Mload). The third model tested loading and intercepts invariance between groups (Mload+int). The fourth model tested loadings, intercepts, and regression coefficient invariance between groups (Mload+int+reg). If the chi-square difference (Δχ2) between nested models is not significant (), then the measurement invariance is confirmed [36]. We also calculated the comparative fit index (CFI) of each model, and we used the difference of CFI (ΔCFI) to test measurement invariance. According to Chen [37], a difference of less than 0.01 (in absolute values) indicates measurement invariance.

Table 5 shows the results of the measurement invariance analysis between the MIMIC models. All Δχ2 are not significant, and all the ΔCFI were <|0.01|, confirming the measurement invariance of MIMIC models. Therefore, there is an agreement between mothers and fathers in judging the mental health problems of their children and their own level of parental stress. Figure 1 shows the standardized parameters (loading and regression coefficients), both for mothers and fathers, and the relative significance (z-test) of each parameter.

Among eating habits, only FR resulted to be significantly predictive of mental health (CBCL), and only DD resulted to be significantly predictive of parental stress (PSI). FR is the principal eating behavior connected to presence of behavioral or psychological problems in children, while DD is the eating behavior with the strongest connection with parental stress.

4. Discussion

Our study has shown that parents’ assessments of their children’s eating disorders are consistent and reliable. Their assessments of children’s mental health and parental stress are also reliable.

Hypothesis H1a was partially confirmed. Children’s eating habits are not affected by children’s BMI, except for mothers. Our study shows a higher sensitivity of mothers toward the physical fitness of their children, indicated by the significant correlations between children’s eating habits and parental stress on the BMI index, while fathers seem to neglect this connection. The greater attention of mothers to the physical fitness of their children has also been highlighted by other studies [14, 16]. Probably, it is due to mothers’ and fathers’ different social and cultural attitudes regarding the physical aspect of their children [10, 16, 17]. This seems the most probable explanation because, in our sample, both mothers and fathers apply the same way of feeding their children unlike the evidences of other studies [1416, 22]. Despite these differences, however, mothers and fathers both prove to be good observers of their child’s eating behavior.

Hypothesis H1b was not confirmed. Children’s anxiety (general and separation anxiety) is not correlated to their eating habits.

Hypothesis H2a was confirmed only for FR, and hypothesis H2b was confirmed only for DD. Eating habits have the same predictive value on children’s mental health and parental stress, both for mothers and fathers. In particular, children show better mental health if they like and enjoy eating. The high desire to drink in children is significantly related to a higher level of stress and discomfort in parents. Therefore, according to parents’ assessment, children’s positive attitude toward food is associated with good and well-adapted behavior, while excessive desire in eating or drinking can generate anxiety and stress in parents. There are no differences between mothers and fathers in their cognitive or emotive reactions to their children’s eating habits. Substantially, our data confirm the results of previous studies that found a positive association between children’s mental health and healthy eating habits [57]. Our data also confirm the negative impact of children’s unhealthy eating habits on parents’ psychological well-being [10]. The fact that only some dimensions of children’s eating habits are related to children’s mental health and parental stress can be explained by the low percentage of children with serious eating problems in our sample.

A limit to our research is due to the fact that the subjects are all part of a restricted geographical complex, given that they were selected within a specific scholastic complex. However, it is quite difficult to recruit volunteers for research purposes without the help of an institution operating on the territory such as a school in this case. Another point is that this study should be extended to children who really suffer from severe eating disorders to see how, in this case, parents’ assessment of their children’s eating habits and behavior change in relation to their problems. It must be said, however, that from a preventive perspective, it is necessary to involve also families with healthy eating habits as control groups in programs for the identification of possible forms of unhealthy eating behavior to prevent the appearance of more serious problems, which then require medical and/or psychological intervention or rehabilitation [4, 7, 10].

5. Conclusion

The results of our study show that parents can play an important role in assessing their children’s eating habits. Parents can assist pediatricians, psychiatrists, and clinical psychologists of the developmental age in diagnosing the presence of eating disorders, which can then affect both the behavior and mental state of children and the psychological well-being of parents. The risk of obesity is a real risk in today's society [38] [2]. In particular, in Italy, 32.4% of male children and 30.9% of female children are overweight or obese [39]. A quick and effective intervention by pediatricians, psychiatrists, and psychologists is essential to avoid the appearance and persistence of unhealthy eating habits [4]. Parents can play an important role in preventing obesity, both because they give primal care to their children’s health and because they are constantly in contact with them. Parents, unlike other caregivers or professional figures, have a greater ability to grasp and detect the presence of an eating disorder in their children [10, 16]. In addition, variations in life habits due to exceptional events, such as the COVID-19 pandemic, that forced people to live in isolation, make every single family essential to ensure correct eating habits in children [40].

Other studies are needed to define better intervention programs to prevent childhood obesity, but our study highlights that in these programs, the involvement of parents in the observation and assessment of children’s eating habits is necessary and affordable.

Data Availability

The data sets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Ethical Approval

Research involving them has been performed in accordance with the Declaration of Helsinki. The research was approved and authorized (ref: 185046, 10/04/2019) by the Comprehensive Scholastic Institute “E. Fermi” of Alba Adriatica (TE), in which participants were recruited. . The privacy of participants was guaranteed according to the Italian and European laws (Italian law no. 196/2003 and EU GDPR 679/2016, respectively).

All participants read and signed informed consent before taking part in the study.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

M.T. projected the work and carried on the statistical analyses. F.T. collected the empirical data and handled the collaboration between universities and schools. A.S. collected the empirical data and helped with statistical analyses. A.S. supervised the work and helped with statistical analyses.

Acknowledgments

The authors thank Prof. Nadia Di Gaspare, the head teacher and manager of the Comprehensive Scholastic Institute “E. Fermi” in Alba Adriatica (TE), the 4th and 5th grade classroom teachers of the primary school, and the families (parents and children) attending the school for their kind attention and for their time and commitment.