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.

Report a problem on this page

Is something not working? Is there information outdated? Can't find what you're looking for?

Please contact us and let us know how we can help you.

Privacy notice

Date modified: