International Journal of Health Geographics Open Access Research Geospatial Analysis of Hiv-related Social Stigma: a Study of Tested Females across Mandals of Andhra Pradesh in India

Background: In Geographical Information Systems issues of scale are of an increasing interest in storing health data and using these in policy support. National and international policies on treating HIV (Human Immunodeficiency Virus) positive women in India are based on case counts at Voluntary Counseling and Testing Centers (VCTCs). In this study, carried out in the Indian state of Andhra Pradesh, these centers are located in subdistricts called mandals, serving for both registration and health facility policies. This study hypothesizes that people may move to a mandal different than their place of residence for being tested for reasons of stigma. Counts of a single mandal therefore may include cases from inside and outside a mandal. HIV counts were analyzed on the presence of outside cases and the most likely explanations for movement. Counts of women being tested on a practitioners' referral (REFs) and those directly walking-in at testing centers (DWs) were compared and with counts of pregnant women. Results: At the mandal level incidence among REFs is on the average higher than among DWs. For both groups incidence is higher in the SouthEastern coastal zones, being an area with a dense highway network and active port business. A pattern on the incidence maps was statistically confirmed by a cluster analysis. A spatial regression analysis to explain the differences in incidence among pregnant women and REFs shows a negative relation with the number of facilities and a positive relation with the number of roads in a mandal. Differences in incidence among pregnant women and DWs are explained by the same variables, and by a negative relation with the number of neighboring mandals. Based on the assumption that pregnant women are tested in their home mandal, this provides a clear indication that women move for testing as well as clues for explanations why. Conclusions: The spatial analysis shows that women in India move towards a different mandal for getting tested on HIV. Given the scale of study and different types of movements involved, it is difficult to say where they move to and what the precise effect is on HIV registration. Better recording the addresses of tested women may help to relate HIV incidence to population present within a mandal. This in turn may lead to a better incidence count and therefore add to more reliable policy making, e.g. for locating or expanding health facilities. Background HIV related stigmas are a driving …


Background
HIV related stigmas are a driving force influencing the behavior and location specific testing results of persons seeking HIV testing [1]. Much has been reported about stigmatized behavior, but little has been investigated on the possible movements of persons in general and women in particular seeking anonymity and thus moving from their residence to other places for getting tested. Misinformation about HIV testing attitudes, and HIV stigmatizing beliefs represent potential barriers to testing [2][3][4][5]. Kaplan et al [3] note that our understanding of the mechanisms by which HIV related stigma perpetuates is limited. To plan improved interventions it is necessary to better understand the behavioral pattern of those getting tested. Various population based studies report major differences from sentinel surveillance based estimates [6][7][8].
Hence, obtaining a good insight into the spread of the HIV incidence requires a reliable registration of those infected. Registration has an effect on the official statistics as well, as for example the Indian government recently reported a change in the official incidence value from 4.5 to 1.5% at the national level, similarly to what happened in Kenya (Appendix I). According to Pandey et al [9], the earlier HIV estimates in India were based on HIV sentinel surveillance (HSS) data. It is assumed that prevalence among attendees of antenatal clinics serves as as a proxy for the prevalence in the general population and prevalence among the patients of sexually transmitted diseases as a proxy for the prevalence among populations with high risk behavior. The absence of HIV surveillance among female sex workers and men having sex with men was a weakness of this system. Those two groups were later included in the estimation but sexually transmitted disease clinics were not discarded. This resulted in double counting. In 2006, improved data became available as the sentinel surveillance among ANC women was expanded covering nearly every district in the country allowing better geographical representation with adequate data for each state. Additionally, community based HIV prevalence measured by the National Family Health Survey-3 [10] provided an opportunity to replace earlier assumptions, validate the HSS data and improve the HIV estimate. Calculations and estimates in Pandey et al [9] reverse the number of total HIV prevalence in India. They quote that the current estimate is a revision based on improved data and methodological changes. The difference between the current estimate and previously published estimates does not represent a true decline at the population level.
Oppong [11] while discussing the data problems in HIV research quotes that sample size, nonrepresentative samples, and geographic and testing bias, tend to make seroprevalence estimates defective if generalized beyond its sample population. A change of scale to the facility level improved representativeness and lead to more promising results. Until recently, methodology was developed and applied to data that were only available at the state and the district level. The analysis presented in our study goes one step further and considers the sub-districts or mandals level.
Disease data as analyzed in this study have a clear spatial component. Registration is done at hospitals and clinics that are located in mandals, the spread of the disease is most likely done by roads and transport networks and spatial components are possibly helpful to provide better estimates. To do so, geographical information system was used in this study, providing opportunities to analyze HIV data and related layers of information in a quantitative way by means of readily available spatial statistical tools. The aim of the study is to quantify the degree to which women move to a different place for HIV-testing and to find explanatory variables. To do so, data sets of women being tested on a practitioners' referral and those directly walking-in at testing centers were compared with data sets of pregnant women. HIV data used was collected in 2006 from the Indian state of Andhra Pradesh, where a well-established registration system exists.

Study area
This study area concerns the state of Andhra Pradesh (Insat, Figure 1). For HIV, Andhra Pradesh is the second worst affected state of India, after the state of Manipur. It is in size the fifth largest state of India with an area of 277,000 km 2 , accounting for 8.4% of India's territory. It is divided into 23 districts and 1103 (2001) sub-districts called mandals. Based on the Census 2001, the total population is 76 million, making it the fifth most populous state, with a population density of 277 people km 2 . Population is mainly rural (approximately 72.7%). Andhra Pradesh is predominantly an agriculture-oriented economy, e.g. it is the largest producer of rice in India. The movement of agricultural products and raw and finished material depends on road transport [12]. Because of its large population, good data accessibility, well established e-governance and its comparatively better health infrastructure it is well suited for this study.

Data Used
Population data from the national census in 2001 have been used. Estimation of the annual projected population for 2006 has been based on the population projection for India [13]. Base data to delineate the different mandals and their boundaries were made available from the National Remote Sensing Center. The Eicher Andhra Pradesh Road map © 2008 was used to generate the roads layer for application in the GIS. This study uses HIV data collected by the NACO (National AIDS Control Organisation) in India on indicators supplied by Voluntary Counseling and Testing Centers (VCTC)s. In 2006, there were 190 VCTCs located within 1103 mandals. These data are the most comprehensive one on the population and thus may provide clues to understand the HIV epidemiology [14]. The role of these centers as a convenient and cost effective tool for monitoring the HIV epidemic is well known. Their high coverage within the state and country is a key in the overall success in combating HIV [4,15,16]. As our study is mandal based, the data do not include mobile VCTCs and primary health care centers since these units are district based and hence can not be assigned to a single mandal. The VCTC data represent unbiased samples from the general population. They distinguish two types of clients: self-referred clients or Direct Walk-ins, (DW) and provider-referred clients or Refferals (REF). DWs voluntarily present themselves at a VCTC, whereas REFs are referred for HIV counseling by health-care providers. The decision to undergo an HIV test is voluntary for REFs [17].
Gynecological units are present in different hospitals and clinics to serve as a Prevention Of Parent to Child Transmission (PPTCT) center. Such units facilitate assistance to pregnant women and take measures to control the transmission of infection to the newborns [18]. 167 units were functional in Andhra Pradesh in 2006, covering approximately 10% of the pregnancies in the country. In 2006 PPTCT data were available from January to August that were used in this study.

Incidence Calculation
This analysis is based on the representativeness of PPTCT data. Therefore, the subset of HIV-positive women from the VCTC data was selected. To make group-wise comparisons, incidence is calculated per 10,000 female population. Table 1, extracted from [10], shows that pregnant women mainly belong to the age class 15-39. HIV-positive women belonging to the same age-group were selected from the DW and REF groups in the VCTC data. Data from these three groups represent the numbers of HIV-positives belonging to age range 15-39 in a particular mandal at a given time.  Methods for this study are based on a spatial pattern analysis, outlier detection and establishment of spatial relationships. An outlier is defined as an unusually high or low HIV frequency as compared to the DW and REF data sets for the same mandal and/or values for the same data set for other mandals in the direct vicinity. The first question considered is which women are represented by the three data sets and what is their behavior and spatial distribution? The second question is what types of movement should be linked to HIV testing and which group is expected to have a particular type of movement behavior? The spatial pattern analysis was done based on the understanding of the type of the three diverse groups and what they represent.
• Pregnant women may be used as a proxy for prevalence in the overall population [19]. They are mostly married and they are equally distributed over the population. There are well-known limitations, however, as not all pregnant women may access the antenatal care services or may accept HIV testing. This apparently is of a limited importance, because in this data set 92% of women attending the antenatal clinics accepted the tests.
• Women registered as REFs show a larger diversity than women registered as DWs. Hence more cases are expected in places with more and better facilities for testing, i.e. with the order of the facilities. REFs, being already asked by a practitioner to get tested, are less likely to move. So they may be more inclined go to the nearest testing center.
• Women registered as DWs are more likely to belong to a high risk category, i.e. being involved in sex trade and injecting drug use. Their spatial pattern may reflect the locations of areas conducive for risk activities like regions rich in trade that are well connected with urban setups, such as roads signifying movement. Because DWs get tested on their own, they may move to any facility of choice. Thus DWs have an opportunity of travelling larger distances and thus have a higher probability of being registered at another place. They might therefore seek testing at anonymous sites and hence they will form the group governing the movement. Any difference in the spatial pattern for the three groups can be attributed to the cause of movement. Based on the scale of analysis, socio-psychologic behaviour of women getting tested and the societal setup, we distinguish accessibility movement and hierarchical movement. Accessibility movement relates facilities that are better connected to urban setups with a higher incidence. Since the choice of connection is important for the DWs, this type of movement should be identifiable by the highest correlation of connectivity with higher incidences of DWs. Hierarchical movement relates the order of the facilities, e.g. from a community health center to medical college, to a higher incidence. In particular for REFs higher order facilities should show a higher incidence, as they can be referred to the best facility, usually equated to the highest order. Three other types of movement that are not detectable are distinguished: random movement, movements at a very short distance (women aware of HIV usually select a VCTC within 60 km away or private clinics suggested by friends (Pers. Comm.)) and movements that neutralize or counterbalance movements between mandals. Based on the above the following assumptions were made for the three groups under study: 1. Women to be tested at a VCTC in a particular mandal as REFs (F REF ) comprise both the women (F RL ) from the same mandal and women from other mandals (F RO ) that aim to maintain their anonymity.
2. Women tested as DWs (F DW ) at a VCTC in a particular mandal comprise both the local walking in women (F DL ) and women walking in from other mandals (F DO ).
3. Pregnant women (F P ) getting tested at a PPTCT center in a particular mandal are those belonging to this mandal only. Their main incentive is to receive antenatal care and the HIV test is additional to that. 4. The proportion of local DWs is assumed to be less than the proportion of local REFs in each mandal, since the local DWs have an opportunity to move to other places: It is thus assumed that the local DWs (usually represented by the sex workers) would generally move and hence would be less represented as compared to the REFs who will remain at their place of stay. Movement was analyzed first on the basis of an outlier detection scheme, showing mandals which deviate from the normal behavior. Secondly, a spatial cluster analysis was applied to detect geographic variation patterns and identifying locations having statistically significant higher incidences as compared to their neighbors [14]. Finally, spatial regression was carried out to quantify the observed patterns.

Outlier detection
Mandals at both ends are outliers: lower-end outliers represent mandals with I P comparatively higher or nearly equal to I DW and I REF , whereas higher-end outliers represent mandals where I P is lower than I DW and I REF . Incidence maps were generated within the ArcGIS [20] environment using the I P , I DW and I REF incidences. These maps were used in turn to first yield two difference maps relating the REF and DW groups to the P group These maps were classified into six classes. Lower ranges in the ID 1 and ID 2 correspond to higher I P values. Therefore, smaller class intervals were chosen in the lower range keeping negative values as one class, and then having class ranges of 0, 2, 5 and 10, respectively. These maps could thus identify mandals with strong differences between pregnant women and women from the general population. Values equal to 0 in the ID 3 map identify mandals with equal I DW and I REF . Mandals with a high ID 3 value are outliers, representing an exceptionally high I REF .

Spatial Cluster analysis
Spatial cluster analysis is commonly used in disease surveillance and spatial epidemiology [21]. For this study, SaTscan™ [22] software was used to compare spatial clustering in the data with a Poisson model showing randomness. In total, 1103 -190 = 913 missing values represent mandals without facilities. To account for these, the missing value adjustment parameter was used, assigning a relative risk of zero to mandals without data. In doing so, the analysis ignores those mandals. The results of the spatial clustering in SaTscan™ were imported into the ArcGIS environment where significant clusters were visualized using p-values.

Establishing spatial relationships
On the basis of the above analysis, it can be predicted if people are moving and a trend can be estimated. A hypothesis is established for the following relations to be possibly significant: 1. The effect of facility hierarchy plays a role in a higher I REF value in a mandal. Thus a higher order facilities will draw more REFs and more pregnant women than a lower order facility.
2. Vicinity of roads may increase the I DW because of better connectivity. Therefore, the distance to a major road may have a positive relation to I DW . Also the number of road intersections within a mandal repre- sents a better connectivity that may increase the movement of DWs.
3. Incidence in pregnant women most likely remains unaffected by connectivity, given their status of pregnancy, whereas it may be related with the number of neighboring mandals. If a mandal has more neighboring mandals then it may be attractive to visit, realizing that many mandals do not have their own testing facility. Therefore, effects of the following exploratory variables are investigated: 1 where the coefficients α i are to be estimated. Initially, to decide upon model composition, contribution of each variable was explored by using ordinary linear regression (OLS). Below, after model identification by identifying possible explanatory variables, an autoregressive approach is used to include spatial dependency in making a final estimate of the parameters. A spatial autoregressive modelling (SAR) is done for those variables that show a significant relation. A SAR model consists of a spatially lagged version of the incidence I x as where the matrix W represents neighbour relations, i.e. w ij = 1 if mandals i and j are neighbors, i.e. dist(i, j) < 50 km and w ij = 0 otherwise. The value of 50 km is used as a balance between a sparse neighbourhood pattern and a full inclusion of all the neighboring mandals. Other values have been tested as well but did not show strong differences. The parameter ρ is the autoregressive parameter establishing autocorrelation and the denotes independent noise. Model (8) is equivalent to (7), except for the neighbourhood structure and the autoregressive parameter ρ. The spatial weights matrix W is standardized such that its rows sum to 1 [23]. The 164 mandals having a PPTCT center were selected for the analysis, the other mandals were discarded. The distance of 50 km for neighbourhood definition resulted into 395 neighboring mandals.
All layers have been created in an ArcGIS environment. OLS has been done in SPSS [24], with one variable at a time and I REF , I DW and I P as the response variable, whereas the SAR analysis has been done using the spdep library in the R package [25]. ID 1 , ID 2 and ID 3 maps represent the mandals that explain movement of HIV-positives ( Figure 3). Assuming that incidence for pregnant women I P is generally lower than incidence for the general population, as they are a subsection of the whole female population, it is noted that HIV-positive females are apparently moving from mandals with negative values and values up to 2 to other mandals for getting tested. Such an approximating approach provides a clue in understanding the differences in the incidence in these mandals. In

Spatial Cluster analysis
Cluster analysis is performed to draw regions in the three classes which represent high rates of incidence. Figure 4 shows the results of the cluster maps for I REF , I DW and I P . Such clusters identify the mandals at higher risks as compared to their neighbors, including their statistical significance. The search radius for the moving window was kept at 5% of the population. Cluster analysis for REFs resulted into 14 clusters of which 9 were significant, for DWs into 11 clusters of which 6 were significant and for pregnant women into 14 clusters of which 11 were significant. The DWs are significantly clustered only at the SE coastal zone, a pattern which can also be witnessed in the incidence maps ( Figure 1). As expected, both REFs and pregnant women are spread more equally over the state though in varying proportions.

Establishing spatial relationships
Relations of the spatial pattern of spread were modeled with underlying factors. Layers of the explanatory variables are shown in Figures 5 and 6. First the relationship between the spatial pattern was explored, outliers and the explanatory variables by means of visualisation. The hypothesis is that bigger facilities would attract more HIV-positives, but it is seen that usually the locations with higher incidences have a smaller facilities, such as a CHC. Similar overlays were prepared for the cluster maps with other layers like distance from roads, number of facilities, road density and number of neighbors. The overlay analysis of the incidence maps of the three cate-gories with the cluster maps and the difference maps was done.
Relations between REFs and pregnant women with the type of facilities and the DWs with the roads were explored. it is observed that a relation between the type of facilities and the pregnant women as CHCs usually have higher incidences, although a significant relation between REFs or DWs with higher order facilities was not discovered. The number of neighbors (N N ) seems to affect incidence on the basis of the visual comparison. The distance from roads (D R ) shows a relation to incidences displayed by the difference maps, although, these patterns are far from uniform.
To have statistical evidence, regression analysis was performed and the results are shown in table 2. Relatively low R 2 values ranging from (0 to 0.05), (0 to 0.07) and (0.01 to 0.3) for I REF , I DW and I P as the response variables respectively are observed. The highest R 2 value equal to 0.307 was observed for the relationship of the type of facilities with I P . The corresponding equation equals This means that the incidence increases with 2.306 if the type of facility increases with one unit. All other variables do not significantly contribute to the incidence of any of the three categories. Using the OLS results, the SAR analysis was performed with IP as the response variable and T F . The following linear relation was found: and an estimated ρ parameter equal to 0.0359 (significant at the α = 0.05 level), hence with slightly different coefficients. Use of the conditional autoregressive (CAR) model model did not lead to any substantial change.
Finally relationships were established between ID 1 and ID 2 as a variable measured at the mandal level to the explanatory variables mentioned in section methods, applying a SAR analysis for quantification. The following model was found to be the best describing the variation in ID 1 : where the autoregressive parameter was estimated as 0.073 (p < 0.001) and an AIC value of 646.5. In this equation the contributions of N R is almost significant (p = 0.0897) whereas that of N F is not significant (p > 0.1). It shows that incidence in REFs is larger than in Ps, and, although somewhat weakly, that this difference could be explained by road density, with a higher difference with an increasing road density. It was somewhat surprising, as initially the hypothesis that the most important explanatory variable would have been N F was not confirmed. Its consequences are also relevant for HIV treatment and follow-up. The next model to be the best describing the variation in ID 2 : ID N N F R 1 9 9 0 923 1 719  where the autoregressive parameter ρ was estimated as 0.069 (p < 0.001) and an AIC value equal to 645.5 was obtained. In this equation the contributions of N R is almost significant (p = 0.062) whereas the other contributions are not significant (p > 0.1). It shows a positive relation between road density N R and differences in I DW and I P , as such supporting the initial hypothesis: the difference increases with increasing road density. This increase is larger for REFs than for DWs, in other words: REFs are more inclined to move to another mandal for being tested than DWs.
None of the variables unambiguously explains the behaviour of the type of tested females. Therefore, although it seems that females might be moving one cannot exactly capture the movement and the attributed reasons do not fully explain any of the hypothesized phenomena. The regions where the incidence in pregnant women is higher than the general population can be identified as the zones of movement and similarly those with high DWs; however no significance or a consistent cause could be attributed to this.

Discussion
This empirical study presents a first step to capture the overall pattern of HIV incidence at the state level to address the movement of people for testing on HIV. Its consequences can be relevant for HIV treatment and follow-up.
Trend analysis by means of maps and graphs revealed that incidence in the referrals group, I REF , shows on the average higher values than incidence in the directly walk-ing-in group, I DW . A possible explanation is that in India there is little movement among women. If women do not belong to the high risk groups, then infection occurs through their partners in marriage and they get tested as a REF instead of as a DW. This spatial pattern analysis also shows that I P is lower than I REF and I DW . The most likely explanation is that the number of HIV-positives from PPTCT centers represents only a fraction of the total female population. Several mandals, however, have larger I P than I REF and I DW values. With an underlying assumption that I P should be the lowest, the mandals defying the trend give us a reason to further explore potential causes. A hypothesis that this occurs at random should be tested against the alternative that a definite and clear cause exists, such as the quality of the unit and reported success stories. The current data set did not allow us to do so, however.
The higher rates of I REF and I DW in the South Eastern coastal zones are clearly shown, both by the spatial pattern analysis and by the cluster analysis. This area is marked with a dense highway network and active port business. According to [26], this is also a favourite destination for the female sex workers, most likely explaining the registered incidence in these areas. A clear distinction exists between mandals where people live, and mandals where their HIV status is recorded. Elevated clusters are found for DWs in this region whereas the pattern of REFs is more scattered. The high variation of I P in terms of spatial spread is caused by the fact that pregnant women are a control group which is supposed to reliably represent the underlying population. Also, a high incidence rate is general population is either not getting tested or moves to another place. In particular, the South Eastern coast zone is attractive, being a well connected urban set-up. Other reasons for comparatively lower I REF and I DW values in the rest of Andhra Pradesh might be caused by the low testing rates and lack of adequate and easily accessible facilities. The attempts to relate I REF , I DW and I P with different parameters reveal a few interesting correlations. I P shows a positive significant relation with the type of facilities. This is in accordance with the social behaviour where women using government facilities usually prefer higher order facilities for anti-natal care. Also, based on visual analysis, it is noticed that community health centers have often been associated with higher incidences of REFs and DWs. This means that it is not the hierarchy of facility based on size that plays a role but it is the presence of a facility. Therefore women are likely to get tested if a facility is present, either small or large, and if they are aware of it. Since no significant relationship was observed with the road infrastructure and the proximity, one can infer that it is not governed by the good connectivity whether women move for getting tested. This may also explain the assumption that capturing movement depends on the type of movement and the transportation modes available in a mandal.
The following recommendations are derived from this study: • HIV is a dynamic disease and a good data capturing is the backbone of all the policies. Further analysis in a spatio-temporal domain may be the key to better understand the interplay of various factors.
• The fact that one can only partially, i.e. non-significantly, explain the relations of differences leads us to assume that at the scale of the study and the available data, much of the movement is random and that a more detailed data set should be collected to exactly identify where people are moving and what factors are governing them in their behaviour.
• From a policy point of view, it may be important to increase self-motivation among women specially belonging to the HRGs (High Risk groups) potentially represented by the DWs to get tested because of the rapid progressing of HIV. More focused and better policies are needed to enlighten women so that they do not wait for a reference but visit a VCTC to get tested. In particular Andhra Pradesh needs special attention to let women abstain from behaviour responsible for the spread, and to take special measures not to allow the disease to spread to other states.
• A better insight into the quality of the data may help to improve describing factors determining HIV spread and to support spatial decision making, like positioning new health care facilities.
• Common policy assumption of coincidence of residence and test place is challenged by the present study, and should be reconsidered in future policies. The study was constrained due to some important factors. Different sources of data sets were used; hence interoperability is a major problem. Census data, administrative boundaries, NACO and road data all have different sources and different procurement time which have to be adjusted for each other. This loses the originality of data to an extent and hence affects the results. For this exploratory study the amount of available data was large, but still more could have been measured. Possibly, the use of additional information could lead to a better analysis with a higher amount of explained variation. The available data set, representative at the level of mandals, however, was already quite unique and as far as we know has not been analyzed before.
The aim of this study was to analyze the whole state of Andhra Pradesh, but since facilities are present only in a limited number of mandals the analysis addresses some 20% of the mandals. This is compensated by the fact that an analysis at the district level integrates data from many hospitals. The main point addressed in this study about HIV policy-making, however, is that a change is needed in a basic assumption that place of testing and residence coincide. Consequences of such divergence need to be further explored in future research. Data quality could further improve if a better registration is done. Women should deliver their home address when visiting a VCTC for being tested. Also, motivating information about their preference of choice should be provided.

Conclusions
Some concrete conclusions follow from this study. First, it was hypothesized that higher order facilities would attract more HIV-positives, but the study shows that mandals with higher incidences usually have a lower order facilities, such as a community health center. Therefore a hypothesis for further research could be that anonymity attracts females to a lower order facility for testing. Second, a pattern is observed between the type of facilities and the pregnant women as community health centers usually had higher incidences. However, significant relation between REFs or DWs with higher order facilities could not be discovered. Finally, there is a significant relation between the incidence in pregnant women and the order of the testing center.
Several trends emerge from the present study. The outlier analysis and the cluster analysis show that women move for getting tested. The present dataset did not allow us to say where they move to and what the precise effect is on HIV registration. The assumption that there is a random movement is not traceable at the given scale, also because of the amount of missing data. Alternatively, movement is perhaps an interplay of other interacting socio-economic factors which need to be further addressed.
Further research involving more spatio-temporal data would be helpful. This study relies on the 2006 data, since only those had the detailed PPTCT information. The number of testing centers is increasing with time, and data from 2007, 2008 and 2009, with less missing values, might be used. Comparing different years may provide us with more conclusive inter-relationships. A next step may be to analyze I REF and I DW differences for males, either relating these incidence on the basis of assumptions to those of female incidences, or by using a different benchmarking. It would be interesting to explore the relationships of the male incidence with different variables. This together with the female analysis will give us a larger picture and better understanding of reasons for people to move and in the end more reliable HIV data of a better quality.