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.
Nenhum comentário:
Postar um comentário