1 Introduction

During the period 1950 to 1980, many developing countries had high illiteracy rates, high fertility rates, and substantial differences in fertility and children’s school enrollment ratios between educated and illiterate parents. Conventional views hold that, when illiterate parents have higher fertility than literate parents, this fertility differential reduces average education for the next generation and impedes the demographic transition and economic growth. However, Brazil achieved faster increases in school enrollment ratios and steeper declines in average fertility than Indonesia even though illiterate parents’ fertility was 2.5 times as high as educated parents’ fertility in Brazil but below educated parents’ fertility in Indonesia in 1960–1980. In developed countries with a high human development index (HDI), fertility increases with parental earnings, as in Myrskylä et al. (2009) and Bar et al. (2018). This paper explores fertility rates and school enrollment ratios between educated and illiterate parents (mothers) in a single gender model as well as the role of education subsidies in promoting education and reducing fertility for development.

To motivate the study, this paper first presents some relevant facts in two large developing countries, Brazil and Indonesia, using multiple rounds of census data in 1960–2010 for fertility rates and school enrollment ratios among families (mothers). Previous empirical studies on fertility and school enrollment ratios are mostly based on cross sectional variation.Footnote 1 Using census data with a large sample size and a long sample period helps us observe the puzzling evolution of the fertility distribution and changes in education status across generations. The data suggest a positive relationship between school enrollment ratios and parental education and overall positive trends of school enrollment ratios in all parental groups. The data also show that more than one-third of households enrolled only some rather than all children in school in the 1970s.

The decline in parental illiteracy rates played an important role in reducing average fertility in Brazil because illiterate parents had higher fertility than educated parents. In Indonesia, however, average fertility increased over 1960 to 1980 despite falling illiteracy rates because illiterate parents had lower fertility than educated parents. After the 1980s, educated parents’ fertility fell below illiterate parents’ fertility in Indonesia and average fertility started to decline.

Voluminous studies have explored fertility and child education and their relationships with parental earnings. From the Malthus theory (e.g., Malthus 1872), fertility increases with parental earnings due to subsistence costs for children (e.g., Boyer 1989; de la Croix et al. 2019).Footnote 2 Conversely, from the quantity-quality theory (Becker 1960), fertility decreases with parental earnings/education due to the opportunity costs of time spent on child rearing (e.g., Kremer and Chen 2002; de la Croix and Doepke 2003, 2004; Fan and Zhang 2013; Yasui 2017; Bar et al. 2018; Raute 2019).Footnote 3 Using US data, Bar et al. (2018) show a typical negative correlation between fertility and parental income in 1980 but an N-shaped relationship between fertility and parental income in 2010. They attribute the eventual positive relationship between fertility and parental income for rich households to the marketization of parental time costs at home.

In de la Croix and Doepke (2003), with unequal initial parental human capital, poor families have more children and invest little in education, and this fertility-differential effect accounts for most of the empirical relationship between inequality and growth. When inequality is high, public education yields higher growth by reducing fertility differentials (de la Croix and Doepke 2004). Using micro-data, Vogl (2016) finds that the relationship between fertility and parental socioeconomic status flipped from positive at the early stage of the demographic transition to negative. The reason for the flip in Vogl (2016) is a corner solution in which low-skill parents forgo investment in children as in a dynastic model of Becker et al. (1990). However, even illiterate parents enroll a fraction of their children in school in the Indonesia census data prior to the flip of fertility differential.

In these models, rearing children costs parents time and educating them costs parents money. In addition, Vogl (2016) assumes a living cost for a child and a subsistence constraint on consumption that creates a threshold level for human capital above which parents cease to be subsistence-constrained when making decisions on fertility. A common assumption in previous models is that a parent chooses equal or no education investment in children. This model allows parents constrained by children’s education and living costs to choose the number of children and the fraction of children who receive education, with the other children working.Footnote 4 The earnings of working children weaken the incentive for education and strengthen the incentive for more children, particularly in poor families. This model also assumes a progressive labor-income tax to finance education subsidies that may be higher for children in poor households. The taxes reduce the opportunity cost of time rearing children and the subsidies reduce the private cost of school enrollment. Thus, progressive income taxes and favorable education subsidies for poor households may cause higher fertility for educated parents than for illiterate parents. A partial catch-up rise in education subsidies for educated parents enlarges the decline in their fertility and thus may induce the flip of the fertility differential when illiterate parents also enroll some children in school.

We show that the implications of this mechanism can account for the patterns of enrollment ratios and fertility in relation to parental education attainments in the data and the roles of progressive labor-income taxes and education subsidies in the demographic transition and development. Our first contribution is to account for the positive relationship between school enrollment ratios and parental education or earnings in the census. When the opportunity costs of children’s education (e.g., education spending and child earnings) are high relative to the forgone earnings of rearing children and living costs, parents have as many children as possible and no incentive for children’s education. When income grows relative to the education costs, parents enroll a fraction of children in school. Thus, identical children at birth in a family become educated or illiterate workers. The enrollment ratio increases with the returns to education, forgone earnings of rearing children, and living costs but decreases with the education costs and child earnings.

Our second contribution is to account for the N-shaped relation between fertility and parental earnings in the data. When children’s living costs are high relative to children’s education costs, fertility increases with parental earnings due to a dominant income effect and thus fertility is higher for educated parents than for illiterate parents. When the living costs are low relative to the education costs, fertility decreases with parental earnings due to a dominant substitution effect and thus fertility is lower for educated parents than for illiterate parents. A restriction on child labor can reduce fertility and increase enrollment ratios.

Third, the illiteracy rate of an economy converges to a steady state—a poverty trap for illiterate parents when child earnings and the costs of education and living increase proportionately with parental earnings for stationary education subsidy rates. This motivates a rise in education subsidy rates in the attempt to reduce illiteracy for prosperity.Footnote 5 The decline in illiteracy rates may raise average fertility in the economy when illiterate parents have lower fertility than educated parents as in Indonesia until 1990. Conversely, the rise in education subsidy rates is particularly effective in reducing illiteracy rates and average fertility when illiterate parents have higher fertility than educated parents.

The results help explain the Brazil-Indonesia comparison with their, respectively, high and low education subsidies. In particular, rising education subsidies in favor of poor families financed by progressive labor-income taxes can initially generate higher enrollment ratios and fertility for rich families than for poor families at low average income. Then, a partial catch-up rise in education subsidy rates for other families can enlarge the decline in their fertility and thus yield the flip of the fertility differential as observed in Indonesia. When all children complete education, fertility may relate positively to parental earnings again as in developed countries with high HDI (e.g., Myrskylä et al. 2009; Bar et al. 2018; Vogl 2020).Footnote 6

The rest of this paper proceeds as follows. Section 2 provides facts in the census data. Section 3 proposes the model. Section 4 presents equilibrium analyses. Section 5 provides quantitative results. Section 6 gives the conclusion.

2 Facts in census data from Brazil and Indonesia

We focus on census data from Brazil and Indonesia to motivate our paper for the following reasons. First, both countries have the most extended sample coverage in developing countries, which helps observe declining fertility rates and rising school enrollment ratios.Footnote 7 Second, Brazil has one of the highest ratios of illiterate parents’ fertility to educated parents’ fertility in the world but Indonesia has one of the lowest ratios as noted by Kremer and Chen (2002). Thus, a comparison between these two countries may shed light on the distribution and evolution of fertility and school enrollment ratios across parental groups and over time.

We extract data from the Integrated Public Use Microdata Series (IPUMS). We use actual fertility at the survey time and restrict our sample to focus on parents (women) aged 40 to 49 in the census year because the number of children born over parents’ lifetime can only be accurately constructed for those at the end of their reproductive cycle.Footnote 8 As a result, we observe the fertility and child school enrollment of six birth cohorts in Brazil and seven in Indonesia.Footnote 9 Since one can only match parents with their children if they lived in the same household, we focus on the education outcomes of school-aged children between 7 and 15 to avoid over-representing children’s sample by children who left home at older ages.Footnote 10 We measure child education by school enrollment status because many children in this age group have not completed their education yet.

Table 1 reports the descriptive statistics of the key variables. To check whether changes in school enrollment ratios differ between genders, we report school enrollment ratios for boys and girls separately. Both countries’ demographic transition in the sample period had similar average fertility rates and comparable fertility trends. For example, on average, the sampled parents born in the 1920s had 5.5 children in Brazil and 5.2 children in Indonesia. The respective numbers declined to 2.4 and 3.0 for parents born in the 1960s.

Table 1 Descriptive statistics

Parental education attainments and school enrollment ratios of both boys and girls increased considerably over time. Real GDP per capita in Brazil was more than twice as that in Indonesia and increased substantially in both countries. The 10-year average annual growth rates peaked first at 6% in 1970–1980 and again at 2% and 3.7%, respectively, in 2000–2010 in Brazil and Indonesia, following surging school enrollment ratios in the preceding decades. Public education spending was below 1% of national income (GNI) in Indonesia until 2000, but 4% in Brazil throughout the sample period. At the end of the sample period, public education spending as a share of GNI in Indonesia was comparable to that in Brazil.

To check whether parents enroll a fraction of children in schools, we focus on children aged 7–15 whose school enrollment is still mostly their parents’ decision. Figure 1(a) plots the share of households who only enroll some of their children in schools, and Fig. 1(b) plots the share of households who enroll all their children in schools. The figures show that more than 30% of households only enroll some of their children to schools in the 1970s in both Indonesia and Brazil.

Fig. 1
figure 1

Share of households by children’s enrollment status. Data source: Various years of Brazil and Indonesia census

Figure 1(a) shows that this ratio increased in Indonesia from the 1970s to 1980s and then declined monotonically, suggesting a possibility that the initial increase in school enrollment is because more parents send at least one child in schools. Even at the end of the sample period, a considerable number of households still did not enroll all of their children in schools, particularly in Indonesia. Figure 1(b) shows that while the proportion of households who enrolled all of their children is lower in Brazil than in Indonesia before 1990, it is higher in Brazil afterwards. In 2010, about 93% of Brazilian households enrolled all of their children to schools, which is about 10 percentage points higher than that in Indonesia.

The figure clearly shows that allowing parents to enroll only some of their children could help us understand the demographic transition process. However, the enrollment data in the figure consist of all education levels. When focusing on secondary education that matters to parental choices between schooling and working for children, school enrollment ratios in developing countries would be much lower than the ratios in Fig. 1 (World Bank data, various years). Consequently, the proportion of families that enroll only a fraction of children in secondary school would be much higher than the ratios in the figure.

In Table 2, we further divide our sample by the older sibling’s enrollment status. Column (1) reports the unconditional school enrollment ratio of the youngest child, \(P(S_y=1)\), where \(S_y=1\) if the youngest child is enrolled in school and 0 otherwise. This ratio \(P(S_y=1)\) increased from 50.3% in 1960 to 97.14% in 2010 for Brazilian boys and from 50.2% to 97.7% for Brazilian girls. The increase in school enrollment ratio in Indonesia is in a similar order, from 66.6% in 1971 to 94.1% in 2010 for boys and from 65.6% to 84.8% for girls. Column (2) reports the conditional school enrollment ratio \(P(S_y=1|S_o=1)\), where \(S_o=1\) if the oldest child had been enrolled in school and 0 otherwise. This ratio can only serve as an upper bound because younger children from poor families have more years of schooling and higher school enrollment ratios than their older siblings, as noted by Parish and Willis (1993). For comparison, we also report \(P(S_y=1|S_o=0)\). The steady upward trend in \(P(S_y=1|S_o=0)\) suggests an increasing probability of sending at least one child to school over time.

Table 2 Child school enrollment ratios by the older sibling’s enrollment status

Figure 2 plots school enrollment ratios of children aged 7–15 by parents’ education attainments in four groups: illiterate; incomplete primary education; completed primary education; and completed secondary education. It shows a positive relationship between school enrollment ratios and parental education attainments in both countries. School enrollment ratios were lower for children born to parents with some or completed primary education in Indonesia than in Brazil. School enrollment ratios of illiterate parent increased faster in Brazil than in Indonesia despite a much higher ratio of illiterate parents’ fertility to educated parents’ fertility in Brazil. This may arise from the much higher education subsidies in Brazil than in Indonesia from Table 1. Nevertheless, the data suggest that a higher ratio of illiterate parents’ fertility to educated parents’ fertility does not necessarily imply slower human capital accumulation during the demographic transition.

Fig. 2
figure 2

Child school enrollment ratios by maternal education. Data source: Various years of Brazil and Indonesia census

Table 3 The impact of maternal education on child school enrollment

It should be noted that the upward trend revealed in Fig. 2 might be biased by changes in the age composition of sampled children. Because parents’ age at childbirth might have increased with their education attainments over time, the average age of sampled children could differ systematically across census years. Given that many children dropped out of school after completing primary education, changes in the age profile of the sampled children could cause year to year variations in school enrollment ratios even if age-specific school enrollment ratios do not change. To control for the impact of changes in the age profile of sampled children on school enrollment ratios, we run the following logit regression

$$\begin{aligned} P(E_{it}=1)=F(\gamma _{0}+\gamma _{1}M_{i}+\sum _{a=8}^{a=15}\gamma _{a}D_{ai} + \sum _{j=1}^{j=4}\sum _{t}Z_{j}C_{t}\delta _{jt}), \end{aligned}$$
(1)

where \(E_{it}=1\) if child i is enrolled in school at time t; \(M_{i}=1\) if child i is a male; \(D_{a}\) is a series of age dummies; \(Z_{j}\) is a series of parents’ education attainment dummies; and \(C_{t}\) is a series of census year dummies.

Table 3 reports estimation results. Columns (1)–(4) are the coefficients on the interaction between the parental education dummy and the census year dummy, \(\hat{\delta }_{jt}\), whereas columns (5)–(8) are the marginal effects of these interaction terms. For the sake of brevity, the coefficients on the gender dummy and age dummies are not reported. We use children of illiterate parents in the 1960 Brazilian census (1971 Indonesian census) as the reference group. Consequently, the marginal effect of the interaction term captures the difference in school enrollment ratios between children in the reference group and children in a concerned census year whose parent belonged to the given education group, when evaluated at the sample mean. The vertical reading reveals the time trend of school enrollment ratios for a given level of parental education, while the horizontal reading shows the impact of parental education on school enrollment ratios at any given time.

Table 4 Population shares by maternal education attainment

Parental education has strong positive effects on school enrollment ratios. In Brazil, for instance, the school enrollment ratio of children from parents with secondary education was 54.7 percentage points higher than that of children from illiterate parents in 1960. The difference in school enrollment ratios of children across parental education attainments in Indonesia was smaller than that in Brazil. In 1980, for instance, the school enrollment ratio of children from parents with at least secondary education was 33.6 percentage points higher than that of children from illiterate parents in Indonesia and 59.7 percentage points in Brazil.

The estimation results also show that school enrollment ratios increased over time regardless of parental education even after controlling for children’s age. For children in the 1960 versus 2010 Brazilian censuses, school enrollment ratios increased by 47.1 percentage points for illiterate parents and by 11 percentage points for middle school graduates. As a result, the difference in child school enrollment ratios between illiterate parents and parents with secondary education narrowed from 54.7 percentage points in 1960 to 18.7 percentage points in 2010. Thus, the population share of less educated parents decreased, as shown in Table 4. The least educated parents’ population share decreased from 47.2% for the 1920s cohort to 0.9% for the 1960s cohort in Brazil and from 79.5% to 10.2% in Indonesia.

Table 5 presents fertility rates by parental education with downward trends. The ratio of illiterate parents’ fertility to that of parents who completed secondary education was 2.5 from the cohort born in 1911–1920 to the cohort born in 1951–1960 in Brazil but 0.9 in Indonesia. Such opposite fertility differentials may explain why average parental schooling years increased slightly slower in Brazil than in Indonesia. Yet, school enrollment ratios grew faster and fertility declined more sharply in Brazil than in Indonesia in contrast to the conventional view that economies with higher fertility from illiterate parents than from educated parents tend to have lower education for the next generation and higher average fertility. These challenging facts motivate our theoretical investigation.

Table 5 Fertility rates by maternal education

3 The model

Consider a quality-quantity model with infinite discrete time and two-period lived agents, children and parents. There are two types of parents in period t: illiterate (or less educated) parents with human capital \(h_t^L\) and literate parents with human capital \(h_t^H>h_t^L\). A child who is not enrolled in school has a lower level of human capital \(h^C_t<h_t^L\) than an illiterate parent does. Children are identical at birth within a family but may grow up as different workers with or without education.

The production of a single final good uses effective labor \(l_t^k h_t^k\) for \(k=H,L\):

$$\begin{aligned} y_t^k=l_t^k h_t^k. \end{aligned}$$
(2)

The wage rates of parents are equal to their own human capital \(W_t^k=h_t^k\) and thus \(W_t^H>W_t^L\). Similarly, the wage rates for children are \(W_t^C=h_t^C\).

The cost and level of education are externally structured to parents. The education cost, such as a tuition fee, is proportional to the human capital of literate adults, \(e_t=\xi _t {h}_t^H\) where \(\xi _t\in (0,1)\) as the bulk of the cost is the implicit cost of teachers. The adulthood human capital of children from schooling depends positively on education investment received in period t

$$\begin{aligned} h_{t+1}^H=A(\psi +B\xi _t)^\phi {h}_t^H; \end{aligned}$$
(3)

and the adulthood human capital of children without schooling is

$$\begin{aligned} h_{t+1}^L=A\psi ^\phi {h}_t^H, \end{aligned}$$
(4)

where \(A>0\) is the efficiency parameter; \(B>0\) is the relative advantage of schooling; \(\psi >0\) captures the spillovers from the human capital of literate parents to all children;Footnote 11 and \(0<\phi <1\) measures the degree of returns to human capital accumulation. When compulsory education at the primary level is implemented, it can be represented by \(\psi \) when children at this school age are too young to work; and accordingly parents choose secondary school enrollment for children.

From (3) and (4), the return to education \(\omega _{t+1}\equiv h_{t+1}^H/h_{t+1}^L\) is determined by

$$\begin{aligned} \omega _{t+1}=\frac{h_{t+1}^H}{h_{t+1}^L}=\left( \frac{\psi +B\xi _t}{\psi }\right) ^\phi \equiv \omega (\xi _t)>1. \end{aligned}$$
(5)

Since \(0<\phi <1\), \(\omega _{t+1}\) increases with \(\xi _t\), at a diminishing rate: \(\omega '(\xi _t)>0\) and \(\omega ''(\xi _t)<0\). It is straightforward to extend this case to more education levels.

A parent allocates one unit of time endowment to labor or childcare. As in Becker et al. (1990), rearing a child needs \(v\in (0,1)\) units of time. A type-k parent chooses the number of children, \( n_t^k\in [0, 1/v]\), and a fraction of them for education, \( \lambda _t^k\in [0, 1]\), taking the cost and structure of education as given. Thus, \(v n_t^k\) is the units of time for child rearing and the remaining time \(1-vn_t^k=l_t^k\) is labor. The average adulthood human capital of children in the family is equal to \(\bar{h}_{t+1}^k\equiv \lambda _t^k h_{t+1}^H+(1-\lambda _t^k)h_{t+1}^L.\)

Parents derive utility from their own consumption \(c_t^k\), the number of children \(n_t^k\), and the average human capital of children \(\bar{h}_{t+1}^k\) as follows:

$$\begin{aligned} U_t^k=\frac{\left( c_t^k\right) ^{1-\sigma }\!-\!1}{1\!-\!\sigma }+\beta \frac{\left[ \left( n_t^k\right) ^\alpha \left( \bar{h}_{t+1}^k\right) ^\eta \right] ^{1-\gamma }\!-\!1}{1\!-\!\gamma }, \end{aligned}$$
(6)

where \(1>\beta >0\) indicates the taste for utility from children; \(1>\alpha >0\) or \(1>\eta >0\) indicates the relative role of the number or average human capital of children; \(1/\sigma >0\) or \(1/\gamma >0\) indicates the elasticity of intergenerational substitution.

A type-k parent pays an income tax at rate \(\tau ^{k}_{t}\in [0,1)\) and receives an education subsidy at rate \(s_{t}^{k}\in [0,1)\). As seen in Fig. 2, a common problem in developing countries is low school enrollment in poor families which motivates various public programs to mitigate poverty and low school enrollment in poor families. For instance, Morley and Coady (2003) study targeted education subsidies in the form of cash or food conditional on keeping children in school in many developing countries including Brazil. Given progressive income taxes and targeted education subsidies in developing countries, we expect that \(\tau _{t}^{H}\ge \tau _{t}^{L}\) and \(s_{t}^{H}\le s_{t}^{L}\). From Table 1, the ratios of public education spending to GDP increased in both countries in 1970–2010 with declines only during times of rapid GDP growth, suggesting a positive trend in education subsidy rates for improvements in enrollment ratios.

Child earnings are supplemental income in parental budget constraints as in Doepke and Zilibotti (2005). The living cost per child exceeds a subsistence level \(d_t\ge \tilde{d}\). Both the education cost and child living cost are externally given to families. The budget constraint of a type-k parent is given by

$$\begin{aligned} c_t^k=(1-\tau _t^k)W_t^k(1-v n_t^k)+W_t^C(1-\lambda _t^k)n_t^k-d_t n_t^k-(1-s^k_t)e_t\lambda _t^k n_t^k. \end{aligned}$$
(7)

Since the budget-feasible set is non-convex due to the product of choice variables \(\lambda _t^k n_t^k\), we need the following restriction as in Ehrlich and Lui (1991):

Assumption 1

\(1>\alpha>\eta >0\).

This restriction means that the number of children is more important than the quality of children in parental preference. A parent maximizes utility in (6) by choosing the number of children \(n_t^k\in [0,1/v]\) and school enrollment ratio \(\lambda _t^k\in [0,1]\) subject to (7). Given education costs, if parents cannot afford education for all children, then enrolling a fraction of them in school may be more favorable than enrolling none. This partial enrollment ratio is consistent with the evidence in Parish and Willis (1993) and in Brazil and Indonesia. From the partial school enrollment ratio, inequality arises not only among children across families as in the literature (e.g., Córdoba et al. 2016) but also among identical siblings.

Overall, the assumption of enrolling a fraction of children in school is based on the low enrollment ratios and high illiteracy rates in 1960–1970 in the data. In contrast, the conventional assumption of equal investment in all children implies either full or zero illiteracy rates. In this model high illiteracy means persistent inequality and poverty for illiterate parents due to low school enrollment ratios, thus motivating relatively high education subsidy rates for low income families in developing countries as noted in Morley and Coady (2003).

The first-order condition with respect to the number of children is

$$\begin{aligned} \frac{\partial U_t^k}{\partial n_t^k}\!\!=\! & {} \!-(c_t^k)^{-\sigma }[(1-\tau _t^k)W_t^k v+d_t+(1-s_t^k)e_t\lambda _t^k-W_t^C(1-\lambda _t^k)]\nonumber \\{} & {} +\alpha \beta (n_t^k)^{\alpha (1-\gamma )-1}(\bar{h}_{t+1}^k)^{\eta (1-\gamma )}= 0\;\; \textrm{for}\; n_t^k\in (0,1/v); \nonumber \\{} & {} \frac{\partial U_t^k}{\partial n_t^k}<0\;\textrm{for}\;n_t^k=0; \frac{\partial U_t^k}{\partial n_t^k}>0\;\textrm{for}\;n_t^k=\frac{1}{v}. \end{aligned}$$
(8)

The first term in the net marginal benefit of having a child is the forgone marginal utility through the time cost \((1-\tau _t^k)W_t^k v\), living cost \(d_t\), and education cost \((1-s^k_t)e_t\lambda _t^k\), net of child earnings \(W_t^C(1-\lambda _t^k)\). The second term is the marginal utility gained from having an additional child. If child earnings are no less than the time and living costs, \(W_t^C\ge (1-\tau _t^k)W_t^k v+d_t\), then the net marginal benefit of having a child \(\partial U_t^k/\partial n_t^k\) will be positive at \(\lambda _t^k=0\). This implies that parents may choose maximum fertility \(n_t^k=1/v\) and no school enrollment for children. Since \(W_t^H>W_t^L\), it is possible that illiterate parents choose maximum fertility while literate parents choose lower fertility and enroll a fraction of children in school.

The first-order condition with respect to the school enrollment ratio is

$$\begin{aligned} \frac{\partial U_t^k}{\partial \lambda _t^k}=-&(c_t^k)^{-\sigma }n_t^k[(1-s_t^k)e_t+W_t^C ] \nonumber \\ +&\beta \eta (h_{t+1}^H-h_{t+1}^L) (n_t^k)^{\alpha (1-\gamma )}(\bar{h}_{t+1}^k)^{\eta (1-\gamma )-1}=0\;\;\textrm{for}\; \lambda _t^k\in (0,1);\nonumber \\&\frac{\partial U_t^k}{\partial \lambda _t^k}<0\;\;\textrm{for} \;\lambda _t^k=0;\;\;\frac{\partial U_t^k}{\partial \lambda _t^k}>0\;\;\textrm{for}\;\lambda _t^k=1. \end{aligned}$$
(9)

The first term in the net marginal benefit of enrolling a child in school is the forgone marginal utility of direct education cost \((1-s_t^k)e_t n_t^k\) as well as forgone child earnings \(n_t^k W_t^C \). The second term is the marginal utility from increasing average human capital through education \(h_{t+1}^H-h_{t+1}^L>0\).

Let \(N_{t}^{k}\) be the number of type \(k \ (k=H,L)\) parents in period t. The government runs a balanced budget:

$$\begin{aligned} N_t^H(1-v n_t^H)W_t^H\tau _t^H+N_t^L(1-v n_t^L)W_t^L\tau _t^L=N_t^H\lambda _t^H e_t s_t^H +N_t^L\lambda _t^L e_t s_t^L+G_t, \end{aligned}$$
(10)

where \(G_t\) is government spending exclusive of education subsidies.

4 Equilibrium analysis

In equilibrium, \(e_t=\xi _t h_t^H\), \(W_t^H=h_t^H\), \(W_t^L=h_t^L\), and \(W_t^C=h_t^C\) where earnings are equal to human capital. Substituting a given state \((h_t^H, h_t^L,h_t^C, d_t,\xi _t)\) and \(\omega (\xi _t)=h_{t+1}^H/h_{t+1}^L\) into (8) and (9), the unique interior solution to these two equations, denoted \(\Lambda _t^k\), is as follows:

$$\begin{aligned} \Lambda _t^k\!\equiv \!\frac{\eta (\omega (\xi _t)\!-\!1)[(1\!-\!\tau _t^k)h_t^k v\!+\!d_t\!-\!h_t^C]\!-\!\alpha [(1\!-\!s_t^k)\xi _t h_t^H+\!h_t^C]}{(\alpha -\eta )(\omega (\xi _t)-1)[(1-s_t^k)\xi _t h^H_t+h_t^C]}. \end{aligned}$$

Since \(\omega (\xi _t>1\), the necessary condition for \(\Lambda _t^k>0\) is \((1\!-\!\tau _t^k)h_t^k v\!+\!d_t\!>\! h_t^C\). Combining it with \(0\le \lambda _t^k\le 1\), the equilibrium school enrollment ratio for type-k parents is

$$\begin{aligned} \lambda _t^k\!=\!\left\{ \begin{array}{ll}\min \left\{ 1,\max \left\{ 0,\Lambda ^k_t\right\} \right\} \;&{}\text {if}\ (1\!-\!\tau _t^k)h_t^k v\!+\!d_t\!>\! h_t^C,\\ 0\;\;&{}\text {otherwise.}\end{array}\right. \end{aligned}$$
(11)

In (11), parents enroll a fraction of children in school only if forgone parental earnings for child rearing plus the living cost per child are greater than child earnings, \( (1-\tau _t^k)h_t^k v+d_t>h_t^C\). Otherwise, parents would not enroll any child in school and would choose the maximum number of children \(n_t^k=1/v\) in (8). Notably, the rate of return to education \(\omega (\xi _t)-1=({h_{t+1}^H-h_{t+1}^L})/{h_{t+1}^L}=({W_{t+1}^H-W_{t+1}^L})/{W_{t+1}^L}\) must be high enough, relative to the opportunity costs of education \((1\!-\!s_t^k)\xi _t h_t^H+\!h_t^C\), to motivate parents to enroll a fraction of children in school.

Figure 3(a) illustrates the relationships between school enrollment ratios and parental earnings under the conditions for school enrollment by numerically solving (11). In doing so, we set the living cost as 7.7% of illiterate parents’ earnings in every period, the education cost as 3% of educated parents’ earnings in every period, and child earnings as 4% of illiterate parents’ earnings in the initial period. At a first glance, we also set an equal subsidy 56% on education spending for both types of parents financed by a 20% tax on educated parents’ earnings and a 3% tax on illiterate parents’ earnings in all periods.Footnote 12 We will discuss unequal education subsidies in favor of children of illiterate parents and the parameter values in the calibration for the simulation in Section 5. For illustration, we increase parental earnings by 10% in each period. In Fig. 3, the enrollment ratio for children of educated parents starts low in period 0 and grows to 100% in period 4. By contrast, the enrollment ratio for children of illiterate parents remains zero until period 4, and reaches 100% in period 8. Formally, we link school enrollment ratios to parental earnings and other factors below.

Fig. 3
figure 3

School enrollment ratios and fertility vs parental earnings. Parameters: \(\beta =3\), \(\gamma =1\), \(s^H=s^L=0.56\), \(\tau ^H=0.2\), \(\tau ^L=0.03\), \(H_0^L=1\), \(H^H/H^L=1.6803\), \(d=0.077 H^L\), \(H^C=0.04 H_0^L\). \(H^L\) grows by 10% per period. Other parameter values are in Table 6

Proposition 1

If \(\eta (\omega (\xi _t)-1)[(1-\tau _t^k)h_t^k v+d_t]>\alpha (1-s_t^k)\xi _t h_t^H+[\alpha +\eta (\omega (\xi _t)-1)]h_t^C\), then school enrollment ratios are positive, increasing with returns to education \(\omega (\xi _t)\), forgone earnings \(h_t^k v\) for rearing a child, and child living costs \(d_t\), but decreasing with education costs \((1-s_t^k)\xi _t H_t^H\), income taxes \(\tau _t^k\), and child earnings \(h_t^C\). If the growth rates of parental earnings, child earnings, and the costs of living and education are equal, then school enrollment ratios are constant over time for stationary rates of taxes and subsidies.

Proof

See Appendix A. \(\square \)

These results differ from those in the related literature with equal or no education for children. Subsidizing education increases school enrollment ratios by decreasing direct education costs, thereby raising upward intergenerational mobility for all families. If tax and subsidy rates are equal, then school enrollment ratios increase with parental education or earnings because \(v(h_t^H-h_t^L)>0\). Thus, literate parents have higher school enrollment ratios than illiterate parents for \(\lambda _t^H\in (0, 1)\), due to the difference in their forgone earnings of rearing a child. The positive relationship between school enrollment ratios and maternal education is observed in Fig. 2.Footnote 13

The positive impact of rising child living costs on school enrollment ratios differs from those in the literature such as Becker et al. (1990). The intuition will become clear when we examine the impacts of child living costs on fertility. When child living costs increase proportionately with income, economies with higher income may not have higher school enrollment ratios unless they have higher education subsidies or lower child earnings than others. Unlike education subsidies, child earnings decrease enrollment ratios by increasing the opportunity cost of schooling. When child earnings increase proportionately with parental earnings, their opposite effects on enrollment ratios cancel out. When child earnings remain constant as in Fig. 3, the increase in parental earnings drives up enrollment ratios.

Equations (7), (8) and (11) yield an explicit solution for fertility in a special case with logarithmic utility \(\sigma =\gamma =1\):

$$\begin{aligned} n_t^k\!=\!\left\{ \begin{array}{ll}\frac{1}{v}\;\;\;\!&{}\textrm{if}\;(1\!-\!\tau _t^k)h_t^k v\!+\!d_t\!\le \!h_t^C; \\ \;\;\;\;\;\;\;\;&{}\mathrm{else\; if} \;\sigma =\gamma =1:\\ \frac{\alpha \beta (1-\tau _t^k)h_t^k}{(1+\alpha \beta )[(1-\tau _t^k)h_t^k v+d_t-h_t^C]},&{}\forall \; \lambda _t^k= 0;\\ \frac{(\alpha -\eta )\beta (1-\tau _t^k) h_t^k(\omega (\xi _t)-1)}{(1\!+\!\alpha \beta )\{[(1\!-\!\tau _t^k)h_t^k v \!+\!d_t\!-\!h_t^C](\omega (\xi _t)\!-\!1)\!-\!(1\!-\!s_t^k)\xi _t h_t^H\!-\!h_t^C\}},&{} \forall \lambda _t^k\!\in \!(0,1);\\ \frac{\alpha \beta (1-\tau _t^k) h_t^k}{(1+\alpha \beta )[(1-\tau _t^k)h_t^k v+d_t+(1-s_t^k)\xi _t h_t^H]}, &{}\forall \; \lambda _t^k=1. \end{array}\right. \end{aligned}$$
(12)

Figure 3(b) depicts the relationships between fertility rates and parental earnings in (12) when the living cost grows proportionately with parental earnings. At low income levels, educated parents initially have higher fertility than illiterate parents. As parental earnings and the living cost grow, fertility rates decline. The decline in fertility is relatively slow at zero school enrollment for children of illiterate parents. Thus, at higher levels of parental earnings, the fertility rate of educated parents falls below the fertility rate of illiterate parents. When the enrollment ratio reaches 100% for children of educated parents, the fertility rate of educated parents becomes constant, while the fertility rate of illiterate parents continues to fall as they increase the enrollment ratio for their children. Eventually, when the enrollment ratio reaches 100% for all types of children, the fertility rate is higher again for educated parented than for illiterate parents. Formally, we link fertility to parental earnings and other factors as follows:

Proposition 2

Under \((1-\tau _t^k)h_t^k v+d_t> h_t^C\) with \(\sigma =\gamma =1\), fertility increases with parental earnings if \(d_t>h_t^C\) at \(\lambda _t^k= 0\), or \(d_t>[(1-s_t^k)e_t+h_t^C\omega (\xi _t)]/(\omega (\xi _t)-1)\) at \(\lambda \in (0,1)\), or \(\lambda =1\); however, fertility decreases with parental earnings if \(d_t<[(1-s_t^k)e_t+h_t^C\omega (\xi _t)]/(\omega (\xi _t)-1)\) at \(\lambda \in (0,1)\). Fertility also increases with child earnings at \(\lambda \in [0,1)\), or with education costs at \(\lambda \in (0,1)\), but decreases with child living costs at \(\lambda \in [0,1]\), returns to education at \(\lambda \in (0,1)\), and education costs at \(\lambda =1\). If the growth rates of parental earnings, child earnings, and the costs of living and education are equal, then fertility is constant for stationary rates of taxes and subsidies.

Proof

See Appendix B. \(\square \)

The conditions signing the effects of parental earnings on fertility are complementary to those in the literature. In the extreme case when parental education is so low that school enrollment ratios are zero (\(\lambda ^k= 0\)), fertility increases with parental earnings if child living costs are greater than child earnings as the (Malthusian) income effect dominates the (Beckerian) substitution effect. A rise in child living costs strengthens the income effect, whereas a rise in child earnings weakens it. When these two factors offset \(d_t=h_t^C\), logarithmic utility would imply that the income and substitution effects of parental earnings should cancel out.

When parental earnings are high enough such that enrollment ratios are positive and below 100%, the impact of parental earnings on fertility is positive if child living costs are high relative to education costs and child earnings, \(d_t> [ (1-s_t^k)e_t+h_t^C\omega (\xi _t)]/(\omega (\xi _t)-1)\), and negative otherwise. Thus, this model can generate a hump-shaped relationship between fertility and parental skills or education with positive investment in child education in contrast to the corner solution in Becker et al. (1990) and Vogl (2016) at low human capital or skill.Footnote 14 A rise in education costs tips the quality-quantity tradeoff in favor of the latter, whereas a rise in returns to education tips the tradeoff against the latter. Subsidizing education reduces the education cost for \(\lambda ^k\in (0,1)\), thus reducing fertility.

When parental education is so high that school enrollment ratios are 100% (\(\lambda ^k=1\)), children do not work and thus fertility increases with parental earnings as opposed to the Beckerian theory since now the income effect dominates as in Bar et al. (2018). This N-shaped relationship between fertility and parental education or earnings extends the hump-shaped relation in previous models. When all children are enrolled in school, a rise in education costs reduces fertility and thus subsidizing education may raise fertility by reducing education costs.

In all cases, fertility decreases with child living costs, explaining why school enrollment ratios increase with child living costs. If the growth rates of parental earnings, child earnings, and the costs of living and education are equal, then their effects on fertility completely cancel out with logarithmic utility and thus fertility remains constant. In Fig. 3, due to the constant child earnings, fertility rates decline when the living and the direct education costs grow proportionately with parental earnings. This is because the constant child earnings strength the negative substitution effect of rising parental earnings on fertility.

When the rates of taxes and subsidies are equal for all parents, fertility is higher for literate parents than for illiterate parents under the conditions for a positive relationship between fertility and parental earnings. If these conditions are reversed, fertility is lower for parents with higher education attainments. The mixed results for fertility differentials are consistent with the mixed observations of fertility differentials in Brazil and Indonesia.

The predictions for partial school enrollment ratios are comparable to those observed in Brazil and Indonesia in Tables 2 and 3. In these two countries, the majority of parents were illiterate or had some incomplete primary education in census years in 1960–1990, as shown in Table 4. From Table 5, fertility rates decreased with parental education in Brazil, but increased with parental education in Indonesia except for those who completed secondary education in census years in 1970–1990. The contrasting patterns can be explained by Proposition 2 with \(\lambda ^k\in (0,1)\) if parental earnings were much higher in Brazil than in Indonesia relative to subsistence costs for children in these years.

This argument may be relevant as real GDP per capita adjusted by purchasing power parity in Brazil was more than three times as that in Indonesia in these years in Table 1. The low income in Indonesia means that subsistence costs were high relative to income, especially in the 1960s. When income per capita in Indonesia reached the 1960s level of Brazil after the 1990s, the majority of sampled parents in Indonesia completed primary education and fertility decreased with parental education. As a result, the fertility differential in Indonesia flipped in Table 5 when reaching the 1960 income per capita in Brazil (2,469 dollars at 2005 price).

Propositions 1 and 2 help explain the demographic transition, the change in school enrollment ratios, as well as the evidence on the flip of fertility differentials in Vogl (2016) across nations. In the early development stage when parental earnings are low relative to child living costs, illiterate parents have fewer children and lower school enrollment ratios for children than literate parents. Rising parental earnings during early development may induce higher fertility but little change in school enrollment ratios, unless governments subsidize education. As income grows, parental earnings and education costs are high relative to child living costs, and thus fertility is higher for illiterate parents than for literate parents. Now, a rise in parental earnings decreases fertility and increases school enrollment ratios.

Eventually, when most children complete all levels of education, fertility increases with parental earnings again in contrast to a downward trend to low fertility in models with equal investment in siblings. According to Myrskylä et al. (2009), developed countries with the human development index (HDI) value exceeding 0.9 have observed increases in fertility along with increases in HDI in the previous decade. These countries have very high tertiary education enrollment ratios, e.g., over 80% in Australia and over 90% in the US.

We explore the equilibrium transition for the model economy as follows. The number of type-L parents in generation \(t+1\) is determined by

$$\begin{aligned} N_{t+1}^{L}=N_{t}^{L}n_{t}^{L}(1-\lambda _t^L)+N_{t}^{H}n_{t}^{H}(1-\lambda _t^H). \end{aligned}$$
(13)

Similarly, the number of type-H parents in generation \(t+1\) is determined by

$$\begin{aligned} N_{t+1}^{H}=N_{t}^{L}n_{t}^{L}\lambda _t^L+N_{t}^{H}n_{t}^{H}\lambda _t^H. \end{aligned}$$
(14)

Let \(\rho _{t}\equiv N_{t}^{L}/N_{t}^{H}\) be the ratio of illiterate to literate parents and \(\mu _{t}\equiv n_{t}^{L}/n_{t}^{H}\) be the ratio of illiterate parents’ fertility to literate parents’ fertility. Here, \(\rho _t\) magnifies illiteracy rates to the range \((0,\infty )\) monotonically. Dividing (13) and (14) by the number of children of literate parents \(N_{t}^{H}n_{t}^{H}\) gives

$$\begin{aligned} \frac{N_{t+1}^{L}}{N_{t}^{H}n_{t}^{H}}=\rho _{t}\mu _{t}(1-\lambda _t^L)+(1-\lambda _t^H), \end{aligned}$$
(15)
$$\begin{aligned} \frac{N_{t+1}^{H}}{N_{t}^{H}n_{t}^{H}}=\rho _{t}\mu _{t}\lambda _t^L+\lambda _t^H. \end{aligned}$$
(16)

Thus, the transition equation of the ratio of illiterate to literate parents is

$$\begin{aligned} \rho _{t+1}=\frac{\rho _{t}\mu _{t}(1-\lambda _t^L)+(1-\lambda _t^H)}{\rho _{t}\mu _{t}\lambda _t^L+\lambda _t^H}. \end{aligned}$$
(17)

The evolution of this ratio has the following properties:

Proposition 3

The ratio of illiterate to literate parents in the next period \(\rho _{t+1}\) increases with its current ratio \(\rho _t\) and the ratio of illiterate parents’ fertility to literate parents’ fertility \(\mu _{t}\) but decreases with school enrollment ratios \(\lambda _t^k\).

Proof

See Appendix C. \(\square \)

The positive relationship between the future and current ratios of illiterate to literate parents suggests the difficulty in improving the education status of the population or reducing illiteracy rates for development. High ratios of illiterate parents’ fertility to literate parents’ fertility also impede the reduction in illiteracy rates. Education subsidies can improve parental education status across generations by increasing school enrollment ratios and reducing fertility.

From (13) and (14), the average fertility is determined by

$$\begin{aligned} \frac{N_{t+1}}{N_t}=\frac{n_t^H(1+\rho _t\mu _{t})}{1+\rho _t}. \end{aligned}$$
(18)

The average fertility increases with the fertility rates of all parents or the ratio of illiterate parents’ fertility to literate parents’ fertility; however, its relationship with the illiteracy rate depends on fertility differentials as follows:

Proposition 4

If illiterate parents have higher fertility than literate parents \(\mu _t>1\), then average fertility increases with the ratio of illiterate to literate parents \(\rho _t\); otherwise, average fertility decreases with \(\rho _t\).

Proof

See Appendix D. \(\square \)

This result has novel implications. When illiterate parents have higher fertility than literate parents, falling illiteracy rates will decrease average fertility. Conversely, when literate parents have higher fertility than illiterate parents, falling illiteracy rates will increase average fertility. Thus, education subsidies are particularly effective to reduce illiteracy rates and average fertility in countries with higher fertility from illiterate parents than from literate parents. As noted in Tables 1 and 5, Brazil with relatively high fertility from illiterate parents and a high ratio of public education spending to GDP attained a larger rise in school enrollment ratios and a sharper decline in average fertility than Indonesia did.

The growth rate of average output in the economy from period t to \(t+1\) is

$$\begin{aligned} g_{t+1}\equiv \frac{N_{t+1}^H(1-v n_{t+1}^H)h_{t+1}^H+N_{t+1}^L(1-v n_{t+1}^L)h_{t+1}^L}{N_t^H(1-v n_t^H)h_t^H+N_t^L(1-v n_t^L)h_t^L}. \end{aligned}$$
(19)

Given a state \((N_t^L, N_t^H, h_t^L, h_t^H,h_t^C)\), the growth rate of average output (\(g_{t+1}\)) increases if fertility rates fall over time (lower \(n^k_{t+1}\)) or if school enrollment ratios rise (higher \(N_{t+1}^H\)). The growth rate also increases with current fertility (\(n^k_t\)) since high \(n^k_t\) means low labor (hence low income) given the time cost of child rearing.

From (17) and (19), the long-run equilibrium is as follows:

Proposition 5

For stationary \((\xi ,\tau ^L,\tau ^H,s^L,s^H, d_t/h_t^k, h^C_t/h_t^k)\) and \(\sigma =\gamma =1\), the economy converges to a unique steady-state growth path with

$$\begin{aligned} \rho ^*=\frac{\mu (1-\lambda ^L)-\lambda ^H+\sqrt{[\mu (1-\lambda ^L)-\lambda ^H]^2+4\mu \lambda ^L(1-\lambda ^H)}}{2\mu \lambda ^L}\ge 0, \end{aligned}$$
(20)
$$\begin{aligned} g^*=A(\psi +B\xi )^\phi \left[ \frac{\omega (\rho ^* n^L\lambda ^L+n^H\lambda ^H)+\rho ^* n^L(1-\lambda ^L)+n^H(1-\lambda ^H)}{\omega (1-vn^H)+\rho ^*(1-v n^L)}\right] . \end{aligned}$$
(21)

Proof

See Appendix E. \(\square \)

The convergent and stable illiteracy rate creates a poverty trap, calling for rises in education subsidy rates. A permanent rise in education subsidy rates raises school enrollment ratios of all children, leading to a lower illiteracy rate and a higher growth rate of average output in the long run. However, the steady state applies to a special case with stationary education subsidy rates, stationary ratios of the living cost, education cost and child earnings to parental earnings, and logarithmic utility.

In Fig. 3, the rise in school enrollment ratios and decline in fertility rates in the process of rising parental earnings arise from constant child earnings under the assumptions of stationary education subsidy rates and stationary ratios of the living and education costs to parental earnings. Although the illustrated numerical results match the observed trends of enrollment ratios and fertility, child earnings and education subsidy rates may rise over time, education subsidy rates may favor children of poor parents, and the utility function may be more general than logarithmic. The illustrated numerical results also have limitations: The fertility rates at high enrollment ratios are still higher than the observed levels in the data; the gap in enrollment ratios for children across families is larger than observed in the early stage; the first flip of fertility differentials occurs when poor parents forgo investment in education for children as in the literature. As shown numerically in Greenwood et al. (2005), if \(\gamma >\sigma =1\) (relatively inelastic fertility), then fertility may decline when income increases. We consider this scenario and more realistic education subsidy rates and child earnings for quantitative implications next.

5 Quantitative implications

We now simulate the equilibrium path with \(1/\gamma <1\) for quantitative implications to account for the patterns of fertility rates and school enrollment ratios in the data of Brazil and Indonesia. Table 6 presents the parameter values for simulations in Cases I and II for Brazil and Indonesia, respectively. One period of the model is 20 years. In de la Croix and Doepke (2003) and Vogl (2016) with logarithmic utility from consumption and fertility, the fraction of parental time for rearing a child is set at \(7.5\%\). Vogl (2016) also sets the living cost per child at 7.7% of median full income. Given that 80% of parents are illiterate in 1971 in Indonesia in Table 4 (57% in 1960 in Brazil), we use the time cost in their studies and the living cost per child as 7.7% of an illiterate parent’s earnings in both cases.

Table 6 Parameterization

As noted in Table 1, the ratios of public education spending to GNI were as high as 4.3% in both countries. Thus, we assume the parameter values for the education technology (\(\xi =0.03, A=7.5, B=2.6, \phi =0.9\), and \(\psi =0.1\)) such that the return to education spending is \(\omega (\xi )=1.6803\) and that education spending per child is \(\xi H^H/H^L=\)5.05% of an illiterate parent’s earnings. From Assumption 1, we set the relative roles of the number and average quality of children in the utility function at \(\alpha =0.6\) and \(\eta =0.4\) to permit plausible values for the enrollment ratios in equation (11) in both countries with different tax rates and education subsidies.

From the progressive income tax rates in Indonesia and Brazil in footnote 12, the respective income tax rates for educated and illiterate parents are 20% and 3% in Case II as used for Fig. 3. Since the top tax rate is lower but the bottom rate is higher in Brazil, we set the respective tax rates at 18% and 5% in Case I. As per capita GDP in Indonesia was between a quarter and one-third of that in Brazil in the 1960s in Table 1, the initial levels of human capital or earnings in Cases I and II of Table 6 reflect the gap in initial income levels in these two countries.

The remaining parameter values to account for the enrollment ratios in Fig. 2 are child earnings and education subsidies. From Fig. 2 and Table 1, the enrollment ratios of children in families with illiterate parents were lower in 1960 in Brazil despite a higher ratio of public education spending to GDP than in Indonesia, suggesting that child earnings were higher in Brazil. Thus, we set child earnings as 3.8% of an illiterate parent’s earnings in Case I and 2% in Case II. Since child earnings and the costs of living and education are all proportional to parental earnings in Table 6, school enrollment ratios should remain constant over time by Proposition 1 unless education subsidy rates rise as observed in the data.

As in Fig. 3, equal subsidy rates for all families yield a much larger gap in enrollment ratios across families than the observed gap. Although the census data has no information about education subsidy rates for each type of family, education subsidy rates are higher for poor families in developing countries as noted by Morley and Coady (2003). From Table 1, public spending on education subsidies in Indonesia started at much lower levels (0.767% of GDP) than in Brazil (3.599%) in 1970–1980 but caught up (to 4.328%) in 2010. Both Brazil and Indonesia have been hit by several political or economic shocks relating to education policy or funding. Brazil passed a new federal constitution in 1988 that expanded the rights and resources in the field of education (Weller et al. 2020). In Indonesia, the government doubled the expenditures on regional development in 1973–1980, particularly for school construction or funding, following the oil boom (Duflo 2001). Yet, the lower education subsidies in Indonesia in 1970–1980 supported a higher enrollment ratio for children of illiterate parents in Fig. 2, who accounted for nearly 80% of the parent population in Table 4, than that in Brazil even though Brazil also had higher GDP per capita and lower illiteracy in 1960. The comparisons suggest a possibility that education subsidies were highly in favor of poor families in Indonesia initially and became more equal over time.

Specifically, we assume that education subsidies begin at respective rates 67.5% and 82% for children of educated parents and children of illiterate parents and grow respectively at 1.6% and 1.5% per period in Case I to account for the rises in enrollment ratios in Brazil in Fig. 2. In Case II, the respective subsidy rates begin at low rates 16% and 35%, and grow respectively at 2.3% and 2% in the first two periods (1.5% and 1.8% afterwards) to account for the rises in enrollment ratios in Indonesia.

From Table 1, the fertility rate in 1960 was much higher in Brazil than in Indonesia despite a much higher income level in Brazil, suggesting stronger taste for utility from children in Brazil. Fertility rates also declined dramatically to the replacement level in both countries in 2010, suggesting an inelastic preference for the number of children relative to consumption. The decline in fertility in data is much larger than the illustrated level in Fig. 3 that arises from the constancy of child earnings and education subsidy rates. To match the patterns of fertility rates in different families in Table 5, we set the taste for utility from children’s quantity and quality at \(\beta =26\) in Case I and 10.5 in Case II, the elasticity of substitution with respect to consumption at \(1/\sigma =1\) (logarithmic utility), and the elasticity of substitution with respect to fertility at \(1/\gamma =1/2.8\). In Greenwood et al. (2005), the elasticity of substitution with respect to consumption is \(1/\sigma =1\) and the elasticity of substitution with respect to fertility is \(1/\gamma =1/1.2\) that reduces fertility over time when income increases. The reason for the lower elasticity of substitution for fertility than in their model is rising child earnings in this model that induce parents to have high fertility. If child earnings increase at a lower rate than parental earnings, then fertility can fall as in Fig. 3 and thus the elasticity of substitution for fertility can be higher than 1/2.8 to match the observed decline in fertility. However, the census data has no information about child earnings.

Figures 4 and 5 present simulation results in Cases I and II in 8 periods based on parameter values in Table 6. Each figure has four panels. Figures 4(a) and 5(a) show that educated parents always have higher school enrollment ratios for children than illiterate parents. Both enrollment ratios increase over time until reaching 100% due to rising education subsidy rates when child earnings and the costs of living and education grow proportionately with parental earnings. From the higher education subsidy rates for children of illiterate parents, enrollment ratios are positive for all children in all periods as in the data and the gap in the enrollment ratios across families is much smaller than in Fig. 3 and closer to the observed level.

Fig. 4
figure 4

Fertility, school enrollment and economic growth I

Fig. 5
figure 5

Fertility, school enrollment and economic growth II

In Figs. 4(b) and 5(b), the simulated fertility rates have a downward trend due to rising parental income, rising subsidy rates, and inelastic preference for fertility \(1/\gamma <1\). In Fig. 4, the fertility rate is lower for educated than for illiterate parents in all periods as in Brazil. In Fig. 5, the fertility rate of educated parents starts higher but decreases more rapidly than that of illiterate parents as in Indonesia. Consequently, the fertility differential flips to a lower fertility rate of educated parents than that of illiterate parents in the first three periods as observed in Indonesian. The difference between the two figures are due to higher earnings and education subsidies in Fig. 4 than in Fig. 5.

The flip of the fertility differential in Fig. 5 arises from the difference in education subsidy rates for children in different families at low income levels under progressive taxes. On the one hand, the higher subsidy rate for children of illiterate parents induces the parents to have higher enrollment ratios and lower fertility than they would have otherwise. On the other hand, the partial catch-up rise in education subsidy rates for children of educated parents in the first two periods yields a sharper decline in the fertility of educated parents than that of illiterate parents. The high tax rate for educated parents also induces educated parents for high fertility, particularly in the initial period or early development when income is low. The flip in the fertility differential occurs when the enrollment ratios are positive for children in all types of families in contrast to Fig. 3 where illiterate parents forgo investment in children’s education at the flipping point as in the aforementioned literature.

Figures 4(c) and 5(c) plot the growth rate of average output and average fertility. Average fertility declines from above five to the replacement level of two as observed in Table 1. The growth rate displays large swings due to the uneven paces of changes in school enrollment ratios and fertility rates over time and across different groups of families. By contrast, changes in growth rates are relatively small in models with equal and divisible education investment for all children within families.

In Figs. 4(d) and 5(d), the share of educated parents in parental population or the literacy rate starts from a low level, and initially declines in Case I but rises slowly in Case II due to the opposing effects of rising school enrollment ratios but rapidly falling fertility of educated parents. The opposite fertility differentials in Cases I and II are harmful for the average education attainment in Case I but conducive in Case II, while the larger increases in education subsidy rates in Case II are also conducive to the average education attainment. When the school enrollment ratios increase and fertility rates decline further for all families, the share of educated parents rises rapidly above 80%. Overall, the simulated results capture the patterns of the movements observed in Brazil and Indonesia in Table 1 and Fig. 2. In particular, the rise in education subsidy promotes development by accelerating human capital accumulation and the demographic transition.

When accounting for the stylized facts, both the endogenous dynamics and the exogenous rises in education subsidies play important roles. Education is a driving force for growth in parental earnings that, in turn, decreases fertility in this endogenous growth model with the choices of school enrollment and fertility. However, when child earnings, the cost of living, and the cost of education increase proportionately with parental earnings, enrollment ratios would remain constant for stationary education subsidy rates by Proposition 1. Thus, restricting child labor or raising the education subsidy rate is important to increase enrollment ratios and decrease illiteracy rates as in Figs. 3 and 4.

At the early development stage with low average earnings but high fertility and illiteracy rates, differentiated education subsidies in favor of children of illiterate parents are useful to induce them to enroll at least some children in school. Otherwise, illiterate parents may have little incentive for children’s education because child earnings are particularly important for illiterate parents with low earnings. Sufficiently favorable education subsidies for children of illiterate families can also yield lower fertility for illiterate parents than for educated parents under progressive labor-income tax rates that reduce the opportunity cost of time for rearing children. In this case, a partial catch-up rise in the education subsidy rate for children of educated parents can enlarge the decrease in their fertility for the observed flip of the fertility differential in early development with positive school enrollment ratios in all families as in Fig. 5.

6 Conclusion

This paper explored differential fertility rates and school enrollment ratios among parents with different education attainments. In this model parents constrained by children’s education and living costs choose the number of children, a fraction of children for education, and the remaining children for supplemental income. This distinct mechanism yields rich implications to account for the observed patterns.

When the opportunity costs of enrolling children in school are high relative to the forgone earnings of rearing children and living costs, parents have as many children as possible and no incentive for children’s education. When income grows relative to the subsidized education costs, parents enroll a fraction of children in school. The enrollment ratio increases with the returns to education, forgone earnings of rearing children, and living costs but decreases with the subsidized education costs and child earnings. The positive relationship between enrollment ratios and parental education or earnings matches the pattern in the census data of Brazil and Indonesia.

The model also accounts for the N-shaped relation between fertility and parental earnings in the data. When children’s living costs are high relative to children’s education costs, fertility increases with parental earnings due to a dominant income effect and thus fertility is higher for parents with higher education attainments as in the data of Indonesia before 1990. When the living costs are low relative to the education costs, fertility decreases with parental earnings due to a dominant substitution effect and thus fertility is lower for parents with higher education attainments as observed in Indonesia after 1990 and in other countries.

When child earnings and the costs of education and living increase with parental earnings, school enrollment ratios remain constant and the illiteracy rate of the economy is convergent and stable for stationary education subsidy rates, creating a hurdle for development. An increase in education subsidy rates raises school enrollment ratios and reduces fertility, thus reducing illiteracy rates. However, the decline in illiteracy rates may raise the average fertility in the economy when illiterate parents have lower fertility than educated parents, especially under progressive income taxes and favorable education subsidies for poor families, as observed in Indonesia in 1960–1990. A partial catch-up rise in education subsidies for children of literate parents can enlarge the decrease in their fertility to flip the fertility differential. Conversely, when literate parents have lower fertility than illiterate parents, the decline in illiteracy rates reduces average fertility and an increase in education subsidy rates is particularly effective to reduce average fertility and illiteracy rates for development as in Brazil. When most children complete all education levels, fertility may relate positively with parental earnings again as observed in developed countries with high tertiary enrollment ratios.