[R-br] set.seed

Daniel Marcelino dmarcelino em live.com
Quarta Maio 1 19:48:53 BRT 2013


Caros, estou realizando um experimento e acabei constatando que a semente
que eu distribuir é fundamental para eu obter as amostras binárias.
Eu estou interessado em distribuir 25 unidades para tratamento e o restante
para controle. Há alguma maneira de saber, sem ser por tentativa e erro,
que minha distribuição terá 25 unidades em cada categoria?
 Por exemplo,

> # 50 Bernoulli trials
> set.seed(51)
> T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50))
> table(T)
T
 0  1
22 28

> set.seed(12345)
> T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50))
> table(T)
T
 0  1
23 27

> set.seed(123)
> T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50))
> table(T)
T
 0  1
25 25

> set.seed(2)
> T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50))
> table(T)
T
 0  1
25 25


-- 
\begin{signature}
Daniel Marcelino
☁ dm.silva em umontreal.ca
☎ (514) 343 6111 #3799
Skype: d.marcelino
✎ 3200 Jean Brillant,  Office C5071
Montreal, QC; H3T 1N8
Canada
\end{signature}

"Small steps toward a much better world"
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130501/dcb29d9a/attachment.html>


Mais detalhes sobre a lista de discussão R-br