table(sample(rep(0:1,25),50,rep=FALSE)) 0 1 25 25 Em 01/05/2013 19:48, Daniel Marcelino escreveu:
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 <mailto: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"
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.