[R-br] Veracidade das Loterias

[Elias T. Krainski]

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 em bol.com.br 
> <mailto:andrebvs em 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 em listas.c3sl.ufpr.br <mailto:R-br em 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 em 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.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140912/5f609495/attachment.html>