Articulo de referencia

Algoritmo de puntuación

El algoritmo de puntuación , también conocido como puntuación de Fisher , [ 1 ] es una forma del método de Newton utilizado en estadística para resolver numéricamente ecuaciones...

El algoritmo de puntuación , también conocido como puntuación de Fisher , [ 1 ] es una forma del método de Newton utilizado en estadística para resolver numéricamente ecuaciones de máxima verosimilitud , llamado así en honor a Ronald Fisher .

Esquema de derivación

DejarY1,,Ynorte{\displaystyle Y_{1},\ldots ,Y_{n}}sean variables aleatorias , independientes e idénticamente distribuidas con función de densidad de probabilidad dos veces diferenciableF(y;θ){\displaystyle f(y;\theta )}y deseamos calcular el estimador de máxima verosimilitud (EMV).θ{\displaystyle \theta ^{*}}deθ{\displaystyle \theta }Primero, supongamos que tenemos un punto de partida para nuestro algoritmo.θ0{\displaystyle \theta _{0}}y consideremos una expansión de Taylor de la función de puntuación ,V(θ){\displaystyle V(\theta )}, acerca deθ0{\displaystyle \theta _{0}}:

V(θ)V(θ0)J(θ0)(θθ0),{\displaystyle V(\theta )\approx V(\theta _{0})-{\mathcal {J}}(\theta _{0})(\theta -\theta _{0}),\,}

dónde

J(θ0)=i=1norte|θ=θ0registroF(Yi;θ){\displaystyle {\mathcal {J}}(\theta _{0})=-\sum _{i=1}^{n}\left.\nabla \nabla ^{\top }\right|_{\theta =\theta _{0}}\log f(Y_{i};\theta )}

es la matriz de información observada enθ0{\displaystyle \theta _{0}}Ahora, configurandoθ=θ{\displaystyle \theta =\theta ^{*}}, usando esoV(θ)=0{\displaystyle V(\theta ^{*})=0}y reordenando obtenemos:

θθ0+J1(θ0)V(θ0).{\displaystyle \theta ^{*}\approx \theta _{0}+{\mathcal {J}}^{-1}(\theta _{0})V(\theta _{0}).\,}

Por lo tanto, utilizamos el algoritmo.

θmetro+1=θmetro+J1(θmetro)V(θmetro),{\displaystyle \theta _{m+1}=\theta _{m}+{\mathcal {J}}^{-1}(\theta _{m})V(\theta _{m}),\,}

y bajo ciertas condiciones de regularidad, se puede demostrar queθmetroθ{\displaystyle \theta _{m}\rightarrow \theta ^{*}}.

Puntuación de Fisher

En la práctica,J(θ){\displaystyle {\mathcal {J}}(\theta )}suele ser reemplazado porI(θ)=mi[J(θ)]{\displaystyle {\mathcal {I}}(\theta )=\mathrm {E} [{\mathcal {J}}(\theta )]}, la información de Fisher , lo que nos da el algoritmo de puntuación de Fisher :

θmetro+1=θmetro+I1(θmetro)V(θmetro){\displaystyle \theta _{m+1}=\theta _{m}+{\mathcal {I}}^{-1}(\theta _{m})V(\theta _{m})}...

Bajo ciertas condiciones de regularidad, siθmetro{\displaystyle \theta _{m}}es un estimador consistente , entoncesθmetro+1{\displaystyle \theta _{m+1}}(la corrección después de un solo paso) es "óptima" en el sentido de que su distribución de error es asintóticamente idéntica a la de la verdadera estimación de máxima verosimilitud. [ 2 ]

Véase también

Referencias

  1. Longford, Nicholas T. (1987). "Un algoritmo de puntuación rápido para la estimación de máxima verosimilitud en modelos mixtos desequilibrados con efectos aleatorios anidados". Biometrika . 74 (4): 817– 827. doi : 10.1093/biomet/74.4.817 .
  2. Li, Bing; Babu, G. Jogesh (2019), "Inferencia bayesiana" , Springer Texts in Statistics , Nueva York, NY: Springer New York, Teorema 9.4, doi : 10.1007/978-1-4939-9761-9_6 , ISBN 978-1-4939-9759-6, S2CID 239322258 , consultado el 3 de enero de 2023 

Lecturas adicionales

  • Jennrich, RI y Sampson, PF (1976). "Newton-Raphson y algoritmos relacionados para la estimación de componentes de varianza de máxima verosimilitud" . Technometrics . 18 (1): 11– 17. doi : 10.1080/00401706.1976.10489395 (inactivo el 12 de julio de 2025). JSTOR 1267911 . {{cite journal}}: CS1 maint: DOI inactivo desde julio de 2025 ( enlace )