Estimating the false negatives due to blocking in record linkage
Section 6. Empirical study

The empirical study is based on staffing data from the Public Service Resourcing System (PSRS), which is used by applicants to the federal public service in Canada. A given user may open many accounts and apply to many jobs using the same account; each account being associated with a distinct email address. To fulfill its mandate, the Public Service Commission needs to identify all accounts from a given applicant. However this is a challenge because there is no unique identifier except for a minority of applicants. Instead, for most records, the linkage must be based on the given name, the surname and the partial birthdate, which are available for all records. The partial birthdate is comprised of the day and month of birth along with the last digit of the birth year.

The empirical study is based on a subset of 126,330 records selected from the PSRS data since 2006. The selection is based on the following criteria.

The selected records represent 63,155 distinct values of the identifier and so many distinct individuals, with two matched records per individual. These records are split into two complete and duplicate-free registers that are linked with the following blocking criteria, and without the unique identifier. A pair is selected if the partial birthdate is the same and the SOUNDEX code (Herzog et al., 2007, Chapter 11) is the same for the given name or the surname. The expected error rates are estimated with the model and compared with the actual values based on the unique identifier.

In Figure 6.1, the histogram shows that the vast majority or records have exactly one neighbour. However 1,659 records have no neighbour, while five records have five neighbours; the maximum number of neighbours of any record.

Figure 6.1 Histogram of the number of neighbours

Description for Figure 6.1

Histogram of the number of neighbours showing that the vast majority or records have exactly one neighbour. However 1,659 records have no neighbour, while five records have five neighbours; the maximum number of neighbours of any record.

Table 6.1 cross-classifies the records by their number of neighbours and linkage errors, in agreement with Table 3.1.


Table 6.1
Number of neighbours and errors
Table summary
This table displays the results of Number of neighbours and errors. The information is grouped by Neighbours (équation) (appearing as row headers), False negatives, False positives and Freq. (appearing as column headers).
Neighbours ( n i ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacPqpw0le9v8qqaqpepeeaY= Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq 0=yq=He9q8qqQ8frFve9Fve9Ff0dmeqabeqadiWaceGabeqabeWade qaeeaakeaaqaaaaaaaaaWdbiaacIcacaWGUbWdamaaBaaaleaapeGa amyAaaWdaeqaaOGaaiykaaaa@3FD3@ False negatives False positives Freq.
0 1 0 1,659
1 1 1 116
1 0 0 53,835
2 1 2 8
2 0 1 6,867
3 1 3 1
3 0 2 602
4 0 3 62
5 0 4 5

The confusion matrix is as follows.


Table 6.2
Confusion matrix
Table summary
This table displays the results of Confusion matrix Link, Non-link and Total (appearing as column headers).
Link Non-link Total
Matched 61,371 1,784 63,155
Unmatched 8,412 3.99E9 3.99E9
Total 69,783 3.988E9 3.989E9

From this matrix, FNR = 1 , 784 / 63 , 155 = 0.0282 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaciOraiaac6eacaGGsbGaeyypa0Ja aGymaiaacYcacaaI3aGaaGioaiaaisdacaGGVaGaaGOnaiaaiodaca GGSaGaaGymaiaaiwdacaaI1aGaeyypa0JaaGimaiaac6cacaaIWaGa aGOmaiaaiIdacaaIYaaaaa@4CF7@ and FPR = 8 , 412 / 3.99 E 9 = 2.11 E 6. MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaaeOraiaabcfacaqGsbGaeyypa0Ja aGioaiaacYcacaaI0aGaaGymaiaaikdacaGGVaGaaG4maiaac6caca aI5aGaaGyoaiaabweacaaI5aGaeyypa0JaaGOmaiaac6cacaaIXaGa aGymaiaabweacqGHsislcaaI2aGaaiOlaaaa@4EB4@ Both measures may be viewed as the estimators E ^ [ FNR ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGabmyra8aagaqca8qacaaMc8Uaai4w aiaabAeacaqGobGaaeOuaiaac2faaaa@4205@ and E ^ [ FPR ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGabmyra8aagaqca8qacaaMc8Uaai4w aiaabAeacaqGqbGaaeOuaiaac2faaaa@4207@ of their respective expectations. Since the false negative rate is the summation of independent and identically distributed random variables, its variance may be estimated by

v a ^ r ( FNR ) = 1 N ( N 1 ) i = 1 N ( 1 n i | M FNR ) 2 , MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaiGacAhaceGGHbGbaKaacaGGYbGaaGPaVlaacIcaqaaa aaaaaaWdbiaabAeacaqGobGaaeOua8aacaGGPaGaaGjbVlaaykW7cq GH9aqpcaaMe8UaaGPaV=qadaWcaaWdaeaapeGaaGymaaWdaeaapeGa amOtaiaaykW7caGGOaGaamOtaiabgkHiTiaaigdacaGGPaaaaiaays W7daaeWbqaaiaaykW7caGGOaGaaGymaiabgkHiTiaad6gapaWaaSba aSqaa8qacaWGPbGaaGjbVlaacYhacaaMc8UaamytaaWdaeqaaOWdbi abgkHiTiaabAeacaqGobGaaeOuaiaacMcadaahaaWcbeqaaiaaikda aaGccaGGSaaaleaacaWGPbGaeyypa0JaaGymaaqaaiaad6eaa0Gaey yeIuoaaaa@6A36@

based on the latent variables n 1 | M , , n m | M , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOBa8aadaWgaaWcbaWdbiaaigda caaMe8UaaiiFaiaaysW7caWGnbaapaqabaGcpeGaaiilaiaaysW7cq WIMaYscaGGSaGaaGjbVlaad6gapaWaaSbaaSqaa8qacaWGTbGaaGjb VlaacYhacaaMe8UaamytaaWdaeqaaOGaaiilaaaa@4FE1@ which are not directly observed in practice. As a result, the estimated FNR variance is v a ^ r ( FNR ) = 4.35 E 7. MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaiGacAhaceGGHbGbaKaacaGGYbGaaGPaVlaacIcacaqG gbGaaeOtaiaabkfacaGGPaaeaaaaaaaaa8qacqGH9aqpcaaI0aGaai OlaiaaiodacaaI1aGaaeyraiabgkHiTiaaiEdacaGGUaaaaa@4AA6@ This means the estimated standard error S ^ E( E ^ [ FNR ] )=6.6E4 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGabmyra8aagaqcaiaabweacaaMe8+a aeWaaeaapeGabmyra8aagaqca8qacaaMe8+aamWaaeaacaqGgbGaae OtaiaabkfaaiaawUfacaGLDbaaa8aacaGLOaGaayzkaaGaaGjbV=qa cqGH9aqpcaaMe8UaaGOnaiaac6cacaaI2aGaaeyraiabgkHiTiaais daaaa@4FF4@ for the estimator E ^ [ FNR], MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGabmyra8aagaqca8qacaaMc8Uaai4w aiaabAeacaqGobGaaeOuaiaab2facaqGSaaaaa@42B3@ and the 95% normal confidence interval E ^ [ FNR ] z α/2 S ^ E( E ^ [ FNR ] )=( 2.82E21.3E3 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGabmyra8aagaqca8qacaaMe8+aamWa aeaacaqGgbGaaeOtaiaabkfaaiaawUfacaGLDbaacaaMe8UaeS4eI0 MaaGjbVlaadQhapaWaaSbaaSqaa8qacqaHXoqycaGGVaGaaGOmaaWd aeqaaOWdbiqadweapaGbaKaacaqGfbGaaGjbVpaabmaabaWdbiqadw eapaGbaKaapeGaaGPaVpaadmaabaGaaeOraiaab6eacaqGsbaacaGL BbGaayzxaaaapaGaayjkaiaawMcaaiaaysW7cqGH9aqpcaaMe8+aae WaaeaapeGaaGOmaiaac6cacaaI4aGaaGOmaiaabweacaaMe8UaeyOe I0IaaGjbVlaaikdacaaMe8UaeS4eI0MaaGjbVlaaigdacaGGUaGaaG 4maiaabweacaaMe8UaeyOeI0IaaGjbVlaaiodaa8aacaGLOaGaayzk aaaaaa@7114@ for the expected FNR, where α = MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaeqySdeMaeyypa0JaaGPaVdaa@3F82@ 0.05 and z α / 2 = MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOEa8aadaWgaaWcbaWdbiabeg7a Hjaac+cacaaIYaaapaqabaGcpeGaeyypa0JaaGPaVdaa@4264@ 1.96. The corresponding 99% confidence interval is ( 2.82 E 2 1.71 E 3 ) . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaaiikaiaaikdacaGGUaGaaGioaiaa ikdacaqGfbGaaGjbVlabgkHiTiaaysW7caaIYaGaaGjbVlabloHiTj aaysW7caaIXaGaaiOlaiaaiEdacaaIXaGaaeyraiaaysW7cqGHsisl caaMe8UaaG4maiaacMcacaGGUaaaaa@5296@ Estimating the FPR variance is more difficult because the FPR involves a second order U statistic (Hoeffding, 1948; Lee, 1990). As a matter of fact, Table 6.1 does not give enough information to estimate this statistic. Estimating the variance of the model-based estimators is also challenging because the n i s MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOBa8aadaWgaaWcbaWdbiaadMga a8aabeaaieaakiaa=LbicaqGZbaaaa@3F50@ are correlated. All the point estimates are given in Table 6.3, where the first row gives the actual FNR and FPR.


Table 6.3
Point estimates
Table summary
This table displays the results of Point estimates (équation), calculated using G = 3, 0.0303 and 2.14E-6 units of measure (appearing as column headers).
E ^ [FNR] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacPqpw0le9v8qqaqpepeeaY= Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq 0=yq=He9q8qqQ8frFve9Fve9Ff0dmeqabeqadiWaceGabeqabeWade qaeeaakeaaqaaaaaaaaaWdbiqadweapaGbaKaapeGaaGPaVlaacUfa caqGgbGaaeOtaiaabkfacaGGDbaaaa@42E8@ E ^ [FPR] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacPqpw0le9v8qqaqpepeeaY= Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq 0=yq=He9q8qqQ8frFve9Fve9Ff0dmeqabeqadiWaceGabeqabeWade qaeeaakeaaqaaaaaaaaaWdbiqadweapaGbaKaapeGaaGPaVlaacUfa caqGgbGaaeiuaiaabkfacaGGDbaaaa@42EA@
Unique id 0.0282 2.11E-6
Model G = 1 0.0301 2.14E-6
G = 2 0.0298 2.13E-6
G = 3 0.0303 2.14E-6

The results show that the model based estimates are very close to the actual FNR and FPR when using one, two or three classes. For the false negative rate, the relative error is 100 × | 0 .0303 0 .0282 | / 0 .0282 = 7 .45 % , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpu0de9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaeaaaaaaaaa8 qacaaIXaGaaGimaiaaicdacaaMe8Uaey41aqRaaGjbVpaalyaabaWa aqWabeaacaaMc8Uaaeimaiaab6cacaqGWaGaae4maiaabcdacaqGZa GaaGjbVlabgkHiTiaaysW7caqGWaGaaeOlaiaabcdacaqGYaGaaeio aiaabkdacaaMc8oacaGLhWUaayjcSdGaaGPaVdqaaiaaykW7caqGWa GaaeOlaiaabcdacaqGYaGaaeioaiaabkdacaaMe8Uaeyypa0JaaGjb VlaabEdacaqGUaGaaeinaiaabwdacaGGLaaaaiaacYcaaaa@6203@ while this relative error is 100 × | 2 .11 2 .14 | / 2 .11 = 1 .42 % MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpu0de9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaeaaaaaaaaa8 qacaaIXaGaaGimaiaaicdacaaMe8Uaey41aqRaaGjbVpaalyaabaWa aqWabeaacaaMc8UaaeOmaiaab6cacaqGXaGaaeymaiaaysW7cqGHsi slcaaMe8UaaeOmaiaab6cacaqGXaGaaeinaiaaykW7aiaawEa7caGL iWoacaaMc8oabaGaaGPaVlaabkdacaqGUaGaaeymaiaabgdaaaGaaG jbVlabg2da9iaaysW7caqGXaGaaeOlaiaabsdacaqGYaGaaiyjaaaa @5D09@ for the false positive rate. The small relative errors are encouraging regarding the accuracy of the proposed estimators, even if the model estimates of the expected FNR lie slightly outside the 95% confidence interval. However, the estimate belongs to the 99% confidence interval when using two classes. Choosing two classes seems optimal because the resulting estimate has the smallest relative error with respect to the actual FNR.


Date modified: