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
? [email protected] <mailto:[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.
_______________________________________________
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.
