Dado un conjunto de eventos {A, B, C, D, E}
que se producen una vez cada mes para n
años:
[A, B, C, C, B, D, A, B, C, C, B, D]
[E, B, C, B, B, D, E, B, C, B, B, D]
[C, B, C, D, E, A, A, D, C, C, B, D] //12 months x 3 years pictured
...
Tengo la función de probabilidad p(w, r) = ...
que calcula la probabilidad de un evento a seguir a otro (por ejemplo, Una después de B), donde w
y r
son parámetros que el modelo de salida de tal manera que cuando correctamente recogidos deben ajustarse a los datos de la muestra.
Mi objetivo final es el de predecir los acontecimientos en un año, pero estoy atascado en la utilización de la EM para determinar w
y r
.
Intuitivamente, lo que hago ahora es:
- Empezar con un total valor aleatorio para
w
yr
- Calcular la probabilidad de cada par en un año (por ejemplo,
A
yB
) el uso dep
función, compararla con la actual distribución de probabilidad extraído de las muestras y, básicamente, obtener la probabilidad (v
) es para la actualw
yr
a de ser para adaptarse a las muestras. - Ahora mi 3er paso sería modificar
w
yr
, de modo quev
deberían converger hacia 1. Esto es donde estoy atascado.
¿Cómo debo usar v
para obtener los nuevos valores de w
y r
, de modo que v
eventualmente convergen hacia 1?
EDITAR:
Me gustaría añadir que p
, básicamente, me da la estimación de la distribución de probabilidad para el caso de la sucesión. Lo que significa que tengo dos superficies 2D: la estimación y las muestras de las distribuciones de probabilidad y quiero usar w
y r
a cambio de escala y la superficie estimada de tal forma que mejor se ajusta a la muestra de uno.
Así que mi problema es, ¿cómo puedo comparar los dos después de un paso y cómo debo conseguir nuevos w
y r
valores para el siguiente paso.
E. g.: Estaba pensando que yo podría usar algún tipo de matriz de la norma, para obtener el grado de similitud entre los dos y, a continuación, utilizar esto para decidir, basado en iteraciones anteriores si debo aumentar w
/r
o disminuir.