Articulo de referencia

Estimador M de dos pasos

Los estimadores M de dos pasos abordan problemas de estimación M que requieren una estimación preliminar para obtener el parámetro de interés. La estimación M de dos pasos difie...

Los estimadores M de dos pasos abordan problemas de estimación M que requieren una estimación preliminar para obtener el parámetro de interés. La estimación M de dos pasos difiere del problema de estimación M habitual porque la distribución asintótica del estimador de segundo paso generalmente depende del estimador de primer paso. Considerar este cambio en la distribución asintótica es importante para una inferencia válida.

Descripción

La clase de estimadores M de dos pasos incluye el estimador de selección de muestra de Heckman , [ 1 ] mínimos cuadrados no lineales ponderados y mínimos cuadrados ordinarios con regresores generados . [ 2 ]

Para fijar ideas, deja {Wi}i=1norteRd{\displaystyle \{W_{i}\}_{i=1}^{n}\subseteq R^{d}}ser una muestra i.i.d.Θ{\displaystyle \Theta }yΓ{\displaystyle \Gamma }son subconjuntos de espacios euclidianosRpag{\displaystyle R^{p}}yRq{\displaystyle R^{q}}, respectivamente. Dada una funciónmetro(;;;):Rd×Θ×ΓR{\displaystyle m(;;;):R^{d}\times \Theta \times \Gamma \rightarrow R}, estimador M de dos pasos θ^{\displaystyle {\sombrero {\theta }}} se define como:

θ^:=argmáximoθΘ1norteimetro(Wi,θ,γ^){\displaystyle {\hat {\theta }}:=\arg \max _{\theta \in \Theta }{\frac {1}{n}}\sum _{i}m{\bigl (}W_{i},\theta ,{\hat {\gamma }}{\bigr )}}

dóndeγ^{\displaystyle {\sombrero {\gamma }}}es una estimación M de un parámetro de perturbación que debe calcularse en el primer paso.

La consistencia de los estimadores M de dos pasos se puede verificar comprobando las condiciones de consistencia para los estimadores M habituales, aunque podría ser necesaria alguna modificación. En la práctica, la condición importante a comprobar es la condición de identificación . [ 2 ] Si γ^γ,{\displaystyle {\hat {\gamma }}\rightarrow \gamma ^{*},}dóndeγ{\displaystyle \gamma ^{*}}es un vector no aleatorio, entonces la condición de identificación es que mi[metro(W1,θ,γ)]{\displaystyle E[m(W_{1},\theta,\gamma ^{*})]}tiene un maximizador único sobreΘ{\displaystyle \Theta }.

Distribución asintótica

Bajo condiciones de regularidad, los estimadores M de dos pasos tienen normalidad asintótica . Un punto importante a tener en cuenta es que la varianza asintótica de un estimador M de dos pasos generalmente no es la misma que la del estimador M usual en el que la estimación del primer paso no es necesaria. [ 3 ] Este hecho es intuitivo porque γ^{\displaystyle {\sombrero {\gamma }}}es un objeto aleatorio y su variabilidad debería influir en la estimación deΘ{\displaystyle \Theta }Sin embargo, existe un caso especial en el que la varianza asintótica del estimador M de dos pasos toma la forma como si no existiera un procedimiento de estimación de primer paso. Este caso especial se produce si:

miθγmetro(W1,θ0,γ)=0{\displaystyle E{\frac {\partial }{\partial \theta \partial \gamma }}m(W_{1},\theta _{0},\gamma ^{*})=0}

dóndeθ0{\displaystyle \theta _{0}}es el verdadero valor deθ{\displaystyle \theta }yγ{\displaystyle \gamma ^{*}} es el límite de probabilidad deγ^{\displaystyle {\sombrero {\gamma }}}. [ 3 ] Para interpretar esta condición, primero observe que bajo condiciones de regularidad, miθmetro(W1,θ0,γ)=0{\displaystyle E{\frac {\partial }{\partial \theta }}m(W_{1},\theta _{0},\gamma ^{*})=0}desdeθ0{\displaystyle \theta _{0}}es el maximizador de mi[metro(W1,θ,γ)]{\displaystyle E[m(W_{1},\theta,\gamma ^{*})]}. Por lo tanto, la condición anterior implica que una pequeña perturbación en γ no tiene impacto en la condición de primer orden . Por lo tanto, en muestras grandes, la variabilidad de γ^{\displaystyle {\sombrero {\gamma }}} Esto no afecta al argmax de la función objetivo, lo que explica la propiedad invariante de la varianza asintótica. Por supuesto, este resultado solo es válido cuando el tamaño de la muestra tiende a infinito, por lo que la propiedad para muestras finitas podría ser bastante diferente.

Con la participación de MLE

Cuando el primer paso es un estimador de máxima verosimilitud , bajo ciertas suposiciones, el estimador M de dos pasos es más asintóticamente eficiente (es decir, tiene una varianza asintótica menor) que el estimador M con un parámetro del primer paso conocido. La consistencia y la normalidad asintótica del estimador se derivan del resultado general sobre los estimadores M de dos pasos. [ 4 ]

Sea {V i ,W i ,Z i } n i=1 una muestra aleatoria y el estimador M de segundo pasoθ^{\displaystyle {\widehat {\theta }}}es lo siguiente:

θ^:=argramomáximoθΘi=1nortemetro(vi,wi,zi:θ,γ^){\displaystyle {\widehat {\theta }}:={\underset {\theta \in \Theta }{\operatorname {arg\max } }}\sum _{i=1}^{n}m(v_{i},w_{i},z_{i}:\theta ,{\widehat {\gamma }})}

dónde γ^{\displaystyle {\widehat {\gamma }}}es el parámetro estimado por máxima verosimilitud en el primer paso. Para el MLE,

γ^:=argramomáximoγΓi=1norteregistroF(vit:zi,γ){\displaystyle {\widehat {\gamma }}:={\underset {\gamma \in \Gamma }{\operatorname {arg\max } }}\sum _{i=1}^{n}\log f(v_{it}:z_{i},\gamma )}

donde f es la densidad condicional de V dado Z. Ahora, supongamos que dado Z , V es condicionalmente independiente de W. Esto se denomina supuesto de independencia condicional o selección en observables. [ 4 ] [ 5 ] Intuitivamente, esta condición significa que Z es un buen predictor de V, de modo que, una vez condicionado a Z, V no tiene dependencia sistemática de W. Bajo el supuesto de independencia condicional, la varianza asintótica del estimador de dos pasos es:

mi[θs(θ0,γ0)]1mi[gramo(θ0,γ0)gramo(θ0,γ0)T]mi[θs(θ0,γ0)]1{\displaystyle \mathrm {E} [\nabla _ {\theta }s(\theta _ {0},\gamma _ {0})]^{-1}\mathrm {E} [g(\theta _ {0},\gamma _ {0})g(\theta _ {0},\gamma _ {0})^{\mathrm {T} }]\mathrm {E} [\nabla _{\theta }s(\theta _{0},\gamma _{0})]^{-1}}

dónde

gramo(θ,γ):=s(θ,γ)mi[γs(θ,γ)]mi[d(γ)d(γ)T]1d(γ)s(θ,γ):=θmetro(V,W,Z:θ,γ)d(γ):=γregistroF(V:Z,γ){\displaystyle {\begin{aligned}g(\theta ,\gamma )&:=s(\theta ,\gamma )-\mathrm {E} [\nabla _ {\gamma }s(\theta ,\gamma )]\mathrm {E} [d(\gamma )d(\gamma )^{\mathrm {T} }]^{-1}d(\gamma )\\s(\theta ,\gamma )&:=\nabla _{\theta }m(V,W,Z:\theta ,\gamma )\\d(\gamma )&:=\nabla _{\gamma }\log f(V:Z,\gamma )\end{alineado}}}

y representa la derivada parcial con respecto a un vector fila . En el caso en que γ 0 es conocido, la varianza asintótica es

mi[θs(θ0,γ0)]1mi[s(θ0,γ0)s(θ0,γ0)T]mi[θs(θ0,γ0)]1{\displaystyle \mathrm {E} [\nabla _ {\theta }s(\theta _ {0},\gamma _ {0})]^{-1}\mathrm {E} [s(\theta _ {0},\gamma _ {0})s(\theta _ {0},\gamma _ {0})^{\mathrm {T} }]\mathrm {E} [\nabla _{\theta }s(\theta _{0},\gamma _{0})]^{-1}}

y por lo tanto, a menos quemi[γs(θ,γ)]=0{\displaystyle \mathrm {E} [\nabla _ {\gamma }s(\theta,\gamma )]=0}El estimador M de dos pasos es más eficiente que el estimador M habitual. Esto sugiere que, incluso cuando γ₀ se conoce a priori, se obtiene una ganancia de eficiencia al estimar γ mediante MLE. Una aplicación de este resultado se encuentra, por ejemplo, en la estimación del efecto del tratamiento. [ 4 ]

Ejemplos

Véase también

Referencias

  1. Heckman, JJ, La estructura común de los modelos estadísticos de truncamiento, selección de muestras y variables dependientes limitadas y un estimador simple para dichos modelos, Annals of Economic and Social Measurement, 5,475-492.
  2. 1 2 Wooldridge, JM, Análisis econométrico de datos de sección transversal y de panel, MIT Press, Cambridge, Mass.
  3. 1 2 Newey, KW y D. McFadden, Estimación de muestras grandes y prueba de hipótesis, en R. Engel y D. McFadden, eds., Manual de econometría, vol. 4, Ámsterdam: North-Holland.
  4. 1 2 3 Wooldridge, JM, Análisis econométrico de datos de sección transversal y de panel, MIT Press, Cambridge, Mass.
  5. Heckman, JJ y R. Robb, 1985, Métodos alternativos para evaluar el impacto de las intervenciones: una visión general, Journal of Econometrics, 30, 239-267.