segunda-feira, 7 de novembro de 2011

Brasileirão 2011 - Previsões - 50 jogos restando

Após a rodada do último fim de semana, restam apenas 50 jogos para o final do campeonato, cartas estão na mesa, tabela é modificada e ... probabilidades mudam!

Optou-se por rodar uma simulação mais longa (50 mil iterações).

Ainda não modifiquei o algoritmo de simulação para calcular corretamente as probabilidades para se chegar à copa toyota libertadores. Contudo, ontem a noite, ao conversar no Twitter com o Henrique Guzella (Pessoa que faz um ótimo trabalho acompanhando o desempenho das equipes disputantes das 4 divisões do Campeonato Brasileiro que recomendo dar uma conferida), surgiu uma idéia que será implementada em meu tempo livre no decorrer desta semana. A idéia consiste em fazer o óbvio: ao final de cada rodada da simulação, armazenar também a classificação final sem a presença de Vasco e Santos, equipes já classificadas para tal competição.

Com certeza, ao se implementar tal solução, criarei um post específico comentando como implementar tanto a classificação quanto a classificação SEM determinadas equipes, utilizando, como sempre, o software R.

seguem os resultados da simulação com 50 jogos restantes para o final da Série A do Campeonato Brasileiro de 2011:


> probabilities(dados)
                 Camp    top5   last4
América/MG    0.00000 0.00000 0.98952
Atlético/GO   0.00000 0.00002 0.00158
Atlético/MG   0.00000 0.00000 0.04864
Atlético/PR   0.00000 0.00000 0.75750
Avaí          0.00000 0.00000 0.98130
Bahia         0.00000 0.00000 0.14464
Botafogo      0.03962 0.61452 0.00000
Ceará         0.00000 0.00000 0.48174
Corinthians   0.39330 0.97714 0.00000
Coritiba      0.00006 0.07714 0.00000
Cruzeiro      0.00000 0.00000 0.56462
Figueirense   0.01978 0.47110 0.00000
Flamengo      0.02462 0.52454 0.00000
Fluminense    0.17970 0.88874 0.00000
Grêmio        0.00000 0.01226 0.00000
Internacional 0.00184 0.24660 0.00000
Palmeiras     0.00000 0.00000 0.03046
Santos        0.00000 0.06338 0.00000
São Paulo     0.00044 0.16908 0.00000
Vasco         0.34064 0.95548 0.00000


respondendo à pergunta que não quer calar feita por 9 entre 10 torcedores de equipes envolvidas na luta para escapar do rebaixamento, com os dados gerados se pode utilizar o comando quantile do software R para responder esta pergunta:

> with(dados[dados$Pos == 17,], quantile(P, c(0.025, 0.5, 0.975, 0.995)))
 2.5%   50% 97.5% 99.5% 
   37    40    43    44 

Observa-se pelo comando acima que a pontuação da última equipe dentro da região que será relegada à Série B em 2012 possui probabilidade de 97,5% de ser igual ou inferior a 43 pontos e em 99,5% igual ou inferior a 44 pontos.

Quando o caso é ficar entre os 5 melhores, o comando:

> with(dados[dados$Pos == 5,], quantile(P, c(0.025, 0.5, 0.975, 0.995)))
 2.5%   50% 97.5% 99.5% 
   59    61    64    64 
nos informa que 64 pontos pode ser um valor mais que suficiente para garantir presença neste seleto grupo.

Quanto ao título, o comando abaixo:

> with(dados[dados$Pos == 1,], quantile(P, c(0.025, 0.5, 0.975, 0.995)))
 2.5%   50% 97.5% 99.5% 
   65    68    73    73 
nos diz que a exigência pelo título permanece na casa dos 73 pontos.