Development and Validation of a Skin Fold Thickness Prediction Equation for Asian Indians Using

Volume 2 | Issue 1 Development and Validation of a Skin Fold Thickness Prediction Equation for Asian Indians Using Hydrodensitometry Shenoy S*1, Shrivastava S1, Sandhu JS1, Malhotra B1 and Gupta SK2 1Department of Sports Medicine and Physiotherapy, Guru Nanak Dev University, Amritsar, India 2Indian Institute of Statistics, New Delhi, India *Corresponding author: Shenoy S, Associate Professor, Department of Sports Medicine and Physiotherapy, Guru Nanak Dev University, Amritsar, India 143001, E-mail: drshweta.sportsmed@gmail.com Research Article Open Access


Background
Considering the global burden of obesity and its co morbidities, emphasis is now shifting towards a preventive approach as well as risk stratification for the identification of high risk individuals and lifestyle changes [6,7].Because obesity is a common denominator in the development of lifestyle disorders, accurately calculating body fat % using field methods has now become a necessity.Several models such as Durnin and Womersley 4 sites, Jackson and Pollock 3 sites and 4 sites equation, BMI body fat %, Garcia equation are widely used for the estimation of body fat %, as convenient methods.Several studies have shown that these models sometimes underestimate [8,9] or overestimate body fat [10] in different populations and these models have not been validated for use in the Indian population.Considering the above factors, we derived a new prediction equation for the estimation of body fat % using skinfold measurement, as well as common kinanthropometric measurements using hydrodensitometry as a gold standard [11,12].We also explored the use of DEXA in the evaluation of body fat, as well bone mineral density (BMD) for evaluation of significance and possible inclusion in the equation.137 healthy subjects 77 males and 60 females were included in the study within the age group of 30-60 years.The study was conducted in the months of May-November 2012.The procedure was explained and an informed consent letter was obtained from the patients before performing the test.The study was approved by the Institutional Ethics Committee.Body weight was measured by a weighing machine in kilograms (kgs).Height was measured by stadiometer without foot wear and BMI was calculated by weight (kg)/height 2 (m 2 ).Skinfolds measurements were taken from the right side using Harpenden skinfold caliper (Harpenden Ltd, England, United Kingdom) and sites involved biceps, triceps, subscapularis, chest, abdominal, suprailliac, thigh, knee, calf and chin.Circumference of waist, hip, thigh, calf and arm were measured with a soft tape measure to the nearest 0.1cm.Breadths (chest, elbow, knee, ankle and wrist) were measured using widespread vernier caliper (Holtain Ltd, Croswell, Crymch, United Kingdom).Each anthropometric site was measured using the techniques described in the Anthropometric Standardization Reference Manual [13].

Methodology
Secondly, the relationship of body density to skin folds may not be rectilinear because of a larger portion of body fat which is deposited subcutaneously with increasing obesity [14].
Durnin and Womersley 1973 [14], explains why the Logarithm of each skinfold thickness equation was taken rather than actual measurement reasons are as follows: Firstly, frequency distribution of most of skin fold measurement in general population is skewed, with a long tail of high readings (although this may well represent a pathological state of common obesity).

Hydrostatic under-water weighing
Body fat percentage was assessed using the hydrostatic underwater weighing machine" Vacumed Turbo fit 5.10" (www.vacumed.com) as described in a previous study from this laboratory [2].Procedure was explained and participants practiced before taking the final readings.An average of three readings was taken as the final reading.The software estimated the residual lung volume using the following equations: Male residual lung volume = Vital Capacity × 0.24 Female residual lung volume = Vital Capacity × 0.28 The % BODY FAT was then estimated by the equation of Brozek et al. [15].DEXA Hologic Discovery Wi was used as a non-invasive technique for whole body and segmental fat-, bone-and lean-mass and bone density [16].The Hologic Discovery Wi with fan beam features (Hologic) and array mode was used for all DXA scans.The characteristics and physical concepts of DXA measurements have been described previously elsewhere [17].The scans were analyzed with the most recent software APEX version 4.0.The scan measurements and analyses were conducted following standard 9 analysis protocol as described in the Hologic User Manual.All scans were subsequently analyzed by a single trained investigator [18].

DEXA
The whole data set consisted of 137 subjects with 77 males and 60 females.The data set was then randomly divided into a prediction and validation set in a ratio of 75:25 (i.e.data of randomly selected 55 males and 49 females were used to predict the model and remaining 33 subjects i.e. 22 males and 11 females were used to validate the model).The descriptive statistics of both males and females along with its significance status by one-way sample t test were used to summarize the continuous variables.Apart from descriptive statistics, ANOVA, Pearson correlation and multiple linear regression was done to create a model.

Results
137 healthy subjects of age group of 30-60 years, both male and female were included in the study.Descriptive statistics summarizing the continuous variables along with its significance are shown in (Table 1 and 2).All the variables which were found to be significant by the application of one way sample t test were then put into the multiple linear regression in SPSS version 20.0 which showed all the variables in our models, it also showed the sets of variables for several models at once, by step wise and backward methods of linear regression.We found 2 best models for both males and females by these methods.Model summary are described in table 3 for both males and females.Table 3 shows that the R Square for males is 0.90 (model 1) and 0.71 (model 2); and for female R Square is 0.90 (model 1) and 0.717 (model 2).The error of estimate for body fat calculation for males is 2.85 (model 1) and 3.24 (model 2) and for females it is 2.24 (model 1) and 2.74 (model 2).A Table 4 indicates the results of analysis of variance (ANOVA).In these tables considering the two sums of squares introduced in class -the regression and residual (or error) sums of squares indicates that the variance of the residuals (or errors) is the value of the mean square error (MSE) which is 8.26 and 10.47 for model 1 and 2 respectively in males and 5.84 and 7.5 for model 1and 2 respectively for females.These values show that the error rate is less in model 1 in both males and females.This table also summarizes results of the F test.This tests the hypothesis that the predictor (here our only predictor) shows no relationship to Y scores.The value of the test for our data is F is Where skf is skin folds, (W*H) is a product of waist circumference and hip circumference in cm, (log 4) includes sum of skin folds of biceps, triceps, subscapular, suprailiac in mm with log base 10, (log 7) includes sum of skin folds of biceps, triceps, subscapular, abdominal, suprailiac, thigh, calf in mm with log base 10, (TRI *SS) is a product of skin fold of triceps and subscapular in mm.(SI * SS) is a product of subscapular and suprailiac in mm.MODEL 2: BODY FAT % = -27.245+ 32.43 (LOG 5), where LOG 5 is a sum of skin folds of biceps, triceps, subscapular, suprailiac and calf in mm with log base 10.Unstandardized coefficient for (suprailiac skf) 3 was very low i.e. -4.94E-005 and was therefore not included in the model.

Comparison of new predicted equations with hydrodensitometry, DEXA and other existing skin fold equations
The results of body fat % calculated from hydro densitometry, DEXA, existing equations and new predicted equations are shown in Table 5A.In the new predicted equations both model 1 (males 33.24 ± 6.8 and females 36.11± 7) and model 2 (males 34.13 ± 5.09 and females 35.57± 3.7) are closely related to hydrodensitometry (males 33.21 ± 5.85 and females 35.99 ± 4.9), which is the gold standard for body fat % estimation.The error differences in model 1 (for males -0.039 ± 5.14 and for females -0.42 ± 4.03) and in model 2 (for males -0.93 ± 3.88 and females -0.121 ± 7.05) are also comparatively lower for both males and females shown in Table 5B.Based on our results, we found that model 1 and 2 are most closely related to hydro densitometry followed by DEXA, DWW, JP 3, JP 4, Garcia, and BMI bf % in males and DWW, BMI bf %, DEXA, JP 4, JP 3 and Garcia in females.
Tables 6A and 6B show the mean values and standard deviations of BMI, BMD, WHR and % BF derived from selected models classified according to BMI categories in males and females respectively.They clearly show that the newly developed models give a higher value of percent body fat than previously used equations.Further in females as percent body fat increases the BMD tends to increase.The assessment of body composition to mark obesity becomes very important to rule out diseases related to body fatness [19].All popular skinfold thickness equation models currently in use, were derived and validated using data obtained from diverse populations, but did not use data primarily of Asian Indians.In this study we evaluated data from several sources including DEXA and kinanthropometric parameters using hydrodensitometry as the gold standard.
The anthropometric parameters used to assess body fat included, skinfold thickness, BMI and waist to hip ratio, an indicator of intra-abdominal obesity [19].The sites of skin folds we used to calculate body fat, included biceps, triceps, subscapular, chest, abdominal, suprailiac, thigh, knee and calf.Circumferences (arm, waist, hip, thigh and calf), bone breadth (elbow, wrist, chest, knee and ankle) ,body density using hydrodensitometry was calculated by mass/volume which is a function of additive densities of four components including fat and those of fat free body consisting of water, mineral and residual components (which consist mostly of protein Though it was not our primary aim to evaluate BMD, we observed a positive relation between BMI and BMD in women.We also found that the BMD of our subjects for comparable age group in Caucasians were much lower [20].Earlier studies have indicated that there was a positive relationship; i.e. as body fat % increased, BMD also increased [21][22][23][24][25][26][27] and in females this observation held true.This was attributed to the fact that with an increase in body weight, the mechanical loading on bones would increase [20].
A limitation of this study was that this study was a single centered trial and its validation in other geographical areas needs to be studied.This could not be accomplished because to date there is no other center in India, which utilizes hydrodensitometry for the estimation of body fat.Another limitation is that we did not take into account menopausal status of women, which would affect body fat accumulation.However, we desired to formulate an equation, which would suit a broad age group of women.Future work in this area could take into account menopausal status, which could be incorporated into the equation.Further it is important to mention here that in our study we had great difficulty in getting volunteers to participate in the process of underwater weighing for cultural as well as other reasons, which included a phobia of staying underwater without breathing, restriction in range of motion at hip and knee.Almost 40 % of those who initially volunteered had to be finally excluded because they were unable to perform complete exhalation and then stay underwater for the required time.

Discussion
Given the difficulties associated with hydrodensitometry many investigators now use DEXA, which is based on the 4-compartment model for the calculation of body fat percentage and bone mineral density, given its non-invasive nature as well as convenience.However in DEXA [28][29][30][31][32] body fat is calculated on the constancy assumption that ≈ 73 % of lean body mass (lean + BMC) is water [29][30][31][32].However the tissue water content in adipose tissue is highly variable e.g.ranging from ± 17 % to 50 % in humans [30][31][32], which can thus affect the accuracy of the result.Despite these limitations, the relatively better accuracy of the equation developed by us demonstrates the need for the development of population specific equations for the assessment of body composition.
We developed a new prediction equation for assessing the total body fat by using hydrodensitometry as the gold standard method to assess body composition with anthropometric measurements.Further validation and cross validation revealed that the new model was more significant, as its error difference is less as compare to the existing models.Further validation of this equation is recommended.

Conclusion
The authors would like to acknowledge and profoundly thank Prof S K Neogy, Head SQCN and OR Unit, Indian Statistical Institute, New Delhi for his whole hearted assistance in the statistical analysis, and the volunteers for their co-operation and participation in the study.
Submit your next manuscript to Annex Publishers and benefit from: Submit your manuscript at http://www.annexpublishers.com/paper-submission.php→ Easy online submission process → Rapid peer review process → Open access: articles available free online → Online article availability soon after acceptance for Publication → Better discount on subsequent article submission → More accessibility of the articles to the readers/researchers within the field

Bone Breadth Wrist (In Cm)
Shows that P value is <0.001 which is highly significant for all the variables.

Table 1 :
Descriptive statistics for both males and females 6.39 and 63.21 for model 1 and 2 respectively for males and 6.4 and 27.87 for model 1 and 2 respectively for female 's.The table shows us this is significant (p < .001).As the value of F is both large and small, we determine that our predictor of body fat % outcome is related to body fat % in our population.Shows that P value is <0.001 which is highly significant for all the variables.The comparison is within a gender i.e. between males of the prediction and validation set and females of the prediction and validation set

Table 2 :
Descriptive data of prediction set and validation set for both males and females

Table 3 :
Model summary for both males and females using correlation between Model 1 and Model 2

Table 4 :
Analysis of variance for predicted models with hydrodensitometry as a dependent variable for both males and females

Table 5B :
Standard Error Differences between Hydrodensitometry, DEXA, Existing Equations and New Predicted Equations

Table 6A :
Mean values and standard deviations of BMI, BMD, WHR and % BF derived from selected models classified according to BMI categories in males

Table 6B :
Mean values and standard deviations of BMI, BMD, WHR and % BF derived from selected models classified according to BMI categories in females ).Using key variables, which were shown to be significantly associated with body fat % derived from hydrodensitometry, with the help of linear regression models, we developed a prediction equation in a wide age group of Asian Indians i.e. from 30-60 years of age.The newly developed skin fold equation is important for 2 reasons: firstly, it was validated and cross validated using hydrodensitometry as a gold standard and secondly, no previous studies have been conducted in a wide age group of Asian Indians.It is clear from our results that model 1 is very accurate and is best for deriving a very close estimation of body fat %.If investigators desire a highly precise body fat estimation with its R 2 of 0.90 and 0.917 for males and females respectively, they can use model 1.However this clearly entails measurements from 20 physical characteristics.Considering economy of time and convenience, model 2, which entails the use of only 3 and 5 physical characteristics, in males and females respectively, is recommended.This model 2, with an accuracy of R 2 of 0.71 and R 2 of 0.84 for males and females respectively, in fact is greater than the R 2 values we obtained for our population using DEXA and other existing skin fold equations.After cross validating model 2, we compared the values of body fat % obtained using this equation, with popular skin fold prediction equations, the Durnin and Womersley equation, Jackson and Pollock 3 site and 4 site equations, Garcia equation, BMI body fat % equations , body fat % calculated by DEXA.We found that the R value of our model 2 with hydrodensitometry was more than all existing skin fold equation and DEXA.