Have a personal or library account? Click to login
Correlation of Different Anthropometric Methods and Bioelectric Impedance in Assessing Body Fat Percentage of Professional Male Athletes Cover

Correlation of Different Anthropometric Methods and Bioelectric Impedance in Assessing Body Fat Percentage of Professional Male Athletes

Open Access
|May 2021

Full Article

INTRODUCTION

Body composition component in form of fat tissue is very important for sports performance (1,2,3,4,5,6,7) It is known that body fat percentage affects the cardiorespiratory ability of male athletes (8), as well that both high and low levels of body fat can put athletes in an unwanted situation (9). Today we are familiar with a large number of methods that can be applied in the analysis of body composition. More precise methods for body composition analysis such as magnetic resonance imaging (MRI), computed tomography (CT), dual-energy X-ray absorptiometry (DXA), plethysmography (ADP) and other second level validity methods are rarely available to us due to financial constraints and the need for trained staff. Because of this, alternative, third level validity methods such as anthropometry and bioelectrical impedance analysis are far more often used on a daily basis for general population, i.e. for non-athletes, people who do recreational sports and athletes. Besides a far more affordable price, the advantages of these alternative methods are also reflected in non-invasiveness, portability and in their simple and fast application in body composition analysis.

On the other hand, one of the biggest disadvantages of these methods is that they are considered less accurate than the aforementioned second level validity methods. In practice we can find numerous anthropometric methods, but the most common ones are those using measured values of skinfolds thickness at defined points of the body and specific mathematical formulas to estimate the body fat percentage. The most common bioelectrical impedances (BIA) are portable models, which today can be seen more and more often in fitness and wellness centers, but also in private homes. Both BIA and anthropometric methods are considered to have a measurement error in estimating body fat of about ± 3.5%, assuming that all preconditions for measurement and analysis are fulfilled (10). However, in situations where athletes have to travel, anthropometric methods still show certain portable advantages. Anthropometric equipment is more practical for transporting compared to BIA analyzers being bigger in size, robust and carry the risk of electronic component damage or failure. Bioelectrical impedance also requires an athlete preparation protocol starting 48 hours prior to measurement which is not practical when going on a trip, while with anthropometry this is not the case. There is also the possibility of temporal disruption of physiological systems caused by air travel or any longer land travel that may lead to somewhat reduced accuracy of BIA, which is not the case with anthropometry where the measurement is performed anatomically. This can influence coaches, sports experts or sports physicians to use anthropometry, especially in situations where athletes often change location and travel to away games, competitions and preparation camps. When it comes to male athletes, Forsythe and Sinning already investigated back in 1973. and pointed out that anthropometric equations developed for the general non-athletic population are not accurate enough in estimating body fat in male athletic population. (11). Throughout decades, numerous anthropometric methods based on sport specificity have been developed through research (12,13,14) and also based on general samples of athletes from different sports, by determining the correlation of these methods with one of the second level validity reference methods (11,15,16,17,18,19) in evaluating body fat of athletes. A smaller number of studies have been conducted in examining the correlation between anthropometric methods and BIA method in general population of male athletes or in sport-specific male population, where second level validity reference method was not included (20,21). However, these studies often examined the correlation of anthropometric methods and BIA method by using only one sports anthropometric method and equation in the study. Determining whether some of the already existing and represented sports anthropometric methods (sports-specific or general sports methods) can replace BIA method in estimating body fat percentage of general male athlete population could be one of the solutions for athletes, coaches and sports experts to avoid confusion and encourage them to use anthropometry in the field of practice, either for traveling needs or using it as personal preference.

In accordance with the previous, the aim of this study is to examine the correlation of different anthropometric methods and BIA methods in estimating body fat percentage of male professional athletes from different sports.

MATERIALS AND METHODS
Participants

The study included eighty five professional athletes (N=85), 17 – 33 years of age, recruited from four different sports, wrestling (n=28), football (n=28), boxing (n=15) and basketball (n=14). All recruited athletes were members of FIT IN health-related fitness club in Belgrade, where they performed part of their conditioning procedures. The inclusion criteria included athletes who were competing in sports for more than 3 years and not having any long training breaks or any rest caused by an injury or any other factor within last six months.

Procedures

Athletes were examined in the laboratory for sports medicine and exercise therapy at the Institute of Physiology "Richard Burian" in Belgrade. Anthropometric measurement procedures was conducted in accordance with the guidelines and recommendations of International Standards for Anthropometric Assessment (ISAK) (22,23). Procedure for BIA analysis was conducted using BIA pre-test guidelines (24). Athletes were not pressured to participate in the study and after being well informed they voluntarily signed the consent form for participating in the study.

Table 1.

Athlete characteristics

SportBasketball (n=14)Football (n=28)Boxing (n=15)Wrestling (n=28)Total (n=85)
VariablesX ± SD rangeX ± SD rangeX ± SD rangeX ± SD rangeX ± SD range
Age (years)25.4 ± 4.223.9 ± 4.223.687 ± 4.221.4 ± 6.423.7 ± 4.3
17.2 – 3117.5 – 3317.4 – 32.217.5 – 32.917.2 – 33

Height (cm)198.3 ± 10.8181.9 ± 5.1186.4 ± 7.2176.7 ± 7.7183.7 ± 10.4
171.2 – 211169.0 – 190.5175 – 198155.0 – 191.5155 – 211

Weight (kg)96.1 ± 19.277.2 ± 5.687.7 ± 17.580.6 ± 11.783.3 ± 14.4
43.7 – 120.667.3 – 86.559.8 – 123.459.5 – 105.243.7 – 123.4

BMI (kg/m2)24.1 ± 323.3 ± 1.125.1 ± 3.425.7 ± 2.524.6 ± 2.6
14.9 – 27.321.2 – 25.319.3 – 32.322.4 – 31.214.9 – 32.3

FFMI (kg/m2)21.7 ± 2.521.2 ± 0.921.9 ± 2.222.4 ± 1.721.8 ± 1.8
14 – 25.419.7 – 22.917.6 – 26.420.1 – 26.114 – 26.4

BIA-BF (%)10 ± 3.28.8 ± 2.112.2 ± 3.612.6 ± 2.810.8 ± 3.3
5.5 – 15.4169.0 – 190.57 – 19.27.4 – 18.44.9 – 19.2

X - mean; SD - standard deviation; n - number of subjects within a group; N - total number of subjects; BMI - body mass index; FFMI - fat-free mass index; BIA-BF - body fat percentage estimated with bioelectrical impedance.

Equipment

Skinfolds measurement was performed using Harpenden caliper (model HSB-BI, produced by company HaB Direct, United Kingdom) with measuring range from 0 to 80 mm (caliper needle can go four full circles around a dial scale with gradation from 0 to 20 mm), measuring pressure on a lifted skinfold of 10 gr/mm2 and reading accuracy of 0.2 mm. Athletes height was measured with roll-up measuring tape with wall attachment (model SECA 206, produced by SECA, Germany), providing measuring range of 0 – 220 cm with 1 mm gradation. Body weight measurement and estimation of body fat were conducted using hand-to-foot type BIA (model InBody 230, produced by InBody, Republic of Korea). All additional needed equipment was prepared and calibrated before measurement took place.

Anthropometric and BIA measurement

Anthropometric measurement of skinfolds was conducted using ISAK guidelines and recommendations (22,23). Bioelectrical impedance analysis was conducted according to the manufacturer's instructions for model InBody 230. Complete testing of every individual athlete was conducted on the same day. Upon testing arrival, athletes performed BIA analysis first, having their body fat percentage assessed and their body weight measured. After BIA analysis was completed anthropometric measurements of body height and skinfolds thickness were conducted. An anthropometrist with a decade of experience in practice was recruited for anthropometric measurement.

Anthropometric methods

Eleven anthropometric methods designed for different male athlete populations were selected for the study and then a total of sixteen skinfold equations (formulas) were selected from these methods (Table 2). The inclusion criterion for anthropometric equations was met with equations being developed for a specific sport population or for general male sport population, by using only measurements of skinfolds or measurements of skinfolds combined with some of the basic anthropometric or descriptive characteristics such as body height and weight, body mass index and age. Also, only equations developed through regression analysis with the highest multiple correlation coefficients between a dependent variable and a group of independent variables (R) or that explained the largest variance in dependent variable by using independent variables (R2) (depending on what was reported in an individual study, R or R2), when these anthropometric equations were correlated with referent methods within a study. Anthropometric measurement included the following ten skinfolds: subscapular, midaxillar, chest (pectoral), abdominal, biceps, triceps, suprailiac, supraspinale, quadriceps and medial calf. Locating and measuring skinfolds was conducted accordance with ISAK guidelines and recommendations (22). Additionally, Siri equation was applied to convert body density to body fat percentage for anthropometric equations that only estimated only body density (29). Determining the test-retest reliability of anthropometric measurement was performed using the method of technical measurement error (TEM) of a measurer, where a deviation of up to 7.5% for skinfolds and up to 1.5% for other anthropometric measures is considered acceptable. The calculation of the TEM was carried out according to the recommendations by Norton (23).

Table 2.

Selected anthropometric methods and equations developed for assessing body fat in different male athlete populations.

Author(s)/methodAnthropometric equation
Yuhasz (25)Equation using 6 skinfolds. %BF = 3.64 + (0.097 (Ch + Tr + Sb + Si + Ab + Th))
Faulkner (26)Equation using 4 skinfolds, today considered a modified Yuhasz method (27). %BF = 5.783 + (0.153 (Tr + Sb + Si + Ab))
Forsyth & Sinning 1 (11)Equation using 2 skinfolds (equation no. 2a). BD = 1.103 − (0.00168 x Sb) − (0.00127 x Ab)
Forsyth & Sinning 2 (11)Equation using 4 skinfolds (equation no. 2b). BD = 1.10647 − (0.00162 x Sb) − (0.00144 x Ab) − (0.00077 x Tr) + (0.00071 x Ma)
Forsyth & Sinning 3 (11)Equation using 2 skinfolds and height (equation no. 3a). BD = 1.02415 − (0.00169 x Sb) + (0.00444 x Ht) − (0.00130 x Ab)
Forsyth & Sinning 4 (11)Equation using 4 skinfolds and height (equation no. 3b). BD = 1.03316 − (0.00164 x Sb) + (0.00410 x Ht) − (0.00144 x Ab) − (0.00069 x Tr) + (0.00062 x Ma)
White et al. (12)Equation using 2 skinfolds. BD = 1.0958 − (0.00088 x Si) − (0.0006 x Th)
Thorland et al. 1 (19)Equation using 7 skinfolds. BD = 1.1091 − (0.00052 (Tr + Sb + Ma + Si + Ab + Th + Ca)) + (0.00000032 (Tr + Sb + Ma + Si + Ab + Th + Ca)2)
Thorland et al. 2 (19)Equation using 3 skinfolds. BD = 1.1136 − (0.00154 (Tr + Sb + Ma)) + (0.00000516 (Tr + Sb + Ma)2)
Withers et al. (28)Equation using 7 skinfolds, not fully published in the original 1987 paper by Withers et al. (28), but can be found in Reilly et al. study (13) derived from Withers et al. data. BD = 1.0988 − (0.0004 (Tr + Sb + Bc + Sp + Ab + Th + Ca))
Evans et al. 1 (18)Equation using 7 skinfolds, gender and race. %BF = 10.566 + (0.12077 (Sb + Tr + Ch + Ma + Si + Ab + Th)) − (8.057 x gender) − (2.545 x race)
Evans et al. 2 (18)Equation using 3 skinfolds, gender and race. %BF = 8.997 + (0.24658 (Ab + Th + Tr)) − (6.343 x gender) − (1.998 x race)
Oliver et al. (14)Equation using 7 skinfolds (equation model number 3). %BF = 3.53 + (0.132 (Ch + Tr + Sb + Ma + Si + Ab + Th))
Reilly et al. (13)Equation using 4 skinfolds. %BF = 5.174 + (0.124 x Th) + (0.147 x Ab) + (0.196 x Tr) + (0.13 x Ca)
Civar et al. (16)Equation using 3 skinfolds and weight. BF% = (0.432 x Tr) + (0.193 x Ab) + (0.364 x Bc) + (0.077 x Wt) − 0.891
Stewart & Hannan (17)Equation using 2 skinfolds and weight. This equation estimates body fat in grams, which are then converted into body fat percentage for BIA comparison. BFM = (331.5 x Ab) + (356.2 x Th) + (111.9 x Wt) − 9108

Ht - height; Wt - weight; BD - body density; BF% - body fat percentage; BFM - body fat mass in grams; Tr - triceps skinfold; Ma - midaxilar skinfold; Sb - subscapular skinfold; Ab - abdominal skinfold; Si - suprailiac skinfold; Sp - supraspinale skinfold; Th - quadriceps skinfold; Ca - calf skinfold (medial calf); Ch - chest skinfold; Bc - biceps skinfold; gender - men = 1; race - African American = 1, Caucasian = 0.

Statistical analysis

Statistical analysis was conducted using SPSS statistical program, package version 25 (30). Linearity for model validity, outliers and data normality distribution were checked using scatter plot graph, Q-Q plot, histogram, skewness and kurtosis and Kolmogorov-Smirnov test. Based on this, examining the correlation of anthropometric methods and BIA method was conducted using Spearman's rang correlation (rs), where values of rs = 0.0 – 0.09) were considered trivial, rs = 0.10 – 0,29 small, rs = 0.30 – 0.49 moderate, rs = 0.50 – 0.69 large, 0.70 – 0.89 very large, 0.90 – 0.99 almost perfect and rs = 1 perfect correlation (31). Descriptive data was described through means and standard deviations (mean ± SD). Statistical significance was set at 0,05. Confidence interval was set at 95%.

RESULTS

Ten anthropometric equations revealed a very large correlation with BIA method (rs = 0.70 – 0.89) while 6 anthropometric equations revealed a large correlation with BIA method (rs = 0.50 – 0.69). All sixteen anthropometric equations showed correlation statistical significance of p < 0.001. Out of the ten anthropometric equations with a very large correlation coefficient, six equations slightly isolated and showed a correlation of 0.760 or larger (Yuhasz, rs = 0.769; Faulkner, rs = 0.792; White et al., rs = 0.761; Thorland et al. 1, rs = 0.760; Evans et al. 1, rs = 0.761; Oliver et al., rs = 0.761). Out of these six equations, along with their highest correlation coefficients, the narrowest bias-corrected and accelerated confidence interval was revealed in equations of Yuhasz with interval length of 0.191 (BCa CI 95% (0.649 – 0.840)), Faulkner with interval length of 0.137 (BCa CI 95% (0.712 – 0.849)) and White et al. with interval length of 0,174 (BCa CI 95% (0.665 – 0.837)), isolating these three equations furthermore from the rest.

Table 3.

Correlation of different anthropometric methods and BIA method in assessing body fat percentage of professional male athletes.

MethodNrsBCa CI 95% lower - upper
Yuhasz (25)vs BIA850.769 ***0.649 – 0.840
Faulkner (26)vs BIA850.792***0.712 – 0.849
Forsyth & Sinning 1 (11)vs BIA850.743 ***0.619 – 0.823
Forsyth & Sinning 2 (11)vs BIA850.738 ***0.627 – 0.019
Forsyth & Sinning 3 (11)vs BIA850.676 ***0.536 – 0.775
Forsyth & Sinning 4 (11)vs BIA850.685 ***0.545 – 0.786
White et al. (12)vs BIA850.761 ***0.665 – 0.837
Thorland et al. 1 (19)vs BIA850.760 ***0.650 – 0.848
Thorland et al. 2 (19)vs BIA850.736 ***0.613 – 0.830
Whiters et al. (13)vs BIA850.674 ***0.502 – 0.796
Evans et al. 1 (18)vs BIA850.761 ***0.642 – 0.845
Evans et al. 2 (18)vs BIA850.674 ***0.518 – 0.791
Oliver et al. (14)vs BIA850.761 ***0.614 – 0.843
Reilly et al. (13)vs BIA850.681 ***0.535 – 0.793
Civar et al. (16)vs BIA850.740 ***0.625 – 0.819
Stewart & Hannan (17)vs BIA850.681 ***0.543 – 0.787

BIA - bioelectrical impedance; rs -Spearman's correlation coefficient; p - statistical significance;

***

- p < 0.001; BCa CI 95% - bias-corrected and accelerated confidence interval, set at 95% confidence.

DISCUSSION

The aim of this study was to identify if any of the existing anthropometric methods developed for male athlete population (general or specific) are applicable to general male athlete population as replacement for BIA method. A trend can be noticed of publishing studies in journals across the world that are developing and suggesting new population-specific anthropometric methods and equations for determination body fat percentage. It is also noticeable that the most numerous of these population-specific methods are developed for specific or general athlete populations (13,14,16,17,32) and for nation-specific populations (specific to anthropology or characteristics of a certain nation) (32,33,34,35,36,37). Publishing and suggesting a large number of anthropometric methods and equations for athlete population, which amounts up to tens of equations developed for popular sports, can be confusing and exhausting for coaches and sport experts when they are making a decision in selecting and applying anthropometric method for their athletes. Also, being exposed to numerous anthropometric methods increases the chances for selecting inadequate anthropometric method for their athletes which can lead to significant mismanagement in regulation of body fat percentage. Examples of inadequate selection of anthropometric methods in practice can be noticed when coaches and sport experts apply popular anthropometric methods developed for general or general male non-athlete population such as Jackson and Pollock (38), Durnin & Womersley (39), Sloan (40) or Lohman (41) for their male athlete population.

The specificity of anthropometric methods should cover parameters or factors of gender, race, age, method and protocol used for skinfold measuring, condition level or competition level, nation (anthropology) etc. All these factors should be taken into consideration during the selection of specific or general athlete anthropometric methods and equations. This means that even if coaches or sport experts decide to select an anthropometric method developed on the sample of athletes from the same sport, chances are that even that method won't be precise enough for their athletes unless all aforementioned factors of specificity match with the athletes which comprised the sample in the development of that anthropometric method. Considering everything said, it can be roughly concluded that the only precise anthropometric method for athletes would have to include coaches or sport experts conducting their own anthropometric measurement and applying regression analysis, thereby developing their own anthropometric equations using only their athletes as the sample. These kind of complications and uncertainties can influence coaches and sport experts to avoid using anthropometry and to prefer the use of BIA method for assessing body fat percentage. As aforementioned, anthropometry also has advantages in field practice. For example, anthropometric equipment takes up a lot less space than BIA instruments, it weights less and doesn't possess any electrical components therefore creating almost no risk from damage during transport. These can be significant advantages of anthropometry compared to BIA when athletes are obligated to travel to preparation camps or competitions, especially if we are considering transporting a higher quality and bigger size (therefore more expensive as well) BIA instrument.

Precision of BIA analysis can be affected by physiological oscillations in human body and potential physiological oscillations during airplane travel or longer land travel, and also it requires athletes to deal with a long and unpopular preparation protocol (48 hours) and requires coaches trust that athlete fully complied with the preparation protocol before BIA analysis, while anthropometry is based on anatomical measurements which require a much simpler athlete testing preparation. In conjunction with the previous, authors of this study considered that it would be useful to examine and provide athletes, coaches and other sport experts with an option of applying an existing anthropometric method to general male athlete population that reveals the closest correlation coefficient to BIA method and therefore can be used as BIA replacement in field practice, either by circumstantial need, personal preference, or simply to eliminate uncertainty and potential selection of inadequate anthropometric methods for specific athletes, which are usually lacking specificity factors.

Anthropometric methods and equations today are numbered by the hundreds. For this study, we selected more common and known anthropometric methods developed for athletes that are used in practice up to a few decades back. Results of all sixteen anthropometric equations revealed a strong correlation with BIA method, which is somewhat interesting considering that some of the included anthropometric methods and equations were developed using samples of athletes coming from different sports, using general samples of athletes, using different reference methods (DXA, MRI, CT, ADP etc.), applying different models for assessing body composition (two-component, three-component, four-component, etc.), probably with large technical error of skin-fold measurement caused by different anthropometry measurers (inter-tester TEM) and protocols, and in the time span of more than 40 years (changes in average athletes body composition and constitution over decades), so it would be more logical that some anthropometric methods reveal much more different correlation coefficients than others in comparison with BIA method.

Results of this study showed that six anthropometric equations with Spearman's coefficient of rs = 0.760 or stronger, somewhat distanced themselves in correlation with BIA method compared to other methods, even though this distance was minimal. Out of these six equations, Yuhasz, Faulkner and White et al. equations separate themselves furthermore at the top with correlation coefficients of rs = 0,769, rs = 0.792 i rs = 0.761, and with bias-corrected and accelerated confidence interval length of 0.191 (BCa CI 95% (0.649 – 0.840)), 0.137 (BCa CI 95% (0.712 – 0.849)) and 0.174 (BCa CI 95% (0.665 – 0.837)). Even though it used to be seen mostly in studies examining statistical difference, lately, a trend of calculating and adding confidence intervals as important correlation indicators can been noticed in correlation studies as well. The reason for this is that correlation coefficients attained in correlation studies refer to our observed sample from a specific population, while for example a confidence interval of 95% tells us that in 95% of cases a true correlation coefficient for our targeted population (in this case professional male athletes) should fall within attained confidence interval range of an individual anthropometric method and BIA method. A shorter length (narrower range) of confidence interval leads to the assumption that the correlation coefficient for our observed sample (rs) must be close to the true correlation coefficient for our specific population, or at least closer than those correlation coefficients with wider confidence interval lengths. Besides revealing the highest correlation coefficients, Yuhasz, Faulkner and White et al. anthropometric equations also revealed the shortest confidence intervals ranges, and alongside Thorland et al. 1 and Civar et al. equations, these three were the only equations revealing confidence intervals shorter than 0.2 and therefore indicated correlation coefficients closer to the true correlation coefficient than other anthropometric methods when correlated with BIA method.

Faulkner's anthropometric method or equation revealed the highest correlation coefficient, reaching almost rs of 0.80 when compared with BIA method. Faulkner's equation was originally considered a method developed for a population of swimmers (26), however motivated by the lack of evidence of it's development, in 2007 Pires-Neto and Graner managed to reveal and point out that this method is actually a modified Yuhasz method (27). Furthermore, by answering the question proposed by Pires-Neto and Graner in 2006 about this methods origin, Faulkner reveals that this equation was modified by Yuhasz himself by combining data of two anthropometric equations which were previously published in his dissertation in 1962. (42), but he couldn't remember how exactly Yuhasz performed this modification (27). Pires-Neto and Graner then suggested that Faulkner's equation should be addressed as Yuhasz's unpublished equation, because evidence revealed that Faulkner did not develop this equation by using a sample of swimmers. By further analyzing Faulkner's equation and comparing it with other studies, among other finding Pires-Neto and Graner state da Faulkner's equation is applicable to young trained male population (27). This part of Pires-Neto and Graner's conclusion matches the results obtained in this study. A young trained male population can be compared with a sample of relatively young athletes (23.7 ± 4.3 years of age) from this study, while a sample of young athletes from different sports in this study (wrestling, football, boxing, basketball) can be compared with the sample used for Yuhasz's method (university martial arts athletes and swimmers, basketball players selected for Olympics and others), or in other words, can be compared with the sample from Yuhasz's population from which the Faulkner's equation was developed.

One of the main limitations of this study was the lack of reference method of second level of validity for estimating body fat, by which besides determining the correlation with BIA method we could also determine if one of the selected anthropometric methods for athletes in this study correlates more accurately than BIA method, for our professional male athlete population. Without a reference method, results of this study can only showcase if one of the existing anthropometric methods is applicable as a replacement for BIA method on general male athlete population. By adding a second level of validity reference method, we would also be able develop a new anthropometric equation based on the athlete sample like the one in this study.

CONCLUSION

This study examined weather any of existing and relatively used anthropometric methods reveal sufficiently similar results to BIA method within male population of professional athletes. This can provide coaches and sports experts an easier selection and enable them an option to replace BIA method with one of the already existing sports anthropometric methods, eliminating at the same time the risk of selecting an inadequate method. Faulkner's anthropometric method (unpublished Yuhasz's equation) showed the largest correlation and the narrowest bias-corrected and accelerated confidence interval with BIA method thus suggesting that this correlation coefficient probably deviates least from the professional male athletes population's true correlation coefficient. This further indicates that it would be very useful and interesting to repeat this kind of study or similar one with the addition of one of the second level of validity reference methods such as magnetic resonance, computed tomography or dual-energy X-ray absorptiometry.

DOI: https://doi.org/10.2478/sjecr-2021-0026 | Journal eISSN: 2956-2090 | Journal ISSN: 2956-0454
Language: English
Page range: 127 - 136
Submitted on: Mar 25, 2021
Accepted on: Apr 17, 2021
Published on: May 29, 2021
Published by: University of Kragujevac, Faculty of Medical Sciences
In partnership with: Paradigm Publishing Services
Publication frequency: 4 issues per year

© 2021 Marko Dimitrijevic, Dijana Lalovic, Djordje Milovanov, published by University of Kragujevac, Faculty of Medical Sciences
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.