Quiero maximizar la función de verosimilitud logarítmica (L) que es una función de los parámetros de $\beta_i$$i=1,..,k$$\alpha_1, \alpha_2$. Idealmente, quiero realizar la estimación de todos los parámetros en un solo paso. Por desgracia, no puedo realizar un paso de la estimación debido a la forma en que el modelo que tengo. Pero para algunos revisión de los valores de $\alpha_1, \alpha_2$, no puedo encontrar el de máxima verosimilitud estimaciones de $\beta_i$$i=1,..,k$. Así que lo que hice es que he creado una función como $f(\alpha_1, \alpha_2)$ lo define como $f(\alpha_1, \alpha_2)=L(\alpha_1, \alpha_2,\hat{\beta})$ es decir, el logaritmo de la probabilidad correspondiente a $\alpha_1, \alpha_2$. Entonces yo maximiza esta función $f$ numéricamente con respecto a $\alpha_1, \alpha_2$.
¿Este enfoque resolver la inconsistencia de las dos etapas de la estimación que tengo? Es este un enfoque válido en todo? Si no, ¿hay algún otro método de estimación que puedo usar?