[R-br] set.seed

[Leonard de Assis]

  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 em umontreal.ca <mailto:dm.silva em 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 em 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.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130502/979bdb60/attachment.html>