This study uses longitudinal data from women and facilities collected in 2011, prior to the start of ISSU program activities (baseline data), and 2015 (endline data) in all six cities. Longitudinal data for women allow us to use fixed effects methods that allow each woman to act as her own control [13] . At baseline, in each city, we collected data from a representative sample of women identified using a two-stage sampling design. In the first stage, we used the 2002 General Population and Housing Census's list of census districts (updated in 2009) as our primary sampling units (PSU) to select a random sample of PSU in each city. On average, census districts contain about 150–200 households. In total, we selected 268 PSU: 64 PSU from Dakar; 32 each from Guédiawaye, Pikine, and Mbao; and 54 each from Mbour and Kaolack. Prior to selection, we worked with municipal leaders to classify neighborhoods as poor and non-poor based on five characteristics: type of housing, residential security, neighborhood density, access to water, and access to flush toilets. Municipal leaders also gave an overall classification of the poverty level of the neighborhood. We weighted the five characteristics of neighborhoods (access to water had a weight of 2; toilets a weight of 1.5; and all others a weight of 1) and then summed. We considered those neighborhoods scoring in the lowest 40% as poor and all others as non-poor. Prior to selection, we stratified PSU by poor and non-poor and half of selected PSU were from the poor strata to increase inclusion of poor households and women. In the selected PSU, we listed/mapped all households and randomly selected 21 households for interview with equal probability [14] .
In total, we selected 5628 households. In each selected household, all women of reproductive age (15–49) were eligible for interview. A trained female interviewer approached all eligible women and asked for their signed consent to participate. At baseline, we interviewed 9614 women across the six cities to provide a representative sample for each city at baseline. At endline, we tracked all baseline women who were usual residents (not visitors) of the household (n=9421) and, if found, asked them for consent to be re-interviewed.
The Institutional Review Board at the University of North Carolina at Chapel Hill and the Comité National d'Ethique in Senegal approved all study procedures.
Women interviewed at both baseline and endline are the analysis sample for examining the impact of the program on modern contraceptive use over the four-year follow-up period. We also undertake a sub-sample analysis of “poor” women, identified from the matched sample as those women who were in the two lowest wealth quintiles (poorest and poor women) at both baseline and endline.
At baseline and endline, we also collected data from health facilities in the study cities. We approached all public and private facilities offering reproductive health services (n=269) for interview. At each study facility, an interviewer administered a facility audit, provider interviews (up to four per facility), and exit interviews [15] .
The key dependent variable for this analysis is use of modern contraception. At baseline and endline, trained interviewers asked women if they were using a contraceptive method to delay or avoid childbearing and if yes, interviewers asked what method the woman was using. Modern methods of contraception include male and female sterilization, daily pill, IUD, implants, injectables, male and female condoms, emergency contraception, Standard Days Method, and lactational amenorrhea. We coded women who reported traditional method use (e.g., rhythm method, withdrawal, or folkloric methods) as non-users. At each time point, we coded women who were abstinent or not sexually active as non-users.
This analysis examines the impact of exposure to various ISSU and non-ISSU demand generation activities and the ISSU and non-ISSU supply-side activities on the probability of modern method use. We classified the exposure variables as demand-side and supply-side. At endline, we included detailed questions on exposure to ISSU specific radio and television activities to assess the contribution of ISSU programming in target cities above and beyond the national-level radio and television programming (see Appendix Table A for details of questions used to measure exposure). We coded each of the exposure variables as dichotomous variables categorized as exposed versus unexposed. All exposure measures that were ISSU-specific and did not exist at baseline were coded zero (unexposed) at baseline.
We included four supply-side variables in this analysis. One of these is specific to the ISSU program: exposure to the ISSU Informed Push Model (IPM). To link the IPM variable to the women's data, we created a variable capturing whether the woman lived within 1 km of a facility with IPM at the time of the endline survey. This variable is coded zero at baseline since the IPM program activity did not exist at that time. We created three other supply variables based on this same 1 km buffer around where women live. These were: (a) woman lives within 1 km of a facility with any stock outs in the last 30 days; (b) woman lives within 1 km of a facility with a quality improvement committee; and c) woman lives within 1 km of a facility with observed FP guidelines.
We undertake descriptive analyses as well as multivariate analyses. We apply fixed effects regression to the pooled samples (baseline and endline) to reduce bias associated with self-selection, recall, and program targeting to underserved areas due to time invariant unobservables. Our estimation methods control for the possibility of endogenous attrition in the endline sample to the extent that attrition is due to unobserved fixed characteristics of the individuals. Fixed effect methods do not control for time varying unobservables. However, given the relatively close spacing of observations, time-varying unobservables are likely not a major source of residual bias in this setting. We estimated a classic fixed effects model whereby changes in modern contraceptive use depend on changes in time-varying individual characteristics and program exposure. All models control for marital status, age, education, city, religion, and wealth (contact corresponding author for full models). We perform two models, one including all women and the other for those in the two lowest wealth quintiles at baseline and endline. All analyses adjust for correlation at the community level using Huber-White type sandwich estimators for standard errors.
Benson A., Calhoun L., Corroon M., Gueye A., Guilkey D., Kebede E., Lance P., O'Hara R., Speizer I.S., Stewart J, & Winston J. (2018). The Senegal urban reproductive health initiative: a longitudinal program impact evaluation. Contraception, 97(5), 439-444.