Estimating the false negatives due to blocking in record linkage
Section 5. Estimation procedure

The model parameters may be estimated by maximizing the composite likelihood (Varin, Reid and Firth, 2011) of the sample n 1 , , n m . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOBa8aadaWgaaWcbaWdbiaaigda a8aabeaak8qacaGGSaGaaGjbVlabgAci8kaacYcacaaMe8UaamOBa8 aadaWgaaWcbaWdbiaad2gaa8aabeaakiaac6caaaa@4677@ For brevity, this composite likelihood is subsequently called likelihood. To develop the EM procedure (Dempster, Laird and Rubin, 1977) it is convenient to first derive the maximum likelihood (ML) equations for the complete data, which are comprised of the latent variables n i | M , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOBa8aadaWgaaWcbaWdbiaadMga caaMe8UaaiiFaiaaykW7caaMc8UaamytaaWdaeqaaOGaaiilaaaa@44BC@ n i | U MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOBa8aadaWgaaWcbaWdbiaadMga caaMe8UaaiiFaiaaykW7caaMc8UaamyvaaWdaeqaaaaa@440A@ and ( c i 1 , , c i G ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaaiikaiaadogapaWaaSbaaSqaa8qa caWGPbGaaGymaaWdaeqaaOWdbiaacYcacaaMe8UaeyOjGWRaaiilai aaysW7caWGJbWdamaaBaaaleaapeGaamyAaiaadEeaa8aabeaak8qa caGGPaaaaa@48CE@ for each i ; MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamyAaiaacUdaaaa@3CFF@ c i g MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaam4ya8aadaWgaaWcbaWdbiaadMga caWGNbaapaqabaaaaa@3E6E@ being the indicator that record i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamyAaaaa@3C40@ is from class  g . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaam4zaiaac6caaaa@3CF0@

After some algebra, the ML equations for the complete data are as follows.

p ^ g = i = 1 m c i g n i | M i = 1 m c i g λ ^ g = i = 1 m c i g n i | U i = 1 m c i g α ^ g = 1 m i = 1 m c i g . ( 5.1 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaauaabeqadmaaaeaaqaaaaaaaaaWdbiqadchapaGbaKaa daWgaaWcbaWdbiaadEgaa8aabeaaaOqaa8qacqGH9aqpa8aabaWdbm aalaaapaqaa8qadaaeWaqaaiaaykW7caWGJbWdamaaBaaaleaapeGa amyAaiaadEgaa8aabeaak8qacaWGUbWdamaaBaaaleaapeGaamyAai aaysW7caGG8bGaaGPaVlaad2eaa8aabeaaa8qabaGaamyAaiabg2da 9iaaigdaaeaacaWGTbaaniabggHiLdaak8aabaWdbmaaqadabaGaaG PaVlaadogapaWaaSbaaSqaa8qacaWGPbGaam4zaaWdaeqaaaWdbeaa caWGPbGaeyypa0JaaGymaaqaaiaad2gaa0GaeyyeIuoaaaaak8aaba WdbiqbeU7aS9aagaqcamaaBaaaleaapeGaam4zaaWdaeqaaaGcbaWd biabg2da9aWdaeaapeWaaSaaa8aabaWdbmaaqadabaGaaGPaVlaado gapaWaaSbaaSqaa8qacaWGPbGaam4zaaWdaeqaaOWdbiaad6gapaWa aSbaaSqaa8qacaWGPbGaaGjbVlaacYhacaaMc8UaamyvaaWdaeqaaa WdbeaacaWGPbGaeyypa0JaaGymaaqaaiaad2gaa0GaeyyeIuoaaOWd aeaapeWaaabmaeaacaaMc8Uaam4ya8aadaWgaaWcbaWdbiaadMgaca WGNbaapaqabaaapeqaaiaadMgacqGH9aqpcaaIXaaabaGaamyBaaqd cqGHris5aaaaaOWdaeaapeGafqySde2dayaajaWaaSbaaSqaa8qaca WGNbaapaqabaaakeaapeGaeyypa0dapaqaa8qadaWcaaWdaeaapeGa aGymaaWdaeaapeGaamyBaaaacaaMc8+aaabCaeaacaaMc8Uaam4ya8 aadaWgaaWcbaWdbiaadMgacaWGNbaapaqabaaapeqaaiaadMgacqGH 9aqpcaaIXaaabaGaamyBaaqdcqGHris5aOGaaiOlaaaapaGaaGzbVl aaywW7caaMf8UaaGzbVlaacIcacaaI1aGaaiOlaiaaigdacaGGPaaa aa@98D0@

Consequently the ML equations for the observed data (the n i s ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOBa8aadaWgaaWcbaWdbiaadMga a8aabeaaieaakiaa=LbicaqGZbGaaiykaaaa@3FFD@ are as follows.

p ^ g = i = 1 m E [ c i g n i | M | n i ; ψ ] i = 1 m E [ c i g | n i ; ψ ] λ ^ g = i = 1 m E [ c i g n i | U | n i ; ψ ] i = 1 m E [ c i g | n i ; ψ ] α ^ g = 1 m i = 1 m E [ c i g | n i ; ψ ] . ( 5.2 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaauaabaqadiaaaeaaqaaaaaaaaaWdbiqadchapaGbaKaa daWgaaWcbaWdbiaadEgaa8aabeaaaOqaaiabg2da98qadaWcaaWdae aapeWaaabmaeaacaWGfbGaaGPaVlaacUfacaWGJbWdamaaBaaaleaa peGaamyAaiaadEgaa8aabeaak8qacaWGUbWdamaaBaaaleaapeGaam yAaiaaysW7caGG8bGaaGPaVlaad2eaa8aabeaak8qadaabbaWdaeaa peGaaGPaVlaad6gapaWaaSbaaSqaa8qacaWGPbaapaqabaGcpeGaai 4oaiabeI8a5bGaay5bSdGaaiyxaaWcbaGaamyAaiabg2da9iaaigda aeaacaWGTbaaniabggHiLdaak8aabaWdbmaaqadabaGaamyraiaayk W7caGGBbGaam4ya8aadaWgaaWcbaWdbiaadMgacaWGNbaapaqabaGc peWaaqqaa8aabaWdbiaaykW7caWGUbWdamaaBaaaleaapeGaamyAaa WdaeqaaOWdbiaacUdacqaHipqEaiaawEa7aiaac2faaSqaaiaadMga cqGH9aqpcaaIXaaabaGaamyBaaqdcqGHris5aaaaaOWdaeaapeGafq 4UdW2dayaajaWaaSbaaSqaa8qacaWGNbaapaqabaaakeaacqGH9aqp peWaaSaaa8aabaWdbmaaqadabaGaamyraiaaykW7caGGBbGaam4ya8 aadaWgaaWcbaWdbiaadMgacaWGNbaapaqabaGcpeGaamOBa8aadaWg aaWcbaWdbiaadMgacaaMe8UaaiiFaiaaykW7caWGvbaapaqabaGcpe Waaqqaa8aabaWdbiaaykW7caWGUbWdamaaBaaaleaapeGaamyAaaWd aeqaaOWdbiaacUdacqaHipqEaiaawEa7aiaac2faaSqaaiaadMgacq GH9aqpcaaIXaaabaGaamyBaaqdcqGHris5aaGcpaqaa8qadaaeWaqa aiaadweacaaMc8Uaai4waiaadogapaWaaSbaaSqaa8qacaWGPbGaam 4zaaWdaeqaaOWdbmaaeeaapaqaa8qacaaMc8UaamOBa8aadaWgaaWc baWdbiaadMgaa8aabeaak8qacaGG7aGaeqiYdKhacaGLhWoacaGGDb aaleaacaWGPbGaeyypa0JaaGymaaqaaiaad2gaa0GaeyyeIuoaaaaa k8aabaWdbiqbeg7aH9aagaqcamaaBaaaleaapeGaam4zaaWdaeqaaa GcbaGaeyypa0Zdbmaalaaapaqaa8qacaaIXaaapaqaa8qacaWGTbaa aiaaykW7daaeWbqaaiaadweacaaMc8Uaai4waiaadogapaWaaSbaaS qaa8qacaWGPbGaam4zaaWdaeqaaOWdbmaaeeaapaqaa8qacaaMc8Ua amOBa8aadaWgaaWcbaWdbiaadMgaa8aabeaaaOWdbiaawEa7aiaacU dacqaHipqEcaGGDbaaleaacaWGPbGaeyypa0JaaGymaaqaaiaad2ga a0GaeyyeIuoakiaaykW7caGGUaaaa8aacaaMf8UaaGzbVlaaywW7ca aMf8UaaiikaiaaiwdacaGGUaGaaGOmaiaacMcaaaa@CFB6@

The EM procedure alternates between the M-step given by Equation (5.2) and the E-step equations in Appendix A.

The above procedure may produce consistent point estimators even if it treats the sample n 1 , , n m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOBa8aadaWgaaWcbaWdbiaaigda a8aabeaak8qacaGGSaGaaGjbVlabgAci8kaacYcacaaMe8UaamOBa8 aadaWgaaWcbaWdbiaad2gaa8aabeaaaaa@45BB@ as if it were independent and identically distributed. However this is likely to generate some bias when estimating the variance and the critical levels of hypothesis tests.


Date modified: