5. Adjusted maximum likelihood

Isabel Molina, J.N.K. Rao and Gauri Sankar Datta

Previous | Next

The estimation methods for A MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiaadgeaaa a@39C3@  described in Section 2 might produce zero estimates. In this case, the EBLUPs will give zero weight to the direct estimators in all areas, regardless of the efficiency of the direct estimator in each area. On the other hand, survey sampling practitioners often prefer to give always a strictly positive weight to direct estimators because they are based on the area-specific unit level data for the variable of interest without the assumption of any regression model. For this situation, Li and Lahiri (2010) proposed the AML estimator that delivers a strictly positive estimator of A . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiaadgeaca GGUaaaaa@3A75@  This estimator, denoted here A ^ AML , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiqadgeaga qcamaaBaaaleaacaqGbbGaaeytaiaabYeaaeqaaOGaaiilaaaa@3D1C@  is obtained by maximizing the adjusted likelihood defined as

L AML ( A ) = A × L P ( A ) . MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiaadYeada WgaaWcbaGaaeyqaiaab2eacaqGmbaabeaakmaabmqabaGaamyqaaGa ayjkaiaawMcaaiabg2da9iaadgeacqGHxdaTcaWGmbWaaSbaaSqaai aadcfaaeqaaOWaaeWabeaacaWGbbaacaGLOaGaayzkaaGaaGOlaaaa @477D@

The EBLUP given in (2.6) with A ^ = A ^ AML MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiqadgeaga qcaiabg2da9iqadgeagaqcamaaBaaaleaacaqGbbGaaeytaiaabYea aeqaaaaa@3E3E@  will be denoted hereafter as θ ^ AML = ( θ ^ AML ,1 , , θ ^ AML , m ) . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiqahI7aga qcamaaBaaaleaacaqGbbGaaeytaiaabYeaaeqaaOGaeyypa0ZaaeWa beaacuaH4oqCgaqcamaaBaaaleaacaqGbbGaaeytaiaabYeacaaISa GaaGymaaqabaGccaaISaGaeSOjGSKaaGilaiqbeI7aXzaajaWaaSba aSqaaiaabgeacaqGnbGaaeitaiaaiYcacaWGTbaabeaaaOGaayjkai aawMcaamaaCaaaleqabaGccWaGGBOmGikaaiaac6caaaa@51B0@  Note that θ ^ AML MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiqahI7aga qcamaaBaaaleaacaqGbbGaaeytaiaabYeaaeqaaaaa@3CE0@  assigns strictly positive weights to direct estimators.

Li and Lahiri (2010) proposed a second order unbiased MSE estimator of θ ^ AML , i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiqbeI7aXz aajaWaaSbaaSqaaiaabgeacaqGnbGaaeitaiaaiYcacaWGPbaabeaa aaa@3EF6@  given by

mse ( θ ^ AML , i ) = g 1 i ( A ^ AML ) + g 2 i ( A ^ AML ) + 2 g 3 i ( A ^ AML ) B i 2 ( A ^ AML ) b AML ( A ^ AML ) , ( 5.1 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrVipeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaauaabaqacm aaaeaacaqGTbGaae4CaiaabwgadaqadeqaaiqbeI7aXzaajaWaaSba aSqaaiaabgeacaqGnbGaaeitaiaaiYcacaWGPbaabeaaaOGaayjkai aawMcaaaqaaiabg2da9aqaaiaadEgadaWgaaWcbaGaaGymaiaadMga aeqaaOWaaeWabeaaceWGbbGbaKaadaWgaaWcbaGaaeyqaiaab2eaca qGmbaabeaaaOGaayjkaiaawMcaaiabgUcaRiaadEgadaWgaaWcbaGa aGOmaiaadMgaaeqaaOWaaeWabeaaceWGbbGbaKaadaWgaaWcbaGaae yqaiaab2eacaqGmbaabeaaaOGaayjkaiaawMcaaiabgUcaRiaaikda caWGNbWaaSbaaSqaaiaaiodacaWGPbaabeaakmaabmqabaGabmyqay aajaWaaSbaaSqaaiaabgeacaqGnbGaaeitaaqabaaakiaawIcacaGL PaaaaeaaaeaacqGHsislaeaacaWGcbWaa0baaSqaaiaadMgaaeaaca aIYaaaaOWaaeWabeaaceWGbbGbaKaadaWgaaWcbaGaaeyqaiaab2ea caqGmbaabeaaaOGaayjkaiaawMcaaiaadkgadaWgaaWcbaGaaeyqai aab2eacaqGmbaabeaakmaabmqabaGabmyqayaajaWaaSbaaSqaaiaa bgeacaqGnbGaaeitaaqabaaakiaawIcacaGLPaaacaaISaaaaiaayw W7caaMf8UaaGzbVlaaywW7caaMf8UaaiikaiaaiwdacaGGUaGaaGym aiaacMcaaaa@7B53@

where b AML ( A ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiaadkgada WgaaWcbaGaaeyqaiaab2eacaqGmbaabeaakmaabmqabaGaamyqaaGa ayjkaiaawMcaaaaa@3ECD@  is the bias of A ^ AML MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiqadgeaga qcamaaBaaaleaacaqGbbGaaeytaiaabYeaaeqaaaaa@3C62@  and it is given by

b AML ( A ) = trace { P ( A ) Σ 1 ( A ) } + 2 / A trace { Σ 2 ( A ) } . MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiaadkgada WgaaWcbaGaaeyqaiaab2eacaqGmbaabeaakmaabmqabaGaamyqaaGa ayjkaiaawMcaaiabg2da9maalaaabaGaaeiDaiaabkhacaqGHbGaae 4yaiaabwgadaGadeqaaiaahcfadaqadeqaaiaadgeaaiaawIcacaGL PaaacqGHsislcaWHJoWaaWbaaSqabeaacqGHsislcaaIXaaaaOWaae WabeaacaWGbbaacaGLOaGaayzkaaaacaGL7bGaayzFaaGaey4kaSYa aSGbaeaacaaIYaaabaGaamyqaaaaaeaacaqG0bGaaeOCaiaabggaca qGJbGaaeyzamaacmqabaGaaC4OdmaaCaaaleqabaGaeyOeI0IaaGOm aaaakmaabmqabaGaamyqaaGaayjkaiaawMcaaaGaay5Eaiaaw2haaa aacaaIUaaaaa@5F87@

Previous | Next

Date modified: