90 likes | 214 Views
Modelos de máxima verosimilitud para secuencias de ADN. Eliseo Martínez, Catherine Lizama. denota el valor de la base en un determinado sitio del ADN. C. A. C. A. T. G. T. A. C. A. C. C. T. G. T. A. probabilidad de que una base mute a la base i.
E N D
Modelos de máxima verosimilitud para secuencias de ADN Eliseo Martínez, Catherine Lizama
denota el valor de la base en un determinado sitio del ADN C A C A T G T A C A C C T G T A probabilidad de que una base mute a la base i Probabilidad de que la base i cambie a la base j luego de transcurrido el tiempo T (varias generaciones) Hipótesis fundamental: la tasa total de mutación, m, es constante, esto es la razón de sustitución de cambio de base de generación a generación es constante.
C A C A T G T A C A C A T G T A La probabilidad de que no ocurra mutación en ningún sitio hasta la generación T es … y en consecuencia la probabilidad de que ocurra (al menos una) mutación hasta la generación T es
C A C i T G T A C A C j T G T A La probabilidad de mutación después del tiempo T … y que el cambio haya sido por j (desde i) es entonces …
C A C i T G T A C A C i T G T A Ahora calcularemos la probabilidad de que el sitio, estando en la base i, al cabo de T generaciones siga con la base i. Esto es… … la probabilidad de no mutación después del tiempo T… … o que haya mutado al mismo estado (estando inicialmente en i) es decir…
Si ocurre que La probabilidad de cambio de base después del tiempo T está dada por
Ojo: recuerde que El modelo anterior considera el tiempo discreto (de generación en generación, digamos en intervalos de 20 años, por ejemplo). Un tratamiento más general considera el tiempo continuo, y entonces se usan probabilidades “instantáneas” de transición de una base a otra. Probabilidad de cambio de una base a otra (de i a j) en un tiempo de duración Dt … y en consecuencia la probabilidad de no cambiar en un tiempo de longitud Dt es
En este poco tiempo debe mutar desde k hasta j ¡cuidado, hay 16 ecuaciones! … y una de ellas es esta…
Afortunadamente estas 16 ecuaciones se pueden escribir en forma matricial (con matrices de 4 x 4) de la siguiente forma: