Descriptive statistics were computed for all variables. Normality of the variables was examined using boxplots, Kolmogorov-Smirnov, and Shapiro-Wilks tests. The linearity of nominal and ordinal data were assessed using χ2-based measures. BMI was defined into four categories—(1) normal (18.5–24.9 kg/m2); (2) overweight (25–29.9 kg/m2); and (3) obese (>30 kg/m2)—and was also analyzed as a continuous variable.
The relationship between DR and body stature measurements was examined using a multivariable logistic regression model adjusting for potential confounders. For each body stature measurement (BMI, WHR, WC, head circumference, neck circumference, and skinfold thickness), we constructed two models: Model 1 adjusted for age, sex, income, education, smoking, and general health status. Model 2 adjusted for all variables in Model 1 plus duration of diabetes, SBP, insulin use, cholesterol, HDL-cholesterol, fasting glucose level and HbA1c. All pertinent variables were examined for correlations and multicollinearity, using the Pearson product moment correlation. Because multicollinearity was present among the examined predictor variables, factor analysis regression (FAR) was used to estimate the regression coefficients of the factor-predicted score. We also examined the association between body stature measurements and severity of DR (no DR, NPDR, and PDR) using an ordinal regression model (ORM). The proportional odds assumption of ordinal regression was checked by an approximate-likelihood ratio test. A two-tailed P < 0.05 was considered statistically significant (all analyses with Stata Corp, College Station, TX).