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
☁ [email protected]
☎ (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
[email protected]
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.
