Identification of Socio-Demographic, Behavioral Patterns and Their Relationship with HIV Status among Female Sex Workers ()
1. Introduction
India, the second-most populous country in the world, has been experiencing a highly varied HIV epidemic, in the past two decades. HIV in India is concentrated mainly among “high-risk” populations such as Female Sex Workers (FSWs), Men having Sex with Men (MSM) and Injecting Drug Users (IDUs). Targeted interventions (TIs) among these populations have been a major aspect of India’s National AIDS Control Programme (NACP). The recent HIV estimates showed a reduction in the adult HIV prevalence from 0.41% in the year 2000, through 0.36% in 2006 to 0.31% in 2009. The estimated number of new annual HIV infections declined from 2.7 lakh in 2000 to 1.2 lakh in 2009 [1] .
New infections reduced by 50% during this period. Therefore sustained focus on prevention is envisaged by NACO.
Andhra Pradesh (AP) has been categorised as “A” with HIV prevalence among FSW > 1% and the state remains to be critical in the spread & control of the epidemic. The state has the second highest estimated adult HIV prevalence of 0.90% after Manipur (at 1.40%). While India has an estimated 24 lakh people living with HIV/AIDS in 2009, AP tops the list of states with about 5 lakh people living with HIV/AIDS; followed by 4.2 lakh in Maharashtra, 2.5 lakh in Karnataka and 1.5 lakh in Tamil Nadu. Of the estimated people living with HIV/AIDS, 9.3 lakh (39%) are women [1] .
HIV prevalence among Female Sex Workers (FSW) in Andhra Pradesh increased from 7.3% in 2006 to 9.7% in 2007. Sentinel surveillance data of 2007 also depicted Andhra Pradesh as the only state with >1% prevalence among ANC attendees [2] . (HIV Sentinel Surveillance and HIV Estimation in India 2007. A Technical Brief- NACO.)
Targeted interventions by the Andhra Pradesh AIDS Control Society (APSACS) and Avahan by Bill and Melinda Gates Foundation (BMGF) are the two major HIV/AIDS prevention programs currently in operation in the state [3] . The state government (in partnership with NACO) has been implementing prevention, care and treatment programs since 1999, which have resulted in significant gains. However, the epidemic continues to be of serious concern, as recent surveillance rounds reveal that, across Andhra Pradesh, prevalence levels continue to display an upward trend, pointing to the need for intensifying preventive intervention programmes.
In this backdrop, the present paper utilizes the “Integrated Behavioral and Biological Assessment” (IBBA) data to analyze and identify the demographic, behavioral and biological patterns among FSWs that contribute or relate to the HIV prevalence in the State of Andhra Pradesh [4] .
Objective
The objective of the present communication is to identify behavioral patterns of FSWs and study their relationship with HIV status.
2. Methodology
2.1. Sampling & Sample Size
Considering consistent use of condom by 50% FSWs with the clients during round 1 with 15% increase during round 2, with 95% CI, 90% power and 1.7 as design effect the required sample size is 385, rounded-off to 400, which adds up to 3200 FSWs for 8 districts [4] . Sampling frame consisting of different sites of solicitation was developed for each district, a fortnight before commencing the actual survey. Probability sampling method with varied sampling approaches was used to select the site as Primary Sampling Units (PSUs) at the first stage. Conventional Cluster Sampling (CCS) as well as Time-Location Cluster Sampling (TLCS) approaches were adopted, based on the characteristics of the different FSW sub-groups in various districts [5] . FSWs who operated from home were grouped as “home based”, while those who operated from other places/fixed locations such as brothel, dabha, lodge, or rented room were grouped as “brothel based”. The sex workers who were mobile and operated in public places such as street, bus stand, railway station and others were grouped as public place-based sex workers. The CCS approach was used for selecting brothel and home-based FSWs, while TLCS method was used to select public place-based sex workers [5] .
2.2. Coverage
The districts of Hyderabad, Karimnagar and Warangal from Telangana Region, Visakhapatnam, East Godavari, Guntur and Prakasam from Costal Region, Chittoor District from Rayalseema Region were covered in the survey, which were under the Avahan interventions, having a high number of Key populations (KPs), diversity among different socio-cultural strata in the region and with a goal of modeling Avahan’s impact in different types of settings. A total of 3083 FSWs were covered from 8 districts in the State.
For the study, FSWs were operationally defined as “any female, ≥18 years, either brothel based (working/ living/operating in red light/brothel areas), home based or public place based (soliciting male clients on the street or in other non-brothel settings), who sold sex in exchange for cash at least once in the last one month”.
The potential risk variables considered in the present analysis included demographic variables such as 1) age 2) literacy, 3) marital status 4) age of initiation of sex work and duration in sex work, the behavioral variables such as 1) having regular clients, 2) having regular unpaid partners, 3) having other regular unpaid partners, 4) having occasional clients, 5) wanted to use condom but could not use due to several reasons, 6) dependency on sex work, 7) had sex during travel, 8) client volume per week, 9) carrying condom at the time of interview, 10) exposed to any interventions, 11) correct knowledge of HIV/AIDS, 12) obtained condom from peer educator and the biological parameters such as 1) history of having any STI during past 12 months and suffering from any STI (syphilis or NG or CT) at the time of interview.
The study was approved by the Health Ministry Screening Committee, Government of India, Family Health International’s Protection of Human Subjects Committee and Andhra Pradesh State AIDS Control Society, the Government of Andhra Pradesh. In addition, the Scientific Advisory Committee of National AIDS Research Institute (NARI) approved the IBBA protocol, and National Institute of Nutrition (NIN) Ethics Committees reviewed and approved the consent forms, study instruments, Standard Operating Procedures (SOPs) and field manuals.
2.3. Statistical Analysis
The percentage distributions of different demographic, behavioral and biological variables were calculated. Behavioral patterns were obtained by exploratory factor analysis for 20 potential risk variables considering simultaneously [6] . The principal component analysis (PCA) technique, followed by factor analysis was used to transform a large set of correlated variables into smaller sets of non-correlated variables, called principal components or factors [7] . The factor analysis of the data objectively points out characteristics those are clustered. The aim of this technique is to identify the underlying structure in data matrix, by summarizing and consigning data to arrive at a systematic measurement of the behavior. To summarize data, factor analysis requires dimension that, when interpreted and understood, describes it in terms of a much smaller number of significant variables than a host of variables.
For the applicability of factor analysis, the uniformity of sample was tested by examining the distribution of variables in a loading plot, contrasting the value observed against those expected in a normal distribution which was verified by Kaiser-Meyer-Olkin (KMO) measurement of adequacy [8] . A KMO value of more than 0.50 was considered acceptable. The presence of correlations between demographic, behavioral and biological variables was tested using the Bartlett test of sphericity (homogeneity of variance). The Bartlett test statistic is approximately distributed with chi-square and was accepted at a minimum significance level of p < 0.05. Principle Component Analysis was used for extraction of factors and orthogonal rotation (varimax option) to derive non-correlated factors, to minimize the number of indicators that have high loading on one factor. The first factor extracted is the one that accounts for the maximum possible variance in the dataset. The second component, independent of the first, will be the one that explains the largest possible share of the remaining variance and so on, without the components being correlated with each other. While deriving behavioral patterns for FSWs dataset, the variable, viz., current HIV status is not included in the analysis as it is intended to study the relationship between behavioral patterns vis-à-vis HIV status.
Kaiser criterion, namely Eigen value of >1.0, is widely-used for choosing the number of factors in the factor analysis. It was also based on the Eigen plot (scree plot), which shows the total variance associated with each other.
Cluster analysis was performed to identify relatively homogeneous groups of cases based on the factor scores derived through factor analysis, using an algorithm that can handle large number of cases [9] . However, number of clusters would depend on the number of cases fall in different clusters. The prevalence of HIV was considered to name and label the cluster. Correctness of labeling of cluster was verified by testing the significant difference between the clusters on the back ground characteristics of the cases.
Discriminant function analysis was used for studying the relationship between the demographic, socioeconomic and behavioral variables with HIV status. The factor scores obtained by the factor analysis for these variables were considered as continuous independent variables, and HIV status of the female sex workers was considered as dichotomous dependent variable. Sensitivity and Specificity of the discriminant function was also calculated. The statistical analysis was performed using SPSS [10] .
3. Results
The distribution of FSWs according to demographic and behavioral characteristics is presented in Table 1. Majority (42%) of FSWs solicit their clients in public places. About 19% FSWs were in the age group of <25 years,
Table 1. Frequency (%) distribution of FSWs according to their demographic, Behavioral and biological profile (n = 3083).
about 47% were in the age group of 25 - 34, while 35% were in the age group of 35+ years. Majority (50%) of FSWs started commercial sex work before 25 years about 50% reported a duration of >5 years of sex work, While half of the FSWs reported that their volume of clients was >10 per week.
About 42% of FSWs were literates and 56% were married and staying with their family. Majority of FSWs had regular unpaid partners (72.8%) and regular clients (83.3%). Majority of the FSWs (81%) reported that they were using condom consistently with occasional clients, about 17.2% reported that they wanted to use condom but could not use due to various reasons and very few (7.7%) reported had sex during last travel. About 47% of FSWs were carrying condom at the time of survey, while 61% reported that they obtained the condom from peer educator. About 33% of FSWs reported the history of having any STI (NG, CT, and Syphilis), while 10.1% were currently suffering from STIs. Only 1% of the FSWs had correct knowledge on HIV. While majority of FSWs (73.5%) reportedly undergone HIV testing, 11.4% of them were HIV positive.
Analysis reveals that nine components for high risk population were extracted by factor analysis using Principal Component Analysis with varimax rotation for demographic, socioeconomic and behavioral variables (Table 2). These nine components (factors) in the initial solution have an Eigen value over 1, and they account for about 62% of the observed variation in the behavioral pattern among female sex workers. The parallel line to horizontal at Eigen value equaling to 1 in scree plot showed that nine factors will be extracted for this population (Figure 1). The first factor, which accounted for 10.9% of the total variance among the female sex workers was labeled as having unpaid regular partners. High factor loading observed for those FSWs who are either married or having a regular unpaid male partner characterized these factors. The second factor explained 9.2% of the total variance and was labeled as exposure to intervention (exposed to any type of intervention and obtained condom from peer educator employed by the NGOs working in that area). The third factor which accounted for 7.8% of the total variance was characterized by the age of the FSW and duration in sex work was labeled as age & duration of sex work. The fourth factor explained 6.8% of the total variance was done labeled as correct knowledge of HIV/AIDS and getting HIV testing done voluntarily. The fifth factor accounted for 6.0% of the total variance and characterized having any STIs for the past 12 months and wanted to use condom but could not use and it was labeled as suffering from STIs. The sixth factor explained 5.4% of the total variance and was labeled as dependency (majorly depend up on sex work and the place of solicitation). The seventh factor accounted for 5.2% of the total variance and was characterized by single variable namely having regular clients who pay every
Table 2. Rotated component matrix for socio-demographic, behavioral and biological parameters.
Extraction method: principal component analysis. Rotation method: varimax with kaiser normalization. aRotation converged in 16 iterations.
Figure 1. Socio-demographic, behavioral and biological characteristics of FSWs.
time which was labeled as regular clients. The eighth factor also explained 5.2% of the total variance and was labeled as having occasional clients (having occasional clients and the volume of these occasional clients per week). The ninth and the last factor explained 5.0% of the total variance and was labeled as positive for STIs (at the time of survey positive).
In the present analysis, two clusters were formed with sizeable numbers in each of the cluster for further analysis. Cluster 1 was formed with 2668 subjects and cluster 2 with 415 subjects. The distributions of the FSWs according to their back ground characteristics in these two clusters were tabulated and look for significant difference in the proportions are presented in Table 3. The data revealed significant (p < 0.01) differences in all the variables studied except dependency on sex work, having other regular unpaid partners, ever taken HIV test and having correct knowledge regarding HIV/AIDS between clusters (formed on the basis of behavioral patterns). The FSWs belongs to cluster 2 had higher age, about 35% had 10+ years of sex work, 24% had started at early age (<20 years). The literacy status of FSWs of cluster 2 was higher (48% vs 42%) had sex during travel (1% vs 57%), having regular clients (82% vs 94%) and consistent condom use with occasional clients (83% vs 70%) compared to cluster 1. Twenty eight percent of FSWs of cluster 2 reported that they want to use condom but could not use due to various reasons during last sex act, 54% were carrying condom at the time of interview and 51% of respondents were obtained condom from peer educator compared to 62% which is lower than that of FSWs belongs to cluster 1. Exposed to intervention by the FSWs belong to cluster 2 was less (68% against 77%) of cluster 1. Significantly high proportion of FSWs of cluster 2 against of cluster 1 had the history of STIs (43% vs 31%) as well as currently suffering 28% vs 7%). The overall prevalence of HIV was 11.4% and the prevalence was significantly high (15.9%) among FSWs belongs to cluster 2 compared to cluster 1 (10.6%). These findings indicate that the FSWs of cluster-2 had different socio-demographic, behavioral characteristics which are likely be the risk factors for STIs and HIV positivity and can be named as relatively bad cluster.
The factor scores obtained for these populations through factor analysis were used for the prediction and classification of subjects having HIV based on behavioral patterns. The present data showed a considerable relationship between the socio demographic, behavioral and biological patterns and the HIV status. About 68% of FSWs were correctly classified as they were positive or negative for HIV test with sensitivity of 51.4% and specificity of 70% (Table 4). The behavioral patterns that are extracted from the data obtained in the study may not be applicable to other populations with altogether different behavioral characteristics.
Table 3. Distribution (%) of FSWs by their background characteristics by clusters 1 & 2.
c2. **Sig p < 0.01.
Table 4. Classification of subjects with HIV and without HIV using observed HIV and predicted HIV using scores derived from factor analysis for demographic & behavioral characteristics.
4. Discussion
In India, Andhra Pradesh is one of the high prevalence States with an estimated 16 per cent prevalence among commercial FSWs, ranging from 8 - 41 per cent. The present study reported 11.4% prevalence of HIV varying from district to district. Half of the FSWs in the present study had their first sexual debut at a very young age (<25 years). Similar observations were made in a survey carried out in 2006 [11] .
The extent of condom use at with occasional clients reported in the present study was similar to an earlier report [11] . With HIV prevalence being high among FSWs, it is important to promote condom use in all types of sexual relationships, including regular, occasional partners and non commercial partners, since even non commercial partners might have risky sexual behaviours such as unprotected sex with multiple partners.
Analysis of behavioral patterns may offer benefit of summarizing behavior using a small number of factors. Yet, these may be potentially of greater relevance to educate the high risk population on more healthy behavior. Both PCA and cluster analysis are useful approaches for the assessment of behavioral patterns. A commonly-cited criticism of the two techniques is that these involve several subjective―but important―decisions, such as group- ing of demographic and behaviors and possible transformations of variables. PCA involves decisions about the number of components to be retained and their subsequent labeling. Similarly, cluster analysis requires choices about the method of clustering and labeling of clusters. Another disadvantage of these techniques is that they generate patterns based on variation in the variables, with no guarantee that these patterns will be predictive of a particular health outcome. However, the techniques have the advantage that they are empirically derived and are, therefore, not limited by mere a priori knowledge.
In developing countries, studies on identifying behavioral patterns and their relationship with HIV status are scarce. Nine patterns among the variables considered explaining 62% of the total variability.
The first factor, labeled as having unpaid regular partners who accounted for 10.9% of the total variance followed by exposure to intervention 9.2%, age & duration of sex work 7.8%, knowledge of HIV/AIDS and getting HIV testing done voluntarily 6.0% and about 5.0% each by the other factors.
The Socio demographic, behavioral and the biological parameters are highly correlated, and the classification methods based on univariate analysis may lead to flawed estimates [12] . Therefore, multivariate methods, such as factor analysis, represent an alternative approach to the evaluation of individual characteristics for the identification of patterns allows us to examine the effect of behavioral factors as whole and helps describe the association with HIV status. Moreover, it should be kept in minds that individuals behave significantly different, which are influenced by many factors, such as, socioeconomic, demographic and behavioral characteristics. Describing socio demographic, behavioral and biological patterns may be useful in developing community programmes for these high risk populations. Rather than changing the behavior, such programmes can be aimed at changing the behavior that is readily recognized by the target group. It is possible to find a smaller number of measures using factor analysis that would permit the identification of risk behaviors. The present data showed a considerable relationship between the socio demographic, behavioral and biological patterns and the HIV status. About 68% of FSWs were correctly classified as they were positive or negative for HIV test with sensitivity of 51.4% and specificity of 70%. The behavioral patterns that are extracted from the data obtained in the study may not be applicable to other populations with altogether different behavioral characteristics.
5. Conclusion
The data allowed identification of patterns defined by factor analysis based on the data collected from FSWs. There exists a strong relationship between specific behavior patterns and HIV status among the high risk populations. These results will be useful for identifying people in the community with risk behaviors, related to HIV/AIDS, which will help the planners to develop suitable interventions for control and prevention of HIV/AIDS.
Acknowledgements
The authors thank the Director, National Institute of Nutrition for his support throughout the study and acknowledge the hard work and commitment of the IBBA Team at, Field Division and Microbiology Division. Financial support was provided by the Bill and Melinda Gates Foundation through Avahan, its India AIDS initiative.
NOTES
*Corresponding author.