<div dir="ltr"><div><div>Vou checar a referência Benilton, obrigado. Fiquei curioso em saber se é possível estabelecer padrões de semente para se obter resultados semelhantes dado o vetor de probabilidades.<br><br></div>Abs<br>
</div>Daniel<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/5/1 Benilton Carvalho <span dir="ltr"><<a href="mailto:beniltoncarvalho@gmail.com" target="_blank">beniltoncarvalho@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Daniel,<br>
<br>
eu to, na real, confuso....<br>
<br>
Pela descricao do seu problema, vc possui 50 unidades... Destas, 25<br>
devem ser alocadas a um tratamento... e as remanescentes, a um segundo<br>
tratamento.<br>
<br>
Se for este o problema, ele nao e' modelado com uma Bernoulli...<br>
(lembre das hipoteses da Bernoulli: a) dois resultados; b)<br>
probabilidade de sucesso constante; c) ensaios independentes)<br>
<br>
Sendo como descrevi acima, tudo o q vc quer e' aleatorizar os labels<br>
"trat1" e "trat2" (25 de cada) nos individuos disponiveis, como<br>
sugerido pelo Gledson.<br>
<br>
sample(rep(0:1, 25))<br>
<br>
Agora, se vc quiser "mmmmeeeeiiiiixxxxxmo" entender o lance da semente<br>
q faz uma Bernoulli distribuir 25/25, entao e' ler:<br>
<br>
Matsumoto, M. and Nishimura, T. (1998) Mersenne Twister: A<br>
623-dimensionally equidistributed uniform pseudo-random number<br>
generator, _ACM Transactions on Modeling and Computer Simulation_,<br>
*8*, 3-30.<br>
<br>
abs, b<br>
<br>
Em 1 de maio de 2013 22:10, Daniel Marcelino <<a href="mailto:dmarcelino@live.com">dmarcelino@live.com</a>> escreveu:<br>
<div class="HOEnZb"><div class="h5">> Valeu Gledson, a solução eu sabia. Eu estou procurando mais informações<br>
> sobre como descobrir o padrão da semente para que isso aconteça, como nos<br>
> dois últimos exemplos que eu enviei.<br>
><br>
> Daniel<br>
><br>
><br>
><br>
> 2013/5/1 gledson picharski <<a href="mailto:gledson.picharski@gmail.com">gledson.picharski@gmail.com</a>><br>
>><br>
>> sorteia sem reposição um conjunto que tenha 25 zeros e 25 uns.<br>
>><br>
>> ## solução 1<br>
>> set.seed(111)<br>
>> sample(rep(c(0,1),25))<br>
>><br>
>> ## solução 2<br>
>> set.seed(111)<br>
>> rep(c(0,1),25)[sample(1:50)]<br>
>><br>
>><br>
>> _________________<br>
>> Gledson Picharski<br>
>><br>
>><br>
>><br>
>><br>
>><br>
>> Em 1 de maio de 2013 19:48, Daniel Marcelino <<a href="mailto:dmarcelino@live.com">dmarcelino@live.com</a>><br>
>> escreveu:<br>
>>><br>
>>> Caros, estou realizando um experimento e acabei constatando que a semente<br>
>>> que eu distribuir é fundamental para eu obter as amostras binárias.<br>
>>> Eu estou interessado em distribuir 25 unidades para tratamento e o<br>
>>> restante para controle. Há alguma maneira de saber, sem ser por tentativa e<br>
>>> erro, que minha distribuição terá 25 unidades em cada categoria?<br>
>>> Por exemplo,<br>
>>><br>
>>> > # 50 Bernoulli trials<br>
>>> > set.seed(51)<br>
>>> > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50))<br>
>>> > table(T)<br>
>>> T<br>
>>> 0 1<br>
>>> 22 28<br>
>>><br>
>>> > set.seed(12345)<br>
>>> > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50))<br>
>>> > table(T)<br>
>>> T<br>
>>> 0 1<br>
>>> 23 27<br>
>>><br>
>>> > set.seed(123)<br>
>>> > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50))<br>
>>> > table(T)<br>
>>> T<br>
>>> 0 1<br>
>>> 25 25<br>
>>><br>
>>> > set.seed(2)<br>
>>> > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50))<br>
>>> > table(T)<br>
>>> T<br>
>>> 0 1<br>
>>> 25 25<br>
>>><br>
>>><br>
>>> --<br>
>>> \begin{signature}<br>
>>> Daniel Marcelino<br>
>>> ☁ <a href="mailto:dm.silva@umontreal.ca">dm.silva@umontreal.ca</a><br>
>>> ☎ <a href="tel:%28514%29%20343%206111%20%233799" value="+15143436111">(514) 343 6111 #3799</a><br>
>>> Skype: d.marcelino<br>
>>> ✎ 3200 Jean Brillant, Office C5071<br>
>>> Montreal, QC; H3T 1N8<br>
>>> Canada<br>
>>> \end{signature}<br>
>>><br>
>>> "Small steps toward a much better world"<br>
>>><br>
>>> _______________________________________________<br>
>>> R-br mailing list<br>
>>> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
>>> <a 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 href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
>>> código mínimo reproduzível.<br>
>><br>
>><br>
>><br>
>> _______________________________________________<br>
>> R-br mailing list<br>
>> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
>> <a 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 href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
>> código mínimo reproduzível.<br>
><br>
><br>
><br>
><br>
> --<br>
> \begin{signature}<br>
> Daniel Marcelino<br>
> ☁ <a href="mailto:dm.silva@umontreal.ca">dm.silva@umontreal.ca</a><br>
> ☎ <a href="tel:%28514%29%20343%206111%20%233799" value="+15143436111">(514) 343 6111 #3799</a><br>
> Skype: d.marcelino<br>
> ✎ 3200 Jean Brillant, Office C5071<br>
> Montreal, QC; H3T 1N8<br>
> Canada<br>
> \end{signature}<br>
><br>
> "Small steps toward a much better world"<br>
><br>
> _______________________________________________<br>
> R-br mailing list<br>
> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
> <a 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 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<br>
> mínimo reproduzível.<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a 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 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.</div></div></blockquote></div><br><br clear="all"><br>-- <br>
<div style="font-family:Calibri,sans-serif;font-size:16px"><font face="Verdana"><span style="font-family:arial;font-size:10pt">\begin{signature}</span><br style="font-family:arial;font-size:small"><span style="font-family:arial;font-size:10pt">Daniel Marcelino</span><br style="font-family:arial;font-size:small">
<span style="font-family:arial;font-size:10pt">☁ </span><a href="mailto:dm.silva@umontreal.ca" style="font-family:arial;font-size:small" target="_blank">dm.silva@umontreal.ca</a><span style="font-family:arial;font-size:10pt"> </span><br style="font-family:arial;font-size:small">
<span style="font-family:arial;font-size:10pt">☎ (514) 343 6111 #3799</span></font></div><div style="font-family:Calibri,sans-serif;font-size:16px"><font face="Verdana"><span style="font-family:arial;font-size:13px">Skype: d.marcelino </span></font></div>
<div style="font-family:Calibri,sans-serif;font-size:16px"><font face="Verdana"><span style="font-family:arial;font-size:10pt">✎ 3200 Jean Brillant, Office C5071</span><br style="font-family:arial;font-size:small"><span style="font-family:arial;font-size:10pt">Montreal, QC; H3T 1N8</span><br style="font-family:arial;font-size:small">
<span style="font-family:arial;font-size:10pt">Canada</span><br style="font-family:arial;font-size:small"><span style="font-family:arial;font-size:10pt">\end{signature}</span></font></div><div style="font-family:Calibri,sans-serif;font-size:16px">
<font face="Verdana"><span style="font-family:arial;font-size:13px"><br></span></font></div><div style="font-family:Calibri,sans-serif;font-size:16px"><font face="Verdana"><span style="font-family:arial;font-size:13px">"Small steps toward a much better world"</span></font></div>
</div>