Script para estimar a distribuicao do numero de apostas (entre 200 milhoes) ser igual a sorteada.

### definindo uma das possibilidades
sorte <- sample(1:choose(60, 6), 1)

### sortea 200 milhoes de combinacoes (apostas)
i <- sample(1:choose(60,6), 2e8, replace=TRUE)

### conta quantas sao iguais a combinacao definida
sum(i==sorte)

### repetindo o prodecimento 100 vezes
###  para estimar a distribuicao de frequencias
repl <- replicate(100, sum(i==sample(1:choose(60,6), 1)))

### resultado
prop.table(table(repl))

##   0    1    2    3    4    5    6    7    8    9   10   11
##0.02 0.07 0.15 0.19 0.21 0.12 0.06 0.09 0.05 0.02 0.01 0.01
### mais uma vez
##   1    2    3    4    5    6    7    8    9   10
##0.05 0.15 0.18 0.27 0.14 0.11 0.05 0.03 0.01 0.01
### mais uma vez
##   0    1    2    3    4    5    6    7    8    9
##0.01 0.11 0.16 0.22 0.16 0.14 0.11 0.05 0.02 0.02

ou seja, a probabilidade de nenhum acertar, nessas 300 simulacoes, esta em 2,67% 9, isto é, (2+5+1)/300. Com as ultimas 100 simulacoes, a média do número de acertos foi 3.78.

aT+,
Elias

On 12/09/14 04:24, Marcelo Laia wrote:

Bem interessante! Tenho curiodidade no desfecho dessa sua postagem.

Laia

Em 11/09/2014 23:20, <andrebvs@bol.com.br> escreveu:
Alguém aqui já chegou a se questionar sobre a veracidade das loterias, em especial a MegaSena?

Então, como pode haver milhões de apostas realizadas no país inteiro e apenas um ou duas pessoas ganharem? insuficiente isso não?

Tentei realizar no R todas as 50.063.860 combinações possíveis de 6 números, mas, o programa não suportou.

Gostaria de verificar se por exemplo: de 175 milhões de apostas simuladas aleatoriamente no R, quantas dessas, poderiam conter 6 números definidos antes da simulação.

Pela abordagem subjetiva de probabilidade, acredito que no mínimo 5 cartões seriam sorteados. O que os colegas acham disso, teria como testar essa hipótese no R?

Att.
André

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.


_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.