Skip to main content

Bayesian spatial modelling of intimate partner violence and associated factors among adult women and men: evidence from 2019/2020 Rwanda Demographic and Health Survey

Abstract

Background

Intimate partner violence (IPV) remains a global public health concern for both men and women. Spatial mapping and clustering analysis can reveal subtle patterns in IPV occurrences but are yet to be explored in Rwanda, especially at a lower small-area scale. This study seeks to examine the spatial distribution, patterns, and associated factors of IPV among men and women in Rwanda.

Methods

This was a secondary data analysis of the 2019/2020 Rwanda Demographic and Health Survey (RDHS) individual-level data set for 1947 women aged 15–49 years and 1371 men aged 15–59 years. A spatially structured additive logistic regression model was used to assess risk factors for IPV while adjusting for spatial effects. The district-level spatial model was adjusted for fixed covariate effects and was implemented using a fully Bayesian inference within the generalized additive mixed effects framework.

Results

IPV prevalence amongst women was 45.9% (95% Confidence interval (CI): 43.4–48.5%) while that for men was 18.4% (95% CI: 16.2–20.9%). Using a bivariate choropleth, IPV perpetrated against women was higher in the North-Western districts of Rwanda whereas for men it was shown to be more prevalent in the Southern districts. A few districts presented high IPV for both men and women. The spatial structured additive logistic model revealed higher odds for IPV against women mainly in the North-western districts and the spatial effects were dominated by spatially structured effects contributing 64%. Higher odds of IPV were observed for men in the Southern districts of Rwanda and spatial effects were dominated by district heterogeneity accounting for 62%. There were no statistically significant district clusters for IPV in both men or women. Women with partners who consume alcohol, and with controlling partners were at significantly higher odds of IPV while those in rich households and making financial decisions together with partners were at lower odds of experiencing IPV.

Conclusion

Campaigns against IPV should be strengthened, especially in the North-Western and Southern parts of Rwanda. In addition, the promotion of girl-child education and empowerment of women can potentially reduce IPV against women and girls. Furthermore, couples should be trained on making financial decisions together. In conclusion, the implementation of policies and interventions that discourage alcohol consumption and control behaviour, especially among men, should be rolled out.

Peer Review reports

Background

Intimate partner violence (IPV) refers to all forms of behaviour within an intimate relationship that causes physical, sexual or psychological harm, including acts of physical aggression, sexual coercion, psychological abuse, and controlling behaviours [1, 2]. It has been recognized globally as an important public health problem and a huge violation of human rights [3] which has been associated with significant morbidity and mortality as well as health and socioeconomic effects among both women and men [4, 5].

Intimate partner violence is an important public health problem and is associated with several reproductive health issues including increased risk of human immunodeficiency virus (HIV) and sexually transmitted infections (STIs), mental health, unintended pregnancy and pregnancy loss [6,7,8]. Intimate partner violence against women has an estimated global prevalence of 30% [3]. The prevalence of IPV however varies across societies, socioeconomic groups, races and sexes depending on various factors [9, 10]. Previously described risk factors associated with IPV against women include a history of exposure to violence (parental history of spousal violence), marital discord or dissatisfaction, difficulties in communicating between partners and male controlling behaviours between partners, lower levels of education, and a history of exposure to child maltreatment [3, 11]. A prevalence study conducted by the WHO in ten countries with varying income levels uncovered a substantial disparity in the occurrence of IPV. The research revealed that approximately 15%–71% of women had experienced IPV at least once in their lifetime [6]. The same study showed prevalence was especially high in low-income countries. On the other hand, IPV against men is not uncommon but is usually neglected, hence, limited research on IPV against men [12,13,14,15,16,17,18].

The prevalence of IPV is particularly high in Sub-Saharan African countries, especially where there are significant disparities between men and women in terms of income, land ownership, and legal rights. Most countries in the Sub-Saharan African region have assimilated IPV practice into their religious and socio-cultural traditions [19,20,21,22]. Interventions that could aid in the reduction of IPV among women include community initiatives such as the promotion of sports for women [23] and equitable land-owning policies for both men and women [24]. These would assist in achieving the fifth United Nations Sustainable Development Goal (SDG) which aims to achieve gender equality and empower all women and girls (UN-women) [25, 26]. The second target of this SDG alludes to ending all forms of violence against women and girls by 2030 (UN-women), which includes IPV.

Prevalence and factors associated with IPV against women in Rwanda have been well described [27], however, there is insufficient data on IPV against men. There has been growing interest in spatial analysis techniques in recent years as a tool for an in-depth understanding of public health problems including identifying IPV hotspots, spatial distribution, patterns and effects. Intimate partner violence has been shown in some studies to be spatially distributed [1, 27, 28], however none of the studies have used flexible structured additive regression models to estimate district spatial effects while adjusting for other covariates in studying IPV against men and women. A recent study [13] using the 2019/2020 Rwandan Demographic and Health Survey (RDHS) dataset looked at the factors associated with IPV against women and men but did not consider spatial heterogeneity of IPV. Health and social issues are known to have spatial dependency; ignoring to take into account such random effects may result in biased estimates [29].

This study estimates the prevalence of IPV against women aged 15–49 years and men aged 15–59 years as well as the associated factors by applying Bayesian structured additive logistic regression modelling to the 2019/2020 Rwanda DHS dataset. The spatial heterogeneity of IPV against men and women is determined and the contrasts between the two are highlighted.

Methods

Study design and setting

This study was a secondary data analysis of the Rwanda Demographic and Health Survey (RDHS) cross-sectional study. Rwanda is a low-income country with a population of about 13.7 million [30]. It shares borders with Burundi, Uganda, the Democratic Republic of Congo (DRC) and Tanzania. Rwanda has two major ethnic groups, that is, the Hutu and Tutsi who account for more than four-fifths of the population, highlighting a more homogeneous population compared to other African countries. About 75% of this population live in rural areas. The country is divided administratively into 5 regions with the capital in Kigali and 30 districts nested in these regions. The country’s economy relies on agriculture, tourism, mining and donor support [31, 32].

Source of data and sample

The study population was women 15–49 years and men 15–59 years who lived in Rwanda at the time of the survey. The detailed sampling methods are described [33] elsewhere. Briefly, a stratified, two stage cluster random sampling of households was done. Among all the respondents in the survey, 1947 women and 1371 men consented and responded to the domestic violence survey questions and these were included in our study.

The dataset was obtained from the 2019/2020 Rwandan Demographic and Health Survey. A request to use the dataset was made to the Demographic and Health Survey (DHS) team (ICF International-DHS) and permission was granted. The RDHS was reviewed and approved by the Rwanda National Ethics Committee (RNEC) and the ICF Institutional Review Board. The DHS is a nationally representative cross-sectional survey that among other things monitors domestic violence and other public health indicators. The domestic violence module in the DHS, from which our data is derived, used a shortened and modified conflict tactics scale (CTS) [34] to measure different forms of IPV [7, 35] and domestic violence in general for both women 15–49 years and men 15–59 years.

Measurements and covariates

In this study, three forms of spousal violence were considered, that is emotional, physical and sexual violence to generate the outcome. The outcome was a binary variable, which measures IPV in three dimensions, from the questions; (i) Ever experienced physical violence? (ii) Ever experienced sexual violence? and (iii) Ever experienced emotional violence? We also generated other variables to capture; (1) the couple’s financial decision-making strategies i.e. husband making decisions alone or with a partner, and (2) a controlling behaviour. Background characteristic variables such as region, place of residence, age, respondent’s level of education, partner’s educational level, respondent’s age and wealth level were considered as covariates. In addition, age difference, from the respective ages of partners/couples in the dataset was generated from existing variables. Since IPV may be associated with the location of women’s or men’s residences, the spatial heterogeneity was estimated at the district level and these spatial effects were modelled simultaneously with linear and nonlinear effects [1, 27, 36, 37]. The district was selected to provide a low administrative IPV spatial estimate that allows the implementation of target interventions for both men and women to reduce IPV rates. The analysis was performed separately for both women and men.

Bayesian structured additive logistic regression model

Let \({y}_{ij}\) be the IPV status for a woman or man \(i\) in district \(j\). \({y}_{ij}=1\) if the woman or man \(i\) in district \(j\) experienced some form of IPV and \({y}_{ij}=0\) otherwise. A vector \({X}_{ij}=({x}_{ij1},{x}_{ij2},\dots ,{x}_{ijp}{)}^{{\prime }}\) contains \(p\) continuous covariate random variables and \({Z}_{ij}=({z}_{ij1},{z}_{ij2},\dots ,{z}_{ijr}{)}^{{\prime }}\) contains some r categorical variables. In our study, \(p=1\) and \(r=5\). This study assumed that the dependent variable, \({y}_{ij}\) is a Bernoulli-distributed random variable with \({y}_{ij}|{p}_{ij}\sim Bernoulli({p}_{ij})\) with an unknown \(E\left({y}_{ij}\right)={p}_{ij}\), being related to the covariates through the logit link function

$$g\left({p}_{ij}\right)={X}_{ij}^{{\prime }}\beta +{Z}_{ij}^{{\prime }}\theta \hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\left(1\right)$$

where \(\beta\) is the \(p\) dimensional vector of coefficients for the continuous random variables, and \(\theta\) is an \(r\) dimensional vector of coefficients for categorical random variables. To assess both the non-linear effects of continuous random variables and spatial autocorrelation, we employed a semi-parametric model which utilizes a penalized regression approach [38]. The penalized regression approach is a non-parametric method of ordinary least squares (OLS) which relaxes the highly restrictive linear predictor for a versatile semi-parametric predictor [38, 39]. The flexible semi-parametric predictor is defined by:

$$g\left({p}_{ij}\right)=\sum _{v=1}^{p}{f}_{v}\left({x}_{ijv}\right)+{f}_{spat}\left({s}_{j}\right)+{Z}_{ij}^{{\prime }}\theta \hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\left(2\right)$$

where \({f}_{v}(.)\) represents the non-linear twice differentiable smooth function for the continuous covariates and \({f}_{spat}\left({s}_{j}\right)\) is the variable that denotes the spatial effects for each region. In our study, as in [38], we consider a convolution approach to the spatial effects. The assumption is that the spatial effects can be decomposed into two pure components, that is, spatially structured and spatially unstructured effects given as \({f}_{spat}\left({s}_{j}\right)={f}_{str}\left({s}_{j}\right)+{f}_{unstr\left({s}_{j}\right)}\). The final model for our study then becomes:

$$g\left({p}_{ij}\right)=\sum _{v=1}^{p}{f}_{v}\left({x}_{ijv}\right)+{f}_{str}\left({s}_{j}\right)+{f}_{unstr}\left({s}_{j}\right)+{Z}_{ij}^{{\prime }}\theta \hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\hspace{0.25em}\left(3\right)$$

Statistical methods

Statistical analysis was performed in STATA version 17 (StataCorp, College Station, TX, USA) and R version 4.1.0. The data were survey set to adjust analysis results for design and other survey effects. Weighted percentage frequencies were used to estimate the prevalence of IPV and associated 95% confidence intervals (CI). Cross-tabulations and chi-square tests were used for the bivariate analysis of categorical variables. Mean and standard deviation or median and interquartile ranges were used to summarize continuous variables depending on the distribution. Univariate and multivariable logistic regression models (conventional methods) were used to determine factors associated with IPV in both women and men. To fit the Bayesian structured additive logistic regression model, variables to be included in the model were selected using nominal p-value of 0.2 or less in the conventional multivariable logistic regression and some were literature informed. Significant factors in the multivariable model were adjusted for in the spatial model. Spatial distribution and patterns including Global Moran I and local Moran I for assessing autocorrelation and local clustering were assessed using R v 4.1.0. A conditional autoregressive (CAR) generalized structured additive logistic regression model with the binomial link was fitted using the BayesX R package accounting for spatial effects. In addition, non-linear effects of experiencing IPV for some continuous covariates were also assessed. All tests were two-sided and a p-value of less or equal to 0.05 was considered to indicate statistical significance. The 95% credible intervals are reported with adjusted odds ratios (aOR) for the full Bayesian inference.

Ethical considerations

We sought permission to use DHS data from the DHS program via their website and agreed to all standards and laws applicable in accessing and utilizing DHS data. The Rwandan DHS was ethically approved by the Rwandan Health Research Committee, Institutional Review Board of ICF Macro, and Centre for Disease and Control (CDC) in Atlanta, GA, USA, and Prevention IRB [40].

Results

Baseline characteristics

A total of 1947 women were questioned on domestic violence. At least half of the women were aged 35–49 years, while most of the women resided in rural areas (83.5%). More than 60% of the women had primary education and most of them were from the eastern region of Rwanda with the least coming from Kigali (13.4%). Utmost 20% were from the richest households and more than 80% of the women were employed in the past 12 months. For marital status, a majority of the women were not living together with their partner and more than 60% of the study participants had partners who completed primary education. Almost 50% of the women stated that their partners had a controlling behaviour and that they did not make financial decisions together with their partners (68.5%). Most of the women had partners who were older than them (79.6%) and partners that consumed alcohol (63.2%). Majority of the participants were not currently pregnant (91.6%), see Table 1.

Table 1 Baseline characteristics of women and men who participated in the domestic violence module in the survey, Rwanda 2019/2020

A total of 1371 men were interviewed about domestic violence with a majority aged between 35 and 49 years. Most of the men (80%) resided in rural areas and had completed their primary education. Almost 25% of the men were from the eastern region with the least from the Kigali region. Less than 20% were from the richest households, while a majority were employed in the past 12 months. More than 60% of the men were not living with a partner and the majority made financial decisions together with their partners. Most of the men were older than their partners (71.9%). Utmost 70% of the study participants had partners who consumed alcohol. The majority of the men had partners who were not currently pregnant (89.8%) and almost all the men had one partner (92.3%), see Table 1.

Out of the 1947 women interviewed on domestic violence, 668 (34.6%) reported experience of emotional violence, 711 (36.3%) physical violence and 297 (15.6%) sexual violence. Overall, 895 (45.9%) experienced all forms of violence. For the men, 216 (16.7%) experienced emotional violence, 114 (8.7%) physical violence and 14 (1.1%) sexual violence. In total, 245 (18.4%) of the men experienced all forms of violence, see Table 2.

Table 2 All forms of violence with subitems [women & men]

Spatial patterns and distribution of IPV

Figure 1 shows the spatial distribution of IPV prevalence for both men and women in Rwanda. High IPV prevalence among women is noted in Gakenge and Gicumbi districts in the Northern province; Gisagara, Kamonyi, Nyanza and Nyaruguru districts in the Southern province; Kirehe and Nyagatare districts in the Eastern province and Karongi, Ngororero, Nyabihu, Nyamasheke and Rutsiro districts in the western province. Kigali districts had lower prevalence of IPV perpetrated to either women or men. For men, high IPV was noted in Gisagara, Huye, Kamonyi, and Muhanga districts in the Southern province; Bugesera, Kirehe, Ngoma and Nyagatare in the Eastern province; and Karongi, Nyamasheke, Rubavu and Rusizi in the Western province. Of the 30 districts in Rwanda, 13 districts had IPV prevalence of more than 50% for women and four districts reported high IPV in both women and men. Districts with high prevalence in the range of 54.2%-61% against women are scattered around Rwanda suggesting there may not be a global systematic spatial pattern. The negative global moran I value (Moran’s I = -0.1124, p = 0.76) for women IPV and (Moran’s I = -0.022, p = 0.456) for men indicates that IPV was not clustered together. However, the local indicators of spatial autocorrelation indicated some local clustering.

Fig. 1
figure 1

Women and Men IPV prevalence in Rwanda

Tables S1 and S2 were used to inform the variables for inclusion in the adjusted generalized structured additive regression model. From the adjusted generalized structured additive regression model, the model with spatial random effects (DIC = 1788) was compared to the model without a spatial component (DIC = 1803) showing the relevance of adding the spatial component although the global Moran’s I was not significant. Table 3 shows the adjusted odds of IPV with corresponding credible intervals (CI) for the selected model.

Table 3 Bayesian Structured Additive Logistic Regression Model

Women making financial decisions together with their partner (aOR = 0.52; CI: 0.40–0.67), women in richer households (aOR = 0.49; CI: 0.39–0.70) and women in richest households (aOR = 0.48, CI: 0.33–0.70) are significantly at lower odds of experiencing IPV. In addition, women living with controlling partners (aOR = 5.80; CI: 4.55–7.39) and women with partners who consume alcohol (aOR = 3.11; CI: 2.44–3.96) are at significantly higher odds of experiencing IPV. In particular, women with controlling partners have about 6 times higher odds of experiencing IPV while those with partners that drink alcohol have about 3 times higher odds of experiencing IPV (Table 3).

Men who make financial decisions together with their partner (aOR = 0.61; CI: 0.41–0.93) and men in the richest households (aOR = 0.54; CI: 0.30–0.94)) are at significantly higher odds of experiencing IPV. In addition, men living with controlling partners (aOR = 7.78; CI: 5.35–11.52) and men with partners who consume alcohol (aOR = 2.31; CI: 1.63–3.27) are at significantly higher odds of experiencing IPV. In particular, men with controlling partners have about 8 times higher odds of experiencing IPV while those with partners that drink alcohol have about 2 times higher odds of experiencing IPV. Spatial correlation explains 64% of the residual spatial variation amongst women who experience IPV while in men experiencing IPV, 62% of the spatial variation is explained by district heterogeneity (Table 3).

Figure 2 (left) shows elevated odds of IPV among women in Nyagatare, Gicumbi, Burera, Musanze, Nyabihu, Gakenke, Rutsiro, Nyarugenge, Gasabo and Kirehe districts which are mainly in the Northern province and Kigali after adjusting for selected individual, partner and household level covariates.

Fig. 2
figure 2

The predicted prevalence of IPV after adjusting for spatial as well as individual, partner and household level covariates among women (left) and men (right) in Rwanda

There are some pronounced IPV district-associated effects after adjusting for individual-level, partner as well as household-level covariates among women. The district-associated adjusted odds of IPV range from 0.56 to 1.40 with the structured spatial effects dominating the unstructured spatial heterogeneity effect as shown by the ratio (\(\phi =\frac{{f}_{str}\left({s}_{j}\right)}{{f}_{str}\left({s}_{j}\right)+{f}_{unstr}\left({s}_{j}\right)}=\frac{0.132}{0.132+0.074}=64.1\%\)).

Figure 2 (right) shows elevated odds of IPV among men in Gakenke, Rubavu, Karongi, Nyamasheke, Ruhango, Nyanza, Gisagara, Bugesera, Kamonyi, Ngoma, Kirehe and Kayonza districts which are mainly in the Eastern and Western provinces. The district-associated adjusted odds of IPV for men range from 0.52 to 1.60 with the spatial heterogeneity effects dominating the structured spatial effects as shown by the ratio (\(\phi =\frac{{f}_{unstr}\left({s}_{j}\right)}{{f}_{str}\left({s}_{j}\right)+{f}_{unstr}\left({s}_{j}\right)}=\frac{0.211}{0.130+0.211}=62.0\%\)). This indicates a general spatial variability in women perpetrated IPV against men in Rwanda with the Southern, Eastern and Western districts having elevated odds for men to experience IPV. Women experienced IPV show some level of structure than that experienced by men.

Discussion

This study explores the spatial distribution, patterns, and associated factors of IPV against men and women using the 2019/2020 RDHS survey data, by fitting Bayesian structured additive logistic regression models. The study also uses districts administrative areas as spatial units, which is beneficial for policy decision-making especially for IPV.

We found a higher IPV prevalence for women (46%) than men (18%.) Our study is one of the very few studies that have reported men as victims of IPV [41]. Though the IPV prevalence in men is lower than females, it is relatively high and raises concern on the magnitude of IPV among men. The lower IPV prevalence in men can be explained by under-reporting of the information in the survey. Men may shy away from reporting as a ways for preserving their muscularity in the society. Moreover, women tend to have formal institutions which protect them; hence, they are normally better informed and feel more comfortable in reporting violence than men. There are many support services for women than men, therefore, it is important to involve the men in IPV campaigns and offer men appropriate support they may need so that they get the help they need like women. Though the recent SDG5.2.1 focuses on eliminating IPV among women and girls, there is silence on the men since this study has found that men are also affected with IPV regardless of the lower rates.

Our results are similar to those found in other studies [13], however our study highlights the specific districts in Rwanda where IPV prevalence is high for both women and men. Using a bivariate choropleth, IPV perpetrated against women was shown to be higher in the North-Western districts of Rwanda whereas for men it was shown to be more prevalent in the Southern districts. A few districts presented high IPV for both men and women. The study also found that reported IPV against women was almost twice that experienced by men. The lower prevalence in men can be explained by the African social norms. Men tend not to report being victims of IPV, especially if perpetrated by women as this may be construed as being weak [12,13,14,15,16,17,18, 42].

For IPV against women, physical violence was the most common (36.3%), followed by emotional violence (34.6%). Whereas, for men emotional violence (16.7%) was the most common. In line with our findings, a study in West Africa reported psychological or emotional violence being more commonly experienced among men [43]. This can be explained by the fact that men are more muscular and therefore do not usually suffer from physical violence perpetrated by their women intimate partners, whereas women tend to talk about their problems more than men and therefore are less likely to experience emotional violence.

Accounting for variability due to respondent age, how the couples made financial decisions, household wealth levels, partner controlling behavior and partner’s alcohol consumption, our study highlights districts that have elevated IPV prevalence. These districts are hotspots for IPV and the unobserved structural effects may need to be interrogated in order to mitigate against high IPV prevalence. We found higher odds of IPV against women mainly in the North-western districts and the spatial effects were dominated by spatially structured effects contributing 64%. On the other hand, higher odds of IPV were observed for men in the Southern districts of Rwanda and spatial effects were dominated by district heterogeneity accounting for 62%. There were no statistically significant district clusters for IPV in both men or women. Contrary to our findings, a study on IPV in Namibia [1], found a strong and significant clustering (structured component of IPV) mainly driven by poor economic conditions and cultural beliefs in the severely affected regions.

Rwanda is a more homogenous country with mostly Hutus and Tutsis constituting over 80% of the population, for which 75% reside in rural areas. Our results, adjusting for spatial effects, reveal that being in a higher socioeconomic status is protective to experiencing IPV for both women and men. This is consistent with results by [13, 43,44,45,46] who showed similar results using 2015 and 2020 DHS data. In addition, we also show that after adjusting for spatial effects, couples making financial decisions together are more protected against experiencing IPV. Several studies [44, 45] that have been conducted in Rwanda support this finding, including a clustered randomized trial study by [46] which highlighted that couples that discussed and completed tasks together had reduced experiences of IPV as a result of fostered higher relationship qualities. [46] also indicate that relationship quality is a key pathway for healthy behaviours for couples.

Our results show that women with higher education had lower odds of experiencing IPV and so were those married to partners with higher education. Making financial decisions together with a partner and living in the southern region were associated with lower odds of experiencing IPV in women. Women whose partners were of the same age, with controlling behaviour as well as alcohol consumers were at significantly higher odds of experiencing IPV. Similar findings were observed in men where higher education, those who were wealthier, and who made financial decisions together with their wives were associated with reduced risk of IPV. Our findings on the factors associated with IPV are in line with findings made by [13, 47,48,49].

The study did not fall short of some limitations, there was a lack of temporality between the factors and IPV, hence, we cannot ascertain causality of these factors on IPV. Not all potential factors were considered in the analysis, particularly in the Bayesian structured additive logistic regression model which is usually computationally intensive to attain convergency. Different variables were used for the male and female models which made it difficult to compare the effects between the two sub-populations. In addition, due to the nature of IPV being self-reported, men may have under-reported IPV perpetrated against them.

Conclusion

This study found significant spatial variations in IPV against both men and women in Rwanda. Interventions that improve relationship quality within couples should be encouraged in the severely affected districts as well as across the country. The findings from our study show that making financial decisions together between partners helps in reducing intimate partner violence. Consumption of alcohol and unfavourable controlling behaviour by a partner may contribute to increase in IPV. There is a need to intensify campaigns that fight GBV and instill defense mechanism skills to individuals for future defense against GBV by their partners. Such campaigns and interventions should target both women and men if the SDG5 goal to eliminate IPV is to be achieved. Men cannot remain a blind spot in the fight of IPV.

Further research involving qualitative interviews with victims of IPV, both women and men, can provide a deeper understanding of drivers of IPV in Rwanda and similar African settings.

Data Availability

The dataset generated and analysed during the current study are not publicly available since we received a data access letter from the DHS team https://dhsprogram.com/ specific to our project but are available from the DHS team upon request.

References

  1. Ojifinni O, Maposa I, Ibisomi L. Bayesian semi-parametric spatial modelling of intimate partner Violence in Namibia using 2013 Demographic Health Survey Data. BMC Womens Health. 2021;21(1):1–10.

    Article  Google Scholar 

  2. Krug EG, Mercy JA, Dahlberg LL, Zwi AB. The world report on Violence and health. The Lancet. 2002;360(9339):1083–8.

    Article  Google Scholar 

  3. Violence. against women [https://www.who.int/news-room/fact-sheets/detail/violence-against-women].

  4. Devries KM, Mak JY, Garcia-Moreno C, Petzold M, Child JC, Falder G, Lim S, Bacchus LJ, Engell RE, Rosenfeld L, et al. Global health. The global prevalence of intimate partner violence against women. Science. 2013;340(6140):1527–8.

    Article  PubMed  CAS  Google Scholar 

  5. Fulu E, Jewkes R, Roselli T, Garcia-Moreno C, Men UNM. Violence research t: prevalence of and factors associated with male perpetration of intimate partner violence: findings from the UN multi-country cross-sectional study on men and Violence in Asia and the Pacific. Lancet Glob Health. 2013;1(4):e187–207.

    Article  PubMed  Google Scholar 

  6. Hindin MJ, Kishor S, Ansara DL. Intimate partner violence among couples in 10 DHS countries: predictors and health outcomes. Macro International Incorporated; 2008.

  7. MacQuarrie K, Winter R, Kishor S. Spousal Violence and HIV: exploring the linkages in five sub-saharan African countries. ICF International; 2013.

  8. Verduin F, Engelhard EA, Rutayisire T, Stronks K, Scholte WF. Intimate partner violence in Rwanda: the mental health of victims and perpetrators. J Interpers Violence. 2013;28(9):1839–58.

    Article  PubMed  Google Scholar 

  9. Organization WH. The world health report 2002: reducing risks, promoting healthy life. World Health Organization; 2002.

  10. Okenwa LE, Lawoko S, Jansson B. Exposure to intimate partner violence amongst women of reproductive age in Lagos, Nigeria: prevalence and predictors. J Fam Violence. 2009;24(7):517–30.

    Article  Google Scholar 

  11. Abramsky T, Lees S, Stöckl H, Harvey S, Kapinga I, Ranganathan M, Mshana G, Kapiga S. Women’s income and risk of intimate partner violence: secondary findings from the MAISHA Cluster randomised trial in North-Western Tanzania. BMC Public Health. 2019;19(1):1–15.

    Article  Google Scholar 

  12. Ayodele JO. The Socio-cultural causes of male victimisation in domestic contexts in Lagos, Nigeria: a qualitative analysis. Int J Criminal Justice Sci 2017, 12(2).

  13. Bahati C, Izabayo J, Munezero P, Niyonsenga J, Mutesa L. Trends and correlates of intimate partner violence (IPV) victimization in Rwanda: results from the 2015 and 2020 Rwanda Demographic Health Survey (RDHS 2015 and 2020). BMC Womens Health. 2022;22(1):368.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Safariolyaei N, Amiri M. Impact of intimate partner violence against men. Int J Law Educ Social Sports Stud. 2017;4(1):43–7.

    Google Scholar 

  15. Moore T. Suggestions to improve outcomes for male victims of domestic abuse: a review of the literature. SN Social Sciences. 2021;1(10):1–18.

    Article  Google Scholar 

  16. Gubi D, Wandera SO. Prevalence and correlates of intimate partner violence among ever-married men in Uganda: a cross-sectional survey. BMC Public Health. 2022;22(1):1–9.

    Article  Google Scholar 

  17. Lanre AO, Olugbenga OO, Adenike Iyanuoluwa O, Olusegun BJ, Olalekan AW. Assessment of prevalence and forms of violence against married men in Olorunda Local Government of Osun State, Nigeria. Int J Soc Behav Sci. 2014;2(1):1–10.

    Google Scholar 

  18. Ameh N, Shittu S, Abdul M, Bature S, Oyefabi O. Burden and characteristics of Domestic Violence among males in a sub saharan African setting. Niger J Med. 2012;21(4):412–5.

    PubMed  CAS  Google Scholar 

  19. Ahinkorah BO. Polygyny and intimate partner violence in sub-saharan Africa: evidence from 16 cross-sectional demographic and health surveys. SSM-Population Health. 2021;13:100729.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Demeke MG, Shibeshi ET. Intimate partner Violence against women of reproductive age and associated factors during COVID-19 pandemic in Northern Ethiopia, 2021: a community-based cross-sectional study. Front Global Women’s Health. 2023;3:977153.

    Article  Google Scholar 

  21. Tiruye TY, Harris ML, Chojenta C, Holliday E, Loxton D. Determinants of intimate partner violence against women in Ethiopia: a multi-level analysis. PLoS One. 2020;15(4):e0232217.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Stiller M, Bärnighausen T, Wilson ML. Intimate partner violence among pregnant women in Kenya: forms, perpetrators and associations. BMC Womens Health. 2022;22(1):1–25.

    Article  Google Scholar 

  23. Barchi F, Winter SC, Mbogo D, Thomas B, Ammerman B. Exploring the relationship between participation in an adult-women’s Soccer League and Intimate Partner Violence in Rural Kenya. J Interpers Violence. 2022;37(9–10):NP7907–31.

    Article  PubMed  Google Scholar 

  24. Hilliard S, Bukusi E, Grabe S, Lu T, Hatcher AM, Kwena Z, Mwaura-Muiru E, Dworkin SL. Perceived impact of a land and property rights program on Violence against women in rural Kenya: a qualitative investigation. Violence against Women. 2016;22(14):1682–703.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Citaristi I. United Nations Entity for Gender Equality and the empowerment of Women—UN women. The Europa Directory of International Organizations 2022. edn.: Routledge; 2022: 189–92.

  26. Transforming our world. : the 2030 Agenda for Sustainable Development [https://sustainabledevelopment.un.org/post2015/transformingourworld].

  27. Habyarimana F, Zewotir T, Ramroop S. Structured spatial modeling and Mapping of Domestic Violence against Women of Reproductive Age in Rwanda. J Interpers Violence. 2021;36(5–6):2430–54.

    Article  PubMed  Google Scholar 

  28. Seid AM, Mishra GD, Dobson AJ. The association between childhood sexual abuse and historical intimate partner violence with body mass index and diabetes: evidence from the Australian longitudinal study on women’s Health. Prev Med. 2022;161:107134.

    Article  PubMed  Google Scholar 

  29. Julian MW. The consequences of ignoring multilevel data structures in nonhierarchical covariance modeling. Struct Equ Model. 2001;8(3):325–52.

    Article  Google Scholar 

  30. Population size and Population characteristics. [https://www.statistics.gov.rw/statistical-publications/subject/population-size-and-population-characteristics].

  31. Enterprises AFOH. Rwanda Economic Update.

  32. Behuria P. Committing to self-reliance and negotiating vulnerability: understanding the developmental change in Rwanda. SOAS University of London; 2015.

  33. National Institute of Statistics of Rwanda MoH, ICF. Rwanda Demographic and Health Survey 2019-20 final report. In.: NISR and ICF Maryland, USA; 2021.

  34. Straus MA. Measuring intrafamily conflict and Violence: the conflict tactics (CT) scales. Physical Violence in American families. edn.: Routledge; 2017: 29–48.

  35. MacQuarrie KL, Mallick L. Intimate partner violence and contraceptive discontinuation in Kenya. Afr J Reprod Health. 2021;25(2):17–27.

    PubMed  Google Scholar 

  36. Fahrmeir L, Sagerer F, Sussmann G. Geoadditive regression for analyzing small-scale geographical variability in car insurance. Blätter Der DGVFM. 2007;28(1):47–65.

    Article  Google Scholar 

  37. Kazembe LN, Chirwa TF, Simbeye JS, Namangale JJ. Applications of bayesian approach in modelling risk of malaria-related hospital mortality. BMC Med Res Methodol. 2008;8(1):1–14.

    Article  Google Scholar 

  38. Ngesa O, Mwambi H, Achia T. Bayesian spatial semi-parametric modeling of HIV variation in Kenya. PLoS ONE. 2014;9(7):e103299.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Demaret L, Dyn N, Iske A. Image compression by linear splines over adaptive triangulations. Sig Process. 2006;86(7):1604–16.

    Article  Google Scholar 

  40. National Institute of Statistics of Rwanda (NISR) [Rwanda] MoHMR, and ICF International Agency for Research on Cancer (IARC). Rwanda Demographic and Health Survey 2014-15. In. Rockville. Maryland, USA: NISR, MOH, and ICF International; 2015.

    Google Scholar 

  41. Wörmann X, Wilmes S, Seifert D, Anders S. Males as victims of intimate partner Violence—results from a clinical-forensic medical examination centre. Int J Legal Med. 2021;135(5):2107–15.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Taylor JC, Bates EA, Colosi A, Creer AJ. Barriers to men’s help seeking for intimate partner Violence. J Interpers Violence. 2022;37(19–20):NP18417–44.

    Article  PubMed  Google Scholar 

  43. Oyediran K, Spencer CM, Stith SM. Men as victims of intimate partner Violence in Cameroon and Sierra Leone. J Interpers Violence 2022:08862605221120883.

  44. Zegenhagen S, Ranganathan M, Buller AM. Household decision-making and its association with intimate partner Violence: examining differences in men’s and women’s perceptions in Uganda. SSM Popul Health. 2019;8:100442.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Raj A, Silverman JG, Klugman J, Saggurti N, Donta B, Shakya HB. Longitudinal analysis of the impact of economic empowerment on risk for intimate partner Violence among married women in rural Maharashtra, India. Soc Sci Med. 2018;196:197–203.

    Article  PubMed  Google Scholar 

  46. Dunkle K, Stern E, Chatterji S, Heise L. Effective prevention of intimate partner violence through couples training: a randomised controlled trial of Indashyikirwa in Rwanda. BMJ Global Health. 2020;5(12):e002439.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Semahegn A, Mengistie B. Domestic Violence against women and associated factors in Ethiopia; systematic review. Reproductive Health. 2015;12(1):1–12.

    Article  Google Scholar 

  48. Tadesse AW, Tarekegn SM, Wagaw GB, Muluneh MD, Kassa AM. Prevalence and associated factors of intimate partner violence among married women during COVID-19 pandemic restrictions: a community-based study. J Interpers Violence. 2022;37(11–12):NP8632–50.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Regassa N. Intimate partners violence in Southern Ethiopia: examining the prevalence and risk factors in the Sidama zone. Int J Sociol Anthropol. 2012;4(1):13–22.

    Article  Google Scholar 

Download references

Acknowledgements

We thank the DHS Program for providing access to the Rwanda Demographic and Health Survey data set.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

IM conceptualized the topic with assistance from all co-authors who contributed equally to the data analysis, writing, research, critical review, and proofreading of the manuscript. All authors read and approved the final manuscript. AM contributed as a senior author, overseeing the writing of the whole manuscript.

Corresponding author

Correspondence to Halima S. Twabi.

Ethics declarations

Ethical approval and consent to participate

The study is a secondary data analysis of publicly available and non-identifiable data. The IRB approved procedures for Demographic and Health Surveys are carried out in accordance with relevant guidelines and regulations (e.g. Declaration of Helsinki). In particular, the ICF IRB ensures that the survey complies with the U.S. Department of Health and Human Services regulations for the protection of human subjects (45 CFR 46), while the host country IRB ensures that the survey complies with laws and norms of the nation (https://dhsprogram.com/Methodology/Protecting-the-Privacy-of-DHS-SurveyRespondents.cfm). The procedures do not, in any way, allow respondents, households, or sample communities to be identified. There are no names of individuals or household addresses in the data files. The geographic identifiers only go down to an enumeration area and each enumeration area has a primary sampling unit number in the data file but the numbers do not have any labels to indicate their names or locations. Before each interview is conducted, an informed consent statement is read to the respondent, who may accept or decline to participate. More details on ethical approval for DHS datasets may be accessed from https://dhsprogram.com/Data/terms-of-use. cfm. Furthermore, we received authorization to use the Rwandan dataset from the DHS Program.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Maposa, I., Twabi, H.S., Matsena-Zingoni, Z. et al. Bayesian spatial modelling of intimate partner violence and associated factors among adult women and men: evidence from 2019/2020 Rwanda Demographic and Health Survey. BMC Public Health 23, 2061 (2023). https://doi.org/10.1186/s12889-023-16988-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-023-16988-8

Keywords