- Research article
- Open Access
A field tool for prediction of body fat in Sri Lankan women: skinfold thickness equation
Journal of Health, Population and Nutritionvolume 35, Article number: 31 (2016)
Valid skinfold thickness (SFT) equations for the prediction of body fat are currently unavailable for South Asian women and would be a potentially robust field tool. Our aim was to assess the validity of existing SFT equations against deuterium (2H2O) dilution and, if invalid, to develop and validate an SFT equation for % fat mass (%FM) in Sri Lankan women.
H2O dilution was used with Fourier transform infrared (FTIR) spectroscopy as the criterion method for the assessment of %FM in urban Sri Lankan women (30–45 years). This data was used to assess the validity of available SFT equations and to generate and validate a new SFT equation for the prediction of %FM against the criterion method. Women (n = 164) were divided into validation and cross-validation groups for the development and validation of the new equation. The level of agreement between the %FM calculated by the final derived prediction equation and the %FM obtained by 2H2O dilution was assessed using Pearson’s correlation coefficient (R) and Bland Altman plots. Student’s t test was used to assess over- or underestimation, and significance was set at p < 0.05.
Existing equations significantly (p < 0.001) underestimated %FM compared with the 2H2O dilution method. The final equation obtained was %FM = 19.621 + (0.237*weight) + (0.259*triceps). When compared with 2H2O dilution, %FM by the equation was not significantly different. There was a significant (p < 0.001) correlation between %FM by the reference method and %FM by the equation. The limit of agreement by Bland Altman plot was narrow with a small mean positive bias.
Existing SFT equations were not applicable to this population. The new equation derived was valid. We report a new SFT equation to predict %FM in women of South Asian ancestry suitable for field use.
Obesity is a major contributor to the global burden of disease with the prevalence of overweight and obesity increasing manyfold in recent years in Asia . The projection for the next decade is that the largest proportional increase in overweight and obesity prevalence would be in low income and middle income countries including Asia . Body mass index (BMI) has been the primary method used to determine body adiposity and is currently used to categorize overweight and obesity among adults . Despite its limitations in both the assessment of body fat as opposed to lean body mass as well as in screening for risk of cardiovascular disease, BMI continues to be in widespread use. The reasons for this are twofold; the ease of use in field settings and the lack of validated, accurate field methods to assess body fat. The evidence for body fat related morbidity and mortality and the rise in non-communicable disease prevalence in Asia is strong, emphasizing the requirement for such body fat assessment methods to be developed and used over BMI . Asian populations are at risk of disease at lower BMI values than global cutoff values  and Asians have three to 5 % higher total body fat when compared with Europeans of the same BMI .
Subcutaneous fat is a valid indicator of both regional and total fat mass (FM) and also insulin resistance  with associations with cardiometabolic risk being demonstrated . Skinfold thickness (SFT) is a valid method in the measurement of body density, FM, fat free mass (FFM), and their changes with time, when carried out by a trained observer using a standard protocol . Skinfold thickness, due to its simplicity and low cost requiring only a skinfold caliper and training, makes it suitable for field use. The most widely used prediction equations for the determination of body fat in women are Siri equation  and Brozek et al. equation. . Since Siri and Brozek et al. equations require a prior body density prediction, Durnin and Wormersley  density equation was used. Deuterium dilution method is considered the gold standard method of estimation of total body water (TBW). Fourier transform infrared (FTIR) spectroscopy has previously been shown to be a robust estimation method of detection of deuterium in saliva. Hence, deuterium dilution determined by FTIR can be used in validation of SFT equations for prediction of %FM .
There is much evidence that existing SFT prediction equations are population specific [13, 14]. Presently, there are no published SFT equations for the prediction of %FM in South Asian adults, with existing prediction equations being developed in Western populations. Our aim was to validate existing equations against the gold standard method of deuterium dilution and, if found invalid, to develop and validate a field tool for body fat assessment using skinfold measurements and anthropometry in Sri Lankan women.
Body fat data of 164 women aged 30–45 years were used for the validation of existing equations and for the development and validation of a new SFT equation for the assessment of %FM. The women were part of a larger study on body composition and glycaemic status which recruited 617, urban Sri Lankan women without a previous history of non-communicable diseases, from the community.
The study design of this larger study, published elsewhere, used a random cluster sampling method using the smallest administrative unit within the Colombo Municipal Council area as a cluster .
All anthropometric measurements were measured by one trained researcher using the same equipment in order to minimize variability, under standard conditions according to the International Society for Advancement of Kineanthropometry (ISAK) protocol . All measurements were taken from the left side of the body. Height and weight were recorded while the subjects were standing barefoot. Weight was measured to the nearest 0.1 kg using a calibrated electronic scale (Seca 813). Height was measured to the nearest 0.1 cm using a Stadiometer (Seca 225, telescopic height measurement), after placement of heel, buttocks, back of shoulder and occiput in the vertical plane and head in horizontal Frankfurt plane. Waist circumference was measured in standing position at the end of normal expiration, in the horizontal plane at the level of the narrowest point between the lower costal border and the iliac crest, to the nearest 0.1 cm using a non-stretchable measuring tape (Seca 200). The World Health Organization (WHO) recommended waist circumference cutoff for high risk of metabolic complications (>80 cm) was used to categorize women . All measurements were taken in duplicate and the mean was calculated. Skinfold thickness (SFT) was measured using a Harpenden Caliper (graduation 0.2 mm, range 80 mm, Model: HSK-BI). Biceps, triceps, subscapular, and suprailiac (also known as iliocristale) SFTs were measured to the closest 1 mm. The skin at the appropriate site was pinched to raise a double layer of skin and the underlying adipose tissue but not the muscle. The calipers were then applied 1 cm below and at right angles to the pinch and a reading in millimeters (mm) taken 2 s later. The mean of two measurements was calculated. If the two measurements differed by 0.5 mm, a third measurement was taken and the median value was used .
2H2O dilution technique for body composition analysis
2H2O doses (Cambridge Isotope Laboratories Inc., MA, USA) were previously prepared (30 g each weighed to the nearest 0.01 g) and stored in wide mouthed leak-proof screw capped bottles at 4 °C. Bottles were not reused. Following an overnight fast and on an empty bladder, all measurements were carried out between 8.00 am and 12.00 noon at room temperature, with the participants at rest. The baseline saliva sample was collected to assess basal 2H2O concentration in the body prior to oral administration of the 2H2O dose. Administration was supervised and ingestion of the entire dose was ensured. A second and a third saliva sample was obtained at 3 and 3.5 h post administration, allowing equilibration with body water and mean enrichment was used for further calculations . Entire procedure and quality control were done according to the International Atomic Energy Agency protocol .
Saliva samples were collected, transported at 4 °C and stored at –20 °C pending analysis. All saliva samples were centrifuged at room temperature prior to measuring enrichment, with each sample being measured in duplicate. Quality control was in accordance with International Atomic Energy Agency (IAEA) guidelines. The equation proposed by IAEA for use as a quality control measure is TBW (kg) = 7.4 × height3 (m3), where TBW from each post-dose sample is used with the height of each participant. Since this equation was not derived specifically for our population, the range incorporating the 95 % confidence interval for the above relationship (<5.7 × height 3or >9.6 ×height 3) which has also been proposed by the IAEA for assessment of accuracy was used. Using this criterion, all samples fell within the range and were therefore included. Isotope analysis was carried out in duplicate for each sample using FTIR (8400S; Shimadzu, Vienna, Austria) at St John’s Research Institute, Bangalore, India.
Data of 164 women were included in the analysis. Statistical analysis was carried out using Statistical Package for Social Sciences (SPSS) version 20.
Applicability of selected SFT equations to predict %FM of healthy adult women
From published equations for females, three equations were selected, with one being for South Asian females. These selected equations were tested by comparing %FM obtained through each equation with %FM obtained by the gold standard 2H2O dilution method. Initially density was calculated using Durnin and Womersley equation from the skinfold values of the 164 women (Table 1). These density values were applied independently to Siri and Brozek et al equations to derive %FM. Percentage FM was also calculated by the equation developed by de Lanerolle-Dias et al using the skinfold thickness data from the women (Table 1). Pearson’s rank correlation coefficients were used to assess the association between %FM obtained from 2H2O dilution and %FM obtained from the above equations. Mean values for %FM by 2H2O dilution were compared with %FM values by the equations using Student t test for over or underestimation. The FM% cut off of >35 % for Asian women was used .
Derivation of a population-specific SFT equation to predict %FM using 2H2O dilution as the reference method
Data from the 164 women were randomly divided into validation (60 %; n = 100) and cross-validation (40 %; n = 64) groups by the statistical package (SPSS) for derivation and validation of SFT equations [13, 20]. The preliminary prediction equation for %FM was derived for the validation group, by linear regression analysis, using forward likelihood ratio, with %FM by 2H2O dilution as the dependant variable and weight, waist circumference (WC), SFT-biceps, SFT-triceps, SFT-subscapular, and SFT-suprailiac as independent variables. Predictability of the preliminary equation was evaluated using the cross-validation group. The final prediction equation was derived after combining both validation and cross-validation groups, by linear regression analysis using the enter method [13, 14, 17]. The independent variables for the final equation were determined by the preliminary equation and %FM by 2H2O dilution was used as the dependent variable. Level of agreement between the final derived prediction equation and the relevant 2H2O dilution value was assessed using Pearson’s correlation coefficient (R) and Bland Altman plots. Student’s t test was used to assess over or underestimation and significance was set at p < 0.05.
The mean age of women was 37.4 ± 3.5 years and their body composition parameters are given in Table 2.
Applicability of selected published SFT equations to predict %FM in urban adult women
Selected published prediction equations were tested in the 164 women (Table 1) to validate the %FM calculated by these equations against the reference method. Density obtained by Durnin and Womersley equation was applied to Siri and Brozek et al equations to derive %FM. Percentage FM derived by these two equations and the equation developed by de Lanerolle-Dias et al. significantly (p < 0.001) underestimated %FM by the 2H2O dilution in this population (Table 3).
Derivation and validation of a SFT prediction equation for %FM in normoglycaemics
In deriving a new population-specific SFT equation for the prediction of %FM, FM data by the 2H2O dilution method was used as the reference. Skinfold thickness measurements and correlations with %FM from the criterion method are given below (Table 4).
Using the sample of 164 women for the validation and cross validation, the following analysis yielded the derivation of the SFT equation. The validation (n = 100) and the cross-validation (n = 64) groups were not significantly different in age, anthropometry, and body composition data by the reference method. Percentage FM assessed by the reference method was significantly correlated with weight (R = 0.590, p < 0.001), SFT-biceps (R = 0.471, p < 0.001), SFT-triceps (R = 0.518, p < 0.001), SFT-subscapular (R = 0.580, p < 0.001), and SFT-suprailiac (R = 0.559, p < 0.001). Using %FM by the 2H2O dilution method as the dependant variable and weight, WC, SFT-biceps, SFT-triceps, SFT-subscapular, and SFT-suprailiac as the independent variables, a preliminary prediction equation was derived for %FM using the validation group. Age was not included as an independent variable as age range was narrow in this population. Of the two models obtained by linear regression analysis, the second model had a higher R 2 (model 1: R 2 = 0.557; model 2: R 2 = 0.609) and a lower standard error of estimate (SEE) (model 1, 4.21; model 2, 4.05) and was selected. While SFT-biceps, subscapular and suprailiac were not significant predictors of %FM in this population, weight (P = <0.001) and SFT-triceps (P = 0.003) were reliable independent variables. The preliminary equation obtained was %FM = 20.037 + (0.227 * weight) + (0.257 * triceps). When the preliminary equation was applied to the cross-validation group (n = 64), it was significantly (p < 0.001) correlated (r = 0.693) with %FM by the 2H2O dilution method (Table 5). There was no significant difference between %FM by the 2H2O dilution method and %FM by the preliminary equation (Table 5).
The final equation was derived following the combination of both the validation and cross-validation groups (n = 164). Linear regression analysis was carried out using %FM by the 2H2O dilution method as the dependant variable and weight and SFT-triceps (identified as the most predictive variables from the preliminary equation) as the independent variables. The final equation obtained was %FM = 19.621 + (0.237 * weight) + (0.259 * triceps). There was no significant difference between %FM by the 2H2O dilution method and %FM by the final equation, and %FM by the final equation was significantly (p < 0.001) correlated (r = 0.647) with %FM by the 2H2O dilution method (Table 5).
Visual comparison by scatter plots of the final equation with the reference method indicates high correlation (Fig. 1a). The Bland Altman plot between the final prediction equation and %FM by the reference method for plot difference is shown in Fig. 1b. The limits of agreement was narrow (+7.9 − (−7.9) %FM), and the equation resulted in a small mean positive bias of +0.0128 ± 3.9.
We report a field tool to predict body composition in women of South Asian ancestry, which uses triceps skinfold and weight measurements. The validity of this equation to predict %FM both for field use and for research purposes in women is indicated by the high correlation between the calculated %FM from the equation and %FM by the reference method (2H2O dilution). The strength of the new equation was the small bias and narrow limit of agreement between the actual and predicted values. Our values for the limit of agreement (+6.8 %FM) are better than those reported elsewhere .
The importance of this equation is in the finding that existing SFT prediction equations were not valid for estimating %FM in Sri Lankan women. Since there were no SFT equations for prediction of FM% developed in South Asian adult female populations, existing equations relevant to our population selected to test for validity were Siri , Brozek et al , and Durnin and Womersley equations . The calculation of %FM by both Siri and Brozek et al. equations requires the calculation of body density prior to FM calculation. These published equations correlated well (R = 0.562–0.624, p < 0.001) with the reference method, yet, significantly underestimated %FM in Sri Lankan urban women, making them less suited for use in this population. The risk of using older equations developed on populations whose characteristics may be different to that in the new population has been demonstrated by Leahy et al.  who applied Durnin and Womersley equations to a current population and found that it was less accurate at higher body fat levels. Durnin and Womersley developed their equations in an era when body fat levels were lower on average in a given population. Such differences in body fat are also seen in different ethnic groups. Further, previous validation studies among Sri Lankan adolescents and children reported disagreement between actual and predicted body fat values, necessitating the development of new equations [13, 14].].
Differences in body proportion in Europeans and South Asians may account for underestimation observed in the present study when using existing equations. Durnin and Womersley equation  was developed in an age and gender stratified sample which included participants between 16 and 72 years who were moderately sedentary males and females with some participants being recruited from obesity clinics. The ethnic difference as well as the body composition differences in the populations could account for its inapplicability in our population. Due to ethnic similarity, a published equation for Sri Lankan adolescent girls was tested in this sample . Here, the underestimation could be due to comparatively higher age, mean weight, BMI, and SFTs in adult women in this study.
The need to accurately assess body composition including body fat in South Asians is due to Asians having higher FM for a given BMI compared to Caucasians and the rising prevalence of non-communicable diseases including diabetes mellitus and overweight in Asia, and its link with body fat. Hence, the negative consequence of underestimation of FM during screening of a population is the resultant loss of access to preventive care. The significance of our new equation is in filling this gap of the requirements for accurate and precise measurement of FM in the general population. Our equation being developed in a current population has addressed the concerns of issues arising due to secular trends in body fat levels in populations.
While the simplicity of our equation is due to there being only a single skinfold measurement and weight, the strength of our equation is that triceps and weight measures were retained following the inclusion of four skinfold thickness measurements in generating the equation using regression. Further, as Freedmen highlights in his discussion of Bergman’s body fat index, studying body fat in men and women separately is important due to the difference in distribution of body fat [8, 23]. The new equations for women that were developed in Leahy et al.’s population included abdominal, biceps, triceps, mid-axillary, and medial calf . These equations were developed using a wide range of SFT and girth measurements, and one could expect very high accuracy. Our limits of agreement for %FM compare well with these recent equations of Leahy and co-workers [21, 22]. Thus, for developing countries for use as a field tool, our equation has the simplicity of only having one SFT measurement and weight while preserving accuracy.
The narrow age range of participants may reduce generalizability of the results to women of reproductive age, without further validation and can be considered a limitation. However, the mean BMI of our sample is close to the mean BMI (22.8 ± 4.5 kg/m2) of a nationally representative study population . Further, the BMI range of our population was 15.2–35.9 kg/m2 which reflects BMI observed in Sri Lanka, where morbid obesity is rare. Populations in other South Asian countries also show similar mean BMI values as reported here . Hence, our sample reflects adequate variability, allowing applicability to such populations.
Percentage FM by our SFT prediction equation significantly correlated with the reference method and provided excellent agreement with the reference measure. The feasibility for field use of this equation is enhanced by the fact that it uses only two parameters, namely triceps skinfold and weight measurements. The equation is suitable for use among urban Sri Lankan women and is possibly appropriate to be used in South Asian females, including women of South Asian origin living in Western countries.
Body mass index
Fat free mass
Fourier transform infrared
International Society for Advancement of Kineanthropometry
Standard error of estimate
Statistical Package for Social Sciences
Total body water
World Health Organization
Ramachandran A, Snehalatha C. Rising burden of obesity in Asia. J Obes. 2010;2010:8. Article ID 868573.
Malik VS, Willett WC, Hu FB. Global obesity: trends, risk factors and policy implications. Nature Rev Endo. 2013;9(1):13–27.
Tan KCB. Appropriate body-mass index for Asian populations and its implications for policy and intervention strategies. Lancet. 2004;363(9403):157–63.
Misra A. Revisions of cutoffs of body mass index to define overweight and obesity are needed for the Asian-ethnic groups. Int J Obes Relat Metab Dis. 2003;27:1294–6.
Consultation, WHO Expert. Waist circumference and waist-hip ratio. 2011.
Deurenberg P, Deurenberg-Yap M, Schouten FJM. Validity of total and segmental impedance measurements for prediction of body composition across ethnic population groups. Eur J Clin Nutr. 2002;56(3):214–20.
Addo OY, Himes JH. Reference curves for triceps and subscapular skinfold thicknesses in US children and adolescents. Am J Clin Nutr. 2010;91(3):635–42.
Freedman DS, Blanck HM, Dietz WH, DasMahapatra P, Srinivasan SR, Berenson GS. Is the body adiposity index (hip circumference/height 1 · 5) more strongly related to skinfold thicknesses and risk factor levels than is BMI? The Bogalusa Heart Study. Br J Nutr. 2013;109(02):338–45.
Reilly JJ, Wilson J, Durnin JV. Determination of body composition from skinfold thickness: a validation study. Arch Dis childhood. 1995;73(4):305–10.
Siri WE. Body composition from fluid spaces and density: analysis of methods. In: Brozek J, Henscel A, editors. Techniques of measuring body composition. Washington DC: National Academy of Sciences; 1961. p. 223–4.
Brozek J, Grande F, Anderson JT, Keys A. Densitometric analysis of body composition: revision of some quantitative assumptions. Ann N Y Acad Sci. 1963;110:113–40.
Durnin JVGA, Womersley J. Body fat assessed from total body density and its estimation from skinfold thickness: measurements on 481 men and women aged from 16 to 72 years. Br J Nutr. 1974;32(01):77–97.
de Lanerolle-Dias M, De Silva A, Lanerolle P, Arambepola C, Atukorala S. Body fat assessment in Sri Lankan adolescent girls; development of a simple field tool. Ann Hum Biol. 2011;38(3):330–6.
Wickramasinghe VP, Lamabadusuriya SP, Cleghorn GJ, Davies PS. Use of skin-fold thickness in Sri Lankan children: comparison of several prediction equations. Ind J Peds. 2008;75(12):1237–42.
Waidyatilaka I, Lanerolle P, Wickremasinghe R, Atukorala S, Somasundaram N, de Silva A. Sedentary behaviour and physical activity in South Asian women: time to review current recommendations? PLoS One. 2013;8(3):e58328.
Stewart A, Marfell-Jones M, Olds T, de Ridder H. International standards for anthropometric assessment. Lower Hutt: International Society for Advancement of Kineanthropometry; 2011.
Waidyatilaka I, Lanerolle P, de Lanerolle-Dias M, Atukorala S, de Silva A. Body composition in urban South Asian women; development of a bioelectrical impedance analysis prediction equation. Ann Hum Biol. 2013;40(4):360–7.
International Atomic Energy Agency. Assessment of body composition and total energy expenditure in humans using stable isotope techniques. IAEA Human Health Series No 3. 2009.
Gallagher D, Heymsfield SB, Heo M, Jebb SA, Murgatroyd PR, Sakamoto Y. Healthy percentage body fat ranges: an approach for developing guidelines based on body mass index. Am J Clin Nutr. 2000;72(3):694–701.
Sun SS, Chumlea WC, Heymsfeild S, Lukaski H, Schoeller D, Friedl K, Kuczmarski RJ, Flegal KM, Johnson CL, Hubbard VS. Development of bioelectrical impedence analysis equations for body composition with the use of a multicomponent model for use in epidemiological surveys. Am J Clin Nutr. 2003;77:331–40.
Peterson MJ, Czerwinski SA, Siervogel RM. Development and validation of skinfold-thickness prediction equations with a 4-compartment model. Am J Clin Nutr. 2003;77(5):1186–91.
Leahy S, O’Neill C, Sohun R, Toomey C, Jakeman P. Generalized equations for the prediction of percentage body fat by anthropometry in adult men and women aged 18–81 years. Br J Nutr. 2013;109(04):678–85.
Bergman RN, Stefanovski D, Buchanan TA, Sumner AE, Reynolds JC, Sebring NG, Xiang AH, Watanabe RM. A better index of body adiposity. Obesity. 2011;19(5):1083–9.
Katulanda P, Jayawardena MAR, Sheriff MHR, Constantine GR, Matthews DR. Prevalence of overweight and obesity in Sri Lankan adults. Obes Rev. 2010;11(11):751–6.
Dudeja V, Misra A, Pandey RM, Devina G, Kumar G, Vikram NK. BMI does not accurately predict overweight in Asian Indians in northern India. Br J Nutr. 2011;86(01):105–12.
The authors would like to thank Professor Anura Kurpad, St John’s Research Institute, Bangalore, India, for arranging sample analysis and Professor Rajitha Wickremasinghe, Faculty of Medicine, University of Kelaniya, Sri Lanka, for assisting in study design.
This research was funded by the International Atomic Energy Agency (IAEA)—RC 15920. IAEA had no role in the design, analysis, or writing of this article.
The authors are solely responsible for the design and conduct of this study. IW, PL, AdeS, and SA participated in the design of the study. IW, PL, AdeS, and SA conducted the study, collected the data, and wrote the manuscript. IW, PL, AdeS, and MdeL contributed to the analysis. IW, PL, AdeS, and MdeL reviewed and edited the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
The study protocol was approved by the Ethics Review Committee of the Faculty of Medicine of the University of Colombo, Sri Lanka. Procedures followed were in accordance with the Ethics standards of the responsible committee on human experimentation. All participants were informed about the study both verbally/orally and in writing and written consent obtained.