Bilateral international migration flow estimates for 200 countries

Data on stocks and flows of international migration are necessary to understand migrant patterns and trends and to monitor and evaluate migration-relevant international development agendas. Many countries do not publish data on bilateral migration flows. At least six methods have been proposed recently to estimate bilateral migration flows between all origin-destination country pairs based on migrant stock data published by the World Bank and United Nations. We apply each of these methods to the latest available stock data to provide six estimates of five-year bilateral migration flows between 1990 and 2015. To assess the resulting estimates, we correlate estimates of six migration measures from each method with equivalent reported data where possible. Such systematic efforts at validation have largely been neglected thus far. We show that the correlation between the reported data and the estimates varies widely among different migration measures, over space, and over time. We find that the two methods using a closed demographic accounting approach perform consistently better than the four other estimation approaches.


Background & Summary
International migration is becoming an ever more important component of population growth 1 , a driver for socio-economic change 2,3 and a topic for policy debate 4 in many countries. Good data on international migration are becoming more important as a means to monitor international development agendas and agreements such as the Sustainable Development Goals and Global Compact for Safe, Orderly and Regular Migration 5 . The estimates of migration flows are of potential use to those studying global or regional migration systems, such as economists and sociologists who study migration drivers and networks, demographers and climate change scientists for base line data for projections, and epidemiologists who study the role of human movements in the transmission of infectious diseases.
However, reliable international migration flow data, to measure the movement of people between countries over a given period, are scarce 6 . In many countries, migration flow data are not collected because data collection systems are expensive. Countries that collect migration flow data often use widely varying criteria for defining migrants (such as their length of stay and purpose for moving). These definitional differences prevent detailed cross-national comparisons. Combined, missing and non-comparable data prohibit accurate measurement of patterns and trends in global migration flows.
The lack of accurate migration flow data has motivated recent efforts to develop synthetic estimates for bilateral flows between all pairs of countries [7][8][9][10][11] . Each of these methods relies on bilateral migrant stock data. Migrant stock data measure the number of migrants at a particular point in time residing in each country by their country of birth. Such measures are easier to collect than flows, requiring only birthplace questions from a country's census, population register or administrative data collection system.
Currently proposed methods to estimate migration flows from changes in migrant stocks fall into three groups. The first set of methods uses the differences in successive bilateral stocks for a given pair of countries to estimate the corresponding migration flows in each direction. These approaches are commonly used by economists and regional scientists who use these migration flow estimates as a dependent variable in their explanatory models. The second group of methods uses migrant stock data to estimate migration flow rates, which are then multiplied by a population at risk. The third group of methods frames the changes in migrant stocks as the residuals in a global demographic account. These so-called "demographic accounting" approaches estimate migration flows to match increases or decreases in the reported bilateral stocks of migrants, and births and deaths during the period.
Thus far, the different estimation methods have used different migrant stock data as inputs, prohibiting a systematic comparison. Many have not used newly published migrant stock and demographic data. Further, systematic comparisons of different sets of estimates are lacking. Previous applications of the estimation methods typically do not compare outputs with either results from alternative methods or reported data. Where comparisons are undertaken, they are performed for only a small fraction of pairs of countries. This lack of comparison has made it difficult for users to know which estimation method produces values most relevant to their application.
To supply migration flow estimates based on the most recent migrant stock data and systematic comparisons of different methods, in this paper we apply multiple methods for estimating migration flows from stocks to produce estimates of the bilateral international migration flows between all pairs of 200 countries for five five-year periods from mid-year (July 1) 1990 to mid-year (June 30) 2015, using the most recent set of bilateral migration stocks published by the United Nations. We also propose a set of validation exercises based on a number of commonly used migration measures and reported migration flow data. These comparisons highlight the strengths and weaknesses of different sets of estimates.
An overview of our study design is shown in Fig. 1. The following sections outline the estimates of bilateral international migration between all pairs of 200 countries for five-year periods during 1990-2015 based on six estimation methods.

Methods
In this section we outline six methods. We use the same notation and a simplified hypothetical data set, introduced next, to allow for easier comparison between the methods. Details on each estimation method are provided towards the end of this section, before an overview of the code used to perform our calculations. contingency tables. Bilateral migration data are commonly represented in square contingency tables with off-diagonal entries (the cell in the intersection of row i and column j, i ≠ j) containing a migration measure from origin i to destination j. Consider two hypothetical migrant stock tables in consecutive years (t and t + 1) in Table 1. Regions A to D represent places of birth in the rows, and place of residence in the columns. The off-diagonal entries represent the number of foreign-born migrants in each area of residence, while the diagonal entries contain the number of native-born residents.
Stock differencing. Stock difference, drop negative. One way to estimate bilateral migration flows from migrant stock tables is to find the difference in the size of migrant stocks at the beginning and end of a given period. The simplest form of differencing to estimate the number of migrant transitions y ij from origin i to destination j during a unit time interval subtracts the migrant stock totals in each migration corridor (a one-directional flow from i to j between a pair i, j of countries), setting all negative differences to zero: where s ij t and + s ij t 1 refer to the number of migrants born in i and residing in j at time t and t + 1 respectively and i = 1,2, …, R and j = 1,2, …, R for R places of birth and residence. This method has been used to create a dependent variable for econometric models to explain migration flows 12,13 . The left side of Table 2 illustrates the migration flows estimated from applying this method to the hypothetical migrant stock data in Table 1.
Stock difference, reverse negative. Beine and Parsons (2015) 14 The estimated migration flows following this method for the hypothetical migrant stock data in Table 1 are shown on the right side of Table 2.
Migration rates. Dennett    www.nature.com/scientificdata www.nature.com/scientificdata/ As the total M is unknown, Dennett (2016) suggested approximating it using the sum of absolute net migration flows for the time interval. The hypothetical data in Table 1 are constructed on the assumption that births minus deaths equal zero in each country over the time interval, hence in both tables the global sum of migrant stocks is the same. Further, the sum of absolute net migration flows for the time interval can be derived as the difference in the population totals (the column totals) in each place of residence. Hence,  Table 3. Table 1  Demographic accounting. Demographic accounting methods to estimate bilateral migration flows have two distinguishing characteristics. First, they use the approach of Abel (2013) to rearrange the tables of bilateral migrant stocks, such as in Table 1, into an array of birthplace-specific migration flow data with missing cell values but known marginal values 7 . Second, they control for changes from births and deaths in migrant stocks over the unit interval. We discuss each of these in turn.
Demographic accounting, minimisation. Migration flows can be detected from the bilateral stocks in Table 1 from differences in the number of people moving to different places of residence, i.e. moving across columns. Moves across rows are not possible as people's birthplace (row) remains fixed. Thus, when bilateral migrant stock data from two successive points in time have matching row totals they can be presented as birthplace-specific origin-destination migration flow tables with known margins as shown in the bold typeface in Table 4. In each birthplace-specific flow table, the marginal sums correspond to stock data; the row (origin) sums are from the distribution of migrant stocks at the beginning of the period, and the column (destination) sums are from the distribution of migrant stocks at the end of the period. The remaining cells of each table are unknown, but are imputed using a variety of methods discussed later. For example, in Table 4, in the second sub-table for place of birth B, the last column on the right, headed "Sum" is exactly the transpose of row B from Table 1, left side, showing data at time t; and, in Table 4, in the second sub-table for place of birth B, the last row at the bottom, labelled "Sum" exactly equals row B from Table 1, right side, showing data at time t + 1.
Unlike our dummy data in Table 1, in which the column headed "Sum" in the left table is identical to the column headed "Sum" on the right, reported bilateral migrant stock data do not have matching row totals at different times. Apart from measurement errors, the differences are predominantly due to demographic changes from births and deaths over the period. Commonly the numbers of births and deaths of the total population (natives and migrants) are known for each place of residence. These counts can be used to modify the migrant stock tables to reduce differences in row sums when estimating bilateral migration flows. Typically, each migrant stock at the beginning of the period is adjusted by proportionally reducing each stock by the total number of deaths in each place of residence. Death rates are assumed to be the same for foreign-born and native populations in the absence of global measures of mortality rates by nativity that could enable estimation methods to adjust for  A  0  5  5  10  A  5  5  30  40   B  0  0  0  0  B  5  10  45  60   C  25  0  10  35  C  25  5  10  40   D  0  0  15  15  D  0  35  15  50   Sum  25  0  20  15  60  Sum  30  45  30  85  190   Table 2. Estimated bilateral migration flows from stock differencing. On the left side are the estimated migration flows when negative differences are set to zero. On the right side are the estimated migration flows when negative stock differences are reversed. Diagonal elements are all zero.  Table 3. Estimated bilateral migration flows from the approach of Dennett (2016) based on rates and migrant stock data from Table 1 and an assumption of 70 total migration flows during the period.

Origin (t)
www.nature.com/scientificdata www.nature.com/scientificdata/ known heterogeneity. The migrant stocks at the end of the period are adjusted by subtracting the number of births to the native-born populations 7 in each place of residence from the diagonal elements.
In greater detail, native born populations -the diagonal cells in migrant stock tables -may increase over the period from return migration and new babies. New babies can come from native or foreign-born mothers (i.e. from a migrant stock that is off-diagonal in the stock table). Regardless of their mother's nativity, new babies increase the native population and hence are counted in the diagonal cell of the stock table. If we do not control for the number of births over the period (by subtracting the total number of births from the diagonal elements), then when we estimate flows from the stock differences (differences in both native and foreign-born) we obtain much higher levels of return migration than expected.
While these controls for demographic changes over the period ensure modified bilateral stock tables with smaller differences in row totals, they do not ensure a perfect match. A match is required to allow successive bilateral migrant stocks to have the same row totals, which permits the migration flows to be estimated. To close the small difference, two approaches have been proposed to further adjust the bilateral migrant stocks. Abel (2013) proposed estimating the number of migrants from each place of birth, in each place of residence, that leave to places of residence outside the migration system 7 . For example in the hypothetical data of Table 1, this would involve estimated counts of migrants to places of residence outside of the four (A to D) where data are available. The values of these estimated flows to places of residence outside of the migration system are derived so that they sum to the residuals of the difference in row sums between the migrant stock tables at times t and t + 1. This approach creates an open demographic accounting system where persons can move to or from countries beyond the set of those in the input bilateral migrant stock tables.
By contrast, Abel and Sander (2014) proposed scaling each migrant stock table so that the row sums matched the mid-point of their difference, whilst the column totals remained unchanged 8 . This approach creates a closed demographic accounting system where all persons either move, do not move, are born or die in the same set of countries.
Our estimates apply both the open and closed demographic accounting methods to produce two different sets of estimates. Both of these estimates use the following minimisation approach to estimate the missing bilateral migration flows, given the marginal totals from the different demographic accounting systems.
Given matching row sums in two successive bilateral migrant stock tables, the cell value m ijk represents the migration flows from origin i to destination j for birthplace k, where i, j, k = 1, 2, …, R, with known margins  www.nature.com/scientificdata www.nature.com/scientificdata/ The diagonal cells' values are set to their highest possible value, which is the minimum of the corresponding row and column marginal sums. These imputed diagonal values are shown in italic font in Table 4. The off-diagonal cells are estimated based on a quasi-independent log-linear model, where α ik , β jk and γ ijk are estimated using an iterative proportional fitting procedure 15 given the row and column sums and diagonal elements. The final estimates from the algorithm are shown in off-diagonal cells of Table 4. Given the full set of m ijk estimates, origin-destination migration flows can be obtained by summing over birthplaces: The estimates from the dummy data are shown on the left side of Table 5.
Demographic accounting, pseudo Bayesian. As the diagonal elements of m ijk represent the maximum number of stayers in each country (people with the same origin and destination), the estimated off-diagonal elements are the minimum migration flows required to match the changes in migrant stocks. Azose and Raftery (2018) 11 addressed this possibly undesirable feature by averaging the minimum flow estimate and an alternative array of migration flows, z ijk , based on an independent log-linear model, without terms for diagonal cells of the array: The α ik and β jk terms are estimated using an iterative proportional algorithm given the known marginal totals for each birthplace-specific flow table. Estimates from this model are shown in Table 6, where the known marginal constraints are shown in bold.
The estimated origin-destination migration flows are a weighted average of the estimates of minimum migration flows and those under the independence model: Azose and Raftery (2018) estimated w = 0.87 using a pseudo-Bayesian approach to minimise a cost function of the logged differences between estimated outward migration flow rates with equivalent estimates of European migration flows between 2002 and 2008 from the IMEM project 16 . Their application of the method used the closed demographic accounting system to obtain marginal totals for the birthplace-specific flow table of m ijk and z ijk . Estimates using the same value of w for the dummy data are shown on the right side of Table 5.

Estimates of Bilateral international Migration Between all pairs of 200 countries
To estimate bilateral migration flows from changes in bilateral migrant stocks for all countries between 1990 and 2015, we obtained and cleaned input data (Fig. 1). Demographic and bilateral migrant stock measures came from the United Nations Population Division (UNPD) as a sequence of bilateral migrant stock data covering all pairs of 232 countries for the mid-years of 1990, 1995, 2000, 2005 and 2015 17 . At the time of writing (February 2019), this was most up-to-date set of bilateral migrant stocks available for all countries. Data were primarily based on birthplaces reported in censuses and population registers provided by national statistical institutes. The UNPD adjusted data to include available refugee statistics. When recorded data did not align to the mid-years at five-year intervals, the UNPD extrapolated values based on the change in the overall population size. For countries or areas without data, the UNPD used a similar country or group of countries to estimate missing bilateral stocks.
The demographic accounting flow estimation methods require population totals at the beginning and end of each time period, and births and deaths over each interval. Further, the migration rates method requires net migration estimates during each period. Estimates of each of the measures were obtained from the 2017 version of the UNPD World Population Prospects (WPP2017) 18 .
Of the 232 countries in the migrant stock data, only 200 had complete estimates of other demographic measures such as births and deaths. A full list of the countries used is given in Supplementary File 1. We excluded all countries that lacked complete information. All of the 32 excluded countries had populations below 100,000.   www.nature.com/scientificdata www.nature.com/scientificdata/ G.J.A. developed the migest R package 19 to calculate each of the estimation methods outlined in the previous section (namely, a. stock differencing approach based on dropping zeros, b. stock differencing approach reversing negative estimate flows, c. migration rates approach, d. demographic accounting using an open system and minimisation approach to estimate missing flows, e. demographic accounting using a closed system and minimisation approach to estimate missing flows and f. demographic accounting using a closed system and pseudo-Bayesian approach to estimate missing flows). These routines were applied to clean sets of complete data for 200 countries five times to estimate bilateral migration flows for each five-year period between 1990 and 2015. The R code for cleaning the data and estimation is briefly described in the code availability section.
The total migration flows and their corresponding crude migration rates, obtained by dividing the total migration flows (summed over all pairs of countries in both directions) by the total population of all countries at the beginning of the period, are shown in Fig. 2. Summary statistics by estimation method are shown in Table 7. The bilateral migration patterns during the most recent period (2010-2015), aggregated to the United Nations regions, are shown in Fig. 3. To facilitate comparisons, the margins of each regional sector of the chord diagram plot are set to their maximum values over all time periods and estimation methods. Plots were produced using the circlize R package 20 . Animations of the changes in regional bilateral patterns for each estimation method, plotted with both fixed and non-fixed margins, can be found on Figshare 21 .

Data records
The international migration flow estimates described in this article, between all pairs of 200 countries listed in Supplementary File 1, are publicly and freely available through Figshare 21 . The estimates are stored in a single file with separate rows for each migration corridor -period combination (200 origins × 200 destinations × 5 periods = 200,000) and separate columns for estimates based on each method described above: sd_drop_neg, sd_rev_neg for the estimates using the stock differencing approaches; mig_rate for estimates using the migration rates approach; and da_min_open, da_min_closed and da_pb_closed for estimates using the demographic accounting approaches.

Technical Validation
We compared the estimated bilateral migration flows with two data sets: bilateral and total migration flows from selected countries reported by the United Nations 22 ; and net migration flows for all countries from WPP2017. The first data set is based on numbers of migrants sent and received by 45 countries using a variety of definitions of migration events and different concepts to determine their origin and destination. We used the reported migration flows between 1990 and 2015 based on previous or next place of residence for the total population, for which there are 32,698 one-directional migration corridors-period combinations that correspond with our estimated bilateral migration flows.
For each of our 200 countries, for each five-year period, for each of the six methods of estimation, we computed an estimated net migration as the difference between the sum of immigrants from all other countries, minus the sum of all emigrants to all other countries (among the 200 countries in our study). A positive net  www.nature.com/scientificdata www.nature.com/scientificdata/ migration therefore meant more immigration than emigration. We compared each estimated net migration with the corresponding net migration flows from the most recent WPP. The WPP2017 net migration data are based on a combination of net residuals from population accounting and official migration statistics, where they exist. WPP2017 estimates reflect demographic accounting residuals, net migration counts (where available), accumulated statistical errors from population, fertility and mortality estimates, and adjustments the UNPD carried out to ensure that global net migration is zero in each five-year period. The WPP2017 net migration estimates are available for all pairs of 200 countries and periods where we have migration flow estimates. calculation of equivalent reported flow. Reported migration flow data collected by the UN are for single years, whereas the estimated flows are based on the migration flows over a five-year period. Thus, for a given five-year period in a migration corridor there are potentially five reported migration flows to compare with a single estimated value. To provide an approximate guide to the level of reported migration flows in each five-year period we multiplied the average annual migration flow in each corridor and period by five.
We believe this modified reported flow should be highly correlated with a good estimated flow. However, we do not believe that the two numbers will match precisely as there are a range of measurement issues that prohibit an exact match. For example, national statistical institutes use different criteria to define how long a new entrant to the country must stay to qualify as a migrant. In some countries, a new entrant must declare an intention to settle permanently. In others, a stay of one year, six months, three months, or one month qualifies a new entrant as a migrant. These definitions may boost or dampen counts of migration flows along specific corridors.    www.nature.com/scientificdata www.nature.com/scientificdata/ We do not adjust the UNPD WPP2017 net migration data for the validation, as both the time period and the measurement of migration match those of the estimated flows from all methods. Explicitly, the UN net estimates and the estimates from our application of any of the methods have corresponding time periods (mid-year to mid-year over 5-year intervals) and both count migrant transitions (based on migrants' location at the beginning and end of the interval) rather than the number of moves during the interval.

Validation measures.
Due to the problems in comparing the reported annual migration flows with estimated five-year flows, we do not measure validity by accuracy compared to a 'true' observed value, as would be measured for example by the mean absolute error or mean squared error. Instead, we use Pearson's correlation between the estimated five-year flows and the modified reported flows to measure the validity of estimates. Hence, we avoid the need for an accurate multiplier for the average of the annual reported flows because the correlation between two positive variables is invariant when either variable is multiplied by a positive constant, or both are.
We assess the correlation between the estimated migration flows and available reported data by using six common measures of migration (Table 8). These six measures reflect possible applications of migration flow estimates to monitor trends (e.g. counts), explore patterns and migration drivers in explanatory models (e.g. logarithm of counts), as weights to detect migrant network structures (proportions), or baseline measures for inputs into population projection models (immigration and emigration rates or net migration counts).
The first measure is the count (number of individuals per five years) of a migration flow. Estimates that are statistically independent of the reported flows will have a sample correlation not far from zero, and estimates that are only weakly related to reported flows will have a low correlation.
The second measure is the natural logarithm of the count. Estimates that are randomly variable percentages of the reported flows will have low correlation.
The third measure is the proportion of migrants, where the denominator is either the total inflow or outflow over the period depending on which reported data source (receiving or sending data) is used. These three measures can be calculated for 32,698 corresponding equivalent reported bilateral flows. Of these, 5,134 observations are based on 2,567 origin-destination-period combinations where both sending and receiving data are available.
The fourth and fifth measures are the emigration rate from a source country and the immigration rate to a destination country, respectively. They divide the flows out of (emigration) and into (immigration) a country by the population of either the origin or the destination. The population sizes in the denominator are the corresponding population sizes in WPP2017 at the beginning of the five-year period.
The sixth measure is the net migration count. The WPP2017 reported net migration counts for our 200 countries in all five periods, enabling us to compare 1,000 estimates and reports. www.nature.com/scientificdata www.nature.com/scientificdata/ overall agreement. For each of these six measures, we calculated the Pearson correlation between the estimated flows and the equivalent reported flows (Fig. 4). For all bilateral measures (count, logarithm of count, proportions), the estimates based on Azose and Raftery (2018) method have the highest correlations (bottom row of Fig. 4). The correlation for the logarithm of count measure is at least as strong as or stronger than the non-transformed counts over all estimation methods (correlations in column 2 of Fig. 4 are greater than or equal to correlations in column 1).
The correlations of the estimated and reported emigration rates (column 4) are lower for all methods than the correlations of all other migration measures. The two stock differencing and migration rates methods have close to zero correlation between their emigration rate estimates and the reported emigration rates. Correlations of the immigration rates are higher for the closed demographic accounting method of Abel (2018) and Azose and Raftery (2018). These two methods produce estimates that are perfectly correlated with reported net migration counts in all countries because both methods use the scaling technique of Abel and Sander (2014). This technique implicitly constrains the estimated bilateral migration flows sum to net migration by adjusting migrant stock totals so that their differences over the total population follow the demographic accounting equation. Other estimation methods that do not use this technique, except for migration rates method (row 3), have higher correlations (column 6) between their net migration counts and the reported data than correlations using other migration measures.
For greater temporal resolution, we calculated the correlation for each estimation method and each validation measure in each time period (Fig. 5). For many of the methods, the migration count correlations appear to decline in more recent five-year periods, while the correlations in the logarithm of counts are relatively stable or slightly increasing. Proportion correlations for most methods fall to lows during 2000-2005 and then increase.   www.nature.com/scientificdata www.nature.com/scientificdata/ Correlations in the emigration rates are generally lower and vary more over time than the correlations for the immigration rates. The highest correlations between the reported and estimated immigration rates occur in later periods based on closed demographic accounting methods (bottom right corner of the panel for Immigration Rate in Fig. 5). The correlations between the reported net migration counts and the estimates peak during 2005-2010, with the exception of the estimates from the closed demographic accounting methods that perfectly correlated in every period (next-to-last column of the panel for Net Count in Fig. 5).
We also analysed the effects of the level of development (more or less developed) on the correlations of estimated flows with the reported data for each validation measure using the development level of the origin and destination in the United Nations World Population Policy Data 23 (Fig. 6). The development category for each country is listed in the un_group column of the country list in Supplementary File 1. Over all periods, there are  www.nature.com/scientificdata www.nature.com/scientificdata/ 471 equivalent reported migration flows between less developed countries, 6,667 between more developed countries, 10,580 from more to less developed countries and 12,413 from less to more developed countries. Only five countries (Armenia, Azerbaijan, Cyprus, Kazakhstan and Kyrgyzstan) reported flows between less developed countries. The correlations for net migration counts were based on all countries (155 less developed, 45 more developed) because the UNPD generated net migration figures regardless of the availability of supporting data.
The correlations of bilateral migration measures with equivalent reported data are highest for migration flows from less to more developed countries for all estimation methods, with the exception of the counts based on the drop negative stock differencing approach (where the correlations between reported and estimated flows to and from less developed countries are highest). The correlations of estimated emigration rates with equivalent reported flows are higher for migration flows from more developed countries for all estimation methods except Stock Difference Drop Negative. There was less uniformity for the immigration rates, where higher correlations were found for flows into more developed countries for only the estimates from the migration rates and pseudo-Bayesian approach. For the four estimation methods without the implicit net migration constraint, all had higher correlation of net migration counts with reported data for less developed countries.
These validation exercises provide new insights into global migration flow estimates. First, the correlation between the reported data and the estimates varies widely among different migration measures. For example, correlations of emigration rate estimates with reported data are lower than correlations of immigration rate estimates with reported data. Second, the correlations between the estimates and reported data vary widely in both time and space. Third, over most migration measures, estimates based on closed demographic accounting methods, especially those using the pseudo-Bayesian method, have the highest correlations.

Usage Notes
The estimates in this paper apply six estimation methods to the most recently available bilateral migrant stock data. The estimation methods based on migration rates and demographic accounting require additional information on demographic measures. Consequently, estimates that require these additional data are susceptible to errors propagated from inaccuracies in the stock and demographic data. As Rees (1980) suggested, demographic accounting methods such as those used here help users to identify consistencies, check for inadequacies, and match data with a conceptual model 24 .
The six methods of estimation used here share some limitations. All methods estimate that at least half of the bilateral migration flows are zero (Table 7). In addition, all estimates generate the number of migration transitions (i.e. the number of people in country i at the beginning of the period and country j at the end of the period); they do not take account of additional migration events to third countries or return moves during the period. Estimates based on stock differences and the demographic accounting methods by Abel (2013) and Abel (2018) produce minimum estimates of migration flows. Estimates based on stock differencing or migration rates ignore changes in stocks that might have occurred from deaths of migrants outside their country of birth. In cases with large elderly migrant stocks, deaths of migrants outside their country of birth could lead to overestimating migration flows.
Our validations of the bilateral migration and rate measures were carried out only where migrant flow data were reported, mainly from rich western countries; if and when migrant flow data become available from less developed countries, the validity of these measures may be different.

code Availability
All estimation routines to produce the open-access estimated five-year migration flows for all pairs of 200 countries between 1990 and 2015 are in the migest R package available on CRAN 19 . The formatted data and R code for the application of the functions are publicly available on Figshare 21 . They are 1) a tidy version of the United Nations bilateral migrant stock data, demographic period measures for net migration, births and deaths and demographic population totals; 2) an R script that given the three imported data sets a) cleans each to find the countries with complete input data; b) calculates the native population for the diagonal elements of the migrant stock data required for some estimation methods; and c) estimates bilateral migration flows for each method and period combination.