The assessment of body size in adult populations has been established by several methods using different parameters 123. For example, somatotype of an individual, which provides information of body shape and size, is derived from selected anthropometric traits. The use of methods like, body mass index (BMI), waist hip ratio (WHR), fat mass index (FMI) and conicity index (CI), derived from anthropometric variables are mainly used to assess the nutritional status of a population. However, these techniques are perhaps neither comprehensive to assess the body size nor an appropriate measure of assessing nutritional status of individuals 456. Another study demonstrated using three techniques namely, mid upper arm circumference (MUAC), BMI (construct only two variables) and confirmatory factor analysis (CFA) (construct with more than two variables) in assessing the under nutrition and found that the statistical model, namely, CFA as the best measure compared to the other two 78.

The use of statistical modelling in determining the best model fit of the data, association and distribution pattern among the parameters has come to use during the end of the 20^th century 91011121314 highlighted statistical model for predicting the human head shape and 15 developed a statistical model for estimating the relationship between different anthropometric measures and standing height.

In India, attempts were made towards statistical modelling, concerned with the human body. For example, 1617 proposed a measure of group divergence and generalized distance (D²); 181920 developed theories of statistics using anthropometric traits; 21 used univariate statistical analysis to describe the body composition and distribution of body fat between high and low altitude people of Himalayas. Moreover, 22 showed the relationship among the anthropometric variables at various high altitude populations. Besides these studies, a number of researchers 232425 found the association of anthropometric traits with various parameters of biological and demographic aspects. However, none of the aforementioned studies attempted to provide an accurate statistical model to estimate different types of body size among individuals using anthropometric traits. Thus, it becomes imperative to develop a relatively improved model to capture the maximum variance among the parameters so that the model can be applied to every individual in a community or a population. Such a model needs to be developed with two important purposes (1) models should help in interpreting observed data using the measurement variables in the fitted model to the data and (2) models should be proposed to study the interaction of underlying variables (latent variables) in the future events.

To the best of my knowledge, in the Indian subcontinent, there is hardly any research which used anthropometric traits in developing a composite body size score to determine the body size.

In view of the above, the aim of the present study was to develop a composite size score using anthropometric traits to estimate body sizes and to assess the contribution of each of these anthropometric traits in determining body size.

The study was conducted on a group of adult individuals belonging to a fishermen community inhabiting the coastal villages of East Midnapore district, West Bengal, India.

Anthropometric data were collected as a part of an intramural research project approved by the Technical Advisory Committee of Biological Science Division, Indian Statistical Institute, Kolkata in India entitled “Health and Disease among Population Inhabiting Contrasting Ecological Niches”. It was not feasible to conduct a sampling technique for selecting a random sample due to the lapse of the effects of the major biological factors regarding pedigree on inheritance. Purposive or representative sampling technique was used in order to involve maximum participants in the study. The bias of the data was previously estimated and was found negligible by 26.

At first multiple linear regression analysis was performed using Weight (Wt) as a dependent variable and the rest of the anthropometric variables (measurement variables) as independent variables. The logic is to examine the contribution of each of the anthropometric traits to weight. Standardized predicted weight was estimated for each individual from the regression equation. Standardized weight could be positive or negative.

In the next step, I assessed the maximum standardized predicted weight (say, X) among the individuals and find its corresponding CBSS (say, Y). Hence Y is a first cut off point between robust and normal body size and was designated as the upper limit of normal body size or lower limit of robust body size. The lowest positive CBSS (say, Z) was designated as the lower limit of normal body size. Therefore, individuals having CBSS between a lower limit (Z) and upper limit (Y) were considered as normal body size score and those fell in this range (Z-Y) indicates the normal body size.

Individuals with positive CBSS and greater than the upper limit of normal CBSS (Y) were designated as robust body size.

Finally, individuals with negative CBSS i.e., < 0 (second cut off point) were categorized as lean body size.

Therefore, the study participants were categorized under Robust, Normal and Lean body sizes based on the following criteria (Figure 4).

Robust body size: Body size score > Y (positive value)

Normal body size: Body size score (Z – Y) (positive value)

Lean body size: Body size score < 0 (negative value)

(Table 1) reveals the results of both first order and second order confirmatory factor analyses. This includes dependent variables, their path, latent variables with estimated unstandardized coefficients, standard errors, critical ratio, significant level and standardized coefficient. All the coefficients were highly significant. The standardized estimates allow for the evaluation of the relative contribution of each anthropometric trait to each latent variable. The result shows that for linear latent variable, transverse chest (Tvc) has the maximum factor loading of 0.89; for circular latent variable, chest circumference (Ccn) has the maximum factor loading of 0.87; and for skinfold latent variable, triceps skinfold (Skt) has the maximum factor loading of 0.92. The standard errors are estimates of the errors of unstandardized coefficients to be expected because of the sampling error.

The first-order main sub-construct latent variables (linear and circular) have a significant loading with 6.05 and 0.13 respectively on second-order main constructs latent variable CBSS. These loading values show that linear latent variable has more influence in determining CBSS compared to other two latent variables. Table 1 further mentions the reference point 1.00 (standardized factor loading), which indicates fixing the variance of the corresponding variables within each group; thereby freely estimating the other concerning with the reference variable.

(Table 2) demonstrates the variance of the default model. The estimated variance explained for each of the measurement errors may be an observed or a latent variable associated with the measurement. The result shows that all the variances were highly significant.

(Table 3) represents the fit indices to evaluate the goodness of the model. SEM was used for modification indices to meet the accuracy and precision of structured prediction equations since the fit of the model was not adequate.

The Maximum Likelihood Estimation method is widely used to maximize the fit of the model and also to estimate all model parameters simultaneously. Modification of the model is the general procedure to improve the fit. In this model, I applied the modification indices to maintain the correlation- covariance between the measurement errors derived from the model that provided maximum improvement in fit and the process was continued until to reach an adequate fit. The determination of model fit in SEM technique produce Chi-square statistic which is the test of absolute fit of the model and the Chi-square value should reach p > 0.05. From the modified model, result shows that the overall model fitting was highly significant (likelihood ratio, χ²= 1006.22, p ≤ 0.01) and do not allow for the consistent estimation of parameters in the model involving latent variables. However, 32 has mentioned that the Chi-square statistic is too sensitive to the large size of the sample for it to be interpreted as a significant test. Therefore, the Chi-square statistics was not used as a fit index for this study. Instead, we have turned to other relative fit indices that are comparatively less sensitive to large data set. Therefore, we used modern approaches in relative fit indices such as Comparative Fit Index (CFI), Tucker-Lewis Index (TLI) and Root Mean Square Error of Application (RMSEA) to fit the model. The values of these relative fit indices for the final model are found to be 0.91 (Comparative Fit Index), 0.85 (Tucker-Lewis Index), and 0.19 (Root Mean Square Error of Application). The values were close to the specified ranges 3334 which indicates that the higher-order three factors model provided a moderate fit to the data.

The fit indices expressed as the χ², CMIN/DF, Akaike’s Information Criteria (AIC), Bayesian Information Criteria (BIC ), Comparative Fit Index (CFI), Tucker-Lewis Index (TLI) and Root Mean Square Error of Application ( RMSEA) values were as 1006.22 (p<0.01), 13.97, 1166.22, 1170.15, 0.91, 0.85 and 0.13 respectively.

(Figure 1) shows the standardized visual diagrams of second-order confirmatory factor analysis as well as of structural equation model (SEM). The figure consists of four latent variables- linear, circular and skinfold, and composite body size score, each representing different executive functions. The first three latent variables were measured using sixteen anthropometric traits derived from the first order confirmatory factor analysis. Later CBSS was estimated by the first three latent components derived from second-order confirmatory factor analysis. The figure represents the CBSS as the main construct (called second-order) while linear, circular, and skinfold are three sub-constructs (called first-order). The single-headed arrows that originate from the latent components (first-order) and terminate in the measured variable represent direct relationships from the latent variable to the measurement variables. In second-order confirmatory factor analysis, three arrows originated from the latent variable CBSS and terminated in the first-order three latent variables. In the figures, values shown from bottom to top are error variances (besides error), first-order factor loadings (beside arrows), and error variances of the first-order factor (in the circle), second-order factor loadings (beside arrows) respectively.

(Figure 2) represents the scatter plot on CBSS with standardized weight (measurement variable), where X- axis represents the standardized weight and Y-axis is the CBSS. The line in the scatter plot based on the origin of CBSS indicates a cut-off line (first cut off point). Individuals having scores below the line are either of normal or robust body size. In the scatter plot, at the top of the right corner, at location ‘A’ is the second cut-off point between normal and robust body size. The coordinate of location ‘A’ is (2.88, 3.09) indicates intercepts of maximum predicted standardized weight and corresponding CBSS. Hence, the area of the scatter plot is divided into three parts, the scores that fall in the area of the upper part represent robust body size, and those in the middle part represents normal body size and those in the lower part represents lean body size for this study population.

Based on the cut off points of CBSS, initially, the highest standardized weight (2.88) among the individuals and its corresponding CBSS (3.09) were noted. This value of CBSS (3.09) was designated as the upper limit of normal body size. The lowest positive CBSS (0.001) remained the lower limit of normal body size score. Individuals with CBSS greater than the upper limit of normal body size score (3.09) was considered as robust body size scores and those having CBSS values below 0.001 were considered with lean body size.

On the basis of the proposed cut-off points of CBSS, it has been found that 50.6%, 48.8% and 0.6% of the study participants fell respectively in lean, normal and robust body size categories. It is remarkable that females frequently exhibit lean body size (73.8%) compared to the males (26.2%) in the population as only lean body size scores are take into consideration. Furthermore, result indicates that female participants are more under nutrition than male.

(Figure 3) also reveals the scatter plot on CBSS with age, where X- axis represents the CBSS and Y-axis the age. Since the plotting points were spread all over the plot areas, it indicates that the two variables were independent fitted lines. The figure further reveals that the values of R² were near to zero on both sexes, indicating that the sex is independent of CBSS.

(Figure 4) represents the bar diagram of mean values of all anthropometric traits separately on three types of body size. It appears from the figures that means of all traits of the linear and circular and skinfold components are high to robust body size, followed by normal and lean body sizes.

Since CBSS is independent of age and sex (shown in Figure 3) it can be considered as a unique method to classify the distinct body size in any population. The accuracy of the classification of body size using this statistical model may be considered better and satisfactory due to the consideration of a number of anthropometric traits in the analysis 8. As a result, the higher order three-factor model was preferred for determining body size. As a general assumption, results focused on mean values of most of the measurement variables maintained increasing higher order from Lean body size to robust body size is the validity of the model.

The high value of coefficient of several anthropometric traits such as, transverse chest (Tvc), chest circumference (Ccn), skinfold triceps (Skt) reflects the high contribution of these respective latent component in the first order model. and the latent component linear was found to have contributed the most in second-order model.

Thus, the derived statistical model suggests that anthropometric traits can be used to determine the CBSS based on the hypothesis of the study.

It may be concluded that the proposed models may be used as a precise estimator to assess the body size of a population and in determining the nutritional status. The advantages of this method lie in terms of assessment of body size are: reliability and feasibility and cost effectiveness.

The model would have been more précised with the increase in sample size. The method is generalized concerning the cut-off points of different categories of body size but may remain inconsistence across populations. An adequate number of community-based studies can give better insights into population variation in body size cut-off values.

The author is indebted to the study participants for their unhesitating involvement and to the authority of the Indian Statistical Institute for providing financial support for this research. The author is also immensely grateful to Prof. Subha Roy and Dr. Indrajit Dasgupta for their valuable comments and throughout editing the manuscript. Thanks to Mr. Akash Mallick for his technical support and revision of the manuscript. Prof. Barun Mukhopadhyay and Prof. Chrales Weitz, are also acknowledged for the crucial grammatical revisions and suggestions.

The dataset used and/or analysed during the current study are available from the author on reasonable request.

The Author solely contributed to the set design, data analysis and interpretation of the results as well as preparing the drafts of the manuscript.