✎ Dinâmica de soma zero

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.

<WRAP tabs>
  * [[en:ecovirt:roteiro:math:zerosumrcmdr|{{:ecovirt:logorcmdr01.png?20|}}]]
  * [[en:ecovirt:roteiro:math:zerosumr|{{:ecovirt:rlogo.png?20|}}]]
</WRAP>



====== Dinâmica de soma zero - Roteiro em R======

{{http://imgs.xkcd.com/comics/ultimate_game.png}}

O conceito de [[http://en.wikipedia.org/wiki/Zero-sum_game|de soma zero]] vem da [[http://en.wikipedia.org/wiki/Game_theory|teoria dos jogos]], e descreve a divisão de uma quantia fixa entre participantes, de modo que só se ganha o que outros perdem.  

Se os ganhos e consequentes perdas acontecem com uma certa probabilidade, o jogo torna-se uma dinâmica estocástica, como na [[https://pup.princeton.edu/chapters/s7105.html|teoria neutra da biodiversidade]]. Seu criador, [[http://shubbell.eeb.ucla.edu/people.php|Stephen Hubbell]], assumiu que as comunidades estão saturadas, de modo que um novo indivíduo só se estabelece se outro morre. A sucessão ao acaso de mortes, nascimentos e chegada de migrantes criaria então uma dinâmica de zoma zero, que explicaria vários padrões das comunidades.

Neste tutorial, simulamos uma dinâmica estocástica de soma zero muito simples no ambiente R. Depois disso, você pode estudar a aplicação desse modelo no roteiro sobre [[en:ecovirt:roteiro:neutr:neutrar|teoria neutra da biodiversidade]].

<box red>
Para prosseguir você deve ter o ambiente R instalado. Se você não tem e não sabe como ter, consulte a página de  [[en:ecovirt:roteiro:soft:instalacaor|Instalação]].
</box>

===== Um Joguinho Besta =====

Vamos imaginar um jogo de apostas entre dois jogadores, sem empates. A cada rodada o perdedor da aposta paga uma quantia fixa ao ganhador. Os dois jogadores têm a mesma probabilidade de ganhar a cada rodada. Esse é [[http://en.wikipedia.org/wiki/Zero-sum|um jogo de soma zero]], pois o valor total em jogo não se altera. O que muda é apenas a fração deste total em poder de cada jogador. Note que a propriedade de soma zero é possível mesmo que os jogadores tenham chances diferentes de vencer, ou que a quantia que o perdedor paga também seja sorteada.

Em nossa simulação, [[http://pt.wikiquote.org/wiki/Vicente_Matheus|o jogo só termina quando acaba]], ou seja, quando um dos dois jogadores perde todo o dinheiro((ou quando você quiser parar o jogo, se achar que está demorando muito)). Vamos simular esta situação com uma função em R. Copie o código da função abaixo para a linha de comando do R:

<code>
jogo <- function(aposta=1,total=20){
  X <- total/2
  results <- X
  while(X>0&X<total){
    X <- X+sample(c(aposta,-1*aposta),1)
    results <- c(results,X)
  }
  plot(1:length(results),results, type="l", col="blue",ylim=c(0,total), xlab="N de rodadas", ylab="Valor")
  lines(1:length(results),total-results, col="red")
  abline(h=c(0,total),lty=2)
}
</code>

Esta função tem apenas dois argumentos:

  * ''aposta'': que é o valor pago pelo perdedor a cada aposta
  * ''total'': o valor total de dinheiro em jogo. No início este total é dividido igualmente entre os jogadores.

A simulação transcorre até o final do jogo. Para simular um jogo com apostas de uma unidade e total em jogo de cem unidades execute o comando no R:

<code>
jogo(aposta=1, total=100)
</code>

==== O que afeta o tempo de jogo? ====

A simulação transcorre até o final do jogo, ou até o tempo máximo se esgotar.  Varie o total em jogo e o valor da aposta e avalie seu efeito na duração do jogo. Sugestões de valores:


  * total = 100, aposta = 1
  * total = 100, aposta = 10
  * total = 1000, aposta = 1
  * total = 1000, aposta = 10



<WRAP center round tip 60%>
Em dinâmicas estocásticas o resultado varia a cada vez, mesmo que os parâmetros sejam os mesmos. Por isso repita cada simulação algumas vezes para assegurar-se dos resultados.
</WRAP>



==== Perguntas ====

  - Qual o efeito do aumento do total em jogo e do tamanho da aposta sobre o tempo para que o jogo acabe?
  - Este jogo também é um processo de [[en:ecovirt:roteiro:math:bebador|caminhada aleatória em uma dimensão]]. Explique porque.
User Tools

Site Tools

Page Tools
