Estimating the false negatives due to blocking in record linkage
Section 4. Finite mixture model

The linkage of the two sources is of interest when it is a viable option even N MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOtaaaa@3C25@ is very large. To capture the essence of such situations, the following two regularity conditions are assumed.

  1. Two matched records are neighbours with a probability that is bounded away from 0 regardless of  N . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOtaiaac6caaaa@3CD7@
  2. Two unmatched records are accidental neighbours with a probability of O ( 1 / N ) . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaam4taiaaykW7caGGOaGaaGymaiaa c+cacaWGobGaaiykaiaac6caaaa@41FD@

These assumptions imply that each record has a bounded expected number of neighbours and that O ( N ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaam4taiaaykW7caGGOaGaamOtaiaa cMcaaaa@3FDD@ pairs (instead of O ( m N ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaam4taiaaykW7caGGOaGaamyBaiaa d6eacaGGPaaaaa@40CF@ pairs and even O ( N 2 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaam4taiaaykW7caGGOaGaamOta8aa daahaaWcbeqaa8qacaaIYaaaaOGaaiykaaaa@40EF@ pairs if m = O ( N ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamyBaiabg2da9iaad+eacaaMe8Ua aiikaiaad6eacaGGPaaaaa@41D7@ ) are selected by the blocking criteria. They also imply that there is enough linkage information to identify matched records with a success probability, which is bounded away from zero, regardless of the population size. The above assumptions further imply a particular limiting distribution for the number of neighbours n i . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOBa8aadaWgaaWcbaWdbiaadMga a8aabeaakiaac6caaaa@3E49@ Indeed, let n i = n i | M + n i | U , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOBa8aadaWgaaWcbaWdbiaadMga a8aabeaak8qacqGH9aqpcaWGUbWdamaaBaaaleaapeGaamyAaiaays W7caGG8bGaaGPaVlaad2eaa8aabeaak8qacqGHRaWkcaWGUbWdamaa BaaaleaapeGaamyAaiaaysW7caGG8bGaaGPaVlaadwfaa8aabeaaki aacYcaaaa@4EB5@ where n i | M MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOBa8aadaWgaaWcbaWdbiaadMga caaMe8UaaiiFaiaaysW7caWGnbaapaqabaaaaa@4279@ is the number of matched neighbours and n i | U MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOBa8aadaWgaaWcbaWdbiaadMga caaMe8UaaiiFaiaaykW7caWGvbaapaqabaaaaa@427F@ is the number of unmatched neighbours. Note that these latter variables are not directly observed expect when n i = 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOBa8aadaWgaaWcbaWdbiaadMga a8aabeaak8qacqGH9aqpcaaIWaaaaa@3F67@ or n i = N MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOBa8aadaWgaaWcbaWdbiaadMga a8aabeaak8qacqGH9aqpcaWGobaaaa@3F80@ (see Table 3.1). They are also conditionally independent given v i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamODa8aadaWgaaWcbaWdbiaadMga a8aabeaaaaa@3D95@ and such that n i | M | v i ~ Bernoulli ( p ( v i )), MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeWaaqGaa8aabaWdbiaad6gapaWaaSba aSqaa8qacaWGPbGaaGjbVlaacYhacaaMc8UaamytaaWdaeqaaaGcpe GaayjcSdGaaGjbVlaadAhapaWaaSbaaSqaa8qacaWGPbaapaqabaGc peGaaiOFaiaabkeacaqGLbGaaeOCaiaab6gacaqGVbGaaeyDaiaabY gacaqGSbGaaeyAaiaaysW7caqGOaGaaGPaVlaadchacaaMe8Uaaeik aiaadAhapaWaaSbaaSqaa8qacaWGPbaapaqabaGcpeGaaeykaiaabM cacaqGSaaaaa@5CCD@ n i | U | v i ~ Binomial ( N 1 , λ ( v i ) / ( N 1 ) ), MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeWaaqGaa8aabaWdbiaad6gapaWaaSba aSqaa8qacaWGPbGaaGjbVlaacYhacaaMc8UaamyvaaWdaeqaaaGcpe GaayjcSdGaaGjbVlaadAhapaWaaSbaaSqaa8qacaWGPbaapaqabaGc peGaaiOFaiaabkeacaqGPbGaaeOBaiaab+gacaqGTbGaaeyAaiaabg gacaqGSbGaaGjbVlaabIcacaWGobGaeyOeI0IaaGymaiaacYcacaaM e8Uaeq4UdWMaaGjbVlaacIcacaWG2bWdamaaBaaaleaapeGaamyAaa WdaeqaaOWdbiaacMcacaGGVaGaaiikaiaad6eacqGHsislcaaIXaGa aiykaiaabMcacaqGSaaaaa@6446@ if an unmatched record is a neighbour with the probability λ ( v i ) / ( N 1 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaeq4UdWMaaGPaVlaacIcacaWG2bWd amaaBaaaleaapeGaamyAaaWdaeqaaOWdbiaacMcacaGGVaGaaiikai aad6eacqGHsislcaaIXaGaaiykaaaa@46CE@ independently of the other unmatched records. When the functions p ( . ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamiCaiaaykW7caGGOaGaaiOlaiaa cMcaaaa@3FDD@ and λ ( . ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaeq4UdWMaaGPaVlaacIcacaGGUaGa aiykaaaa@409C@ do not depend on N MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOtaaaa@3C25@ and N MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOtaaaa@3C25@ is large, we have n i | U | v i ~ ˙ Poisson ( λ ( v i )) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeWaaqGaa8aabaWdbiaad6gapaWaaSba aSqaa8qacaWGPbGaaGjbVlaacYhacaaMc8UaamyvaaWdaeqaaaGcpe GaayjcSdGaaGjbVlaadAhapaWaaSbaaSqaa8qacaWGPbaapaqabaGc peGabiOFa8aagaGaa8qacaqGqbGaae4BaiaabMgacaqGZbGaae4Cai aab+gacaqGUbGaaGPaVlaabIcacqaH7oaBcaaMc8UaaeikaiaadAha paWaaSbaaSqaa8qacaWGPbaapaqabaGcpeGaaeykaiaabMcaaaa@59B7@ (Billingsley, 1995), where  ~ ˙ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGabiOFa8aagaGaaaaa@3C6C@  means approximately distributed as. Hence, n i | v i ~ ˙ Bernoulli ( p ( v i )) * Poisson ( λ ( v i ) ) , MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeWaaqGaa8aabaWdbiaad6gapaWaaSba aSqaa8qacaWGPbaapaqabaGccaaMe8oapeGaayjcSdGaaGjbVlaadA hapaWaaSbaaSqaa8qacaWGPbaapaqabaGcpeGabiOFa8aagaGaa8qa caqGcbGaaeyzaiaabkhacaqGUbGaae4BaiaabwhacaqGSbGaaeiBai aabMgacaaMc8UaaeikaiaaykW7caWGWbGaaGjbVlaabIcacaWG2bWd amaaBaaaleaapeGaamyAaaWdaeqaaOWdbiaabMcacaqGPaGaaGjbVl aabQcacaaMe8Uaaeiuaiaab+gacaqGPbGaae4CaiaabohacaqGVbGa aeOBaiaaykW7caGGOaGaeq4UdWMaaGjbVlaacIcacaWG2bWdamaaBa aaleaapeGaamyAaaWdaeqaaOWdbiaacMcacaGGPaGaaiilaaaa@6DB8@ where * MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaaiOkaaaa@3C00@ is the convolution operator. Note that, in general, the functions p ( . ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamiCaiaaykW7caGGOaGaaiOlaiaa cMcaaaa@3FDD@ and λ ( . ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaeq4UdWMaaGPaVlaacIcacaGGUaGa aiykaaaa@409C@ are unknown high-dimensional parameters. To simplify, further assume that ( p ( . ) , λ ( . ) ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaaiikaiaaykW7caWGWbGaaGPaVlaa cIcacaGGUaGaaiykaiaacYcacaaMe8Uaeq4UdWMaaGPaVlaacIcaca GGUaGaaiykaiaacMcaaaa@4A48@ is (well approximated by) a piecewise constant function with G MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaam4raaaa@3C1E@ levels, such that we have the finite mixture model n i ~ g = 1 G α g ( Bernoulli ( p g ) * Poisson ( λ g ) ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOBa8aadaWgaaWcbaWdbiaadMga a8aabeaak8qacaGG+bWaaabmaeaacaaMc8UaeqySde2damaaBaaale aapeGaam4zaaWdaeqaaaWdbeaacaWGNbGaeyypa0JaaGymaaqaaiaa dEeaa0GaeyyeIuoakiaaykW7caGGOaGaaeOqaiaabwgacaqGYbGaae OBaiaab+gacaqG1bGaaeiBaiaabYgacaqGPbGaaGPaVlaacIcacaWG WbWdamaaBaaaleaapeGaam4zaaWdaeqaaOWdbiaacMcacaaMe8Uaai OkaiaaysW7caqGqbGaae4BaiaabMgacaqGZbGaae4Caiaab+gacaqG UbGaaGPaVlaacIcacqaH7oaBpaWaaSbaaSqaa8qacaWGNbaapaqaba GcpeGaaiykaiaacMcaaaa@6947@ holds approximately. When G MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaam4raaaa@3C1E@ is fixed, the unknown model parameters are given by the vector ψ = [ ( α g , p g , λ g ) ] 1 g G MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaeqiYdKNaeyypa0ZaamWaa8aabaWd biaacIcacqaHXoqypaWaaSbaaSqaa8qacaWGNbaapaqabaGcpeGaai ilaiaaysW7caWGWbWdamaaBaaaleaapeGaam4zaaWdaeqaaOWdbiaa cYcacaaMe8Uaeq4UdW2damaaBaaaleaapeGaam4zaaWdaeqaaOWdbi aacMcaaiaawUfacaGLDbaapaWaaSbaaSqaa8qacaaIXaGaaGPaVlab gsMiJkaaykW7caWGNbGaaGPaVlabgsMiJkaaykW7caWGhbaapaqaba aaaa@5AD5@ that may be estimated with the Expectation-Maximization (EM) procedure in the next section.

The connection between the error rates and model parameters is made by first noting that the FNR and FPR definitions imply

FNR = 1 m i = 1 m ( 1 n i | M ) , ( N 1 ) FPR = 1 m i = 1 m n i | U . ( 4.1 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaauaabaqaciaaaeaaqaaaaaaaaaWdbiaaywW7caaMf8Ua aGjbVlaaysW7caaMe8UaaeOraiaab6eacaqGsbaapaqaa8qacqGH9a qpcaaMe8UaaGPaVpaalaaapaqaa8qacaaIXaaapaqaa8qacaWGTbaa aiaaysW7daaeWbqaaiaaysW7caGGOaGaaGymaiaaysW7cqGHsislca aMe8UaamOBa8aadaWgaaWcbaWdbiaadMgacaaMe8UaaiiFaiaaykW7 caWGnbaapaqabaGcpeGaaiykaaWcbaGaamyAaiabg2da9iaaigdaae aacaWGTbaaniabggHiLdGccaGGSaaapaqaa8qacaGGOaGaamOtaiaa ysW7cqGHsislcaaMe8UaaGymaiaacMcacaaMe8UaaeOraiaabcfaca qGsbaapaqaa8qacqGH9aqpcaaMe8UaaGPaVpaalaaapaqaa8qacaaI Xaaapaqaa8qacaWGTbaaaiaaysW7daaeWbqaaiaaysW7caWGUbWdam aaBaaaleaapeGaamyAaiaaysW7caGG8bGaaGPaVlaadwfaa8aabeaa a8qabaGaamyAaiabg2da9iaaigdaaeaacaWGTbaaniabggHiLdGcca GGUaaaa8aacaaMf8UaaGzbVlaaywW7caaMf8UaaiikaiaaisdacaGG UaGaaGymaiaacMcaaaa@8ED4@

When m = N MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamyBaiabg2da9iaad6eaaaa@3E1D@ almost surely, the above equations imply that

E [ FNR ] = 1 E [ n i | M ] , = 1 E [ p ( v i ) ] , ( N 1 ) E [ FPR] = E [ n i | U ] = E [ λ ( v i ) ] , ( 4.2 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaauaabaqaeiaaaaqaaabaaaaaaaaapeGaaGzbVlaaywW7 caaMe8UaaGjbVlaaykW7caaMc8UaamyraiaaykW7caGGBbGaaeOrai aab6eacaqGsbGaaiyxaaWdaeaacqGH9aqpcaaMe8UaaGPaV=qacaaI XaGaaGjbVlabgkHiTiaaysW7caWGfbGaaGPaVlaacUfacaWGUbWdam aaBaaaleaapeGaamyAaiaaysW7caGG8bGaaGPaVlaad2eaa8aabeaa k8qacaGGDbGaaiilaaWdaeaaaeaacqGH9aqpcaaMe8UaaGPaV=qaca aIXaGaaGjbVlabgkHiTiaaysW7caWGfbGaaGPaVlaacUfacaWGWbGa aGPaVlaacIcacaWG2bWdamaaBaaaleaapeGaamyAaaWdaeqaaOWdbi aacMcacaGGDbGaaiilaaWdaeaapeGaaiikaiaad6eacqGHsislcaaI XaGaaiykaiaaysW7caWGfbGaaGPaVlaacUfacaqGgbGaaeiuaiaabk facaqGDbaapaqaaiabg2da9iaaysW7caaMc8+dbiaadweacaaMc8Ua ai4waiaad6gapaWaaSbaaSqaa8qacaWGPbGaaGjbVlaacYhacaaMc8 UaamyvaaWdaeqaaOWdbiaac2faa8aabaaabaGaeyypa0JaaGjbVlaa ykW7peGaamyraiaaykW7caGGBbGaeq4UdWMaaGPaVlaacIcacaWG2b WdamaaBaaaleaapeGaamyAaaWdaeqaaOWdbiaacMcacaGGDbGaaiil aaaapaGaaGzbVlaaywW7caaMf8UaaGzbVlaacIcacaaI0aGaaiOlai aaikdacaGGPaaaaa@A994@

where E [ p ( v i ) ] = g = 1 G α g p g MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamyraiaaykW7caGGBbGaamiCaiaa ykW7caGGOaGaamODa8aadaWgaaWcbaWdbiaadMgaa8aabeaak8qaca GGPaGaaiyxaiabg2da9maaqadabaGaaGPaVlabeg7aH9aadaWgaaWc baWdbiaadEgaa8aabeaak8qacaWGWbWdamaaBaaaleaapeGaam4zaa WdaeqaaaWdbeaacaWGNbGaeyypa0JaaGymaaqaaiaadEeaa0Gaeyye Iuoaaaa@52E8@ and E [ λ ( v i ) ] = g = 1 G α g λ g MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamyraiaacUfacqaH7oaBcaaMc8Ua aiikaiaadAhapaWaaSbaaSqaa8qacaWGPbaapaqabaGcpeGaaiykai aac2facqGH9aqpdaaeWaqaaiaaykW7cqaHXoqypaWaaSbaaSqaa8qa caWGNbaapaqabaGcpeGaeq4UdW2damaaBaaaleaapeGaam4zaaWdae qaaaWdbeaacaWGNbGaeyypa0JaaGymaaqaaiaadEeaa0GaeyyeIuoa aaa@52DB@ with the finite mixture model. When m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamyBaaaa@3C44@ is random and such that

1 m i = 1 m n i | M p E [ n i | M ] , 1 m i = 1 m n i | U p E [ n i | U ] , ( 4.3 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaauaabeqacmaaaeaaqaaaaaaaaaWdbmaalaaapaqaa8qa caaIXaaapaqaa8qacaWGTbaaaiaaykW7daaeWbqaaiaaykW7caWGUb WdamaaBaaaleaapeGaamyAaiaaysW7caGG8bGaaGPaVlaad2eaa8aa beaaa8qabaGaamyAaiabg2da9iaaigdaaeaacaWGTbaaniabggHiLd aak8aabaWaaCbiaeaapeGaeyOKH4kal8aabeqaa8qacaWGWbaaaaGc paqaa8qacaWGfbGaaGPaVlaacUfacaWGUbWdamaaBaaaleaapeGaam yAaiaaysW7caGG8bGaaGPaVlaad2eaa8aabeaak8qacaGGDbGaaiil aaWdaeaapeWaaSaaa8aabaWdbiaaigdaa8aabaWdbiaad2gaaaGaaG PaVpaaqahabaGaaGPaVlaad6gapaWaaSbaaSqaa8qacaWGPbGaaGjb VlaacYhacaaMc8UaamyvaaWdaeqaaaWdbeaacaWGPbGaeyypa0JaaG ymaaqaaiaad2gaa0GaeyyeIuoaaOWdaeaadaWfGaqaa8qacqGHsgIR aSWdaeqabaWdbiaadchaaaaak8aabaWdbiaadweacaaMc8Uaai4wai aad6gapaWaaSbaaSqaa8qacaWGPbGaaGjbVlaacYhacaaMc8Uaamyv aaWdaeqaaOGaaiyxa8qacaGGSaaaa8aacaaMf8UaaGzbVlaaywW7ca aMf8UaaiikaiaaisdacaGGUaGaaG4maiaacMcaaaa@88DC@

as N , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaaabaaaaaaaaapeGaamOtaiabgkziUkabg6HiLkaacYca aaa@4033@ the error rates and the model parameters are related as follows

FNR p 1 E [ p ( v i ) ] , ( N 1 ) FPR p E [ λ ( v i ) ] . ( 4.4 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBamXvP5wqonvsaeHbmv3yPrwyGmuy SXwANjxyWHwEaebbnrfifHhDYfgasaacH8rrps0lbbf9q8qqaqpepe c8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXd ar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadi qaaqaaaOqaauaabiqacmaaaeaaqaaaaaaaaaWdbiaabAeacaqGobGa aeOuaaWdaeaadaWfGaqaa8qacqGHsgIRaSWdaeqabaWdbiaadchaaa aak8aabaWdbiaaigdacqGHsislcaWGfbGaaGPaVlaacUfacaWGWbGa aGPaVlaacIcacaWG2bWdamaaBaaaleaapeGaamyAaaWdaeqaaOWdbi aacMcacaGGDbGaaiilaaWdaeaapeGaaiikaiaad6eacqGHsislcaaI XaGaaiykaiaaysW7caqGgbGaaeiuaiaabkfaa8aabaWaaCbiaeaape GaeyOKH4kal8aabeqaa8qacaWGWbaaaaGcpaqaa8qacaWGfbGaaGPa VlaacUfacqaH7oaBcaaMc8UaaiikaiaadAhapaWaaSbaaSqaa8qaca WGPbaapaqabaGcpeGaaiykaiaac2facaGGUaaaa8aacaaMf8UaaGzb VlaaywW7caaMf8UaaiikaiaaisdacaGGUaGaaGinaiaacMcaaaa@6F16@


Date modified: