Este no es un problema de muestras finitas. Está relacionada con la complejidad de la función de error y la función de enlace. Como el tiempo es cuadrática no habrá múltiples maxima. Si no lo más probable es que a más de uno. Por ejemplo, si usted está bajo la hipótesis de que los datos tienen una distribución normal, entonces usted tiene:
$N(f(y|x;\theta),\sigma)$
donde f es la función de enlace. Ahora si se considera el pdf:
$f(y) = \frac{1}{(2\pi)^{1/2}\sigma} e^{-\frac{(y-f(x;\theta))^2}{2\sigma^2}}$
La combinación de la función, la función de enlace enchufado en la distribución, en este escenario es cuadrática (no te preocupes por el exponencial a medida que desaparece cuando se trabaja con su logaritmo) si $f(y|x;\theta)$ es una función lineal. Esto sucede en la regresión lineal con errores de Gauss, por ejemplo.
En cualquier caso, si su distribución no tiene una forma cuadrática o cuando se combina con la función de enlace termina no cuadrática tendrá el problema de las múltiples maxima.
No hay una buena solución para el problema como la búsqueda del máximo global es NP-duro.
http://web.maths.unsw.edu.au/~rsw/lgopt.pdf
páginas 20-21
Hay algoritmos que intentan cubrir todo el espacio de parámetros, pero son muy lentos como uno espera. Así que no es un problema de la MLE. Si se conecta la función de Probabilidad para una más sofisticada optimizador como el recocido simulado, búsqueda de la solución global.
http://en.wikipedia.org/wiki/Simulated_annealing