2.1 Study population
The National Center for Health Statistics (NCHS) of the Centers for Disease Control and Prevention conducts a thorough cross-sectional survey every two years, from which the analytic data were taken. This survey assesses the nutritional and health status of Americans through questionnaires and physical examinations [17]. Each participant in the NHANES study signed the informed consent form, approved by the NCHS Research Ethics Review Board.
Our study specifically utilized data from the 2001 to 2004 dataset, as information on ED was exclusively available during this timeframe. Among the 21161 participants, due to the ED questionnaire being administered solely to adult males aged 20 and older, we excluded 10860 female participants, and 5347 participants younger than 20. Furthermore, the participants with incomplete ED assessment(n=838), a former diagnosis of prostate cancer(n=120), missing BMI(n=101), triglyceride(n=2167), fasting glucose, and serum insulin(n=8) data were excluded. As shown in Figure 1, we ultimately included 1720 participants in the final analysis.
2.2 Measurement of IR surrogates
Blood samples for triglycerides, low-density lipoprotein (LDL), insulin, and fasting glucose were taken in the morning after fasting. The TyG index was generated as Ln[fasting triglycerides (mg/dL) × fasting glucose (mg/dL)/2][8]. The TyG-BMI index was generated as Ln[fasting triglycerides (mg/dL) × fasting glucose (mg/dL)/2] × BMI (kg/m2)[10]. The HOMA-IR index was generated as fasting glucose (mmol/L) × fasting insulin (µU/mL)/22.5[18]. IR was characterized by a threshold of HOMA-IR index ≥ 2.2. [19].
2.3 Assessment of ED
ED was evaluated as an outcome variable using a single question adapted from the Massachusetts Male Aging Study [20]. Participants were asked, “How would you describe your ability to achieve and maintain an erection adequate for Satisfying sexual behavior?” Options for responses included “almost always or nearly always,” “usually,” “sometimes,” and “never.” In the present study, ED was defined as participants who responded “never able” or “sometimes able” to keep the erection going.
2.4 Covariates of interest
The list of covariates encompassed age, race, marital status, education level, the family poverty income ratio (PIR), smoking status, drinking status, hypertension, diabetes, physical activity, CVDs, and LDL. Marital status was classified as Living alone or Married/living with a partner. Educational attainment was grouped as less than high school, high school, and more than high school. Family PIR served as a proxy for household income and was divided into three categories: "< 1.3" "≥ 1.3 and < 3.5" and " ≥3.5"[21]. Smoking status was determined by questionnaire (having smoked a minimum of 100 tobacco products throughout their lifetime; Positive replies classified respondents as smokers). Drinking status was assessed using a questionnaire (Consumption of at least 12 units of alcoholic beverages in the past year; Positive replies identified individuals as drinkers). Hypertension, diabetes, the level of physical activity, and coronary heart diseases were also determined through questionnaire responses. CVDs meant participants previously diagnosed with heart attack, angina, or heart failure.
2.5 Statistical analysis
Based on the HOMA-IR index, participants in this study were divided into two groups: the "IR group" and the "Non-IR group." Whereas continuous variables were displayed as means and standard deviations, categorical data were represented as numbers (n) and percentages (%). Logistic regression models were employed to examine the association between ED and different IR surrogates, presenting results as OR with corresponding 95% CI. Following STROBE guidelines, three regression models were established: Model 1 without any covariate adjustments, Model 2 adjusting for race and age, and Model 3 adjusting for race, age, education level, marital status, family PIR, hypertension, diabetes, smoking status, drinking status, physical activity, CVDs, and LDL. On the basis of these analyses, we assessed discrepancies further in ED risk across different TyG-BMI quartiles. Additionally, in Model 3, smooth curve fitting was employed to investigate potential curvilinear associations with ED and the TyG-BMI index. Threshold effects analysis was also conducted to identify associations and inflection points between the ED and TyG-BMI index.
Furthermore, we conducted interaction and stratified analysis dependent on race, age., education level, marital status, family PIR, smoking, drinking, hypertension, diabetes, moderate activity, and CVDs.
Finally, we contrasted the predictive performance of these three IR surrogates for IR and ED using ROC curves and corresponding areas under the curve (AUC).
The statistical analyses were carried out using R studio, release 3.4.3, and EmpowerStats, release 2.0. (http://www.empowerstats.com, X&Y Solutions, Inc.). P-value <0.05 was regarded as significant.