Maximum entropy classification for record linkage
Section 2. Problems with the classical approach

Suppose that we have two data files A MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyqaaaa@36AF@  and B MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOqaaaa@36B0@  that are believed to have many common entities but no duplicates within each file. Any record in A MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyqaaaa@36AF@  and another one in B MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOqaaaa@36B0@  may or may not refer to the same entity. Our goal is to find the true matches among all possible pairs of the two data files. Let the bipartite comparison space Ω=A×B=MU MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuyQdCLaaG jbVlaai2dacaaMe8UaamyqaiaaysW7cqGHxdaTcaaMe8UaamOqaiaa ysW7caaI9aGaaGjbVlaad2eacaaMe8UaeyOkIGSaaGjbVlaadwfaaa a@4C5D@  consist of matches M MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamytaaaa@36BB@  and non-matches U MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyvaaaa@36C3@  between the records in files A MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyqaaaa@36AF@  and B. MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOqaiaac6 caaaa@3762@  For any pair of records ( a,b )Ω, MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaeWabeaaca WGHbGaaGilaiaaysW7caWGIbaacaGLOaGaayzkaaGaaGjbVlabgIGi olaaysW7cqqHPoWvcaGGSaaaaa@425F@  let γ ab MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaC4SdmaaBa aaleaacaWGHbGaamOyaaqabaaaaa@3921@  be the comparison vector between a set of key variables associated with aA MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyaiaays W7cqGHiiIZcaaMe8Uaamyqaaaa@3C33@  and bB, MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOyaiaays W7cqGHiiIZcaaMe8UaamOqaiaacYcaaaa@3CE5@  respectively, such as name, sex, date of birth. The key variables and the comparison vector γ ab MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaC4SdmaaBa aaleaacaWGHbGaamOyaaqabaaaaa@3921@  are fully observed over Ω. MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuyQdCLaai Olaaaa@3829@  In cases where the key variables may be affected by errors, a match ( a,b ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaeWabeaaca WGHbGaaGilaiaaysW7caWGIbaacaGLOaGaayzkaaaaaa@3B83@  may not have complete agreement in terms of γ ab , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaC4SdmaaBa aaleaacaWGHbGaamOyaaqabaGccaGGSaaaaa@39DB@  and a non-match ( a,b ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaeWabeaaca WGHbGaaGilaiaaysW7caWGIbaacaGLOaGaayzkaaaaaa@3B83@  can nevertheless agree on some (even all) of the key variables.

In the classical approach of Fellegi and Sunter (1969), one recognizes the probabilistic nature of γ ab MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaC4SdmaaBa aaleaacaWGHbGaamOyaaqabaaaaa@3921@  due to the perturbations that cause key-variable errors. The related methods are referred to as probabilistic record linkage. To explain the probabilistic record linkage method of Fellegi and Sunter (1969), let m( γ ab )=f( γ ab | ( a,b )M ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyBaiaayk W7daqadeqaaiaaho7adaWgaaWcbaGaamyyaiaadkgaaeqaaaGccaGL OaGaayzkaaGaaGjbVlaai2dacaaMe8UaamOzaiaaykW7daqadeqaai aaho7adaWgaaWcbaGaamyyaiaadkgaaeqaaOGaaGPaVpaaeeqabaGa aGPaVpaabmqabaGaamyyaiaaiYcacaaMe8UaamOyaaGaayjkaiaawM caaaGaay5bSdGaaGjbVlabgIGiolaaysW7caWGnbaacaGLOaGaayzk aaaaaa@580A@  be the probability mass function of the discrete values γ ab MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaC4SdmaaBa aaleaacaWGHbGaamOyaaqabaaaaa@3921@  can take given ( a,b )M. MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaeWabeaaca WGHbGaaGilaiaaysW7caWGIbaacaGLOaGaayzkaaGaaGjbVlabgIGi olaaysW7caWGnbGaaiOlaaaa@41A5@  Similarly, we can define u( γ ab )=f( γ ab | ( a,b )U ). MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyDaiaayk W7daqadeqaaiaaho7adaWgaaWcbaGaamyyaiaadkgaaeqaaaGccaGL OaGaayzkaaGaaGjbVlaai2dacaaMe8UaamOzaiaaykW7daqadeqaai aaho7adaWgaaWcbaGaamyyaiaadkgaaeqaaOGaaGPaVpaaeeqabaGa aGPaVpaabmqabaGaamyyaiaaiYcacaWGIbaacaGLOaGaayzkaaaaca GLhWoacaaMe8UaeyicI4SaaGjbVlaadwfaaiaawIcacaGLPaaacaGG Uaaaaa@573F@  The ratio

r ab = m( γ ab ) u( γ ab ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOCamaaBa aaleaacaWGHbGaamOyaaqabaGccaaMe8UaaGypaiaaysW7daWcaaqa aiaad2gacaaMc8+aaeWabeaacaWHZoWaaSbaaSqaaiaadggacaWGIb aabeaaaOGaayjkaiaawMcaaaqaaiaadwhacaaMc8+aaeWabeaacaWH ZoWaaSbaaSqaaiaadggacaWGIbaabeaaaOGaayjkaiaawMcaaaaaaa a@4B6D@

is then the basis of the likelihood ratio test (LRT) for H 0 :( a,b )M MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamisamaaBa aaleaacaaIWaaabeaakiaayIW7caaI6aGaaGjbVlaaykW7daqadeqa aiaadggacaaISaGaaGjbVlaadkgaaiaawIcacaGLPaaacaaMe8Uaey icI4SaaGjbVlaad2eaaaa@481D@  vs. H 1 :( a,b )U. MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamisamaaBa aaleaacaaIXaaabeaakiaayIW7caaI6aGaaGjbVlaaykW7daqadeqa aiaadggacaaISaGaaGjbVlaadkgaaiaawIcacaGLPaaacaaMe8Uaey icI4SaaGjbVlaadwfacaGGUaaaaa@48D8@  Let M * ={ ( a,b ): r ab > c M } MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamytamaaCa aaleqabaGaaiOkaaaakiaaysW7caaI9aGaaGjbVpaacmqabaWaaeWa beaacaWGHbGaaGilaiaaysW7caWGIbaacaGLOaGaayzkaaGaaGPaVl aaiQdacaaMe8UaaGjbVlaadkhadaWgaaWcbaGaamyyaiaadkgaaeqa aOGaaGjbVlabg6da+iaaysW7caWGJbWaaSbaaSqaaiaad2eaaeqaaa GccaGL7bGaayzFaaaaaa@51C2@  be the pairs classified as matches and U * ={ ( a,b ): r ab < c U } MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyvamaaCa aaleqabaGaaiOkaaaakiaaysW7caaI9aGaaGjbVpaacmqabaWaaeWa beaacaWGHbGaaGilaiaaysW7caWGIbaacaGLOaGaayzkaaGaaGPaVl aaiQdacaaMe8UaaGPaVlaadkhadaWgaaWcbaGaamyyaiaadkgaaeqa aOGaaGjbVlabgYda8iaaysW7caWGJbWaaSbaaSqaaiaadwfaaeqaaa GccaGL7bGaayzFaaaaaa@51CC@  the non-matches, the remaining pairs are classified by clerical review, where ( c M , c U ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaeWabeaaca WGJbWaaSbaaSqaaiaad2eaaeqaaOGaaGilaiaaysW7caWGJbWaaSba aSqaaiaadwfaaeqaaaGccaGLOaGaayzkaaaaaa@3D9E@  are the thresholds related to the probabilities of false links (of pairs in U) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyvaiaacM caaaa@3770@  and false non-links (of pairs in M), MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamytaiaacM cacaGGSaaaaa@3818@  respectively, defined as

μ= γ u( γ )δ( M * ;γ ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiVd0MaaG jbVlaaykW7caaI9aGaaGjbVlaaykW7daaeqbqaaiaaykW7caWG1bGa aGPaVpaabmqabaGaaC4SdaGaayjkaiaawMcaaiaaysW7cqaH0oazca aMc8+aaeWabeaacaWGnbWaaWbaaSqabeaacaGGQaaaaOGaai4oaiaa ysW7caaMc8UaaC4SdaGaayjkaiaawMcaaaWcbaGaaC4Sdaqab0Gaey yeIuoaaaa@55E4@    and    λ= γ m( γ )δ( U * ;γ ), (2.1) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4UdWMaaG jbVlaaykW7caaI9aGaaGjbVlaaykW7daaeqbqaaiaaykW7caWGTbGa aGPaVpaabmqabaGaaC4SdaGaayjkaiaawMcaaiaaysW7cqaH0oazca aMc8+aaeWabeaacaWGvbWaaWbaaSqabeaacaGGQaaaaOGaaG4oaiaa ysW7caaMc8UaaC4SdaGaayjkaiaawMcaaiaacYcaaSqaaiaaho7aae qaniabggHiLdGccaaMf8UaaGzbVlaaywW7caaMf8UaaGzbVlaacIca caaIYaGaaiOlaiaaigdacaGGPaaaaa@61EA@

where δ( M * ;γ )=1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiTdqMaaG PaVpaabmqabaGaamytamaaCaaaleqabaGaaiOkaaaakiaaiUdacaaM e8UaaGPaVlaaho7aaiaawIcacaGLPaaacaaMe8UaaGypaiaaysW7ca aIXaaaaa@4612@  if γ ab =γ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaC4SdmaaBa aaleaacaWGHbGaamOyaaqabaGccaaMe8UaaGypaiaaysW7caWHZoaa aa@3E4B@  means ( a,b ) M * MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaeWabeaaca WGHbGaaGilaiaaysW7caWGIbaacaGLOaGaayzkaaGaaGjbVlabgIGi olaaysW7caWGnbWaaWbaaSqabeaacaGGQaaaaaaa@41CE@  and 0 otherwise, similarly for δ( U * ;γ ). MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiTdqMaaG PaVpaabmqabaGaamyvamaaCaaaleqabaGaaiOkaaaakiaaiUdacaaM e8UaaGPaVlaaho7aaiaawIcacaGLPaaacaGGUaaaaa@4230@

In practice the probabilities m( γ ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyBaiaayk W7daqadeqaaiaaho7aaiaawIcacaGLPaaaaaa@3B2F@  and u( γ ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyDaiaayk W7daqadeqaaiaaho7aaiaawIcacaGLPaaaaaa@3B37@  are unknown. Neither is the prevalence of true matches, given by π= | M |/ | Ω | := n M /n . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiWdaNaaG jbVlaai2dacaaMe8+aaSGbaeaadaabdeqaaiaaykW7caWGnbGaaGPa VdGaay5bSlaawIa7aiaaykW7aeaacaaMc8+aaqWabeaacaaMc8Uaeu yQdCLaaGPaVdGaay5bSlaawIa7aaaacaaMe8UaaGPaVlaaiQdacaaI 9aGaaGjbVlaaykW7daWcgaqaaiaad6gadaWgaaWcbaGaamytaaqaba aakeaacaWGUbaaaiaac6caaaa@58F6@  Let η MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaC4Tdaaa@372C@  be the set containing π MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiWdahaaa@37A6@  and the unknown parameters of m( γ ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyBaiaayk W7daqadeqaaiaaho7aaiaawIcacaGLPaaaaaa@3B2F@  and u( γ ). MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyDaiaayk W7daqadeqaaiaaho7aaiaawIcacaGLPaaacaGGUaaaaa@3BE9@  Let g ab =1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4zamaaBa aaleaacaWGHbGaamOyaaqabaGccaaMe8UaaGypaiaaysW7caaIXaaa aa@3D74@  if ( a,b )M MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaeWabeaaca WGHbGaaGilaiaaysW7caWGIbaacaGLOaGaayzkaaGaaGjbVlabgIGi olaaysW7caWGnbaaaa@40F3@  and 0 if ( a,b )U. MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaeWabeaaca WGHbGaaGilaiaaysW7caWGIbaacaGLOaGaayzkaaGaaGjbVlabgIGi olaaysW7caWGvbGaaiOlaaaa@41AD@  Given the complete data { ( g ab , γ ab ):( a,b )Ω }, MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaiWabeaada qadeqaaiaadEgadaWgaaWcbaGaamyyaiaadkgaaeqaaOGaaGilaiaa ysW7caWHZoWaaSbaaSqaaiaadggacaWGIbaabeaaaOGaayjkaiaawM caaiaaiQdacaaMe8UaaGPaVpaabmqabaGaamyyaiaaiYcacaaMe8Ua amOyaaGaayjkaiaawMcaaiaaysW7cqGHiiIZcaaMe8UaeuyQdCfaca GL7bGaayzFaaGaaiilaaaa@526B@  Winkler (1988) and Jaro (1989) assume the log-likelihood to be

h( η )= ( a,b )Ω g ab log( πm( γ ab ) )+ ( a,b )Ω ( 1 g ab ) log( ( 1π )u( γ ab ) ).(2.2) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiAaiaayk W7daqadeqaaiaahE7aaiaawIcacaGLPaaacaaMe8UaaGypaiaaysW7 daaeqbqaaiaaykW7caWGNbWaaSbaaSqaaiaadggacaWGIbaabeaaae aadaqadeqaaiaadggacaaISaGaaGPaVlaadkgaaiaawIcacaGLPaaa caaMe8UaeyicI4SaaGjbVlabfM6axbqab0GaeyyeIuoakiGacYgaca GGVbGaai4zaiaaykW7daqadeqaaiabec8aWjaad2gacaaMc8+aaeWa beaacaWHZoWaaSbaaSqaaiaadggacaWGIbaabeaaaOGaayjkaiaawM caaaGaayjkaiaawMcaaiaaysW7cqGHRaWkcaaMe8+aaabuaeaacaaM c8+aaeWabeaacaaIXaGaaGjbVlabgkHiTiaaysW7caWGNbWaaSbaaS qaaiaadggacaWGIbaabeaaaOGaayjkaiaawMcaaaWcbaWaaeWabeaa caWGHbGaaGilaiaaykW7caWGIbaacaGLOaGaayzkaaGaaGjbVlabgI GiolaaysW7cqqHPoWvaeqaniabggHiLdGccaaMe8UaciiBaiaac+ga caGGNbGaaGPaVpaabmqabaWaaeWabeaacaaIXaGaaGjbVlabgkHiTi aaysW7cqaHapaCaiaawIcacaGLPaaacaaMe8UaamyDaiaaykW7daqa deqaaiaaho7adaWgaaWcbaGaamyyaiaadkgaaeqaaaGccaGLOaGaay zkaaaacaGLOaGaayzkaaGaaiOlaiaaywW7caaMf8UaaGzbVlaaywW7 caaMf8UaaiikaiaaikdacaGGUaGaaGOmaiaacMcaaaa@A156@

An EM-algorithm follows by treating g Ω ={ g ab :( a,b )Ω } MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbf9F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4zamaaBa aaleaacqqHPoWvaeqaaOGaaGjbVlaai2dacaaMe8+aaiWabeaacaWG NbWaaSbaaSqaaiaadggacaWGIbaabeaakiaayIW7caaI6aGaaGjbVl aaykW7daqadeqaaiaadggacaaISaGaaGjbVlaadkgaaiaawIcacaGL PaaacaaMe8UaeyicI4SaaGjbVlabfM6axbGaay5Eaiaaw2haaaaa@52CE@  as the missing data.


Date modified: