A short note on quantile and expectile estimation in unequal probability samples 3. Expectile estimation

An alternative to quantiles are expectiles. The expectile function M ( α ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamytamaabm aabaGaeqySdegacaGLOaGaayzkaaaaaa@3B06@  is thereby defined by replacing the L 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamitamaaBa aaleaacaaIXaaabeaaaaa@38C4@  loss in (2.1) by the L 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamitamaaBa aaleaacaaIYaaabeaaaaa@38C5@  loss leading to

M ( α ) = arg min m { i = 1 N w α ( Y i m ) ( Y i m ) 2 } . ( 3.1 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamytamaabm aabaGaeqySdegacaGLOaGaayzkaaGaaGypaiGacggacaGGYbGaai4z amaaxababaGaciyBaiaacMgacaGGUbaaleaacaWGTbaabeaakmaacm aabaWaaabCaeqaleaacaWGPbGaaGypaiaaigdaaeaacaWGobaaniab ggHiLdGccaaMc8Uaam4DamaaBaaaleaacqaHXoqyaeqaaOWaaeWaae aacaWGzbWaaSbaaSqaaiaadMgaaeqaaOGaeyOeI0IaamyBaaGaayjk aiaawMcaamaabmaabaGaamywamaaBaaaleaacaWGPbaabeaakiabgk HiTiaad2gaaiaawIcacaGLPaaadaahaaWcbeqaaiaaikdaaaaakiaa wUhacaGL9baacaaIUaGaaGzbVlaaywW7caaMf8UaaGzbVlaaywW7ca GGOaGaaG4maiaac6cacaaIXaGaaiykaaaa@6680@

Note that M ( α ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamytamaabm aabaGaeqySdegacaGLOaGaayzkaaaaaa@3B06@ is continuous in α MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaeqySdegaaa@38AB@ even for finite populations. Moreover M ( 0.5 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamytamaabm aabaGaaGimaiaai6cacaaI1aaacaGLOaGaayzkaaaaaa@3B98@ equals the mean value Y ¯ = i = 1 N Y i / N . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGabmywayaara GaaGypamaaqadabeWcbaGaamyAaiaai2dacaaIXaaabaGaamOtaaqd cqGHris5aOWaaSGbaeaacaWGzbWaaSbaaSqaaiaadMgaaeqaaaGcba GaamOtaaaacaGGUaaaaa@41B6@ Using the sample y 1 , , y n MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyEamaaBa aaleaacaaIXaaabeaakiaaiYcacqWIMaYscaaISaGaamyEamaaBaaa leaacaWGUbaabeaaaaa@3DA6@ with inclusion probabilities π 1 , , π n MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaeqiWda3aaS baaSqaaiaaigdaaeqaaOGaaGilaiablAciljaaiYcacqaHapaCdaWg aaWcbaGaamOBaaqabaaaaa@3F24@ we can estimate M ( α ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamytamaabm aabaGaeqySdegacaGLOaGaayzkaaaaaa@3B06@ by replacing the sum in (2.2) by its sample version, i.e.,

M ^ ( α ) = arg min m { j = 1 n 1 π j w α , j ( y j m ) 2 } MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGabmytayaaja WaaeWaaeaacqaHXoqyaiaawIcacaGLPaaacaaI9aGaciyyaiaackha caGGNbWaaCbeaeaaciGGTbGaaiyAaiaac6gaaSqaaiaad2gaaeqaaO WaaiWaaeaadaaeWbqabSqaaiaadQgacaaI9aGaaGymaaqaaiaad6ga a0GaeyyeIuoakmaalaaabaGaaGymaaqaaiabec8aWnaaBaaaleaaca WGQbaabeaaaaGccaWG3bWaaSbaaSqaaiabeg7aHjaaiYcacaWGQbaa beaakmaabmaabaGaamyEamaaBaaaleaacaWGQbaabeaakiabgkHiTi aad2gaaiaawIcacaGLPaaadaahaaWcbeqaaiaaikdaaaaakiaawUha caGL9baaaaa@592E@

with w α , j MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaam4DamaaBa aaleaacqaHXoqycaaISaGaamOAaaqabaaaaa@3B78@ as defined above. It is easy to see that the sum in M ^ ( α ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGabmytayaaja WaaeWaaeaacqaHXoqyaiaawIcacaGLPaaaaaa@3B16@ is a design-unbiased estimate for the sum in M ( α ) . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamytamaabm aabaGaeqySdegacaGLOaGaayzkaaGaaiOlaaaa@3BB8@ The estimate itself is however not design-unbiased like for the quantile function above. However the same arguments as for Q N ( α ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqWqpipeea0xe9Lq=Je9 vqaqFeFr0xbbG8FaYPYRWFb9fi0FXxbbf9Ff0dfrpm0dXdHqVu0=vr 0=vr0=fdbaqaaeGacaGaaiaabeqaamaabaabaaGcbaGaamyuamaaBa aaleaacaWGobaabeaakmaabmaabaGaeqySdegacaGLOaGaayzkaaaa aa@3C13@ in (2.2) may be used to establish design-consistency.

Date modified: