<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Script para estimar a distribuicao do numero de apostas (entre 200
milhoes) ser igual a sorteada.<br>
<br>
### definindo uma das possibilidades <br>
sorte <- sample(1:choose(60, 6), 1)<br>
<br>
### sortea 200 milhoes de combinacoes (apostas)<br>
i <- sample(1:choose(60,6), 2e8, replace=TRUE)<br>
<br>
### conta quantas sao iguais a combinacao definida <br>
sum(i==sorte)<br>
<br>
### repetindo o prodecimento 100 vezes<br>
### para estimar a distribuicao de frequencias<br>
repl <- replicate(100, sum(i==sample(1:choose(60,6), 1)))<br>
<br>
### resultado<br>
prop.table(table(repl))<br>
<font face="Courier New, Courier, monospace"><br>
## 0 1 2 3 4 5 6 7 8 9 10 11 <br>
##0.02 0.07 0.15 0.19 0.21 0.12 0.06 0.09 0.05 0.02 0.01 0.01 <br>
### mais uma vez<br>
## 1 2 3 4 5 6 7 8 9 10 <br>
##0.05 0.15 0.18 0.27 0.14 0.11 0.05 0.03 0.01 0.01 </font><br>
<font face="Courier New, Courier, monospace">### mais uma vez<br>
## 0 1 2 3 4 5 6 7 8 9 <br>
##0.01 0.11 0.16 0.22 0.16 0.14 0.11 0.05 0.02 0.02<br>
</font><br>
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.<br>
<br>
aT+,<br>
Elias<br>
<br>
<div class="moz-cite-prefix">On 12/09/14 04:24, Marcelo Laia wrote:<br>
</div>
<blockquote
cite="mid:CAEEYVUChBYL9AT_PRGb5btNHN5KsrALPN7o9TgMYuP7-mX5vFw@mail.gmail.com"
type="cite">
<p dir="ltr">Bem interessante! Tenho curiodidade no desfecho dessa
sua postagem.</p>
<p dir="ltr">Laia</p>
<div class="gmail_quote">Em 11/09/2014 23:20, <<a
moz-do-not-send="true" href="mailto:andrebvs@bol.com.br">andrebvs@bol.com.br</a>>
escreveu:<br type="attribution">
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>Alguém aqui já chegou a se questionar sobre a veracidade
das loterias, em especial a MegaSena?<br>
<br>
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?<br>
<br>
Tentei realizar no R todas as 50.063.860 combinações
possíveis de 6 números, mas, o programa não suportou.<br>
<br>
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.<br>
<br>
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? <br>
<br>
<em>Att.</em><br>
<em>André</em></div>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a moz-do-not-send="true"
href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a moz-do-not-send="true"
href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br"
target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a moz-do-not-send="true"
href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>)
e forneça código mínimo reproduzível.<br>
</blockquote>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
</body>
</html>