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@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"