The economics of missionary expansion: evidence from Africa and implications for development

How did Christianity expand in Africa to become the continent’s dominant religion? Using annual panel census data on Christian missions from 1751 to 1932 in Ghana, and pre-1924 data on missions for 43 sub-Saharan African countries, we estimate causal effects of malaria, railroads and cash crops on mission location. We find that missions were established in healthier, more accessible, and richer places before expanding to economically less developed places. We argue that the endogeneity of missionary expansion may have been underestimated, thus questioning the link between missions and economic development for Africa. We find the endogeneity problem exacerbated when mission data is sourced from Christian missionary atlases that disproportionately report a selection of prominent missions that were also established early.


Introduction
One of the most powerful cultural transformations in modern history has been the rapid expansion of Christianity to regions outside Europe. Conversion has been markedly rapid in sub-Saharan Africa. While Christians made up about 5% of the population in 1900, their share has grown to 57% today (Johnson and Grim 2020), making Africa the continent with the highest number of Christians (Todd et al. 2018). At current trends, Africans will comprise 42% of the global Christian population by 2060. According to the World Values Survey, sub-Saharan Africa is also home to the world's most observant Christians in terms of church attendance, making it "the future of the world's most popular religion" (Economist 2015).
The process of African mass-conversion was facilitated by vast Christian missionary efforts. Throughout the colonial era missions provided the bulk of formal education (Frankema 2012;Cogneau and Moradi 2014) and health care (Doyle et al. 2020). A recent and rapidly growing literature explores the long-term effects of missions in colonial Africa, as well as in India and Latin America. Web Appx. Table A1 summarizes 50 studies that find strong effects on local economic development (Bai and Kung 2015;Castelló-Climent et al. 2018;Valencia Caicedo 2019a), education (Gallego and Woodberry 2010;Acemoglu et al. 2014;Waldinger 2017;Baten et al. 2021;Calvi et al. 2021), health (Calvi and Mantovanelli 2018;Cagé and Rueda 2020;Menon and McQueeney 2020), fertility (Guirkinger and Villar 2022), social mobility (Wantchekon et al. 2015;Alesina et al. 2021), culture (Nunn 2010(Nunn , 2014Fenske 2015;Calvi et al. 2022), and political participation (Cagé and Rueda 2016). 1 While a large literature has studied the economic consequences of evangelization, the economics behind the expansion of Christianity remains poorly understood. One narrative describes missionaries as adventurers, with little to no information on local circumstances, crossing political boundaries and whose objective was to save souls no matter the cost (Oliver 1952;Cleall 2009). Their locational choices were highly erratic, but once settled, missions persisted. For example, Wantchekon et al. (2015, p. 714) describe how a series of historical accidents led missionaries to settle left rather than right of a river. An alternative view is that missionaries were following clear expansion strategies laid out by their own mission society (Johnson 1967;Terry 2015). For example, missionary Thauren (1931, pp. 19-20) described the strategy in Togo as follows: "The mission leadership knew that choosing suitable places would be crucial for missionizing the interior. Therefore, no effort was spared to get to know the interior better [...]. The mission society aimed to establish new stations in larger cities, so that few missionaries could spread the gospel to many. In particular, they preferred cities with regular markets. Moreover, the places needed to be centrally located [...]." Missions also aimed for financial independence, which would make them target wealthy locations. To our knowledge, the causal determinants of mission expansion have not been studied. 2 In this paper, we study the determinants and effects of missionary expansion in Ghana and sub-Saharan Africa as a whole. For various reasons, Ghana is an ideal testing ground. First, Ghana is one of the most devout Christian African countries today. 3 Second, Ghana received missionaries since the early 19th century, which allows us to investigate the role of quinine as a groundbreaking treatment for malaria ca. 1850. Missionary expansion in East Africa began later in the 1890s. Third, in Ghana missionaries operated in a comparatively free religious market with little interference by the colonial state (Gallego and Woodberry 2010). 4 Finally, Ghana holds rich historical data. Altogether, this enables us to shed light on a wide range of determinants of mission expansion.
For Ghana, we create a novel annual panel dataset on the locations of missions at a precise spatial level over two centuries: 2091 grid cells of 0.1 × 0.1 degrees (12 × 12 km) in 1751-1932, when the number of missions increased from one to 1,832. Creating new data on the determinants of mission locations, descriptive results suggest that missionaries might have gone to healthier, more accessible, and richer areas, privileging these "better" locations first. They might also have invested more into these missions/locations (with administrative functions, European missionaries and schools). These results are confirmed using identification strategies for malaria, railroads, and cash crops, in Ghana, but also in sub-Saharan Africa.
Once missionary expansion is better understood, we can study the long-term effects of missions. With a few exceptions, existing studies do not analyze the effects of missions on economic development per se, focus on one mechanism at a time, and, more importantly, use data from mission atlases that are marred by omissions (see Web Appx. Table A1). Figure 1a plots the mission locations obtained from the two most commonly used sources: (i) the Atlas of Protestant Missions (Beach 1903), first used by Cagé and Rueda (2016), and (ii) the Ethnographic Survey of Africa (Roome 1925), digitized by Nunn (2010Nunn ( , 2014. 5 Using ecclesiastical census returns and other primary sources instead, we uncover vast discrepancies. Figure 1b shows that atlases omit more than 90% of missions in most African countries. 6 Beach (1903) and Roome (1925) imply congregation sizes of more than 3,000 and 15,000 per mission station respectively. 7 For Ghana, the two atlases miss 91-98% of missions, most of them in the hinterland (see Fig. 2a-b). Our census data suggest an average congregation size of 198 in 1925, which appears far more plausible. 8 We discuss the implications of the observed patterns of missionary expansion and the over-reliance on mission atlases for the study of the long-term effects of missions. We highlight three sources of endogeneity: (i) attenuation bias, (ii) selection bias, and (iii) omitted variable bias.
First, we scrutinize the standard source of mission data. Atlases do not show the distribution of the universe of missions. Classical measurement error would lead to attenuation bias. However, measurement error may be non-classical, e.g. if omissions are not random. 4 In contrast, Belgian, Portuguese and Spanish colonies granted Catholic missions a quasi-monopoly (Meier zu Selhausen 2019). French colonies followed from the 1890s on a more neutral treatment of missionaries in that they restricted most of the activities of both Catholics and Protestants (Sundkler and Steed 2000, p. 430). 5 30 studies use Roome (1925), 6 studies use Beach (1903), 5 use Bartholomew et al. (1911), and 6 use Streit (1913). 6 Weighting the country-specific omission rates in Beach (1903) and Roome (1925) by the countries' population ca. 1920 (Frankema and Jerven 2014), we find for the whole of sub-Saharan Africa that 76% of Protestant missions are missing from Beach (vs. 91% for Ghana) and 91% of them are missing from Roome (vs. 98% for Ghana). 7 Beach (1903) and Roome (1925) report 677 Protestant and 1,212 Christian mission stations respectively. According to the World Christian Database (Johnson et al. 2021), Protestants in 1900 numbered 2.2 million. For the year 1925, there were an estimated 20.2 million Christians in sub-Saharan Africa (excluding Ethiopia).

3
We use our census-based data on Ghana to show that atlases positively select missions, in particular early, main or European residence missions located in more developed areas. This would bias estimates upward.
Next, we examine whether the source of mission data affects the estimated long-term impact of missions in Ghana. Employing the control variables commonly used in the literature on missions in Africa (e.g., Nunn 2010Nunn , 2014Cagé and Rueda 2016), we find strong effects that are twice as large for atlas missions than for census missions. If we use a rich and historically relevant set of control variables, estimated effects are considerably weaker. This also applies to Africa as a whole.
Our interpretation is that mission atlases, because they report "better" missions endogenously located in "better" places, may capture effects that are not due to the missions themselves. Moreover, relying on atlases raises concerns of external validity. Estimated effects are measured for a selected subset of missions that are not representative of the average mission. To measure causal long-term effects of missions, one must properly control for the spatio-historical determinants of missions, or use an identification strategy that can be implemented even with imperfect historical controls. While a handful of studies have relied on innovative identification strategies (e.g., Wantchekon et al. 2015;Cagé and Rueda 2016;Waldinger 2017; Valencia Caicedo 2019a) 9 , many studies use some controls and sometimes add spatial fixed effects.
Our contribution is threefold. First, our paper relates to the literature on the determinants of the adoption and diffusion of religion (Bisin and Verdier 2000;Barro and McCleary 2005;McCleary and Barro 2006;Becker and Woessmann 2013;Becker et al. 2017). There are few quantitative works on religious diffusion, and none for sub-Saharan Africa.  pointed to the role of geography and trade in the diffusion of Islam. Cantoni (2012) linked the spread of Protestantism in 16th century Germany to strategic choices of territorial lords, while Rubin (2014) shows that Protestantism followed cities' adoption of the printing press. We show that in Africa Christianity spread along railroads, with the cultivation of cash crops, and with economic development more generally, before it diffused more broadly.
Second, our contribution is methodological. Non-classical measurement error in historical data and their consequences for the analysis of path dependence are under-investigated. Measurement error in survey data, in contrast, has received a lot of attention (e.g. Bollinger 1996;Mahajan 2006). We emphasize the importance of: (i) reliable mission data, especially given external validity issues; (ii) relevant controls capturing the dynamics and factors behind missionary expansion; and, when needed, (iii) identification strategies that bypass issues related to these measurement issues.
Finally, we add to the literature on the long-term impact of missions. 10 A large body of literature has found positive effects on various proximate determinants of economic growth (Web Appx. Table A1). Missions were the most important human capital promoting institutions to have emerged in developing countries during the colonial era. By shaping attitudes towards factors such as education and fertility, factors that alter the religious orientation of society may have important repercussions for the growth process (Galor 2011;  Roome (1925) (N = 1212). Subfigure 1(b) shows when the data is available the share of missions in census sources that are missing in Beach (1903) (for Protestants only in 1900) and Roome (1925) (for both Protestants andCatholics in 1924). See Web Data Appendix for data sources.  Beach (1903) (N = 24) and in our census data (304). Subfigure 2(b) shows the 1924 missions in Roome (1925) (N = 24) and in our census data (1213). See Web Appendix for data sources Becker et al. 2017). With our improved data we find comparably weaker long-term effects of missions in Ghana and Africa. 11

New data on Ghana and Africa
To examine the determinants of missionary expansion, and draw implications for the analysis of long-term development, we require data on the location of missions over time, the locational factors potentially driving missions' spatial expansion, and local economic development today. The Web Appendix provides more details on sources. Appx. Table A2 shows summary statistics.

Missions in Ghana.
We partition Ghana into 2091 grid cells of 0.1 x 0.1 degrees (12 x 12 km) and construct an annual panel data set from 1751 to 1932 (181 years). Merriam-Webster (2020) defines a mission station as "a place of missionary residence in or from which missionary activity in a given area is carried on." We recreate the history of every mission station (N = 2144) for all mission societies (classified as Presbyterian, Methodist, Catholic and other) and geocode their locations. Our main sources are the ecclesiastical census returns published in the Blue Books of the Gold Coast, 1844-1932 (see Web Appx. Fig. A1 for an example). Each mission society was required to submit annual reports on its activities to the colonial administration, thereby listing all of their stations. Churches received annual grants from the colonial government for their pastoral services, which provided a strong incentive to report. Our data thus represents an exhaustive census of missions. The early origins of mission societies are then well documented by society-specific anniversary reports, and we have no difficulties reconstructing missions before 1844.
Using the same sources, we identify if a mission was a main station or an out-station. Main stations are the principal centers of a "circuit" -a society's administrative unit. Main stations are larger and centrally located (Thauren 1931); they are headed by an ordained missionary. The congregations of out-stations are smaller but still of significant size and taken together have more members than main stations. All mission posts, not just main stations, built churches. 12 We also identify if a mission had a school. We focus on "assisted schools", which followed the government school curriculum and certified quality standards (Williamson 1952). As they received grants-in-aid, they were reported accurately. The existing mission 11 Within the Unified Growth Theory (UGT) literature (see Galor 2005 for surveys), various studies examine how the emergence of human capital promoting institutions has affected the pace of the transition from stagnation to growth . These studies have focused on the expansion of public schooling, such as the creation of a publicly supported secondary school system in the U.K. in 1902 (Galor and Moav 2006) and the U.S. High School Movement in the first half of the 20th century (Moav and Vollrath 2009). The weak long-term economic effects of colonial missions in Africa could be due to African countries having not yet fully entered the process of human capital demanding "industrialization" (see, for example, Gollin et al. 2016) that the UGT has shown to be the main driver behind human capital formation and the demographic transition, and thus the transition to modern economic growth. 12 Typically, the first investment of every mission post was to build a church. For some years in the 1860s and 1870s, the ecclesiastical returns report the type of building, attendance and seat capacity of the mission post, the latter being a measure of how large the church building was. In 1867, 5% and 8% of Wesleyan and Presbyterian mission stations only had a "house" or a "room" respectively. Not only is this a small percentage, it is also of transitory nature. Within 10 years, all those stations had a church building. Seat capacity and attendance are strongly correlated. literature has not distinguished between main, out-stations and schools. We follow this decision in our main analysis. However, we will also study heterogeneous effects with respect to the type of mission.
Missionaries in Ghana. We create a data set of all 338 male European missionaries stationed in Ghana during 1751-1890 from a variety of sources (data not available post-1890). For the mortality analysis in Sect. 3, we reconstructed dates of service and death in Ghana. African missionary careers are less well-documented. From the Blue Books 1846-90, we retrieved the localities where European missionaries were permanently based and which missions they served occasionally. 13 Locational Factors for Ghana. We construct a GIS data set at the same grid resolution: (i) Geography: Historical malaria intensity (based on genetic data) comes from Depetris-Chauvin and Weil (2018). We compute the distance to the coast and obtain ports c. 1850 from Dickson (1969); (ii) Political conditions: Data on large pre-colonial cities before 1800 are from Chandler (1987) and headchief towns in 1901 are from Guggisberg (1908). From Dickson (1969), we derive the boundary of the Gold Coast Colony established by the British c. 1850; (iii) Transportation: We obtain from Dickson (1969) navigable rivers in 1850-1930and trade routes circa 1850. Railroads (1897-1932 and roads (1932) come from Jedwab and Moradi (2016); (iv) Population: From census gazetteers, we compile a GIS database of towns above 1000 inhabitants in 1891, 1901 and 1931. We also collect rural population data for 1901 and 1931. Because all cells have the same area, population levels are equivalent to densities; 14 (v) Economic activities: Slave export and slave market data come from Nunn (2008) and Osei (2014) respectively. We obtain cash crop production areas from Dickson (1969) and total export value of cash crops from Frankema et al. (2018). 15 Mines are from Dickson (1969); and (vi) Other: We obtain area, mean annual rainfall (mm) in 1900-1960, mean altitude (m), ruggedness (standard deviation of altitude), and soil fertility.
Contemporary Data for Ghana. We use satellite data on night lights in 2010 as a proxy for local economic development (NGDC 2015). We use the radiance calibrated version of this data, to avoid issues related to top-coding. 16 Census data on religion, education, fertility, mortality, employment and urbanization in 2000 are from Ghana (2000). To study child mortality, we use seven Demographic and Health Surveys (DHS) between 1993(USAID 2020 (Schott 1879). 14 While we have exhaustive urban data for all census years, we only have georeferenced rural population data for Southern Ghana in 1901. We thus include a dummy if any locality in the cell was surveyed by the 1901 census. 15 We obtain soil suitability for the same cash crops from the 1958 Survey of Ghana Classification Map of Cocoa Soils for Southern Ghana, Survey of Ghana, Accra, as well as Gyasi (1992) and Globcover (2009). 16 This data records levels of luminosity beyond the normal digital number upper bound of 63. 17 Since we only have data for 10% of the population census, the most rural cells of our sample do not have enough observations to correctly estimate these shares. We keep 1895 cells with at least 10 observations. When using the DHS, we keep cells with at least 5 children. Due to the smaller sample size, we must indeed adopt a smaller threshold.
1 3 each station systematically over such a long period. Yearbooks of other colonies are completely silent. We thus use mission location data widely used in the literature. These stem from mission atlases. (Fahs 1925, p. 271) writes that these atlases are based on "hundreds of documents" and "society field reports" and admits "problems of what to include in or to exclude (...)." Keeping this caveat in mind, we use Beach (1903), compiled by Cagé and Rueda (2016), which reports the locations of 677 Protestant missions in 1900, and added the year of foundation ourselves. We then use Roome (1925), digitized by Nunn (2010), which shows the respective locations of Catholic andProtestant missions in 1924 (N = 1212).
The set of locational factors for Africa follows the one for Ghana: (i) Geography: Historical malaria intensity is from Depetris-Chauvin and Weil (2018) and tsetse fly ecology from Alsan (2015). We compute distance to the coast, and 19th century slave ports are from Nunn and Wantchekon (2011); (ii) Political conditions: Data on large pre-colonial cities before 1800 are from Chandler (1987). Data on the capital, largest and 2nd largest cities come from Jedwab and Moradi (2016). The year of colonization for each ethnic group is from Henderson and Whatley (2014). Using the Murdock (1967) map of ethnic boundaries (Nunn 2008), we then assign a year of colonization to each cell. From the same sources, we know if the cell was in an ethnic area with a centralized state before colonization. We compute the distance to historical Muslim centers based on Ajayi and Ade and Michael Crowder, (1974) and Sluglett (2014); (iii) Transportation: We obtain navigable rivers and lakes from Johnston (1915), pre-colonial explorer routes from Nunn and Wantchekon (2011) and railroads from Jedwab and Moradi (2016); (iv) Population: We control for population density c. 1800 and log urban and rural population c. 1900 from HYDE (Goldewijk et al. 2010), and log city population c. 1900 for towns above 10,000 from colonial administrative sources; 18 (v) Economic activities: We know if slavery and polygamy was historically practiced (Murdock 1967). The log number of slaves exported per land area is from Nunn and Wantchekon (2011). Land suitability measures for various cash crops are from FAO (2012). We obtain cash crops' national export value in 1850-1924(Frankema et al. 2018. Mines in 1900 and 1924 come from Mamo et al. (2019); (vi) Other: We obtain area, mean annual rainfall (mm, 1900-60), mean altitude (m), ruggedness, and soil fertility. 19 Finally, we also use satellite data on night lights in 2010 (NGDC 2015).

Background: missionary expansion in Ghana
Colonization. Since the 15th century European powers had established trading posts along the West African coast. By 1821, Britain had established an informal protectorate " the Gold Coast " in the coastal regions of Ghana. In 1874, the British defeated the inland Ashanti Kingdom centered around its capital Kumasi. The ensuing peace treaty of 1875 transformed the protectorate into a formal British colony with the same boundaries, the Gold Coast Colony. Between 1896 and 1900, two wars with Ashanti eventually forced the kingdom into a colony in 1902. British control was extended to the north at this occasion. 18 Goldewijk et al. (2010) do not rely on census data for earlier centuries (there were no censuses then). These population estimates are unreliable. We nonetheless use them when replicating controls used in the literature. 19 We also add a dummy if the main ethnic group in the cell has data in the Murdock (1967) Atlas and a dummy if the underlying anthropological survey used to create this data precedes 1900 or 1924.
Rail construction began in 1897, which helped the British to consolidate their control over Ghana and lowered trade costs. This motivates the choice of five turning points for our descriptive analysis: 1751, 1850, 1875, 1897, and 1932. Missionary Expansion. The first mission station was established in 1751. Thereafter, various mission societies failed to maintain a permanent presence until Presbyterian and Methodist missionaries reached the Gold Coast in 1828 and 1835, respectively. Figure 3 shows the number of Christian missions, main stations, and mission schools from 1840 to 1932. By 1850, 904 Ghanaians had converted and 21 missions existed (Isichei 1995, p. 169;Miller 2003, p. 23). Mass-evangelization did not take off until the 1870s, when 67 Protestant missions served 6,000 Ghanaians. Catholic missions started their conversion efforts from 1880 onward. By 1932, the number of missions had expanded to 1832 with 340,000 followers (9% of the population; ≈ 186 followers per mission). The Christian share has since grown to 80% (USAID 2020). Missions viewed the provision of education as a way to attract new Christians. As such, they provided the bulk of schooling in colonial Ghana (Cogneau and Moradi 2014). As seen in Fig. 3, early missions differed from later missions in that many were main stations and had a school.
Constraints. Missions initially settled along the coast ( Fig. 4a-b). Missionaries shunned away from creating inland stations before they gathered essential intelligence traveling the country (Thauren 1931, pp. 19-21;Engel 1931, p. 14). The Ashanti Kingdom was hostile to Christian proselytizing, restricting missionary activities to the territory of the Gold Coast Colony (Fig. 4b-c). Access to the interior was facilitated by rail and road building since the early 20th century. By 1932, missions covered large parts of Southern Ghana (Fig. 4d).
Malaria was the biggest killer, striking Europeans soon after arrival (Curtin 1961). This changed when quinine was introduced circa 1850 as cure and prophylaxis. 20 The significance was well understood. As Curtin (1973, p. 362) explains, quinine "helped to close an epoch. [It] ... was understood well enough in official and missionary circles to reduce sharply the most serious impediment to any African activity. ... [T]he price in human life was much lower." Figure 5 confirms the high mortality among European missionaries in Ghana prior to the introduction of quinine. 21 Simultaneously with quinine, the presence of European missionary staff expanded considerably, from less than 20 pre-1850 to circa 60 post-1870 (Fig. 5a). However, their numbers always remained below 70 because employing African converts as missionaries was a cost-efficient strategy in both the pre-and postquinine eras. Firstly, Africans acquired immunity to malaria during childhood (Curtin 1973, p. 197). 22 Secondly, their salaries were lower and they spread the gospel in the local vernaculars (Schlatter 1916;Graham 1976;Agbeti 1986, p. 57). 23 20 In 1848, the Director-General of the Medical Department of the British Army sent a circular to all West African Governors, advising the general use of quinine prophylaxis (Curtin 1961, p. 108). Fischer (1991, p. 73;99) explains that from 1845 on, Presbyterian missionaries received medical training before leaving for the Gold Coast and quinine was not missing in any of the missionaries' pharmaceutical orders. However, the proper application of quinine was a gradual process. Missionaries did not apply the correct dose during fever until the late 1850s (Fischer 1991, pp. 74-75). 21 Figure 5b shows that in Ghana the likelihood of European missionaries surviving more than three years during the pre-quinine era was about 30%, whereas in the post-quinine era it was about 80%. 22 African missionary mortality was significantly lower than for Europeans in both eras (see Fig. 5b). 23 By 1890, there were four African missionaries for every European. By 1918, Europeans constituted 2% and 8% of total Methodist and Presbyterian mission staff respectively (Parsons 1963, p. 4;Sundkler and Steed 2000, p. 717). Based on our data, the ratio of total mission stations to European missionaries increased substantially over time.

3
Financing the Mission. Protestant mission societies initially depended on the financial support from congregations and philanthropists in Europe and the US (Miller 2003;Quartey 2007). Cash-strapped mission committees relied on print propaganda, which sensationalized images of tropical missionary benevolence to elicit funding from Western readers (Pietz 1999;Maxwell 2015). Those donations paid for the missionaries' homeland training, the sea journey to Africa and initial set-up costs (Johnson 1967). Metropolitan funding remained limited however. In order to expand, the missionary budget had to be raised from within Ghana. Moreover, the mission societies' declared ultimate goal was to develop selffinancing African churches (Welbourn 1971).
African congregations contributed to the costs in various ways (Schott 1879, pp. 18-19). First, the bulk of the construction and operation of missions was financed by the local community, often in conjunction with local chiefs (Johnson 1967), who donated land, materials and labor to build the church and school (Williamson 1952;Summers 2016). Second, congregations were responsible for providing housing and food to the missionaries (Smith 1966, pp. 156-157); (Debrunner 1967, p. 249). Third, revenues were raised by donations from wealthier church members (Meyer 1999, p. 17), and more generally through Sunday offerings. Furthermore, school fees constituted another substantial part of the mission budget (Frankema 2012). For Africans, these sums were non-trivial, representing in 1926 about 20 days of unskilled wage labor. 24 Missionary expansion also became associated with trade and the cash crop economy: cocoa, kola, palm oil/kernels and rubber (Debrunner 1967, pp. 54, 132, 203). In particular, cocoa farming dramatically increased incomes from the 1890s onwards (Hill 1963;Austin 2003). By 1911, Ghana had become the world's leading cocoa producer. Ghanaians invested their cocoa revenues in their children's education at mission schools (Foster 1965;Meyer 1999). Debrunner (1967, p. 54) made it clear: "Cocoa money helped the African Christians to pay school fees and church taxes and to pay off old debts from the building of schools and chapels". Consequently, "Ghana Churches and the Christians became very dependent on cocoa for their economic support" (Sundkler and Steed 2000, p. 216). This also applies to other parts of Africa. Various Protestant mission societies established   trading companies that exported African cash crop produce and used their profits to sustain missionary activities (Johnson 1967;Gannon 1983). Catholic missions, in contrast, were less constrained as they relied on the financial backing of the Vatican and its missionary associations in Europe (Spitz 1924;Schmidlin 1933, pp. 560-564;Debrunner 1967 t is a dummy equal to one if there is a mission in cell c in year t and X c is the set of locational factors described before. As we control for missions in the first year of the period t-s ( M t−s ), the coefficients t show the long-difference correlation between the factors and missionary expansion in each period. Standard errors account for spatial correlation within 100 km (Conley 1999). 25 Table 1 presents for each period the coefficients of selected variables among all the variables included: 1751-1850 (col. (1)), 1850-1875 (2), 1875-1896 (3) and 1897-1932 (4). In the earliest periods, missions appear to have avoided high-risk malaria areas and settled at their port of entry, in close proximity to the coast (col. (1)-(2)). 26 It was only after the British had defeated the Ashanti Kingdom in 1896 that missionaries expanded northwards beyond the borders of the Gold Coast Colony (col. (4)). Next, while earlier missions expanded along 19th century trade routes (col. (1)-(2)), later missions opened in proximity to railroads (col. (4)). The negative correlations for navigable rivers during the early period (col. (1)-(2)) mirror the correlations for malaria (river floodplains provide breeding grounds for mosquitoes). Missions then concentrated in dense urban areas (col. (1)-(4)). Mission expansion appears to have followed urban population patterns of 1891, 1901 and 1931. Once urban demand was partly satisfied, missions appear to have spread into densely populated rural areas (col. (3)-(4)). Expansion also took place in cash crop growing areas (col. (2)-(4)). Finally, these descriptive results hold if we include 35 ethnic group fixed effects or exclude controls defined ex-post (Jedwab et al. 2019).
Overall, mission societies might have chosen, and might have been better received in, healthier, more accessible, and more developed areas ( R 2 = 0.50-0.61 in col. (2)-(4)). However, these results are not causal. Hence, our focus on malaria, railroads, and cash crops in the next section.

Malaria and missionary expansion
Difference-in-Difference (DiD). Section 3 described the substantial costs of European missionary mortality and how quinine reduced malaria death rates in the 1850s, after which the number of missionaries gradually increased. We test this more formally in Table 2. For 2091 cells c and 98 years t from 1800 to 1897 (N = 204,918), we regress a dummy if there is a mission in cell c and year t on the historical malaria index of cell c interacted with a post-quinine dummy (if year t is after 1850) while simultaneously including cell and year fixed effects. We choose the end of our third period -1897 -as the final year of the post-treatment window (1850-1897). To ensure a pre-treatment window of similar length we choose 1800 as our starting year. 27 Table 2 shows that missions expanded into higher-risk malaria areas after 1850 (col. (1)). The effect of quinine is strong: A one standard deviation in malaria is associated with a 0.18 standard deviation increase in the mission dummy in 1850-1897 (relative to 1800-1849). In col.
(2), we interact malaria with a dummy if year t is between 1800 and 1824. This separates the pre-treatment window into two sub-periods. We find no differential effect for malaria before 1850, implying parallel trends. The effects hold but are smaller when adding district (as of 1931; N = 38)-year fixed effects to compare neighboring cells within the same district over time (col. (3)). 28 Using the baseline DiD specification, we study the intensive margin. Web Appx. Table A3 shows that, conditional on having a mission (i.e. the extensive margin), higher-risk malaria areas had more missions, more main stations and more mission schools per cell by 1897. Moreover, denominational differences confirm that economic considerations mattered for missionary expansion. Mainline Protestants depended more on local contributions and valued entrepreneurship and education (Barro and McCleary 2017). As such, their missionaries had to be better trained, which made the issue of their low tropical life expectancy particularly acute. Consistent with these facts, Web Appx. Table A3 show stronger post-1850 effects for Mainline Protestants than for non-Mainline Protestants. Catholic missions started their conversion efforts only postquinine from 1880 onward. We therefore exclude them from this analysis.
Fuzzy Panel Event Study. Next, we show the effects of the introduction of quinine in a panel event study design. We restrict the analysis to the period 1800-1897, include cell fixed effects and district (1931)-year fixed effects, and aim to capture year-specific effects of historical malaria intensity before and after 1850 (1849 is the omitted year). Subfigure 6(a) shows the coefficients for each of the 15 years before and after 1850 as well as for the years "1835-" (we use a single dummy for all the years before 1835) and "1865+" (we use a single dummy for all the years after 1865). 29 European missionaries entered higher-risk malaria areas after 1855 and further expansion took place in 1857, 1860 and 1865+. There is no trend before 1850. As explained in Sect. 3, quinine had been available in Ghana as early as the late 1840s, but mission societies learned about the proper medication of quinine much later. The panel-event study is thus not sharp but "fuzzy". It also took mission societies time to respond and train and send enough European missionaries to Ghana. Spatial expansion was necessarily gradual. The sudden expansion in missions observed in 1855, the gradual expansion observed in 1855-1865, and the twice higher effects in 1865-1897 when the treatment is not "partial" anymore appear consistent with that.
Of course, as we include more post-treatment years, the panel-event study is less sharp, and the treatment variables may pick up developments other than the diffusion of quinine. There is a trade-off between only considering the partial treatment window and considering a longer window. In our case, we focus on the pre-1897 period, so before the rail, road, and cash crop eras.
Colonization may be a potential confounder. However, district-year fixed effects should take care of any spatial expansion of colonial rule. Furthermore, the boundaries of the Gold 27 For cells c and year t, the specification is as follows (using Conley standard errors): Mission Dummy c,t = + 1 Historical Malaria Index c *Dummy 1850−1897 + c + t + v c,t . 28 Results also hold when adding ethnic group-year fixed effects or cell-specific linear trends (Jedwab et al. 2018). 29 For cells c, districts d, and year t, the specification is as follows (using Conley standard errors): Mission We thus use 1850-1870 and 1871-1897. As seen in col. 6, the baseline effect is larger post-1871, consistent with Fig. 6a. If we only focus on the 1850-1870 period, a one standard deviation in malaria is associated with a 0.08 standard deviation increase in the mission dummy. The results with district-year fixed effects are weaker (col. 7) but the specification may ask too much of the data, removing much of the spatial variation in malaria. More generally, given the fuzziness of the event study, we cannot be sure that the estimated effects are fully causal, especially in the late 19th century. These results should thus be taken with caution.
Missionary Data. We estimate the same DiD model as before but we replace the dependent variable with a dummy indicating whether, for the years 1846-1890, mission stations were permanently inhabited or only monitored by European missionaries. Column (6) of Table 2 confirms a general increase in the number of missionaries in higherrisk malaria regions, which was partly driven by Europeans (col. (7)). Column (8) shows that quinine had a positive effect on the expansion of European permanent residences. Column (9) then shows that quinine had a stronger effect on African missionaries. Once standardized, the effect is about twice larger than for Europeans.
This suggests that the expansion into malaria areas was driven by African missionaries. This result may seem counter-intuitive since African missionaries had acquired natural immunity in some form. However, one needs to take into account specialization within mission societies. European missionaries were engaged in training and supervision activities. They mostly lived in coastal towns from where they would train African staff in church seminaries and routinely visit and supervise African missionaries in the hinterland. Second, in that period ordained priests were overwhelmingly European. Priests performed Christian rites that catechists were not allowed to do. These rites were important services provided at the stations. 30 Third, African catechists were cheaper to train and had a comparative advantage, for example due to their knowledge of local languages.
As such, European missionaries had to routinely visit African-run stations, often traveling to/through high-risk malarial areas. Before quinine, European missionaries would not have been able to expand their hinterland activities (via African personnel). With quinine, more Europeans lived on the coast from where they trained and supervised African catechists, hence leading to a spatial expansion of missionary activities, especially in malarial areas. 31 Finally, quinine likely also helped expansion in less-malarial areas. Thus, our effects may be downward-biased and we may under-estimate the contribution of quinine to missionary expansion.

Railroads and missionary expansion
Once the British had consolidated their control in 1896, they sought to build railroads to permit military domination and boost trade (Gould 1960;Luntinen 1996). By 1932, they had built three lines (see Web Appx. Fig. A3): (i) A western line (1901)(1902)(1903)(1904)(1905)(1906)(1907)(1908)(1909)(1910)(1911), which British capitalists lobbied for, to connect two gold fields in the interior to the port of Sekondi (Fig. 4a maps the cities mentioned here). Construction began in 1897 and its first segment was officially opened in 1901. The line was extended in 1903 to Kumasi, the capital of the annexed Ashanti Kingdom, to facilitate quick dispatch of troops. A small extension was then built in 1911; (ii) An eastern line (1909)(1910)(1911)(1912)(1913)(1914)(1915)(1916)(1917)(1918)(1919)(1920)(1921)(1922)(1923) aimed at connecting the coastal, colonial capital Accra to Kumasi. Other motivations were cited, including agriculture and the exploitation of gold fields; and (iii) A central line (1926)(1927) was built parallel to the coast to connect fertile land and a diamond mine. For the three lines, evangelization was never mentioned as a reason for construction nor missionaries acting as lobbyists.
Cross-Sectional Strategies. Five alternative routes were proposed but never built. We can address concerns regarding endogeneity by using these placebo lines as a placebo check of our identification strategy. Presumably random events such as a war and the retirement or premature death of colonial governors explain why the construction of these routes did not go ahead. 32 We run the same regression as in col. (4) of Table 1. The dependent variable is whether a cell had a mission in 1932. The variable of interest is a dummy for whether the cell is  Table A4 motivates the 0-30 km distance. When including dummies for whether the cell is within 0-10, 10-20, 20-30, 30-40 and 40-50 km from a railroad, we find an effect until 30 km only. The table shows that railroads built after 1897 had no significant positive effect on mission settlement before 1897, thus confirming parallel trends.
Panel A of Table 3, row 1 shows a baseline 0-30 km railroad effect of 0.162***. There is no effect of the 0-30 km rail dummy in the periods before 1897-1932 (rows 2-4). The main result is robust to: (i) Adding 34 ethnic group or 38 district (1931) fixed effects (rows 5-6); (ii) Confining the rail dummy to the more exogenous western line (row 7). Its goal was to connect a port, two mines, and the Ashanti capital Kumasi, without consideration for locations in between. Because we include dummies for whether there is a port, mine and large city as controls, we capture their effects, and identification relies on cells connected by the railroad by chance; (iii) Using cells within 0-30 km of a placebo line, for which no spurious effect is found (row 8); and (iv) Instrumenting the 0-30 km rail dummy by a dummy equal to 1 if the cell is within 30 km from the Euclidean minimum spanning tree between the nodes of the triangular rail network: Sekondi, Kumasi and Accra (see Web Appx. Fig. A3). We drop the nodes, which means that we rely on cells connected by chance (row 9; IV-F. = 39). Overall, the effect is strong: A one standard deviation in the rail dummy is associated with a 0.14-0.15 standard deviation in the mission dummy.
Timing of Rail Building. In Panel B of Table 3, for 2091 cells c in years 1897-1932, we study the effect of the 0-30 km rail dummy for cell c in year t on whether the same cell c has a mission in year t, while adding cell and year fixed effects. 34 Row 10 shows a strong effect (0.179***). Row 11 shows there is no effect when adding one lead of the rail dummy. Row 12 shows that the contemporaneous effect of railroads in t on missions in t is captured by the lag of the rail dummy, suggesting that missions followed railroads. Rows 13-14 show that results hold when adding ethnic group-year or district-year fixed effects, to compare connected and unconnected neighboring cells over time. Row 15 indicates that the baseline effect is unchanged if we include 2091 cell-specific linear trends. Overall, panel estimates are similar to cross-sectional estimates.
Results at the intensive margin point to the same direction (Web Appx. Table A5). When focusing on denomination-specific effects the various specifications show consistently stronger effects for Mainline Protestants than for non-Mainline Protestants and Catholics (see Web Appx. Table A6). This is in line with the fact that economic considerations mattered more for Mainline Protestants.
Panel Event Study. We show the effects of the 30 km rail dummy (based on the year of completion) on mission expansion in a panel event study framework. The first year of completion is 1901. Our last year of data is 1932. We thus restrict the analysis to 1870-1932 so that we observe about 30 years before and after the event. We focus on the 301 cells that were within 30 km from a railroad before 1932. Because dynamic effects cannot be identified without never-treated cells in the sample (Borusyak and Jaravel 2018), we add 69 33 For cells c, the specification is (Conley SEs): Mission Dummy c,1932 = + 0-30 km Rail Dummy c,1932 + X c B + v c . For the cross-sectional analyses on the effects of railroads we omit the controls capturing economic development (e.g., population and cash crops) circa 1932. Existing roads at that time were then purposely built as feeder roads for the railroads (Dickson 1969 placebo cells that were proposed to be connected (= be within 30 km from a railroad) but never were. Our sample includes 370 cells x 63 years = 23,370 cell-years. We include cell and district-year fixed effects. We then show the effects for the 14 years before and after (− 1 is the omitted year). We also include a dummy that captures all the years before year − 15 (incl.; "15-") and a dummy that captures all the years after year 15 (incl.; "15+"). Lastly, in settings with staggered adoption, the estimated effect is a weighted average of the effect of each subtreatment (Athey and Imbens 2022). Therefore, we investigate each line separately interacting each pre-and post-treatment dummy with a dummy for each line (Western, Eastern, Central). 35 Subfig. 6(b) shows the estimated coefficients. Both the Western and Eastern lines had strong effects on mission placement, especially after five years. The Central line does not show any positive effects (there is also a pre-trend), which is plausible. The Central line was built well after the other two lines (1925-1927 vs. 1901-1911 and 1909-1923). By then, mission societies had already established many missions along the Western and Eastern lines. In addition, the Central Line reached locations that were already quite connected. This line was built parallel to the coast, between the other two lines (Web Appx. A3), in areas that already had a few motor roads by the mid-1920s (Luntinen 1996). As such, the line barely reduced trade costs and the line was an economic failure. 36 To some extent, it is reassuring that we do not find an impact for this line.

Cash crops and missionary expansion
Export commodities were an important source of African income during the colonial era (Austin 2003). Ghana experienced various commodity export booms and busts as a result of new crop diffusion and changing world demand (Dickson 1969, pp. 143-178): palm oil (1860-1910s), rubber (1890-1910s), kola (1900-1920s), and cocoa (1900-1930s). 37 We explore the relationship between cash crop cultivation, as a proxy for African incomes, and the expansion of missions. The fact that each export boom took place at different times and in different areas facilitates identification. 38 Panel- Bartik. In the absence of data on annual crop production at the cell level, we study the reduced-form effects of a Bartik-type shift-share instrument predicting labor demand for each crop sector s in cell c and year t. Bartik IVs are used to generate exogenous labor demand shocks by averaging national employment growth across sectors using local sectoral employment shares as weights (Bartik 1991). We use a modified version of these: (i) We know the national export value of crop s (palm, rubber, kola and cocoa) in year t for the 1846-1932 period; (ii) We know in which cells c crop s was produced at 35 For lines r, cells c, districts d and year t, we use (Conley SEs): Mission Dummy c,t = + Σ s=15+ s=15− s,r + c + d,t + v c,t . 36 For example, Luntinen (1996, p. 120) writes: "From the Central Province cocoa was carried on road to the small beach ports, at the distance of fifty or sixty miles, whereas the railway to Takoradi via Huni Valley was more than one hundred miles long. In order to meet this competition the rates for cocoa were again drastically reduced in January 1930, but 'the results were most disappointing'. The railway management tried to explain their plight." 37 Export statistics for the period 1846-1932 (no data available before) confirm this (see Web Appx. Fig. A4). 38 We implement identification strategies due to the concern that missionaries locally promoted the cultivation of cash crops. However, in Ghana cash crop cultivation was essentially an African enterprise (Hill 1963;Austin 2014). any point during 1846-1932; (iii) We know the number of producing cells for crop s; (iv) Assuming that each producing cell was producing an equal amount, we predict the export value of crops s in cell c in year t; (v) Our exogenous measure of crop income in cell s and year t is then log export value of all crops s in cell c and year t; and (vi) When studying its effects on missions, we add cell fixed effects, which capture the time-invariant production dummies, and year fixed effects, which capture changing national export values. 39 Row 1 of Table 4 shows a large positive effect (0.028***) of log predicted cash crop export value at the cell level. In terms of magnitude, the effect is strong as a one standard deviation in the value of cash crops is associated with a 0.20-0.23 standard deviation in the mission dummy.
As explained by Goldsmith et al. (2020), shift-share instruments are not valid if the initial shares used are not exogenous. Row 2 shows that results hold if we construct the Bartik for the most important crops, palm oil and cocoa, using soil suitability dummies instead of production dummies (Globcover 2009). 40 Next, to ensure results are not biased by soil suitability having changed over time as a result of economic development, row 3 shows results hold when we use historical measures of soil suitability for palm oil, rubber and cocoa. 41 Goldsmith et al. (2020) also demonstrate that the Bartik estimator can be decomposed into a weighted combination of the estimates of all sectors, with the weights depending on how much each sector contributes to the identifying variation. Thus, Bartik estimators may hide the fact that they are driven by a few sectors only, and the specific exogeneity of their shares must be discussed. In our case, we rely on four crops -rather than say hundreds of industries -and rows 4-7 show results hold for each crop one by one. Furthermore, we build another Bartik for cocoa for which we have detailed geospatialized data on soil suitability. For each cell, we know the respective shares of moderately, highly and very highly suitable soils, as well as the relative average yields of these different types of soils. As seen in row 8 the estimated effect is similar to the baseline. 42 Moreover, no spurious effects are found when adding one lead of the Bartik (row 9). Instead, the effect of cash crops in t on missions in t is captured by a lag of the Bartik (row 10). This suggests missions followed cash crop incomes. Rows 11-13 show that results hold when adding ethnic group or district fixed effects interacted with year fixed effects or cell-specific linear trends.
Finally, the same Bartik analysis shows additional effects on the number of missions or the opening of main stations once we control for whether the cell had a mission (see Web Appx. Table A7). The same Web Appendix table also confirms that cash crop incomes have stronger effects on Mainline Protestant missions than on Catholic missions or other Protestant missions. This is another piece of evidence that economic considerations strongly mattered for missionary expansion.
Summary. A one standard deviation in malaria, railroads and cash crops is associated with a 0.08-0.18, 0.14-0.15 and 0.20-0.23 standard deviation increase in the mission 39 For cells c and year t, the model is (Conley SE): Mission Dummy c,t = + Ĉ ash Crop Export c,t + c + t + v c,t . 40 No soil suitability maps could be found for kola and rubber. 41 No maps could be found for kola. 42 Very highly suitable soils include class 1 and class 2 ochrosols and produce about 2000 pounds of dry cocoa per acre. Highly suitable soils include class 3 ochrosols and produce 1338 pounds per acre. Suitable soils include other types of cocoa soils and produce 223 pounds per acre. See Web Data Appx. for details on the sources used. dummy, respectively. However, there is a possibility that our estimates for malaria are not causal. Nonetheless, crudely adding these numbers, we obtain 0.42-0.56 (0.34-0.38 excl. malaria). Overall, the variables possibly explain a significant share of missionary expansion over time. The analysis also does not account for the effects of other locational factors or the general equilibrium effects of the three variables.

Dynamics of missionary expansion
This section highlights the dynamics of missionary expansion by documenting the changing locational characteristics in the stock of missions over time. We construct a measure that summarizes how "attractive" a location was to missionaries. More precisely, we regress the mission dummy in 1932, M c,1932 , on all the possible determinants of mission placement X c of Table 1. We then obtain the predicted probability M c,1932 = X cB , or locational score. We distinguish between four groups of cells in 1840-1932: 43 (i) cells with a mission in both t − 1 and t ("remains 1"); (ii) cells with no missions in t − 1 but a mission opening in t ("becomes 1"); (iii) cells with a mission in t − 1 that exits in t ("becomes 0"); and (iv) cells with no missions in both t − 1 and t ("remains 0"). Figure 7a plots a quadratic fit of the average score for those four groups.
The pattern suggests that the best locations received missions first, and that marginally less good locations were increasingly added to the existing stock of mission locations. Indeed, cells with surviving missions ("remains 1") rank consistently higher than cells that gain or lose missions ("becomes 1" or "becomes 0") and their scores significantly exceed those of the "remains 0" group. Scores of all the four groups decrease over time. Scores of the "becomes 1" group decrease, because less and less attractive mission locations are added over time. Scores of the "remains 0" group decrease, because switchers are among the best locations of the cells with no missions. 44 Results hold if (not shown): (i) We use the period 1751-1840 to estimate the coefficient of each factor and study predicted scores post 1840; (ii) The urban share in 1931 is the predicted variable.

Replication of the results for sub-Saharan Africa
We replicate the analysis of the determinants of missionary expansion for sub-Saharan Africa. Obviously, unlike for Ghana, we have to rely on limited and possibly mismeasured data.
Descriptive Analysis. In Table 6, for 203,574 cells in 43 countries, we regress a dummy if there is a mission according to the maps of Beach (1903), supposedly representing the year 1900 (col. (1)), or Roome (1925), supposedly representing the year 1924 (col. (2)), on various locational factors and country fixed effects. The regression is the same as for Ghana. Calculating Conley standard errors is too computationally intensive for that many cells. We therefore cluster standard errors at the district level (as of 2000). With 3284 districts, we get about 62 cells per district. 45 From the year of foundation reported for 83% of Protestant missions in Beach (1903), we construct a quasi-panel. 46 We then study in col.
(3)-(5) how the correlations vary across three periods defined around four turning points: 1792 (first year with a mission), 1850 (date when anti-slavery efforts intensified), 1881 (start of the Scramble for Africa), and 1900 (last year of mission data in Beach 1903).
Missionaries appear to have chosen locations with healthier environments (malaria, tsetse). Especially before 1850, missionaries seem to have avoided large pre-colonial cities, ethnic homelands colonized later, and Muslim centers, three measures of potential local resistance. 47 Transportation possibly played an important role: ports and coastal proximity may have facilitated initial access, while rivers, explorer routes, and railroads possibly enabled internal diffusion. Missionaries seem to have preferred large colonial cities and dense urban areas. We find positive correlations for slavery and cash crops. Finally, these correlations hold if we include 1158 country-ethnic group fixed effects or exclude controls defined ex-post (Jedwab et al. 2019).
Overall, missions were established in better areas. However, the adjusted R 2 are low, at 0.03-0.04 in columns (1)-(2) vs. 0.35-0.61 for Ghana (Table 1). This is due to two reasons. First, the locations of the missions mapped in Beach (1903) and Roome (1925) are mismeasured due to inaccuracies in the georefencing of missions by Cagé and Rueda (2016) and Nunn (2010), respectively. 48 When combining the cells into 3x3 cells, the adjusted R 2 increases to 0.15 (not shown). Second, for Ghana, we compiled a rich data set of controls, but such data do not exist for the whole of Africa.
Regarding causal effects, we do not know when quinine became the regular treatment in each country. However, we implement various identification strategies for railroads and cash crops. Table 5 shows the baseline effects of the 0-30 km rail dummy when including our control variables of Table 6. To avoid any over-controlling problem, we exclude controls that measure local economic development (population, crops, mining) c. 1900 (col. (1)) or 1924 (col. (2)). 49 Results hold if we apply the same cross-sectional identification strategies as for Ghana, whether: (i) adding ethnic (N = 1158) or district (as of 2000; 3284) fixed effects to compare neighboring cells (rows 2-3); (ii) using military or mining lines only (row 4), since their goal was to connect large pre-colonial cities or mines to a port without consideration for locations in between. Given the controls (we re-add mining), we capture the independent effects of locations that mattered 45 The main results (Table 5, Table 6 and panel B of Table 8) hold when using larger grid cells to cluster standard errors, for example 10x10 mega-cells (N = 94 cells per mega-cell) or even 15x15 mega-cells (N = 199 cells) (not shown). 46 To our knowledge, no one has ever used the panel aspect of this data. Panel data does not exist for Roome (1925). 47 Distance to Muslim centers has a negative and almost significant correlation in col. (5) (which captures missionary expansion between 1881 and 1900), suggesting a more nuanced picture for Islam. Large cities in 1400 were mostly located in Muslim areas and their negative coefficients could reflect the fact that they were major Muslim centers. 48 For example, for 109 missions reported in both Beach and Roome and digitized by Cagé and Rueda (2016) and Nunn (2010) respectively, we found an average distance of 2 cells between their georeferenced locations. 49 For cells c and countries n and each separate year t = {1900; 1924} , the specification is as follows (standard errors clustered at the district level): Mission Dummy c,n,t = + t 0-30 km Rail Dummy c,n,t + X c,n,t B t + n + v c,n,t . for military domination or mining, and identification relies on cells connected by chance; 50 and (iii) instrumenting the rail dummy by a dummy if the cell is within 30 km from the Euclidean minimum spanning tree between the capital, largest and second largest cities c. 1900, while simultaneously dropping these cities (row 5; IV-F.=49; 82). We also find no spurious effects when using placebo lines planned c. 1916-1922 but never built (row 6).

Railroads. Row 1 of Panel A in
For the panel analysis of the effects of railroads, we restrict the sample to the 1885-1900 period. Until 1885, only South Africa had railways and rail construction in Africa only expanded from the late 1880s on. Unlike for Ghana, we only know the opening year of the railroad, not the announcement year. In addition, for Africa the exact foundation year of each mission and the year in which each railroad line was completed are both likely to be mismeasured. We therefore use 5 year panel data. 51 Lastly, to reduce the number of observations, we restrict the sample to 103,866 cells in 15 countries with railroads opened at any point during 1885-1900, thus obtaining 103,866 x 4 = 415,464 observations. 52 Panel B shows the baseline effect of row 7 is positive and significant. Row 8 shows there is no effect when adding one lead of the rail dummy. However, row 9 indicates that the contemporaneous effect of railroads in t on missions in t is not captured by the lag of the rail dummy in t − 5 . Next, Rows 10-11 show that point estimates remain similar when adding ethnic group or district fixed interacted with year effects, to compare connected and unconnected neighboring cells over time (effect not significant at 10% with district-year fixed effects). Finally, including 103,866 cell-specific linear trends is unfortunately too computationally intensive. 53 Subfigure 7(b) shows the coefficient of correlation between a dummy if there is an atlas mission in the cell in 1900 in Beach (1903) or 1924 in Roome (1925) and the census mission dummy in year t . 50 Using the raw railroad data from Jedwab and Moradi (2016), we verify that Christianization was never mentioned as a reason for the 200 or so railroad lines that were built in sub-Saharan Africa (not shown). 51 The same reason also prevents us from performing an event study analysis. 52 For cells c, countries n, and years t, the model is (standard errors clustered at the district level): Mission Dummy c,n,t = + 30 km Rail Dummy c,n,t + c + t + n,t + v c,t . The 15 countries include 3 Central African, 5 East African, 3 Southern African and 4 West African countries, respectively. 53 Results generally hold if we add linear trends for 0.75 x 0.75, 0.60 x 0.60, 0.50 x 0.50, 0.40 x 0.40, 0.30 x 0.30 and 0.25 x 0.25 degree mega-cells (not shown, but available upon request). The corresponding number of mega-cells increases from 2110 to 17,209 for railroads and from 1523 to 11,817 for cash crops (see the next analysis). Table 1 Correlates of missionary expansion, long-diff., selected factors For 2091 cells and each period [t-s;t], we regress a dummy if there is a census mission in year t on a dummy if there is a census mission in year t-s and various locational factors. We do not report the coefficients of land area, rainfall, altitude, ruggedness, soil fertility, four dummies if the cell contained a large pre-colonial city circa 1800, the largest or 2nd largest city in 1901, a headchief town in 1901, or the central location of a mine in 1932, four dummies if the cell was within 10 km from a road in 1930, within 50 km from a slave market c. 1800, within 50 km from a rubber plantation in 1900-1936, or Table 2 Malaria and missionary expansion, panel Overall, cross-sectional and panel estimates are similar, and a one standard deviation in the rail dummy is associated with a 0.01-0.09 standard deviation in the mission dummy. The magnitude of the effects is smaller than for Ghana (0.14-0.15). However, mission and railroad openings are mismeasured for the whole of Africa making the comparison with the Ghana results less relevant.
Lastly, using the same cross-sectional or panel identification strategies, we find stronger effects for Mainline Protestants than for Catholics or other Protestants (Web Appx. Tables A8-A9). Thus, economic considerations also mattered more for Protestants when studying the determinants of mission expansion for the whole continent.    54 However, we only have panel data on missions before 1900. In addition, for the seven crops studied for Africa (cocoa, coffee, cotton, groundnut, palm oil, tea and tobacco), we could only find pre-1900 export values for 20 countries (N = 70,546 cells) c. 1850, 1860, 1870, 1875, 1880, 1890 and 1900 (N = 493,822). Using this data, we use the same method as for Ghana to construct the log export value of all seven crops s in cell c and year t. We then regress a dummy if there is a mission in cell c in t on log predicted cash crop export value in cell c in t, while adding cell and year fixed effects as well as country-year fixed effects. 55 Row 1 of col.
(2) in Panel B of Table 5 shows a strong positive effect (0.001***) of log predicted cash crop export value. This effect is driven by one crop, palm oil, as other crops were not important at that time, with the exception of groundnuts that overlap with core Muslim areas. Next, no spurious effects are found when adding one lead of the Bartik (row 8), but the effect of cash crops in t on missions in t is mostly captured by a lag of the Bartik (row 9). Rows 10-11 show results hold when we add ethnic or district fixed effects interacted with year fixed effects.
In terms of magnitude, however, the effect is not that strong: A one standard deviation in cash crop value is associated with a 0.02 standard deviation in the mission dummy Table 4 Cash crops and missionary expansion, panel For 2091 cells c and 87 years t (1846-1932), we regress a dummy if there is a mission in cell c in t on log predicted cash crop export value in cell c in t (cell FE and year FE included). Rows 2-8: Using alternative ways to construct the predicted value of cash crop exports. 100 km Conley SE's (clust. at district (1931) level in row 13): * p<0.10, ** p<0.05, *** p<0.01. See Web Appendix for data sources. There is a concern that missionaries locally promoted crop cultivation in other countries than Ghana. However, in most of colonial Africa, cash crop cultivation was driven by African initiative (Austin 2014)not mission societies. 55 For cells c, countries n, and years t, the specification is (standard errors clustered at the district level): Mission Dummy c,n,t = + Ĉ ash Crop Export Value c,n,t + c + t + n,t + v c,n,t .
(vs. 0.20-0.23 for Ghana). However, both mission and cash crop values are mismeasured for Africa.
Finally, the Bartik analysis shows an additional effect on the number of missions once we control for whether the cell had a mission (see Web Appx. Table A9). The table also confirms that cash crop incomes have stronger effects on Mainline Protestant missions than on other Protestant missions.

Economic development and the adoption of a new religion
To summarize, missions were established in more developed areas. We now discuss why increased missionary supply was met by African demand for Christianity at these locations.  (1)) or 1924 (col. (2)) on the 0-30 km rail dummy and the controls of Table 6 1860, 1870, 1875, 1880, 1890 and 1900 (N = 493,822), we regress a dummy if there is a mission in cell c in t on log predicted cash crop export value in cell c in t (cell FE, year FE, and country-year FE included). See text for details on each regression. Robust SE's clustered at the district (2000) level: * p<0.10, **p<0.05, *** p<0.01. See Web A ppx. for data sources.  First, our results do not exclude the possibility that it was the poorest individuals in the richest places who converted to Christianity (Hastings 1994;Maxwell 2016). By the mid-19th century, Christianity had broadened its appeal among the commercial elite, such as cash-crop farmers and merchants (Debrunner 1967). Indeed, missions required financially capable members to contribute to church activities (see Sect. 3).
Second, Barro and McCleary (2003) argue that if participating in religious activities increases wages, for example because religion and human capital are complements, growth and religiosity go hand in hand. Missions supplied the bulk of formal education (Frankema 2012), which commanded a wage premium and facilitated occupational mobility (Ekechi 1971;Frankema and Van Waijenburg 2019;Meier zu Selhausen et al. 2018). However, the complementarity between Christianity and schooling weakened over time. In Ghana, since the 1870s missions increasingly opened without government approved schools. In 1932, only one out of six missions had a school (Fig. 3). After World War II states expanded the supply of state schools and missions lost their monopoly on schooling. Hence, schooling cannot fully account for the appeal of Christianity.
Third, Christianity disrupted the monopoly, and spread at the expense of, African traditional religions. It has been argued that such religions constrained individual ownership and restricted the pursuit of self-interest (Pauw 1996;Alolo 2007). Christianity, in particular the Protestant denominations, are more capitalist in nature. There are also spiritual needs in a world where established systems of meaning became increasingly disrupted by changing social and economic circumstances, including new technologies (e.g., railroads and steamships). Africans sought a measure of conceptual control over these forces by turning to the new ideas offered by Christianity (Maxwell 2016). Thus, Christianity may have been for converts a more this-worldly religion.

Implications for the study of long-run economic development
Our results have several implications for the study of the long-run effects of colonial missions. Most studies retrieve mission location data from a source different from ours: Historical mission atlases. Using Ghana as an example, we first show that atlases select the most important missions (e.g. early, main or European residence stations). Second, we will examine the long-term economic and non-economic effects of missions when relying on atlas missions or census missions and the standard controls used in the literature or our controls, and compare different types of missions. For 203,574 cells in 43 sub-Saharan African countries, we regress a dummy if there is an atlas mission in Beach in 1900 (col. (1)) or in Roome in 1924 (col. (2)) on various locational factors (country fixed effects always included). We do not report the coefficients of land area, rainfall, altitude, ruggedness, soil fertility, a dummy if there was a large pre-colonial city in the cell circa 1800, a dummy if the cell is within 10 km from a lake, a dummy if the cell is within 50 km from a mine in 1900 or 1924, a dummy if the main ethnic group in the cell before colonization according to Murdock (1959) does not have information on state centralization, slavery and polygamy, a dummy if we know the year of the anthropological study used by Murdock (1959) to create his data, and a dummy if that year precedes 1900year precedes /1924 (3)-(5): For Beach and each period [t-s;t], we regress a dummy if there is a mission in t on the locational factors and a dummy if there is a mission in t-s. Robust SE's clustered at the district (2000 boundaries) level: * p<0.10, ** p<0.0 5, *** p<0.01. See Web Appx. for data sources.

Mission Atlases and non-classical measurement error in missionary activity
In the literature, two mission atlases feature prominently: Beach (1903) and Roome (1925). 56 Yet, we find that atlases significantly underreport missions. For Ghana, atlases show far fewer missions than census returns: 26 vs. 304 (91% are missing) in 1900 (Beach 1903) and 23 vs. 1213 (98%) in 1924 (Roome 1925) (see Fig. 2a-b). For Africa, the extent of misreporting is of similar scale: Beach (1903) and Roome (1925) omitted 93% and 98% of missions (see Fig. 1b). 57 If measurement error is classical, estimates of the contemporary effect of missions will be downward biased. However, this is far from certain. Fahs (1925) pointed out that atlases overwhelmingly plot residence stations of European missionaries. In this case, measurement error is non-classical. 58 We can investigate the measurement error thanks to our unique data on census missions and their characteristics. Examining the geographic distribution of missions in Ghana (see Fig. 2a and b) reveals two stylized facts. First, atlases miss most hinterland missions. Second, Roome (1925) does not capture the exponential growth of missions between 1900 and 1924. We study this further by calculating the coefficients of correlation between a dummy equal to one if there is an atlas mission (1900,1924) and a dummy equal to one if there is a census mission in a cell for each year . Figure 7b shows high correlations for earlier years (0.8 for Beach (1903) and 0.4 for Roome (1925)). 59 Atlases thus best represent missions circa 1850 (that survived) instead of 1900, or 1924. Atlas missions also differ qualitatively. In Table 7 we regress a dummy equal to one if there is an atlas mission on cell-level mission characteristics derived from our data for the years 1900 and 1924. The correlations show that atlases capture not only early missions (col. (2) and (7)), but also main stations as well as missions with schools (col. (3) and (8)) and the residences of European missionaries (col. (4) and (9)). Some characteristics are then correlated with each other (col. (5) and (10)). 60 The R 2 is then 0.64-0.45 (0.79-0.55 with 2x2 cells). Overall, atlases selected "better" missions. But, do atlases at least reliably select the "best" missions that received major investments? The answer is no. If we define the "best" missions as the ones that were main stations, where Europeans resided and which had a school, atlases still miss 68-83% of them. 56 Based on Web Appx. Table A1, about 30 studies rely on either Beach (1903) or Roome (1925). 57 The extent of omissions is relatively similar across movements and denominations: Beach (1903) misses 89% of Methodist missions and 92% of Presbyterian missions whereas Roome (1925) misses 97% of Protestant missions and 99% of Catholic missions. Other atlases compare as follows: 25 vs. 432 Protestant missions (i.e., 94% of missions are missing) in 1908 (Bartholomew et al. 1911) and 6 vs. 49 Catholic missions (88%) in 1913 (Streit 1913). 58 Fahs (1925, p. 271) wrote that "[...] mapping the Christian advance which puts a red underline under some place-name to indicate the residence of a British or Continental or American missionary, but does not indicate where an Azariah or a Kagawa serves his people, fails to give needful perspective." Mission atlases thus ignore the very significant contribution of African missionaries to Christian expansion (Meier zu Selhausen 2019). 59 When combining cells into 2x2 cells to account for possible inaccuracies in the geocoding of atlas missions, the correlation of earlier years is 0.9 for Beach (1903) and 0.7 for Roome (1925) (not shown). 60 The variables still describe different mission characteristics. Correlation coefficients between the eight variables average 0.45 (28 combinations), so much lower than 1.

missions, Ghana
For 2091 cells c, we regress a dummy if there is an Atlas mission in cell c in Beach (1900) or Roome (1924) on the mission characteristics in the same cell the same year (1900 or 1924, unless

Non-classical measurement error and long-term economic effects
We now study the long-term economic effects of missions, depending on whether we rely on census missions or atlas missions. We run a regression correlating a dummy equal to one if the location had a mission ( M c,g ) in 1900 or 1924 and economic development today ( D c,g,today ): Given the lack of data on incomes at the local level, we use log average night light intensity in 2010 -i.e., the total sum of night lights divided by area -as development indicator ( D c,g,today ). To avoid issues related to top-coding, we use the radiance calibrated version of this data (NOAA 2012), which records levels of luminosity beyond the normal digital number upper bound of 63. 61 For both Ghana and Africa, we examine how varies if we include: (i) No controls X c,g ; (ii) The "standard" controls commonly used in the literature on Africa (Std). In particular, we merge the lists of controls from Nunn (2010Nunn ( , 2014 and Cagé and Rueda (2016). Small variations over this set of controls have been used, for example, by Fenske (2015), Okoye (2021), Henn et al. (2021) and Cagé and Rueda (2020); 62 (iii) Our full set of possible determinants from Table 1 (Ghana) or Table 6 (Africa) (Ours); and (iv) Our full set of controls as well as controls for the type of mission. The latter include (only available for Ghana; from Table 7) whether the mission was created in 1751-1850 or 1851-75, is a main station in t, has a school in t, a European missionary lived there in 1846-90, and a European missionary frequently visited the mission in 1846-90. When controlling for the type, the effect is estimated for average missions, which minimizes external validity concerns. Table 8 presents the results. For Ghana, using the census mission dummy (row 1 of Panel A), we find large unconditional effects of 1.05***-0.84*** (col. (1) and (6)). Adding the standard controls reduces the estimated effect to 0.78***-0.73*** (col. (2) and (7)). An early colonial mission thus (relatively) increases night lights by 117-107%, or 0.35-0.49 if expressed in terms of standard deviations. 63 Next, with the imperfectly measured Africa controls (from Table 6), the effect is further reduced by almost 50%, to 0.41***-0.38*** (col. (3) and (8)). If we use our extended set of Ghana controls, the effect is further divided (1) D c,g,today = a ��� + M c,g + X c,g + � g + w c 61 Nightlights have become a standard measure in the stream of literature that studies development at the local level in Africa (Michalopoulos et al. 2013(Michalopoulos et al. , 2014Alesina et al. 2016;. 62 Nunn (201062 Nunn ( , 2014 4) and (9)). 64 Lastly, also controlling for the type of missions, the effect decreases to 0.01-0.17*** (col. (5) and (9)). One standard deviation in the mission dummy is then associated with a 0.00-0.12 increase in the standard deviation in log night lights. 65 Therefore, we find that the average (census-based) missions had possibly no or weak effects on contemporary development in Ghana. 66 However, atlas missions, because they capture important missions, could still have strong long-term economic effects, especially if we do not control for their potentially endogenous placement and type. In row 2 of Panel A, we use the same model as in row 1 but, instead of the census mission dummy, we include an atlas mission dummy and a nonatlas mission dummy (defined as census mission dummy -atlas mission dummy). Without controls, we obtain twice larger point estimates for atlas missions. The effect of atlas missions is still much larger, by 0.74*** (not shown), if we add the standard controls (col. (2) and (7)). Once we add our controls, including the type of missions (col. (5) and (10)), we again find nil or small effects for non-atlas mission (0.01-0.17***). Although atlas missions have a stronger effect, the difference with respect to non-atlas missions is barely significant (0.43*; col. (5)) or insignificant (0.15; col. (10)). 67 More generally, using our full controls rather than the standard controls decreases the effects of atlas missions for Ghana by (1-0.44/1.39)*100 = 68% for Beach and (1-0.33/1.43)*100 = 77% for Roome. 68 In row 3, Panel A of Table 8, we study the effects for Ghana when we only include the atlas mission dummy, in order to mimic the regressions we are running for Africa next. When doing so, the effect of atlas missions is slightly lower than when controlling for the non-atlas missions. However, the same patterns as before emerge. Using our full controls rather than the standard controls decreases the effects of atlas missions for Ghana by 65-85%.
Finally, in Panel B, we examine for 43 African countries the effects of atlas missions on contemporary night lights. The unconditional effect is 0.70*** (col. (1) and (6)). Adding the standard controls reduces the effect to 0.60-0.59*** (col. (2) and (7)). Thus, a colonial mission increases night lights by 82-80%, or 0.10-0.13 if expressed in terms of standard deviations. Adding our Africa controls then lowers the effect to 0.37-0.40*** (we found 0.43-0.59*** for Ghana, see row 3), or 0.07-0.08 if expressed in terms of standard deviations. If we also add two dummies for whether a Beach mission already existed pre-1850 or in 1850-1881 -our only type control for Africa -, the effect for the year 1900 decreases further to 0.28*** (not shown), which is about half the effect with the standard controls. However, the Africa controls are not as precisely measured as those for Ghana. 69 We also just showed how using better controls reduces the effects of atlas missions by 65-85% for Ghana. Were this the case for Africa, the standardized effect would only be 0.02.

Extensions
Measurement Error in the Controls. If the best missions are selected, the long-term effects of missions will be over-estimated. The same is true if only mismeasured controls of mission placement are available. For example, for the Africa regressions, we have no choice but to use the Murdock (1967) map to obtain data on state centralization, slavery and polygamy. As such, the reported estimates when using our full controls likely remain upward-biased. As we just showed, the better measured Ghana controls, which do not rely on continental maps, reduce the effects twice as much as the Africa controls. Now, if we drop the Murdock controls and instead add ethnicity, country-ethnicity, province or district (as of 2000) fixed effects, the Africa effects slightly decrease to about 0.30***-0.35*** (see Panel A of Web Appx. Table A10). However, the mean province is large in Africa, with an area of 228,000 sq km vs. 197,000 on average for U.S. states. The mean ethnic territory (141,000) and the mean district (84,000) are then 1.4 and 2.3 times smaller than the average U.S. state, which is still large. For Ghana, when using the standard controls (see Panel B of Web Appx. Table A10), adding province, ethnicity or district fixed effects only reduces the 1900-based effect from 0.78*** to 0.74***, 0.58*** and 0.51***, respectively, vs. 0.01 with our full controls. The 1924-based effect only decreases from 0.73*** to 0.68***, 0.60*** and 0.50***, respectively, vs. 0.17*** with our full controls. Thus, the standard controls do not capture the fact that missionaries went to more developed areas within these spatial units.
Spillovers. Our main unit of analysis are grid cells of 0.1 x 0.1 degrees (12 x 12 km). The minimal distance from the centroid of the cell to its edge is 6 km, or a one-hour walking distance. Due to the lack of transportation technology, missions were likely to only influence local development during the colonial era. However, spillovers remain a possibility. We use the same model as before but also add the log of the minimal distance to a cell with a mission. When doing so, we find strong spillover effects with the standard controls but no effects when adding our controls (see Panel A of Web Appx. Table A11). Alternatively, if we use dummies if the cell is 0-10, 10-20 or 20-30 km from a cell with a mission, we find strong effects that are decreasing with distance with the standard controls, and no spillover effects for the year 1900 and a barely significant 10-20 km effect for the year 1924 when using our controls including for the type of missions (Panel B).

Effects of Missions in More Developed
Areas. The fact that mission societies disproportionately established mission stations in more developed areas makes it even more surprising that we find nil or small long-term effects of colonial missions. For example, railroads and cash crops might have planted the seeds of future economic development, thus creating complementarities raising the long-term economic effects of missions. To examine how strong such complementarities could be, we use the model where we study the effects of the census mission dummy on log night light intensity in 2010. We then interact the census mission dummy with a dummy if the cell was ever within 30 km from a 1932 railroad as well as a dummy if the cell was potentially producing one of the four dominant cash crops during the colonial era (cocoa, kola, palm oil and rubber). As seen in Web Appx. Table A13, once we include our controls, we find no significant positive differential effects of missions in rail or cash crop areas vs. missions outside these areas.
Effects of Missions Depending on Their Type. Our analysis suggests that the average census mission might have had limited long-term economic effects. Yet, it could be that only census missions of a particular "importance" indeed mattered. We can examine this hypothesis for Ghana thanks to our unique data. We use the model where we study the effects of the census mission dummy on log night light intensity in 2010. However, we now also add variables capturing the type of missions. 70 First, in columns (1) and (5), Table 9 we explore whether contemporary effects vary with the foundation year of the mission. With the standard controls (Panel A), we find strong positive effects for missions created early (in 1751-1850 or 1851-1875). With our controls (Panel B), the effects of early missions are about halved. It might be the case that early missions were able to accumulate growth effects over a longer period of time. Alternatively, early missions may have received more investments, i.e. they could have been main / school / European missions. In columns (2) and (6), we study the effects of the log number of missions -i.e., the intensive margin -and whether there was a main station or a school. With our controls (Panel B), we only find a significant effect of the log number of missions in 1900. In columns (3) and (7), we find strong effects of more European missions. The effects are about halved with our controls. Finally, if we add all characteristics and our controls (Panel B), the only effects that remain significant are for missions 'created 1851-1875' and the log number of missions in 1900. Hence, controlling for early creation and density, we find no effects for main / school / European missions, which is surprising since these are characteristics that one might have expected to matter. 71 Other Strategies. We find similar non-results on long-term economic development if we employ other strategies. In Panel A of Web Appx. Table A16, we compare cells with Protestant missions, Catholic missions, and no missions (Nunn 2014;Cagé and Rueda 2020). With the standard controls, we find higher point estimates for Protestant missions. With our controls, point estimates are much lower, and their difference is reversed. We also compare cells where missions survived with cells where missions were abandoned as in Valencia Caicedo (2019a). For the year 1924, we find strong effects with the standard controls and no effects with our controls (Panel B). Table 9 Long-term effects of missions depending on their type, Ghana For 2069-2091 cells in Ghana (note that the identity of the missionary is missing for 2091 -2069 = 22 obs.), we regress our main measure of local economic development, log mean night light intensity in the cell in 2010, on the census mission dummy and mission characteristics in the cell, all defined in the same year (1900 / 1924 service sectors (Gollin et al. 2016) -(row 5), or skilled occupations 75 (row 6), we find that our controls reduce the long-term effects of missions on these by 5-9, 5-8 and 4 times, respectively. In terms of standardized effects, we find 0.00-0.12 for night lights, 0.02-0.03 for urbanization and 0.05-0.09 for employment. Row 7 of Web Appx. Table A14 uses as the dependent variable the population share of Protestants and Catholics in 2000. Our controls barely change the long-run effect of missions on Christianity. The standardized effects are 0.08-0.15, higher than for economic outcomes. If we use the number of years of education (row 8) or the primary completion rate (row 9) for individuals aged 18 or more in 2000, our controls reduce the effects by about 3 times relative to the standard controls. If we use height-or weight-for-age z-scores for children below 5 in 1993-2017 (source: DHS), our controls reduce the effects by 1.3-2.0 times relative to the standard controls (rows 10-11). Thus, our controls reduce the effects by more than they do for Christianity, and less than they do for economic outcomes. Also, the standardized effects are 0.10-0.17 for education and 0.04-0.09 for health, similar or higher than for economic outcomes. 76 Our controls reduce the effects of missions on completed fertility (row 12) and net completed fertility (row 13) -for 35-49 aged women in 2000 -by 2-4 times and 2.5-5 times relative to the standard controls, respectively. Completed (child) mortality then does not vary with missions (row 14). In terms of standardized effects, with our full controls, we find − 0.03/− 0.08 for net fertility, which is similar in size to economic outcomes.
Overall, we find stronger, but still weak, effects on education, health, and fertility, implying that missions might have affected development outcomes without spurring economic development. 77 Way Forward and External Validity. Our paper proposes a number of potentially important predictors for colonial missions. In principle, those predictors can be used in matching procedures or as instrumental variables for the placement of missions. However, our set of determinants is not exhaustive. Hence, the exogeneity assumption may be violated. In addition, our results are specific to the African context. In other contexts such as in Latin America, missionaries activities and colonial infrastructure and economic investments were often jointly determined by colonial governments and the Catholic Church (Valencia Caicedo 2019a, 2019b). Christian missions were then more likely "proxies" for colonization in general, making it more difficult to isolate the effects of missions themselves. 78 76 In Ghana, missions did not invest into health care. Elsewhere, missions were often the only health care providers during the colonial era (Cagé and Rueda 2020). Still, few missions offered health services. Of all missions only 2% and 4% of Catholic (in 1910) and Protestant (in 1923) missions respectively operated a hospital (Streit 1913;Beach and Fahs 1925). Such investments were costly to missions because they were free of charge for most patients (Doyle et al. 2020). 77 As shown in Web Appx. Table A15, we do not find any significant positive effects of main / school / European missions on educational outcomes today. We also find nil or contradictory effects for health outcomes. Finally, whether main stations or mission schools significantly reduce net fertility depends on whether we consider 1900 or 1924. 78 In Ghana and other British colonies in Africa, the activities of the various mission societies and colonial infrastructure and economic investments were not coordinated. In Belgian, French, Portuguese and Spanish colonies, coordination was also minimal, thus more similar to the Anglo-Saxon model of privatized missionary markets than to the Spanish/Portuguese model of "bundled" colonization/missionization historically observed in Latin America. 75 These include managers, professionals, technicians, and clerks.

Conclusion
In this paper, we studied Christian missions, a major human capital promoting institution that emerged during colonial times. We argued that in order to estimate colonial missions' long-term economic impact, one needs a better understanding of the dynamics and economics of missionary expansion. Focusing on Ghana and Africa, we described how missionaries privileged healthier, more accessible, and richer places first. In particular, we showed how -with quinine -malaria became less of an impediment, and how transport infrastructure and agricultural export commodities attracted missions. The diffusion of Christianity, education and culture depends on spatial patterns of economic development, even at low income levels. We then showed how non-classical measurement error in mission atlases and omitted variable bias are concerns that could lead to an overly optimistic account of the importance of missions for contemporary development.