A generalization of inverse probability weighting
Section 5. The generalized Godambe-Joshi lower bound

For any unbiased estimator θ ^ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbwacuaH4o qCgaqcaaaa@3AFD@  of the population total θ , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbiqaaWsgjugybi abeI7aXjaacYcaaaa@3BE5@  if V p ( θ ^ ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbwacaWGwb GcdaWgaaWcbaGaamiCaaqabaGcdaqadeqaaKqzGfGaaGjcVlqbeI7a XzaajaGaaGjcVdGccaGLOaGaayzkaaaaaa@4292@  is the variance of θ ^ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbwacuaH4o qCgaqcaaaa@3AFD@  under the sampling plan, Godambe and Joshi (1965) have given a lower bound for the value of E ξ V p ( θ ^ ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbwacaWGfb GcdaWgaaWcbaGaeqOVdGhabeaajugybiaadAfakmaaBaaaleaacaWG WbaabeaakmaabmqabaqcLbwacaaMi8UafqiUdeNbaKaacaaMi8oaki aawIcacaGLPaaaaaa@4624@  under the assumption that the variance matrix V ξ ( y ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbiqaaqxajugybi aadAfakmaaBaaaleaacqaH+oaEaeqaaOWaaeWaaeaajugybiaahMha aOGaayjkaiaawMcaaaaa@3FF1@  was diagonal. That lower bound is the sum of the elements of the diagonal matrix ( ( E ( Δ s ) ) 1 I ) V ξ ( y ) . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbiqaaqxadaqada qaamaabmaabaqcLbwacaWGfbGcdaqadaqaceaa0fqcLbwacaWHuoGc daWgaaWcbaGaam4CaaqabaaakiaawIcacaGLPaaaaiaawIcacaGLPa aadaahaaWcbeqaaiabgkHiTiaaigdaaaGccaaMe8EcLbwacqGHsisl caaMe8UaaCysaaGccaGLOaGaayzkaaGaaGjbVNqzGfGaamOvaOWaaS baaSqaaiabe67a4bqabaGcdaqadaqaaKqzGfGaaCyEaaGccaGLOaGa ayzkaaGaaiOlaaaa@539E@  That result is generalized in the following paragraph.

For any linear unbiased total estimator, θ ^ , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbwacuaH4o qCgaqcaiaacYcaaaa@3BAD@ if V ξ ( y ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbiqaaqxajugybi aadAfakmaaBaaaleaacqaH+oaEaeqaaOWaaeWaaeaajugybiaahMha aOGaayjkaiaawMcaaaaa@3FF1@ is positive definite, then E ξ V p ( θ ^ ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbwacaWGfb GcdaWgaaWcbaGaeqOVdGhabeaajugybiaadAfakmaaBaaaleaacaWG WbaabeaakmaabmqabaqcLbwacaaMi8UafqiUdeNbaKaacaaMi8oaki aawIcacaGLPaaaaaa@4624@ is not lower than the sum of the elements of the matrix ( E ( Δ s V ξ ( y ) Δ s ) ) 1 V ξ ( y ) . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbiqaaqxadaqada qaaKqzGfGaamyraiaaykW7kmaabmaabiqaaqxajugybiaahs5akmaa BaaaleaacaWGZbaabeaajugybiaadAfakmaaBaaaleaacqaH+oaEae qaaOWaaeWaaeaajugybiaahMhaaOGaayjkaiaawMcaaiaaysW7jugy biaahs5akmaaBaaaleaacaWGZbaabeaaaOGaayjkaiaawMcaamaaCa aaleqabaGaaiiiGaaaaOGaayjkaiaawMcaamaaCaaaleqabaGaeyOe I0IaaGymaaaakiaaysW7jugybiabgkHiTiaaysW7caWGwbGcdaWgaa WcbaGaeqOVdGhabeaakmaabmaabaqcLbwacaWH5baakiaawIcacaGL PaaacaGGUaaaaa@5D1E@ It is easily verified that the usual Godambe-Joshi lower bound is obtained if V ξ ( y ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbiqaaqxajugybi aadAfakmaaBaaaleaacqaH+oaEaeqaaOWaaeWaaeaajugybiaahMha aOGaayjkaiaawMcaaaaa@3FF1@ is diagonal.

Just as the calibration estimator asymptotically attains the Godambe-Joshi lower bound, the generalized calibration estimator with Σ = V ξ ( y ) , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbwacaWHJo GaaGjbVlabg2da9iaaysW7caWGwbGcdaWgaaWcbaGaeqOVdGhabeaa kmaabmaabaqcLbwacaWH5baakiaawIcacaGLPaaacaGGSaaaaa@4578@ asymptotically attains the generalized Godambe-Joshi lower bound, regardless of the value of the matrices X , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCiwaiaacY caaaa@39F9@ T MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCivaaaa@3945@ and U . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCyvaiaac6 caaaa@39F8@ The link between the value of those three matrices and the value of V ξ ( y ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbiqaaqxajugybi aadAfakmaaBaaaleaacqaH+oaEaeqaaOWaaeWaaeaajugybiaahMha aOGaayjkaiaawMcaaaaa@3FF1@ is not examined in this paper, but the calibration problem stated in Section 3 does clarify the role of each of those matrices. The derivation of the generalized lower bound and the proof of the optimality of the generalized calibration estimator are given in Théberge (2017).

The fact that θ ^ GCAL ( V ξ ( y ) ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbwacuaH4o qCgaqcaOWaaSbaaSqaaiaabEeacaqGdbGaaeyqaiaabYeaaeqaaOWa aeWaaeaajugybiaadAfakmaaBaaaleaacqaH+oaEaeqaaOWaaeWaae aajugybiaahMhaaOGaayjkaiaawMcaaaGaayjkaiaawMcaaaaa@46FA@ asymptotically attains the generalized Godambe-Joshi lower bound shows that the generalized inverse probability estimator performs well when applied to residuals, as it does in (3.1), even though it is not recommended in general. Similarly, the ordinary inverse probability estimator can be inefficient if the sample size is random, but will perform well if applied to residuals.

It should be noted that, contrary to the ordinary Godambe-Joshi lower bound, the generalized lower bound applies only to linear unbiased estimators. In fact, an example with V ξ ( y ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpC0xe9LqFHe9Lq pepeea0xd9q8as0=LqLs=Jirpepeea0=as0Fb9pgea0lrP0xe9Fve9 Fve9qapdbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbwacaWGwb GcdaWgaaWcbaGaeqOVdGhabeaakmaabmaabaqcLbwacaWH5baakiaa wIcacaGLPaaaaaa@3F79@ not diagonal, of a non-linear unbiased estimator which does better than the lower bound is given in Théberge (2017).


Date modified: