Mesure de l’incertitude associée aux estimateurs pour petits domaines basés sur un modèle
Section 3. Estimateurs de l’EQM basée sur le modèle

À la présente section, nous nous concentrons sur l’EQM basée sur le modèle des estimateurs EB sous les modèles de base au niveau du domaine et au niveau de l’unité. Il n’existe aucune expression analytique de l’EQM, sauf dans quelques cas particuliers. Ce problème a fait couler beaucoup d’encre dans la littérature traitant de l’estimation sur petits domaines (EPD), et a mené à des approximations d’ordre deux de l’EQM qui, à leur tour, sont utilisées pour obtenir des estimateurs sans biais d’ordre deux de l’EQM sous les modèles hypothétiques.

3.1  Modèle de base au niveau du domaine

Nous nous concentrons sur les estimateurs du REML des paramètres du modèle, notés β ^ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabCOSdyaaja aaaa@3747@ et σ ^ v 2 . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafq4WdmNbaK aadaqhaaWcbaGaamODaaqaaiaaikdaaaGccaGGUaaaaa@3A6C@ Un estimateur d’ordre deux sans biais inconditionnel de l’EQM sous le modèle de l’estimateur EB est donné par

eqm ( θ ^ i EB ) = g 1 i ( σ ^ v 2 ) + g 2 i ( σ ^ v 2 ) + 2 g 3 i ( σ ^ v 2 ) . ( 3.1 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeyzaiaabg hacaqGTbWaaeWaaeaacuaH4oqCgaqcamaaDaaaleaacaWGPbaabaGa aeyraiaabkeaaaaakiaawIcacaGLPaaacqGH9aqpcaWGNbWaaSbaaS qaaiaaigdacaWGPbaabeaakmaabmaabaGafq4WdmNbaKaadaqhaaWc baGaamODaaqaaiaaikdaaaaakiaawIcacaGLPaaacqGHRaWkcaWGNb WaaSbaaSqaaiaaikdacaWGPbaabeaakmaabmaabaGafq4WdmNbaKaa daqhaaWcbaGaamODaaqaaiaaikdaaaaakiaawIcacaGLPaaacqGHRa WkcaaIYaGaam4zamaaBaaaleaacaaIZaGaamyAaaqabaGcdaqadaqa aiqbeo8aZzaajaWaa0baaSqaaiaadAhaaeaacaaIYaaaaaGccaGLOa GaayzkaaGaaiOlaiaaywW7caaMf8UaaGzbVlaaywW7caaMf8Uaaiik aiaaiodacaGGUaGaaGymaiaacMcaaaa@6688@

Ici, le terme principal dans (3.1) est donné par (2.3) avec σ v 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aa0 baaSqaaiaadAhaaeaacaaIYaaaaaaa@39A0@ remplacé par σ ^ v 2 , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafq4WdmNbaK aadaqhaaWcbaGaamODaaqaaiaaikdaaaGccaGGSaaaaa@3A6A@ et les deux autres termes dans (3.1) sont d’ordre inférieur et rendent compte de l’estimation de β MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOSdaaa@3737@ et σ v 2 , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aa0 baaSqaaiaadAhaaeaacaaIYaaaaOGaaiilaaaa@3A5A@ respectivement (voir Rao et Molina, 2015, chapitre 6 pour des précisions). L’estimateur de l’EQM (3.1) est positif et sans biais à l’ordre deux en ce sens que son biais est d’ordre inférieur à 1 / m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSGbaeaaca aIXaaabaGaamyBaaaaaaa@37BC@ pour m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyBaaaa@36EB@ grand. Des méthodes bootstrap paramétriques ont également été utilisées pour obtenir un estimateur de l’EQM. Cependant, l’estimateur résultant n’est pas sans biais à l’ordre deux et un ajustement supplémentaire du biais est effectué pour s’assurer de l’absence de biais à l’ordre deux. Ces ajustements nécessitent habituellement des méthodes bootstrap doubles et certains des estimateurs bootstrap ajustés de l’EQM peuvent prendre des valeurs négatives; voir Rao et Molina (2015), chapitre 6.

3.2  Modèle de base au niveau de l’unité

Nous nous concentrons de nouveau sur l’estimation du REML des paramètres du modèle dans le modèle au niveau de l’unité (2.5). Un estimateur inconditionnel sans biais à l’ordre deux positif de l’EQM de l’estimateur EB μ ^ i EB MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafqiVd0MbaK aadaqhaaWcbaGaamyAaaqaaiaabweacaqGcbaaaaaa@3A67@ est donné par

eqm ( μ ^ i EB ) = g 1 i ( σ ^ v 2 , σ ^ e 2 ) + g 2 i ( σ ^ v 2 , σ ^ e 2 ) + 2 g 3 i ( σ ^ v 2 , σ ^ e 2 ) , ( 3.2 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeyzaiaabg hacaqGTbWaaeWaaeaacuaH8oqBgaqcamaaDaaaleaacaWGPbaabaGa aeyraiaabkeaaaaakiaawIcacaGLPaaacqGH9aqpcaWGNbWaaSbaaS qaaiaaigdacaWGPbaabeaakmaabmaabaGafq4WdmNbaKaadaqhaaWc baGaamODaaqaaiaaikdaaaGccaGGSaGaaGjbVlqbeo8aZzaajaWaa0 baaSqaaiaadwgaaeaacaaIYaaaaaGccaGLOaGaayzkaaGaey4kaSIa am4zamaaBaaaleaacaaIYaGaamyAaaqabaGcdaqadaqaaiqbeo8aZz aajaWaa0baaSqaaiaadAhaaeaacaaIYaaaaOGaaiilaiaaysW7cuaH dpWCgaqcamaaDaaaleaacaWGLbaabaGaaGOmaaaaaOGaayjkaiaawM caaiabgUcaRiaaikdacaWGNbWaaSbaaSqaaiaaiodacaWGPbaabeaa kmaabmaabaGafq4WdmNbaKaadaqhaaWcbaGaamODaaqaaiaaikdaaa GccaGGSaGaaGjbVlqbeo8aZzaajaWaa0baaSqaaiaadwgaaeaacaaI YaaaaaGccaGLOaGaayzkaaGaaiilaiaaywW7caaMf8UaaGzbVlaayw W7caaMf8UaaiikaiaaiodacaGGUaGaaGOmaiaacMcaaaa@784E@

où le premier terme est le terme principal donné à la section 2.2, le deuxième terme est dû à l’estimation de β MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOSdaaa@3737@ et le dernier terme est dû à l’estimation de σ v 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aa0 baaSqaaiaadAhaaeaacaaIYaaaaaaa@39A0@ et σ e 2 . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aa0 baaSqaaiaadwgaaeaacaaIYaaaaOGaaiOlaaaa@3A4B@ L’estimateur EB μ ^ i EB MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafqiVd0MbaK aadaqhaaWcbaGaamyAaaqaaiaabweacaqGcbaaaaaa@3A67@ et l’estimateur de l’EQM inconditionnelle associée (3.2) sont valides quand la fraction d’échantillonnage f i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpgpC0xc9LqFf0xc9 qqpeuf0xe9q8qiYRWFGCk9vq=dbbf9v8Gq0db9qqpm0dXdHqpq0=vr 0=vr0=edbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzamaaBa aaleaacaWGPbaabeaaaaa@37FE@ est négligeable. Le lecteur est invité à consulter (Rao et Molina, 2015, section 7.2.3) pour l’estimation de l’EQM dans le cas de fractions d’échantillonnage non négligeables.


Date de modification :