한국   대만   중국   일본 
Processo de Moran ? Wikipedia, a enciclopedia livre Saltar para o conteudo

Processo de Moran

Origem: Wikipedia, a enciclopedia livre.

Um processo de Moran ou modelo de Moran e um processo estocastico simples usado em biologia para descrever populacoes finitas . O processo recebe este nome em homenagem a Patrick Moran, que propos o modelo pela primeira vez em 1958. [ 1 ] Pode ser usado para modelar processos que aumentam a variedade, tais como mutacoes , assim como efeitos que diminuem a variedade, tais como a deriva genetica e a selecao natural . O processo pode descrever a dinamica probabilistica em uma populacao finita com tamanho constante na qual dois alelos A e B estao competindo por dominancia. Os dois alelos sao considerados replicadores reais (isto e, entidades que fazem copias de si mesmas).

Em cada momento, um individuo aleatorio (que pode ser tanto do tipo A, como do tipo B) e escolhido para reproducao e um individuo aleatorio e escolhido para morte, garantindo assim que o tamanho da populacao permaneca constante. Para a selecao do modelo, um tipo deve ter uma aptidao maior, tendo por isso maior probabilidade de ser escolhido para reproducao. O mesmo individuo pode ser escolhido para morte e para reproducao no mesmo momento.

Deriva neutra [ editar | editar codigo-fonte ]

A deriva neutra e a ideia de que uma mutacao neutra pode se espalhar por uma populacao, de modo que eventualmente o alelo original seja perdido. Uma mutacao neutra nao traz qualquer vantagem ou desvantagem em aptidao ao portador. O caso simples do processo de Moran pode descrever este fenomeno. [ 2 ]

Se o numero de individuos A for dado por , entao, o processo de Moran e definido no espaco de estados . Ja que o numero de individuos A pode mudar no maximo em uma unidade a cada momento, a transicao existe somente entre o estado e os estados e . Assim, a matriz de transicao do processo estocastico e tridiagonal em forma e as probabilidades de transicao sao

A entrada denota a probabilidade de se passar do estado para o estado . Para entender as formulas para as probabilidades de transicao, e necessario observar a definicao do processo que afirma que sempre um individuo sera escolhido para reproducao e um individuo sera escolhido para a morte. Uma vez que os individuos A tenham morrido, eles nunca mais serao reintroduzidos na populacao, ja que o processo nao modela mutacoes (A nao pode ser reintroduzido na populacao uma vez que tenha morrido e vice-versa), logo, . Pela mesma razao, a populacao de individuos A permanecera sempre , uma vez que tenha alcancado este numero e tomado conta de populacao, logo, . Os estados e sao chamados de absorventes, enquanto os estados sao chamados de transientes. As probabilidades intermediarias de transicao podem ser explicadas considerando o primeiro termo como a probabilidade de escolha do individuo cuja abundancia aumentara em uma unidade e o segundo termo a probabilidade de escolha do outro tipo para a morte. Obviamente, se o mesmo tipo e escolhido para a reproducao e para a morte, entao, a abundancia deste tipo nao muda.

Eventualmente a populacao atingira um dos estados absorventes e nele permanecera para sempre. Nos estados transientes, flutuacoes aleatorias ocorrerao, mas eventualmente a populacao sera extinta ou atingira fixacao. Esta e uma das mais importantes diferencas em relacao a processos deterministicos que nao podem modelar eventos aleatorios. O valor esperado e a variancia do numero de individuos A no momento podem ser computados quando um estado inicial e dado:

Para uma derivacao matematica da equacao acima, clique em "expandir" para que ela seja mostrada

Para o valor esperado, o calculo ocorre como segue abaixo. Escrevendo p = i N ,

Escrevendo e e aplicando a lei da esperanca total, Aplicando o argumento repetidamente, temos ou

Para a variancia, o calculo ocorre como segue abaixo. Escrevendo temos

Para todo t , e sao identicamente distribuidos, logo, suas variancias sao iguais. Escrevendo como antes e e aplicando a lei da variancia total,

Se , obtemos

Rescrevendo a equacao

temos

como desejado.


A probabilidade de que A atinja fixacao e chamada de probabilidade de fixacao. Para o processo de Moran simples, esta probabilidade e

Ja que todos os individuos tem a mesma aptidao, todos terao a mesma chance de se tornar o antecessor de toda a populacao. Esta probabilidade e , logo, a soma de todas as probabilidades (para todos os individuos A) e . O tempo medio ate a absorcao, a partir do estado , e dado por

Para uma derivacao matematica da equacao acima, clique em "expandir" para que ela seja mostrada

O tempo medio gasto no estado j , a partir do estado i , e dado por

Aqui, δ ij denota o delta de Kronecker . Esta equacao recursiva pode ser resolvida usando-se uma nova variavel q i de modo que , logo, > Reescrevendo, temos

A variavel e usada e a equacao se torna

Sabendo que e

podemos calcular :

Portanto

com . Agora, k i , o tempo total do estado i ate a fixacao pode ser calculado


Para um grande, a aproximacao

esta correta.

Selecao [ editar | editar codigo-fonte ]

Se um alelo tem uma vantagem de aptidao comparado a outro alelo, e mais provavel que ele seja escolhido para reproducao. Isto pode ser incorporado ao modelo se individuos com alelo A tiverem a aptidao e os individuos com alelo B tiverem a aptidao g i , em que e o numero de individuos do tipo A, descrevendo assim um processo geral de nascimento e morte. A matriz de transicao do processo estocastico e tridiagonal em forma e as possibilidades de transicao sao [ 3 ]

A entrada denota a probabilidade de se passar do estado para o estado . Para entender as formulas para as probabilidades de transicao, e necessario observar a definicao do processo e ver que a aptidao entra apenas no primeiro termo nas equacoes quanto a reproducao. Assim, a probabilidade de que um individuo A seja escolhido para reproducao nao e mais , mas dependente da aptidao de A, logo,

Tambem neste caso, as probabilidades de fixacao a partir do estado sao definidas pela recorrencia

A forma fechada e dada por

em que por definicao e sera simplesmente para o caso geral.

Para uma derivacao matematica da equacao acima, clique em "expandir" para que ela seja mostrada

Tambem neste caso, as probabilidades de fixacao podem ser computadas, mas as probabilidades de transicao nao sao simetricas. A notacao e e usada. A probabilidade de fixacao pode ser definida recursivamente e uma nova variavel e introduzida.

Agora, duas variaveis obtidas da definicao da variavel y i podem ser usadas para encontrar uma solucao de forma fechada para as probabilidades de fixacao:

Combinando (3) e x N = 1 :

que implica:

Isto, por sua vez, nos da:


Este caso geral em que a aptidao de A e B depende da abundancia de cada tipo e estudado em teoria evolutiva dos jogos.

Resultados menos complexos sao obtidos se uma diferenca constante de aptidoes for assumida. Individuos do tipo A se reproduzem a uma taxa constante e individuos com alelo B se reproduzem a uma taxa 1. Assim, se A tiver uma vantagem de aptidao sobre B, sera maior que um, de outra forma, sera menor que um. Assim, a matriz de transicao do processo estocastico e tridiagonal em forma e as probabilidades de transicao sao

Neste caso, e um fator constante para cada composicao da populacao e, assim, simplificando a probabilidade de fixacao da equacao (1), temos

em que a probabilidade de fixacao de um unico mutante A em uma populacao que, de outra forma, seria inteiramente do tipo B e frequentemente de interesse, sendo denotada por .

Tambem no caso da selecao, o valor esperado e a variancia do numero de individuos A podem ser computados

em que p = i N e r = 1 + s .

Para uma derivacao matematica da equacao acima, clique em "expandir" para que ela seja mostrada

Para o valor esperado, o calculo ocorre como segue abaixo

Para a variancia, o calculo ocorre como segue abaixo, usando a variancia de um unico passo


Taxa de evolucao [ editar | editar codigo-fonte ]

Em uma populacao inteiramente composta por individuos B, um unico mutante A tomara conta de toda a populacao com a probabilidade [ 4 ]

Se a taxa de mutacao (intensidade com que se vai do alelo B ao alelo A) na populacao for , entao, a taxa com a qual um membro da populacao mutacionara a A e dada por e a taxa com a qual toda a populacao ira de inteiramente B a inteiramente A e a taxa a que o alelo mutante A surge vezes a probabilidade de que tomara conta de toda a populacao (probabilidade de fixacao):

Assim, se a mutacao for neutra (isto e, se a probabilidade de fixacao for apenas ), entao, a taxa a que um alelo surge e toma conta da populacao e independente do tamanho da populacao, sendo igual a taxa de mutacao. Este importante resultado e a base da teoria neutralista da evolucao e sugere que o numero de mutacoes pontuais observadas nos genomas de duas especies diferentes seria simplesmente dado pela taxa de mutacao vezes o dobro do tempo decorrido desde a divergencia. Assim, a teoria neutralista oferece um relogio molecular , dado que os pressupostos tenham sido preenchidos, o que pode nao ser o caso na realidade.

Referencias [ editar | editar codigo-fonte ]

  1. Moran, P. a. P. (1 de janeiro de 1958). ≪Random processes in genetics≫ . Mathematical Proceedings of the Cambridge Philosophical Society . 54 (1): 60?71. ISSN   1469-8064 . doi : 10.1017/S0305004100033193  
  2. Nowak, Martin A. (29 de setembro de 2006). Evolutionary Dynamics (em ingles). [S.l.]: Harvard University Press. ISBN   9780674023383  
  3. Moran, Patrick Alfred Pierce (1962). The statistical processes of evolutionary theory (em ingles). [S.l.]: Clarendon Press  
  4. Lieberman, Erez; Hauert, Christoph; Nowak, Martin A. (20 de janeiro de 2005). ≪Evolutionary dynamics on graphs≫ . Nature (em ingles). 433 (7023): 312?316. ISSN   0028-0836 . doi : 10.1038/nature03204