2.1 Official data on Roma

Roma are the largest ethnic minority in Europe. In all the countries where they live, they experience severe social exclusion and poverty. They mainly perform low-skilled jobs, live in segregated areas of the main cities, and do not participate in the political and cultural life [Open Society Institute (2007)]. Their living conditions are often so different from those of the majority population that it is difficult to find official data documenting their situation. For most Central and Eastern European countries where the majority of the Roma population lives, official data on them are scarce and inaccurate. The 2011 Serbian Census counts 147,604 Roma, corresponding to 2.05% of the total Serbian population, while the Open Society Institute (2007) estimates a number between 350,000 and 500,000, approximately 6% of the overall population. In Belgrade, the 2011 Census records 27,325 Roma (1.65% of the population) and the Open Society Institute (2007)'s numbers are three times higher: they are roughly 80,000 (5%).

The UNICEF MICS from 2010 and the Living Standard Measurement Survey (LSMS) from 2007 are valuable sources of information on the living conditions of the Roma population in Serbia. However, the MICS does not report information on where the households interviewed live and therefore cannot be used in our study. Furthermore, the LSMS only interviewed a boosted sample of internally displaced people, making it a very selected sample.

We make use of these sources in order to compare the characteristics of the Roma and non-Roma populations. As reported in Table A in the Appendix, the average Roma household is composed of 5.6 members vs. a national average of 3.5. The average number of children aged 18 or below is 2.4 per Roma household, while the population average is only 0.86. Almost half of the Roma population (43%) is below 18 years old and the average age is 25, whereas the national average is 35. Half of the Roma households are poor: their average consumption is below the absolute poverty line.Footnote ¹¹ While male employment rates are comparable to those of the majority population (56%), female employment remains very low with only one woman out of 10 working vs. a national average of 40%. Only 89% of children from Roma settlements aged 6–15 attend school and among the adults, 29% have not finished primary school.Footnote ¹² Conversely, 99% of non-Roma aged 6–15 are enrolled in school and only 4% of adults have not completed primary school.

2.2 The sample

We use first-hand collected data obtained through a survey conducted with 300 Roma households of Belgrade, originally aimed at evaluating a remedial education program introduced in Serbia in 2009.Footnote ¹³ All surveyed households have at least one child in the lower four grades of primary school in the year of the survey. They were randomly selected among pupils attending primary schools that were involved in the program at different points in time.

Belgrade is divided into 17 municipalities, among which 10 have Roma settlements, as defined by a population of at least 15 households or 100 individuals [Jaksic and Basic (Reference Jaksic and Basic2010)]. Five of these municipalities are included in our survey design, as the five remaining had no schools involved in this initial stage of the remedial program.Footnote ¹⁴ We suspect that schools in those municipalities did not satisfy the eligibility criterion of having a share of Roma pupils between 5% and 40%.Footnote ¹⁵ The municipalities included in the final sample host 59% of all Belgrade schools. Using MICS (2010), we estimate that our survey design amounts to randomly select pupils from schools enrolling 40% of all Roma children in primary school age.Footnote ¹⁶ The survey took place in Fall 2010. The response rate was 93.46%: 321 households had been contacted and 300 answered. Households were not compensated for their participation.

We acknowledge our survey design suffers from some limitations: we consider women with (i) at least one child aged 6–10; (ii) the child needs to be enrolled in a school; and (iii) schools need to have between 5% and 40% of Roma pupils. While we agree that this pool is certainly different in various ways from the universe of all Roma women, we nonetheless think that the concerns relative to the representativeness of the sample can be mitigated. Indeed, in Serbia, enrollment of Roma aged 6–10 into schools is close to universal and there are not any documented gender differences that could affect our results. Dropouts from primary schools usually arise after the first cycle of 4 years, when children are supposed to enter the second cycle of primary education [Open Society Institute (2007); Battaglia and Lebedinski (Reference Battaglia and Lebedinski2015)]. Because schools with less than 5% of Roma pupils are not selected in the program, we are aware that we might have under-sampled families living in settlements with few Roma households. Nonetheless, the lack of available official data on the actual distribution of Roma in Belgrade does not allow us to quantify the resulting bias. MICS (2010) Roma sample itself is formed by excluding all enumeration areas with 17 or less Roma households, sampling thus from a pool of 46% of the all Roma households. We believe though that we most likely underestimate differences associated with segregation as we trim the sample from its most desegregated households.

The design of the sample also implies that all households observed count at least one child. This may distort the representativity of the results in terms of fertility of the whole Roma population as it de facto removes the extensive margin of the fertility decision [or childlessness as put by Gobbi (Reference Gobbi2013) and Baudin et al. (Reference Baudin, de la Croix and Gobbi2015)]. Using Census data for 2011, we obtain that the rate of childlessness of women over 39 in the Belgrade region was of 7% for Roma and 13% for Serbs. Though substantial, the magnitude of this difference is not likely to explain fully the much higher fertility of Roma women. We thus consider that studying only the intensive margin of the fertility decision is an important step per se.

Figure 1 displays a map of Belgrade with the 13 settlements where the survey was carried out. We use the definition of settlements from an NGO called the Society for Improvement of Roma Settlement, which made an inventory of the Roma presence in Belgrade in 2002 and classified it into clusters based on geographic concentration, natural and urban frontiers, as well as origin and time of in-migration. The number of households selected from each settlement is proportional to its size.

Figure 1. Map of Belgrade with settlements.

We classify settlements as composed of either only Roma people, mostly Roma people, or few Roma people. In the survey, we asked respondents whether in their community/neighborhood, defined as the area corresponding to 200 square meters around their house, there were only Roma people or both Roma and non-Roma. In the latter case, we further asked whether Roma were a minority or a majority. Notice that the definition of neighborhoods here does not exactly coincide with our settlements. Indeed, a neighborhood generally does not cover the entirety of a settlement, but conversely, the neighborhood of someone living at the periphery of the settlement may extend besides the settlement's limits. This is why we use the median perception about own neighborhood in the settlement. In almost two-third of the settlements, all households have the same perception. In the remaining third, there are either a few only Roma in a settlement otherwise perceived as mostly Roma, or the other way around. We believe that taking the median perception, because it smooths out potential outliers in individual perceptions, actually gives a relatively accurate and reliable representation of the reality.

We do not use metrics that are more commonly employed in the segregation literature such as the index of dissimilarity because there is no available information on the total Roma population residing in one particular area of Belgrade. The lack of official data on ethnic composition at the municipality level is one of the reasons why our data—though with some limitations—are a useful source of information.

Our data include household members' demographic characteristics, such as education level, religion, language spoken at home, and information on their dwellings. We also have detailed information on their settlement of residence from the database of Roma settlements, such as the number of inhabitants and main current utilities.

Panel A in Table 1 reports households' characteristics for the 13 settlements in our sample in column 1. Overall, they are in line with official data [LSMS (2007); MICS (2010)], as reported in Table B in the Appendix. On average, women in our sample are 32.5 years old, which is slightly older than those of the MICS sample. This is consistent with our conditioning on having at least one child enrolled in the four lower grades of primary school. As a consequence of this age difference, we also observe that women in our sample are slightly less educated (as educational attainment has risen slightly in recent cohorts), that the share of children below 6 is smaller, while that of children between 6 and 14 is larger, and that the number of adults in the household, and particularly older adults, is smaller. However and importantly for our analysis, the age at first birth as well as the total fertility rate in our sample are very similar to those found in the MICS.Footnote ¹⁷ The slightly larger fertility in the MICS may come from the fact that we measure fertility as the number of children alive at the time of the survey. This way, unlike the MICS, we disregard any child who did not survive.

Table 1. Households' and settlements' characteristics

Columns 5, 6, and 7 represent normalized differences, in bold when statistically significant at the 5% level.

^a 3,643 dinars correspond to roughly 35 euros (1 RSD = 0.009626 euros, November 2011).

^b 15,209 dinars correspond to roughly 146 euros (1 RSD = 0.009626 euros, November 2011).

^c The wealth index ranges between −3.135 and 2.865.

^d 1,709 dinars correspond to roughly 17 euros (1 RSD = 0.009626 euros, November 2011).

Wealth is measured by the first component of a principal component analysis on the presence in the household of various durables and utilities.Footnote ¹⁸ Women received an average monthly income of 3,600 Serbian dinars, which is worth around 35US$, while men made over four times that amount, about 15,000 Serbians dinars, or 150US$. Households mainly receive income from labor, either in the formal or informal sector, rather than social transfers.Footnote ¹⁹

Almost all households in the sample are nuclear with on average a little over two adults, although there are a few exceptions with more than four. They are most likely Muslim and never moved from the settlement they are currently living in. In total, 30% of households comprise adults named with only Serbian names.Footnote ²⁰ They expect that one extra year of schooling increases monthly income by roughly 17US$, corresponding to 5% of the minimum wage.Footnote ²¹ Roma people usually do not perform jobs for which high levels of education are required. They mainly work in the informal sector, without written contracts, often self-employed especially in flea markets and more rarely in factories [LSMS (2007)].Footnote ²²

Columns 2, 3, and 4 report separately means for only, mostly, and few Roma settlements. Households are overall comparable in terms of observable characteristics across settlement types: many normalized differences are smaller than one fourth of the combined sample variation, suggesting that linear regression methods are unlikely to be sensitive to specification changes [Imbens and Wooldridge (Reference Imbens and Wooldridge2009)]. The differences in means are not statistically significant in almost all cases between only and mostly. More substantial differences are found between only or mostly and few. Wealth and households with only Serbian names are higher in few Roma, while share of Muslims is lower.

Panel B of Table 1 reports the characteristics of the settlements. Households are equally located in urban and suburban areas, but only Roma settlements are more likely to be located in suburban areas.Footnote ²³ Settlements do also vary in their access to schools and hospitals, few Roma settlements being better connected. While those could influence fertility through their impact on infant mortality, access to contraceptives or family planning services, as well as the cost of child quality, we believe that these differences should not play a role in the gender imbalances in terms of birth timing, therefore not compromising our general conclusions.

On average, Roma women in our sample have 3.2 children currently alive, of which 54% are boys (panel C of Table 1). A preliminary investigation of our outcome of interest shows that the number of children is significantly lower in few Roma settlements than in the other types of settlements.Footnote ²⁴ On average, in few Roma neighborhoods there are 2.8 children per household, while in mostly Roma and in only Roma neighborhoods there are, respectively, 3.2 and 3.6 children per household. The proportion of boys is not significantly different across the three groups, with slightly more boys in families in few Roma settlements.

3.1 OLS results

The summary statistics show that households in only Roma settlements have a higher fertility. Nonetheless, they could be the reflection of different age structures, socio-economic conditions, family arrangements, or returns to education. In this section, we test whether the gap in fertility across more or less-segregated settlements persists once we take into account household and settlement characteristics. To do so, we estimate the following regression equation using ordinary least squares (OLS)Footnote ²⁵:

(1)$$F_{ijs} = \beta _0 + \beta _1mostly{\rm \_}roma_s + \beta _2\ few{\rm \_}roma_s + \gamma _1X_{ijs} + \delta _1S_s + \epsilon _{ijs}$$

where F _ijs stands for the number of living children for woman i in household j in settlement s and $mostly{\rm \_}roma_s$ and $few{\rm \_}roma_s$ are dummies equal to 1 if the household j lives in a mostly or few Roma settlement, respectively, the omitted category being only Roma.

The set of individual and household characteristics X _ijs and the vector of observable settlement characteristics S _s include controls for different mechanisms that could affect fertility. First, mother's age and age squared are present in all specifications, in order to make sure that our results are not driven by differences in the age structure of the female population across settlements. Columns 2 through 6 in all tables test each of the mechanisms mentioned in the Introduction. In column 2, we include mothers' literacy and income, which are our closest proxies for the opportunity cost of female time, husband's income, as more bargaining power to men may translate in a larger family, household wealth, because richer household may afford more children, and the main source of income (either from social benefits, informal, or formal employment).Footnote ²⁶ In column 3 instead, we add expected returns to education, which may influence the way parents allocate resources to quality vs. quantity of children.Footnote ²⁷ In column 4, we use the number of adults in the household, to capture the fact that grandparents may help in taking care of larger cohorts of children. In column 5, we control for whether the settlement is in an urban or suburban area as a proxy for the cost of space. In column 6, we include cultural variables such as religion and whether parents' names are of Serbian origin.Footnote ²⁸ Indeed families who declare being Christian Orthodox or whose first names sound typically Serbian might have been more influenced by Serbian cultural and social norms, among which that of having a small number of children.Footnote ²⁹ Finally in column 7, we keep controls for all mechanisms and test their robustness.

Robust standard errors are clustered at the settlement level with Moulton confidence intervals in the case of linear regressions [Imbens and Kolesár (Reference Imbens and Kolesár2016)]. Results of specification 1 are presented in Table 2.

Table 2. Number of children—OLS

Robust standard errors clustered at the settlement level with Moulton confidence intervals in parentheses: *p < 0.10, **p < 0.05, ***p < 0.01.

The first striking result is that fertility differences documented in the descriptive statistics persist once controlled for individual, household, and settlement characteristics. In particular, women living in mostly Roma settlements seem to have about between a third and a quarter fewer children than those in only Roma locations, but the coefficient is not always precisely estimated. More importantly, women in few Roma settlements have around 0.5–0.8 fewer children than similar women in only Roma locations, and that difference is significant at the 1% level in all specifications. Such a large fertility gap across settlement is surprising as differences in individual characteristics did not seem particularly large (see Table 1) but also because Roma are usually regarded as a culturally homogeneous group.

The second observation is that adding controls in columns 2–6 tends to reduce the fertility gap, while including all controls reduces it by even more, as shown in column 7. Among significant controls, mothers' income and household wealth are strongly negatively associated with fertility, suggesting that the poorest households and households where mothers have a lower opportunity cost of time tend to have larger families. These results confirm previous findings, for instance by Jensen (Reference Jensen2012), who finds that improved labor market opportunities not only substantially increase women's investment in education but as well delay marriage and fertility. The source of income instead does not seem to matter, while husband's income is positively related to a larger family. This suggests that a higher fertility is associated with a more pronounced specialization in the household.

Similarly, number of adults in the household and urban status, which are included to capture, respectively, different types of family arrangements and different housing prices, barely affect fertility. Higher-perceived returns to education instead are strongly associated with a lower fertility, which illustrates the presence of a quality/quantity trade-off. A one standard deviation change in perceived returns to education is associated with a fall in fertility of 0.2 child. Battaglia and Lebedinski (Reference Battaglia and Lebedinski2017) use the same data to show how enhancing people's perceived returns to education can impact positively labor market prospects and educational attainment.

Finally, religious affiliation does not seem to matter much, whereas exposure to the Serbian culture, as measured by whether parents hold Serbian sounding first names, is a strong predictor of a lower fertility. Indeed, shifting the prevalence of Serbian names from that in few Roma settlements down to the one observed in only Roma locations is associated with an increase in fertility of about 0.14 child.

When all controls are included together, the gaps between mostly Roma and only Roma and between only Roma and few Roma decrease in size but remain highly statistically significant. Most of the channels that appear significant taken separately remain significant, except for Serbian names. However, the fact that it is not enough to close the gap indicates that several of these factors are probably confounded and that we do not estimate causal effects.

In Table D in the Appendix, we focus on the subsample of women who reside in their settlement of birth. The rationale for this exercise is to avoid that what we claim is associated with segregation is actually due to recent migration patterns. Mothers who were born in the settlement in which they reside tend to be richer and to have significantly fewer children on average. Examination of the table however shows that the exact same pattern emerges regarding the comparison between only, mostly, or few Roma settlements. We can therefore be reassured that recent migration waves were not responsible for the higher fertility observed in more-segregated areas.

We finally investigate the robustness of these associations using alternative measures of segregation in Table E in the Appendix. Results remain qualitatively consistent using the median perception about whether the neighborhood is only, mostly, or few Roma at the street level rather than at the settlement level.

3.2 Birth timing and son preference

In the previous subsection, we have examined the number of children women have controlling for a second-order polynomial in their age. It implicitly imposes a structured, though quite flexible, relationship between age and the number of children that is common to all women in the sample. Our conclusions are therefore valid for completed fertility if birth timing does not differ significantly across settlement types.

In this subsection instead, we do away with this assumption by analyzing the pace at which women give birth instead of their total number of children only. In addition to confirming the conclusions drawn from the analysis of total number of children, examining birth timing allows us to investigate whether there are gender-specific patterns. Indeed, in a context where sex-selective abortion is not significantly used, looking at the sex ratio at birth is not informative. Instead, son preference can be detected by looking either at the difference in waiting time after having a boy vs. having a girl or at the sex ratio of last birth as originally suggested by Dalla Zuanna and Leone (Reference Dalla Zuanna and Leone2001).

To this end, we start by using a proportional hazard model, which leaves the baseline hazard rate unspecified and assumes that it is shifted multiplicatively by covariates.Footnote ³⁰ We look separately at first births and then all subsequent births. In the first case, we define the at-risk period as starting at age 14, while we make it start 9 months after the previous birth we observe in the latter case. We acknowledge a potential limitation for first births as we observe only children still in the household at the time of the survey, which may differ from the universe of all births. Indeed some of the oldest mothers in the sample may have had an early child who already left the sample and whom we consequently do not observe. Imbalances in the sex ratio of first born children confirm our hypothesis, as reported in Figure B in the Appendix. The at-risk period ends either with a birth or with the woman actually leaving the sample (i.e., when the running variable reaches her age at the time of the survey). Notice that now the unit of observation is the post birth spell and not the mother, hence the increase in the number of observations. Covariates are assumed to affect the baseline hazard multiplicatively, so that the equation we estimate is the following:

(2)$$h( {t\vert x_{\rm c}} ) = h_0( t ) e^{x_{\rm c}\beta _{\rm x}}$$

where x _c includes the same controls as in the previous subsection. Robust standard errors are clustered at the settlement level with Moulton confidence intervals.

The baseline results shown in Table 3 shed some further light on the findings from the previous subsection. Panel A indicates that there are no significant differences in the timing of first births, while panel B adds that women living in few Roma settlements have subsequent children at a significantly slower pace than their mostly or only Roma counterparts. The coefficient −0.449 on few Roma in column 1 corresponds to women in these neighborhoods being 36% less likely to have an extra child than comparable women living in an only Roma settlement. Controlling for potential mechanisms reduces but does not close the gap between settlements with different levels of segregation. As before, we investigate the robustness of the analysis by using an alternative measure of segregation at the street level. Results remain consistent, as reported in Table F in the Appendix.

Table 3. Birth spacing by settlement type

Robust standard errors clustered at the settlement level with Moulton confidence intervals in parentheses: *p < 0.10, **p < 0.05, ***p < 0.01.

Then we turn to whether birth timing differs according to the gender of the previously born child. To this end, we interact our measure of segregation with a dummy indicating whether the previously born child was a male or a female.Footnote ³¹ To make sure that children leaving the household early on do not introduce too much error in the measurement of the gender of the previously born or first born child, we choose to restrict the sample to mothers below 33 years of age.Footnote ³² Results are presented in Table 4.

Table 4. Birth spacing by settlement type and gender of previously born child—sample of women aged 33 or less

Robust standard errors clustered at the settlement level with Moulton confidence intervals in parentheses: *p < 0.10, **p < 0.05, ***p < 0.01.

The main result to highlight in this table is that the difference in birth spacing across settlement types comes mostly from the difference in birth spacing after a boy was born. Indeed, households in only Roma and mostly Roma do not exhibit significantly different spacing patterns after boys or girls. However, households living in a few Roma settlement tend to space substantially more after a boy, which is illustrated by the coefficient on “few Roma—male.” The coefficients in column 1 indicate that the hazard ratio of having an extra child after a male is 23% smaller than after a girl for households in only Roma settlements (but this difference is not statistically different from zero), while it is 58% smaller in the case of few Roma households. This latter difference is significantly different from zero at the 5% level. The same pattern emerges as controls are added: they allow us to close only partially the gap in birth timing across settlement types.

3.3 Comparing son preference across Roma and non-Roma communities in Serbia

In this subsection, we attempt to give further evidence that a cultural transmission process is at work. We use data on Roma and non-Roma populations from the MICS (2010) to analyze how their birth spacing patterns depend on the gender of the previously born child. We illustrate those patterns comparing the Kaplan–Meier survival curves for all births following a male baby to those following a female baby (which requires that we exclude first births) for both Roma and non-Roma. Results are shown in Figure 2. It appears that Roma and non-Roma Serbs exhibit a very similar pattern of shorter spells after a female baby. Table G in the Appendix confirms that non-Roma populations do exhibit a low fertility—1.73 children per woman on average—and a longer spacing after boys than after girls. While this in theory could be driven by common determinants across communities, we believe it points at a phenomenon of cultural assimilation from the Roma minority into the Serbian majority. Indeed, our data reveal that these specific gender patterns are dominantly prevalent in Roma settlements where they are a minority surrounded by non-Roma populations, as shown in Figure 3.

Figure 2. Kaplan–Meier survival curves by gender of the previously born child using MICS (2010) data.

Figure 3. Kaplan–Meier survival curves by gender of the previously born child.

While we are not aware of other studies documenting higher birth spacing after boys in Serbia or in other countries of the Western Balkans, high sex ratios at birth have been recently observed in several countries of the region [UNFPA (2012)].Footnote ³³ According to Guilmoto and Duthé (Reference Guilmoto and Duthé2013), together with the fertility decline and the development of modern healthcare services, the persistence of traditional patriarchal values is central to the son preference observed in these countries. Conversely, traditional Roma societies, especially in rural areas, show a female-biased sex ratio at birth and invest more heavily in daughters since they are more likely than sons to help their parents in taking care of siblings [Bereczkei and Dunbar (Reference Bereczkei and Dunbar1997, Reference Bereczkei and Dunbar2002)].

As a final check that the birth spacing patterns observed do reveal a higher son preference in less-segregated Roma settlements as well as in non-Roma communities, we look at the sex ratio at last birth. Indeed, as originally suggested by Dalla Zuanna and Leone (Reference Dalla Zuanna and Leone2001), the sex ratio at last birth captures gender-biased preferences in a context where couples decide to halt fertility after the birth of a son. To focus on last births, we include women who have at least one child and have not had children in the past 5 years. The sex ratios at last birth are shown in Table H in the Appendix for our data on the left panel and for Roma and non-Roma in the MICS on the right panel.

The sex ratio at last birth appears to be significantly higher than 0.5 for Roma both in our data, at 0.596, and in the MICS, at 0.564. Our data reveal that this skewness is largely driven by a very unbalanced ratio in few Roma settlements at 0.778. The sex ratio at last birth for non-Roma in the MICS is much lower than that, at 0.528, but statistically different from 0.5 and measured on a much larger sample. Ideally, we would like to measure the sex ratio at last birth of non-Roma populations that are in contact with Roma settlements. Unfortunately, this is not possible using MICS data. Ultimately, the sex ratio at last birth is an additional element pointing at a cultural transmission channel from the Serbian majority to the least-segregated Roma communities.

The sex ratio at last birth points at the importance of differential stopping behavior. For comparison purposes, we reproduce Table 4 on a restricted sample limited to all censored inter-birth spells and censored post-birth spells of less than 5 years. The idea is to check whether inter-birth spacing behavior (excluding therefore differences in stopping behavior) also differs across settlement types. Results are shown in Table I in the Appendix. Although it appears that few Roma communities consistently space more after boys than girls, the picture is less clear-cut than in the overall sample, as mostly Roma do so as well to some extent. Additionally, few Roma also space comparatively more after girls than mostly and only Roma. Altogether, these results show that both spacing and stopping happen at a different rate after a son and after a daughter, but differential stopping behavior is key to illustrate the higher son preference in few Roma communities.