our understanding of the mechanisms and patterns of international migration over time are impeded both by the lack of data and by inconsistencies in the measurement and collection of the data that are available .
in fact , it is well known that the patterns of migration vary significantly depending on which country is reporting the data ( kupiszewska and nowok 2008 ; nowok et al .
considering that international migration is the main factor contributing to population growth in europe , this is very unfortunate . in response to the problem of inconsistent migration data
, we have developed a methodology for harmonising the data available to us from countries in europe .
more specifically , we make use of doubly counted information obtained from migrant sending and migrant receiving countries to estimate adjustment factors necessary for producing a consistent set of migration flows .
harmonisation of migration data is required for the development of policies on immigration ( kraier et al .
differences in both the concepts and techniques used to measure migration make any international comparison of migration difficult .
there has been a lot of work on data issues and migration definitions , for example , see champion ( 1994 ) , kelly ( 1987 ) , kraly and gnanasekaran ( 1987 ) , poulain ( 1993 ) , poulain et al .
( 2006 ) , raymer and willekens ( 2008 ) , united nations ( 2002 ) and willekens ( 1994 , 1999 ) . several international institutes such as the international labour organization , the organisation for economic co - operation and development , the united nations and the european commission have all invested heavily in the harmonisation of international migration data , but without much success or progression ( bilsborrow et al . 1997 ; herm 2006a ; fassmann 2009 ) .
in fact , the situation today in terms of migration definitions and measurement is not much better than it was , say , 20  years ago .
recently , some renewed efforts have been made to improve the migration data situation in europe . in 2007 , the european parliament adopted a new regulation on migration statistics .
this regulation provides clear definitions of immigration and emigration ( official journal of the european union 2007 ) , and lists the migration data that must be supplied to eurostat , the statistical office of the european union ( eu ) , by member states .
however , this regulation leaves the member states free to decide how they will provide these data , including the use of estimation methods ( fassmann 2009 ) .
the methodology presented in this paper should help national statistical offices to improve and harmonise the data they currently provide to international organisations , such as eurostat .
the migration definition set out in the 2007 regulation corresponds to the definition recommended by the united nations ( 1998 ) , where an international migrant is defined as  a person who moves to a country
other than that of his or her usual residence for a period of at least a year. one problem affecting the implementation of this definition is that some countries are unable to identify their nationals who have left ( fassmann 2009 ) .
furthermore , many european countries exclude the immigration of nationals from the published statistics , as they are not considered to be  migrants. another important obstacle has to do with the recommended duration of residence in the country of destination .
it may take up to 2  years to identify all persons who have stayed at least 1  year , as they may arrive anytime during the annual time period of interest .
this means that the publication of migration statistics based on the actual duration of stay may be delayed for some time . to provide statistics to the user community in a quicker fashion
, many countries simply count those migrants who have stayed for at least 3  months , which leads to higher numbers than if the 1-year criterion was applied .
other countries use the intended duration of stay as the criterion ( fassmann 2009 ) .
this is mainly caused by the fact that migrants have little incentive to report their move to the administration of the country they have emigrated from .
moreover , it is difficult to count persons leaving the country because they are no longer present in the country collecting the data . in this situation ,
comparisons of sending country data with receiving country data provide important information on the degree of underestimation found in reported emigration flows ( unece 2009 ) .
in fact , the analysis of the so - called  double - entry matrix of migration flows produced by unece since the early 1970s , and more recently by eurostat , has been found to be very useful and informative .
kelly ( 1987 ) and poulain ( 1999 ) , for example , have used the information contained in this matrix to assess the degree of harmonisation amongst reported data . in doing so
, the possibility that very narrow or loose definitions of migration may be used for reported immigration statistics must be taken into account , which results in lower or higher levels of migration flows , respectively , in relation to , say , the united nations recommended 1-year definition ( unece 2009 ) .
the aim of this paper is to illustrate how reliable estimates of harmonised migration statistics may be obtained from a set of origin  destination flows , where two reported flows are available for each particular flow , i.e. , from receiving and sending countries . the new method that we present
is based on earlier efforts by poulain ( 1993 , 1999 ) , and is applied to reported flows between 19 european countries from 2002 to 2007 .
note , however , that this paper does not consider flows outside the 19 country system , or those that are missing .
the reliability of migration statistics can be measured by how well they correspond to a particular country s definition or concept of migration . however , as definitions differ across countries , reliability does not guarantee comparability .
moreover , under - registration , under - coverage and accuracy of the collection system also affect the measurement of migration ( bilsborrow et al .
. this may be the case if the data depend on declarations by the migrants themselves .
the willingness to report changes in places of residence varies both between countries and between groups of migrants .
in general , migrants have more incentive to report their arrival than their departure , as there are usually direct benefits in doing so ( e.g. , access to social services ) .
therefore , immigration statistics are generally considered more reliable than emigration statistics ( thierry et al . 2005 ; unece 2009 ) .
this measurement category refers to the non - inclusion of particular migrant groups . here ,
the differences are most often caused by the absence or inclusion of nationals , students , asylum seekers or irregular ( illegal ) migrants in the data .
in general , asylum seekers are included only when they have been granted refugee status and received a temporary or permanent residence permit .
however , in some instances , they are registered at an earlier stage of the asylum process . in other instances ,
irregular migrants are generally not included in migration statistics , as they are especially difficult to measure ( for obvious reasons ) .
in fact , spain is the only eu country that includes irregular migrants in the official statistics .
finally , data based on sample surveys may be unreliable due to sampling errors . furthermore , unless the sample size is very large , the data are likely to show irregularities in the patterns across ages or in the distribution of origins or destinations over time , as flows of migrants represent a relatively small proportion of the overall population being surveyed .
the main sources of the differences in the definitions used by eu countries to measure migration are the concepts of place of residence and duration of stay ( zlotnik 1987 ; bilsborrow et al .
the de jure ( legal ) approach to residence implies that in order to become a resident , a migrant must comply with certain regulations , which tend to differ between nationals and foreigners , and among foreigners , between eu- and non - eu - nationals .
for example , it is not uncommon for emigrants to be registered in their country of citizenship ( origin ) even after several years of living abroad ( thierry et al .
thus , having a place of residence does not necessary imply a presence in that country .
the de facto ( actual ) approach is connected with physical presence in a country , usually for a specified minimum period of time . to prevent the delay caused by measuring actual duration of stay , most european countries use the intended duration of stay instead ( nowok et al . 2006 ) .
alternatively , the intended duration of stay may be used to provide provisional statistics , which are updated at a later point with the actual duration of stay statistics .
another group of countries measure  permanent change of residence only ( e.g. poland and slovakia ) , which is very restrictive and tends to produce flow levels that are much lower relative to other definitions . the duration of stay criterion used by the majority of eu countries is between 3  months and 1  year .
only three countries ( cyprus , sweden and uk ) apply strictly the 1-year criterion for immigration , as well as for emigration and for both nationals and non - nationals ( thierry et al .
in fact , some countries do not take duration of stay into account at all .
germany is such an example , where everybody taking up a residence is counted as a migrant . because of differences in definition , coverage , registration and accuracy of the collection mechanism , the origin 
destination matrix of migration flows between european countries based on immigration data reported by the countries of origin tends to differ from the matrix reported by the countries of destination . with respect to definitions ,
for example , the german definition is wider than the dutch definition which , in turn , is wider than that of sweden .
in fact , germany reports higher figures than the netherlands , and the figures of the netherlands are higher than those reported by sweden ( kupiszewska and nowok 2008 ) .
a comparison of the size of these reported flows provides information on the effects of differences in definition on the size of migration flows ( bilsborrow et al .
. however , as mentioned above , not all differences can be explained by differences in definition . in some cases ,
furthermore , sudden jumps in observations may be caused by changes in definitions or by changes in the registration method .
data on immigration and emigration flows by country of origin and destination are usually presented in an origin 
destination matrix with off diagonal entries containing the number of people moving from any origin i to any destination j in a given calendar year . for this study , we have collected migration data for the 19 countries set out in table  1 . as each flow can be reported by both sending and receiving countries , two migration tables may be produced .
the average 20022007 values of migration between the 19 european countries set out in table  1 are presented .
table 2 contains flows reported by the countries of destination and table 3 contains the flows reported by the countries of origin
. clearly , there are large differences between the two sets of reported numbers ( see , e.g. , spain to the united kingdom or poland to germany).table  1list of european countries reporting both immigration flows by country of origin and emigration flows by country of destination , 20022007countryabbreviationaustriaatcypruscyczech republicczgermanydedenmarkdkspainesfinlandfiicelandisitalyitlithuanialtluxembourglulatvialvnetherlandsnlnorwaynopolandplswedensesloveniasislovakiaskunited kingdomuktable  2reported migration by country of destination , averages 20022007fromtoatcyczdedkesfiisitltat4131014,2573037741093377417cy22132762325231303cz1,3161189,218262833564267224de15,4473321,3624,00115,98292125512,809490dk20325462,6879643651,41326585es700457114,7031,758644682,044252fi27021382,1734148444523543is31042361,665131503510it1,6084925422,1969869,3202507482lt17935474,4961,0342,27473272378lu67322,28216212350272135lv83104132,1554573008793183175nl7917025513,6818644,7622615590541no9814241,3783,1481,69684536416787pl5,2317521,608136,9272,4368,2771872,2299,045120se48988673,3483,3131,8263,50249237991si5569171,79846136693212sk3,19243214,06411,14814978822456904uk1,2223,17050613,2633,48238,6749462284,553875total31,5045,30618,702256,22124,50287,7258,3975,74133,6952,407fromtolulvnlnoplsesiskuktotalat895591111803071002081,39519,496cy025115761222,5333,087cz4155111164516469794,10918,489de4541669,1822,2682,8763,37429944619,03989,701dk11464752,943345,2643211,87416,721es24183,1017681191,30083614,58140,239fi34337979963,204166849,208is067537311462124173,509it67331,811246309599791095,82943,900lt123630292643574052,50713,380lu216118590516823,897lv21252336264051,2275,511nl272071116397912416,79930,436no224453485,0981241,66715,135pl19455,7444,6023,718327636,759217,977se14546964,91711315203,21322,635si1190142421603,064sk444652381811044,58435,961uk391905,8201,6241,1263,1142211678,969total68291330,00020,9215,11128,7235592,311107,897671,315source : eurostattable  3reported migration by country of origin , averages 20022007fromtoatcyczdedkesfiisitltat189376,665166429231271,022111cy62157619120399cz18613560243528211210de17,7872718,1043,09516,8072,37128731,2352,455dk228241792,6121,6693681,347716655es1559572 68615711091 163120fi9723427584006715320321is132172051,800594810564it58866710,2061491,5081361711lt488541,2691586288723204lu31313911997935191754lv1886302451846551138nl6165029810,4935333,774322541,27854no6917437093,093789855412146108pl538156314,41711134120465056se298731041,6343,1591,3483,40341346348si311314589527411861sk177162925541610420uk1,5934,0602,69212,5791,93233,4316821035,2701,074total22,7584,60013,33966,90514,93361,6498,7582,81842,9144,887fromtolulvnlnoplsesiskuktotalat4542426872,4013884021,77890116,076cy21810211113032371724cz3781165832489,53921911,449de1,6861,4949,2932,122100,8273,9742,0049,45617,233230,499dk1383166022,9478335,25331953,88921,898es871986915939820310453 4309 684fi7127233777633,2164101,1757,842is37294948287247825562324,570it2188531121417199151403,50817,879lt18163116199122233352,6385,975lu4971223735111661,760lv62034266712196987nl191337311,020900451387,95328,482no23692872815,0835611,39513,444pl2335571273032105,21922,306se127625224,74635427293,90520,713si2403055386701,319sk201331583691,235uk3623245,9431,9936,5072,66601,05382,264total3,0622,61919,67614,561114,85423,1172,72422,36452,567499,105source : eurostat list of european countries reporting both immigration flows by country of origin and emigration flows by country of destination , 20022007 reported migration by country of destination , averages 20022007 reported migration by country of origin , averages 20022007
the differences between reported immigration and emigration numbers are useful for improving and harmonising the migration data .
if reported emigration numbers for a given country turn out to be systematically lower than the corresponding immigration numbers reported by the countries of destination , this suggests that the reported emigration numbers are too low . adjusting these numbers in an upward direction
we can estimate one adjustment factor for immigration and one for emigration in such a way that the adjusted immigration and emigration numbers are closer to each other than the reported numbers . to prevent arbitrary judgments biasing the results
, we believe the estimation of adjustment factors for immigration and emigration flows should be estimated simultaneously .
moreover , it should be noted that immigration is not necessarily recorded more accurately than emigration . in some situations , sending country data may be considered better ( nowok et al .
poulain ( 1993 , 1999 ) was the first to develop a method to adjust reported immigration and emigration numbers for the purpose of obtaining a consistent set of migration flows . 
correction factors were estimated by minimising the sum of squares \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \sum\limits_{i , j } { ( \hat{\alpha } _ { j } i_{ij } - } \hat{\beta } _ { i } e_{ij } ) ^{2 } $ $ \end{document } , where iij denotes migration from country i to country j reported by the receiving country j , eij denotes the same flow reported by the sending country i , j is the adjustment factor for all immigration to country j and i is the adjustment factor for all emigration from country i. poulain and dal ( 2008 ) refined this method by dividing the squared differences by the sum of the reported numbers , i.e.,1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \sum\limits_{i , j } { ( \hat{\alpha } _ { j } i_{ij } - } \hat{\beta } _ { i } e_{ij } ) ^{2 } /(i_{ij } + e_{ij } ) .
$ $ \end{document } this refinement prevents flows from ( or to ) large countries from biasing the estimates .
for instance , following the iterative approach to harmonising migration flows suggested by van der erf and van der gaag ( 2007 ) , poulain and dal ( 2008 ) proposed that the estimates should be normalised to swedish immigration data , as they are generally considered to be highly reliable and in agreement with the un recommended measure , as well as with the new eu regulation ( herm 2006b ) .
the parameters j and i may be estimated by solving a system of linear equations , which result from applying the method of lagrange multipliers .
multiplying iij by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\alpha } _ { j } $ $ \end{document } and eij by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\beta } _ { i } $ $ \end{document } produces two sets of migration flow estimates from country i to country j. the final set of estimates are obtained by simply taking the average of the two , i.e. , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{n}_{ij } = ( \hat{\alpha } _ { j } i_{ij } + \hat{\beta } _ { i } e_{ij } ) /2 $ $ \end{document } , where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{n}_{ij } $ $ \end{document } denotes the harmonised migration flows . note , poulain and dal ( 2008 ) applied their correction method first to countries with relatively reliable data to prevent countries with less reliable data influencing the overall patterns . here , the main concern is that the less reliable data have origin  destination patterns that are not consistent with the actual patterns .
thus , less reliable flows were adjusted in a hierarchical fashion , i.e. , by using the harmonised reliable data as a basis .
first , the reported numbers included in the denominator of eq .  1 are known to be incorrect ( abel 2009 ) .
second , the row and column totals of the two estimated matrices are not equal . as a result
, the row and column totals of the average harmonised migration matrix do not correspond to the row and column totals estimated using the adjustment factors .
finally , the method can only be applied to a limited set of countries with reasonably reliable data .
this implies that the estimates of the adjustment factors depend on the selection of countries , which may not reflect the broader patterns of interest . for these reasons ,
first , the row - sums and column - sums of the two estimated matrices are set to be equal .
second , we introduce additional constraints on individual cells in the migration matrices , so that more countries ( with less reliable data ) may be included .
the adjustment factors for our method can be estimated by solving a system of linear equations and imposing a constraint .
if we have a n    n receiving country and an equivalent n    n sending country matrix , the adjustment factors for receiving country , j , and the adjustment factors for sending country data , i , can be estimated by2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \sum\limits_{j } { \hat{\alpha } _ { j } i_{ij } = } \hat{\beta } _ { i } \sum\limits_{j } { e_{ij } } \quad   { \text{for } } i = { 1 } , \ldots , \ , n;\quad i \ ,
\ne \ , j $ $ \end{document}3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\alpha } _ { j } \sum\limits_{i } { i_{ij } = } \sum\limits_{i } { \hat{\beta } _ { i } e_{ij } }   \quad{\text{for } } j = { 1 } , \ldots , \ , n;\quad i \ , \ne \ , j
$ $ \end{document } equation  2 states that for each country the emigration total estimated on the basis of the adjusted matrix of flows reported by receiving countries equals the emigration total estimated on the basis of the adjusted matrix of flows reported by sending countries .
equations  2 and 3 can be written as a homogeneous system of 2n linear equations with 2n unknowns , i.e.,4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \begin{array}{l } \hat{\alpha } _ { 2 } i_{12 } + \hat{\alpha } _ { 3 } i_{13 } + \cdots + \hat{\alpha } _ { n } i_{1n } - \hat{\beta } _ { 1 } \sum\limits_{j } { e_{1j } } = 0 \\ \vdots \\ \hat{\alpha } _ { 1 } i_{n1 } + \hat{\alpha } _ { 2 } i_{n2 } + \cdots + \hat{\alpha } _ { n - 1 } i_{nn - 1 } - \hat{\beta } _ { n } \sum\limits_{j } { e_{nj } } = 0 \\ \hat{\alpha } _ { 1 } \sum\limits_{i } { i_{i1 } } - \hat{\beta } _ { 2 } e_{21 } - \hat{\beta } _ { 3 } e_{31 } - \cdots - \hat{\beta } _ { n } e_{n1 } = 0 \\ \vdots \\ \hat{\alpha } _ { n } \sum\limits_{i } { i_{in } } - \hat{\beta } _ { 1 } e_{n1 } - \hat{\beta } _ { 2 } e_{n2 } - \cdots - \hat{\beta } _ { n - 1 } e_{nn - 1 } = 0 \end{array } $ $ \end{document } this system has an infinite number of solutions for j and i .
for each set of values of j and i that solve this system , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ k\hat{\alpha } _ { j } $ $ \end{document } and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ k\hat{\beta } _ { i } $ $ \end{document } are solutions as well . in order to find a unique solution one restriction needs to be imposed . in accordance with poulain and dal ( 2008 ) , we assume that the adjustment factor for swedish immigration is equal to one , since sweden uses a definition of migration that is consistent with the new eu regulation and the quality of swedish immigration data is considered to be adequate .
this also means that the resulting estimates are harmonised in line with the new european regulation .
the basic assumption underlying our estimation procedure ( as described above ) is that the distributions of reported immigration by country of origin and reported emigration by country of destination correspond to the distribution of actual migration flows under the harmonised definition .
this implies that the reported emigration of country a is x% higher or lower than the actual number ( based on the standard definition ) for all countries of destination .
however , as we find in the next section , the estimated receiving country flows by country of origin and the estimated sending country flows by country of destination are not always consistent with each other . in a number of cases , specific origin 
let us assume that the estimated receiving country migration flow from country p to q , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\alpha } _ { q } i_{pq } $ $ \end{document } , differs substantially from the estimated sending country flow , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\beta } _ { p } e_{pq } $ $ \end{document}. to make them consistent , we can multiply \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\alpha } _ { q } i_{pq } $ $ \end{document } by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\gamma } _ { pq } $ $ \end{document } or \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\beta } _ { p } e_{pq } $ $ \end{document } by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\delta } _ { pq } $ $ \end{document } so that both estimates of migration are equal .
the question whether we should adjust the estimate based on the reported receiving country or the estimate based on the reported sending country depends on our knowledge of the data .
given the estimated values of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\alpha } _ { q } $ $ \end{document}and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\beta } _ { p } $ $ \end{document}we can calculate the value of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\gamma } _ { pq } $ $ \end{document } easily from \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\gamma } _ { pq } = \hat{\beta } _ { p } e_{pq } /\hat{\alpha } _ { q } i_{pq } $ $ \end{document } or the value of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\delta } _ { pq } $ $ \end{document } from \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\delta } _ { pq } = \hat{\alpha } _ { q } i_{pq } /\hat{\beta } _ { p } e_{pq } $ $ \end{document}. however , introducing \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\gamma } _ { pq } $ $ \end{document } or \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\delta } _ { pq } $ $ \end{document } changes the estimates of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\alpha } _ { q } $ $ \end{document } or \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\beta } _ { p } $ $ \end{document}. this also means that the row and column totals of both estimated migration matrices no longer tally .
therefore , we adjust the system of linear eqs .  2 and 3 by adding constraints on individual cells of the matrices .
if we assume that the emigration number reported by country p needs to be adjusted , eqs .  2 and 3 can be rewritten as5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \sum\limits_{j } { \hat{\alpha } _ { j } i_{ij } = \ , } \hat{\beta } _ { i } \sum\limits_{j } { e_{ij } ( 1 + \hat{\delta }
_ { pq}^ { * } d_{ij } ) } \quad   { \text{for } } i = { 1 } , \ldots , \ , n;\quad i \ , \ne \
, j $ $ \end{document}6\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\alpha } _ { j } \sum\limits_{i } { i_{ij } = } \sum\limits_{i } { \hat{\beta } _ { i } e_{ij } ( 1 + \hat{\delta }
_ { pq}^ { * } d_{ij } ) } _ { { } }   \quad{\text{for } } j = { 1 } , \ldots , \ , n;\quad i \ , \ne \ , j $ $ \end{document}where dij  =  1
if i  =  p and j  =  q , dij  =  0 otherwise , and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\delta } _ { pq}^ { * } = \hat{\delta } _ { pq } - 1 $ $ \end{document}. the equations including ipq and epq in the system of eq . 
4 can be rewritten as follows:7\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\alpha } _ { 1 } i_{p1 } + \cdots + \hat{\alpha } _ { q } i_{pq } + \cdots + \hat{\alpha } _ { n } i_{pn } - \hat{\beta } _ { p } e_{p1 } - \cdots - \hat{\delta } _ { pq } \hat{\beta } _ { p } e_{pq } - \cdots - \hat{\beta } _ { p } e_{pn } = 0 $ $ \end{document}\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\alpha } _ { q } \sum\limits_{i } { i_{iq } } - \hat{\beta } _ { 1 } e_{1q } - \cdots - \hat{\delta } _ { pq } \hat{\beta } _ { p } e_{pq } - \cdots - \hat{\beta } _ { n } e_{nq } = 0 $ $ \end{document } in contrast with eq .  4 , these are non - linear equations , because they include the term \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \hat{\delta } _ { pq } \hat{\beta } _ { p } e_{pq } $ $ \end{document}. the values of the coefficients can be estimated by an iterative procedure .
the model can be extended in a straightforward way to include additional constraints . however , for any particular country , the number of constraints should not be too high , as this reduces the available information to estimate  and .
the sending and receiving country migration data have been provided by the national statistical institutes of the eu member states in response to annual rounds of data collection conducted jointly by five international organisations and coordinated by eurostat ( kupiszewska and nowok 2008 ) .
as concerns europe , eurostat processes and disseminates data received from 37 countries on their website ( epp.eurostat.ec.europa.eu ) .
data sources used by eu member states to produce migration statistics are very diverse ( kupiszewska and nowok 2008 ; nowok et al .
the major types of sources are population registration systems , statistical forms , other administrative registers related to foreigners ( such as alien registers , residence permit registers and registers of asylum seekers ) , sample surveys and censuses .
alien registers and residence permit registers are used in seven countries , sometimes in addition to population registers .
cyprus and the uk rely on passenger surveys conducted at the borders , while portugal and ireland rely on household surveys .
some countries derive their emigration statistics from data on residence permits by assuming a migrant has left the country when a residence permit has expired .
moreover , they often assume that the country of next residence is the country of their citizenship .
the result , we believe , is an overestimation of actual emigration to those particular countries .
finally , several countries include in their so - called  administrative corrections emigration that has not been declared , which can not be disaggregated by country of next residence .
data on immigration by country of previous residence or emigration by country of next residence are not always available or complete ( nowok et al .
thus , the sending country and receiving country matrices , when combined into a double - entry matrix may be incomplete .
for some countries , a large share of emigrants have an unknown country of destination : around 75% in slovenia , 40% in luxembourg , 35% in austria , 31% in the netherlands and 39% in spain , for example .
fortunately , the estimation of adjustment factors takes this into account . in the next section
, we present our harmonised estimates of migration between 19 european countries that provide data on both immigration by country of origin and emigration by country of destination for the calendar years 20022007 .
although there are some data for ireland , portugal and romania , these have not been used because they cover only a part of the migration flows ( e.g. only foreigners or nationals ) . for iceland , italy and luxembourg ,
for these countries , the adjustment factors are estimated for averages over the available years .
the results presented in this section are obtained by applying the estimation method described in sect .
3 . table  2 shows the average values of migration between 19 european countries reported by receiving countries for the years 20022007 and table  3 shows the corresponding numbers reported by the sending countries .
the countries listed in the row headings refer to origins and those listed in the column headings refer to destinations .
a comparison of tables  2 and 3 reveals large differences between numbers reported by sending and receiving countries . according to the numbers reported by receiving countries , 671,315 migrants per year moved between these 19 countries , whereas the numbers reported by sending countries total 499,105 . for 11 countries , the reported receiving country immigration totals are higher than the corresponding sending country totals .
for example , germany reported that 256,221 immigrants arrived from the 18 countries in this study , whereas these countries reported that only 66,905 emigrants moved to germany .
poland reported that 22,306 persons emigrated to the other 18 countries which , for their part , reported receiving 217,977 immigrants from poland , suggesting that polish emigration data are around 10 times too low . for 15 of the 19 countries ,
the emigration total reported by the sending country is lower than the corresponding totals reported by receiving countries .
keep in mind that receiving country data should not always be considered better than sending country data .
consider , for example , the flows from poland to germany in tables  2 and 3 . here , germany received an average of 136,927 migrants from poland , whereas poland reported that they only sent an average of 14,417 .
this difference could be explained by the duration criteria used by these countries , with germany having a very loose definition ( instant ) and poland having a very restrictive definition ( permanent ) .
so , in comparison with the harmonised definition of a 1  year period , germany s reported number is considered too high and poland s too low .
we indicated above that in order to estimate the adjustment factors a restriction was introduced , i.e. , the adjustment factor for swedish immigration is set equal to one . for 16 of the 19 countries ,
the eij adjustment factor exceeds one , indicating that sending country numbers tend to be underestimated .
however , table  4 also shows that iij numbers seem to be underestimated in the majority of countries as well .
this may seem contradictory since for 11 of the 19 countries the reported immigration totals exceed the corresponding emigration numbers reported by the sending countries .
this is because the reported receiving country numbers should be compared with the adjusted sending country numbers rather than the reported numbers .
for example , the immigration total reported by the uk ( 107,897 ) exceeds the reported emigration from sending countries to the uk ( 52,567 ) .
however , since the reported emigration from poland is too low ( the adjustment factor equals 18.31 , see table  4 ) the reported emigration from poland to the uk is adjusted from 5,219 to 55,506 . moreover , the adjustment factor for spanish emigration data equals 4.32 , so the reported emigration from spain to uk is adjusted from 3,430 to 16,792 . for several other countries , emigration to the uk is adjusted upwards as well . as a consequence ,
the adjusted emigration numbers to the uk exceed the total of immigration reported by the uk and thus the reported immigration is adjusted upwards as well .
note that the adjustment factors for immigration for most countries are closer to one than the adjustment factors for emigration , which indicates that the reported immigration numbers are more accurate than the emigration numbers.table  4estimates of adjustment factors for immigration and emigration , 20022007immigrationemigrationaustria1.061.74cyprus1.065.29czech republic2.143.33germany1.030.69denmark0.740.80spain0.824.90finland1.261.22iceland0.570.74italy1.422.92lithuania2.332.45luxembourg5.652.43latvia2.926.22netherlands0.971.25norway0.841.19poland17.8510.64sweden1.001.21slovenia5.182.71slovakia18.9043.69united kingdom1.211.18 estimates of adjustment factors for immigration and emigration , 20022007 multiplying the reported numbers in table  2 by the adjustment factors for receiving country data and the reported numbers in table  3 by the adjustment factors for sending country data results in two tables for which the row and column totals are equal ( not presented here for space reasons ) .
the differences between the cells in these two matrices are considerably smaller than those in tables  2 and 3 . in fact , the root mean squared error ( rmse ) is reduced from 8,966 to 2,131 . in other words ,
the differences between the two reported migration flow tables are reduced by 77% . however , we still found some substantial differences in the two estimated migration flow tables .
for example , the migration from poland to germany estimated on the basis of german immigration data equals 141,035 , whereas the estimate based on polish emigration data is equal to 153,399 .
these differences reflect the fact that the distribution of reported polish emigration by country of destination is not consistent with the share of immigration from poland in the total reported immigration numbers of other countries . as a result , the estimate of the migration flow from poland to germany based on polish data exceeds that based on german data , whereas for most other countries , the adjusted polish emigration numbers are lower than the corresponding adjusted immigration numbers .
this means that one substantial inconsistency in the estimates is likely to influence the estimates of other migration flows . to prevent such inconsistencies from affecting the overall estimates , we have added constraints to individual cells ( flows ) in the model .
the introduction of constraints to individual cells in the matrix allows us to consider special cases , such as the poland to germany flow described above . in total
specifically , these flows were poland to germany , poland to uk , germany to poland , germany to uk , czech republic to slovakia and uk to poland . after identifying the flows with large differences , we then had to decide whether the constraint should be applied to the numbers of the receiving country or of the sending country .
since we believe that reported emigration numbers are generally considered to be less reliable than reported immigration numbers , we apply the constraints to the sending country data , except for the germany to poland and uk to poland flows ( i.e. , poland s immigration data is considered to be of lower quality that the corresponding emigration data reported by both germany and the uk ) .
the adjustment factors taking into account the six constraints on individual flows are set out in table  5 .
the coefficients ( lagrange multipliers ) for the poland to germany and poland to uk flows are both equal to 0.42 .
this raises the adjustment factor for emigration from poland from 10.64 ( table  4 ) to 18.31 ( table  5 ) , while at the same time , the adjustment factor for polish emigration to germany and the uk falls to 7.69 ( i.e. , 18.31    0.42 ) . for polish immigration ,
the high adjustment factor for polish receiving data was mainly a consequence of the big difference between the two figures for migration from germany to poland . including a constraint for this flow raises the adjustment factor for poland s reported flow from germany by a factor of 1.74 ( i.e.
in contrast , the adjustment factor for poland s reported flow from the uk falls to 10.40 ( i.e. , 14.25    0.37 ) . for the czech republic ,
the reported emigration numbers are considerably lower than the corresponding reported immigration numbers with one big exception : the number of emigrants reported to slovakia is relatively large .
clearly , the emigration flows from the czech republic to all other countries need to be adjusted by a different factor than the emigration flow to slovakia.table  5estimates of adjustment factors for immigration and emigration , 20022007 , including six additional constraints on individual flowsimmigrationemigrationaustria1.171.35cyprus0.884.71czech republic1.978.92germany0.810.71denmark0.720.74spain0.734.32finland1.181.12iceland0.590.69italy1.482.44lithuania2.162.15luxembourg5.452.08latvia2.785.44netherlands1.041.06norway0.811.10poland14.2518.31sweden1.001.10slovenia4.902.33slovakia8.3439.40united kingdom1.090.91coefficients for additional constraints ( lagrange multipliers )  immigration to poland from germany1.74  immigration to poland from the uk0.37  emigration from poland to germany0.42  emigration from poland to the uk0.42  emigration from germany to the uk1.70  emigration from the czech republic to slovakia0.10 estimates of adjustment factors for immigration and emigration , 20022007 , including six additional constraints on individual flows the adjustment factors in table  5 illustrate how substantial improvements in the estimated adjustment factors can be made by introducing constraints on specific 
for example , the inclusion of a constraint for the migration flow from the czech republic to slovakia lowered the adjustment factor for slovakia s receiving migration data from 18.90 to 8.34 .
this is mainly explained by the reduction of the estimate of polish emigration to germany .
since germany has a wide definition of migration , one would expect the adjustment factor to be below one .
thus , the adjustment factors in table  5 appear more plausible than those set out in table  4 . the harmonised migration tables that used the additional constraints
the introduction of these constraints led to a further strong reduction in the differences between both tables , as indicated by the rmse , which fell from 2131 to 952 or by a further 54% . to obtain a final single set of harmonised flows , we believe it is better to rely on table  6 than on table  7 .
this table gives more weight to the receiving country data , which we consider more reliable .
poulain , on the other hand , advocated taking the average of the two estimated matrices .
destination patterns in the reported sending country data are as reliable as those in the reported receiving country data.table  6estimated migration by country of origin and destination , including constraints on six individual flows , 2002/2007 , based on numbers reported by receiving countriesfromtoatcyczdedkesfiisitltat3661011,526217563129201,14437cy26252231618281447cz1,5471037,453187606662599251de18,1482912,6792,86411,6311,09215218,9201,057dk23821902,172701432840391183es8223913911,8871,259763403,020544fi31719751,7572966142734792is37081911,19295595122it1,8894350017,9457066,78229744177lt21131923,6357411,65586162558lu79251,84511689591631511lv9791261,74232721810355270378nl9306150211,0606193,465309331,33689no11512471,1142,2541,2341,001216246187pl6,1466593,163110,7011,7446,0242211,32413,361259se574771322,7062,3721,3284,150292560197si6538321,4543399754734sk3,75037927,6589,01210757426261,0209uk1,4362,78099610,7222,49328,1451,1211356,7251,886total37,0154,65436,780207,14517,54263,8419,9503,41149,7725,190fromtolulvnlnoplsesiskuktotalat4426584902,5703074881,7301,52221,642cy145312100618182,7643,408cz234353494644164288,1584,48225,200de2,4754629,5861,84371,5143,3741,4623,71520,770172,034dk601284962,3914805,264151762,04416,124es132493,2376241,6981,3004129715,90741,798fi18119395649883,2044517468,819is217783031624626184553,157it365931 8912004,4025993869066,35943,582lt76563167526105741392,73412,858lu416815699024107443,661lv9131189902640381,3385,367nl146555772,323979593437,41730,305no12684736795,09822011,81814,778pl1011255,9973,7393,718122,29840,102199,695se741497273,9961,608741683,50522,690si649412244212903,079sk2412486193254110205,00148,662uk2135286,0771,3205,9073,11410596874,670total3,7152,54031,32117,00093,22228,7232,73519,264117,708751,530table  7estimated migration by country of origin and destination , including constraints on six individual flows , 2002/2007 , based on numbers reported by sending countriesfromtoatcyczdedkesfiisitltat241,2618,973223578311361,376150cy2699266268756018142cz1,6581164,9992133092532199889de12,6161925,7482,19511,9211,68120322,1541,741dk168171321,9241,229271991527482es6693724711,592678475375,020516fi10925478524497546022824is91111421,24341337344it1,4331516424,8773643,6753314126lt102161152,7303391,35218750438lu655261,89520516573403648lv9844301,64024510025027277748nl6565331711,1655674,016342571,36057no7618477803,399867940453160118pl9,8572721,160110,7012,0236,2443698459,247116se327801141,7893,4601,4773,72845350752si7257331,37412631024352sk6,9744624,78410,0281446172001,6420uk1,4463,6852,44411,4181,75330,345619934,784975total37,0154,65436,780207,14517,54263,8419,9503,41149,7725,190fromtolulvnlnoplsesiskuktotalat61575731183,2325235412,3931,21321,642cy108645115205901481,7453,408cz25587211405,203213708,1581,95625,200de1,1951,0606,5911,50571,5142,8181,4216,70720,770172,034dk1012334432,1706133,86822702,86316,124es377833,7496861,7188754219514,80241,798fi8030262873703,6175111,3218,819is25203433360233017391603,157it531201,2942951,015484368968,55143,582lt393512504282635027115,67612,858lu9201264715210223453,661lv331091821413624111,0685,367nl203357771,085958481468,46230,305no26753163095,5875671,53414,778pl4126110,2052,3195,5393418940,102199,695se139685715,19938829314,27722,690si57170111188131643,079sk7204931185843281122,69948,662uk3282945,3951,8095,9072,420095674,670total3,7152,54031,32117,00093,22228,7232,73519,264117,708751,530 estimated migration by country of origin and destination , including constraints on six individual flows , 2002/2007 , based on numbers reported by receiving countries estimated migration by country of origin and destination , including constraints on six individual flows , 2002/2007 , based on numbers reported by sending countries the average adjustment factors estimated for the period 20022007 ( table  5 ) can be applied to the annual reported migration data to create a time series of harmonised flows . in fig .  1 , the estimated total immigration and emigration flows for germany from and to the other 18 countries in this study are compared . as expected ,
the estimated numbers are lower than the reported numbers because the definition for germany is much wider than the harmonised definition .
the figure also shows that estimated emigration increases more gradually over time than the reported numbers . in fig .  2 , the immigration and emigration flows for the uk
, the average levels of the reported and estimated numbers do not differ much , but the estimated flows show a more gradual pattern over time than the reported flows .
one reason for the sharp fluctuations in the reported numbers is that they are based on sample surveys.fig . 
1reported and estimated immigration from and emigration to 18 european countries , germany , 20022007fig . 
2reported and estimated immigration from and emigration to 18 european countries , united kingdom , 20022007 reported and estimated immigration from and emigration to 18 european countries , germany , 20022007 reported and estimated immigration from and emigration to 18 european countries , united kingdom , 20022007
the aim of this paper has been to obtain a reasonable and consistent set of international migration statistics . for this purpose
the method is based on an idea originally proposed by poulain ( 1993 , 1995 ) .
first , we have estimated a set of adjustment factors for receiving and sending country data in a way that ensures consistency in the two sets of marginal totals .
third , instead of calculating the arithmetic average of the two estimated matrices , we believe it is better to use the matrix giving more weight to the reported immigration numbers ( i.e. table  6 ) . in this way we take advantage of the fact that the information on countries of origin in receiving country data tend to be more reliable than the country of destination information in sending country data .
finally , our estimates are consistent with the harmonised migration definition based on an ( intended ) minimum duration of stay of 12  months . due to differences in definition , coverage and registration , the origin 
destination matrix of migration flows between european countries based on receiving country data tends to differ from the matrix based on sending country data .
germany has a wide definition of migration , as it does not include a time constraint and thus the reported number may well include short term migrants .
in contrast , poland has a very narrow definition of migration and , as a consequence , the reported numbers are very low . by comparing
corresponding reported immigration and emigration flows for 19 european countries , we have assessed to what extent german migration statistics are higher than they would be under a harmonised definition and to what extent polish migration statistics are lower . however , the large differences between european countries can not be explained by differences in definitions alone .
first , these differences can not explain why emigration flows are more likely to be underestimated than immigration flows .
second , whereas 11 countries employ a duration limit that is shorter than that of the harmonised definition ( kupiszewska and wisniowski 2009 ) , only five of these countries have an adjustment factor of immigration below one .
the other six countries with durations of 6  months or shorter have adjustment factors for immigration greater than one .
these include austria , czech republic , italy , luxembourg , the netherlands and slovenia .
thus , to an important extent , the differences must also be caused by problems of coverage .
this is confirmed by a study comparing migration statistics between sweden , denmark and belgium which suggests that less than 25% of differences are due to differences in the duration criterion ( nowok et al .
the effects of differences in definition and coverage may offset each other to some extent .
one would expect the under - registration of short term migrants to exceed that of long - term migrants .
a wide definition of migration ( i.e. a short duration of stay ) would lead to a higher reported number of migrants than would be expected on the basis of the harmonised definition . under - registration
this may explain why the adjustment factors for germany are not as low as one might expect from applying the very wide definition .
the main reason for the relatively low numbers reported by sending countries is that emigrants do not have strong incentives to report leaving a country . in particular , this applies to eu citizens who can live in another eu country without asking for a residence permit .
, any person leaving country a would be required to fill in a form to be given to the authorities in country b at arrival .
after country b has determined whether or not the person is an international migrant under a harmonised definition , it would then inform country a of the arrival .
the nordic countries have such a system and their immigration and emigration statistics are mutually consistent ( herm 2006a ) .
however , policy makers tend to be more interested in migrants from outside europe and asylum seekers than intra - european migrants , and therefore such a system is not likely to have a high priority in the future . as long as such a system is lacking , cross - country comparability of migration statistics can only be achieved by comparing statistics from different countries . to the extent that the differences between countries are caused by differences in definitions and coverage ,
the method developed in this paper aims to assess the size of these systematic differences .
table  5 shows that for 10 out of the 19 countries in this study , the adjustment factor for sending country data exceeds two , meaning that reported emigration numbers are underestimated by more than 50% in relation to the 1-year duration definition . as a consequence ,
in addition to  correcting the reported receiving and sending country migration data for differences in definition and coverage , our method contributes to producing estimates that tend to fluctuate less strongly over time .
since the uk uses a general purpose passenger survey , the reported flows fluctuate considerably over time .
moreover , flows to some ( smaller ) countries may not be observed in some years .
we believe our method produces more stable estimates of migration flows for the uk ( and other countries relying on sample data ) .
this implies that the sample survey used for estimating migration to and from the uk provides a reasonably reliable estimate of total migration flows on average , but that the annual estimates are affected by sizeable random fluctuations .
the adjustment factors shown in table  5 can be used to adjust migration numbers to and from countries not included in the matrix , so that total immigration and emigration numbers and total net migration can be estimated for the 19 countries in this study . before doing so
, one - first has to make sure that the share of unknowns in the migration statistics can be distributed evenly across all origins or destinations .
thus , for estimating total immigration and emigration numbers , the adjustment factors should be applied to total migration numbers excluding unknowns .
raymer ( 2008 ) developed a two - step estimation method for countries with missing data ( see also de beer et al .
the first step estimates missing immigration and emigration totals based on harmonised migration flows and covariate information .
destination interaction patterns of the harmonised migration flows and covariate information to estimate the missing interaction patterns .
this estimation step takes into account the fact that migration is relatively high , for example , between neighbouring countries and countries belonging to a similar language group .
finally , work is currently being carried out to integrate harmonisation and estimation of missing data into a single ( bayesian ) model that also includes measures of uncertainty and expert judgements .
the integrated modelling of european migration ( imem ) project recently funded by new opportunities for research funding agency co - operation in europe ( norface ) is expected to develop such a model ( see http://www.norface.org/migration12.html ) over the next couple of years .
we hope this study will provide an important foundation for work such as this , and other projects aiming to improve our knowledge and understanding of the complexity of international migration .
 due to differences in definitions and measurement methods , cross - country comparisons of international migration patterns are difficult and confusing . 
 emigration numbers reported by sending countries tend to differ from the corresponding immigration numbers reported by receiving countries . in this paper , a methodology 
 is presented to achieve harmonised estimates of migration flows benchmarked to a specific definition of duration . 
 this methodology accounts for both differences in definitions and the effects of measurement error due to , for example , under reporting and sampling fluctuations . 
 more specifically , the differences between the two sets of reported data are overcome by estimating a set of adjustment factors for each country s immigration and emigration data . 
 the adjusted data take into account any special cases where the origin  destination patterns do not match the overall patterns . 
 the new method for harmonising migration flows that we present is based on earlier efforts by poulain ( european journal of population , 9(4 ) : 353381 1993 , working paper 12 , joint ece - eurostat work session on migration statistics , geneva , switzerland 1999 ) and is illustrated for movements between 19 european countries from 2002 to 2007 . 
 the results represent a reliable and consistent set of international migration flows that can be used for understanding recent changes in migration patterns , as inputs into population projections and for developing evidence - based migration policies . 
Introduction
Comparability of International Migration Data
Method
Data
Results
Discussion