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"
_______________________________________________
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.