Data collected were checked manually for completeness and consistency, and then cleaned, categorized, coded, and entered into EPI info version 7. For the sake of analysis, data were transferred to SPSS version 24. During the analysis, frequencies and OR with 95% CI were calculated to determine the association of selected variables with the outcome variable. Descriptive, bivariate, and multivariate logistic regression analyses were implemented to determine the association between predictors and outcome variables and to control cofounders. All variables associated with risky sexual behaviors in bivariate analysis were considered as candidate variables with p < 0.25 for the multivariable logistic regression model. A p-value of less than 0.05% and 95% CI that excludes 1 (unity) were considered to be statistically significant. Model fitness was tested using the Hosmer–Lemeshow test at p > 0.05 as a cutoff point. Qualitative data were recorded, transcribed, and translated to English, and finally used to supplement the quantitative findings and analyzed thematically.