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
--
Skype:
d.marcelino
✎ 3200 Jean
Brillant, Office C5071
Montreal,
QC; H3T 1N8
Canada
\end{signature}
"Small steps
toward a much better world"