[R-es] Aleatoriedad

[Jorge I Velez]

Hola Miguel,

Quizas lo siguiente ayude un poco a entender que ocurre:

# para reproducir los resultados
set.seed(123)

# proporcion de 'todos cero'
all0 <- mean(rbinom(1000, 8, .5) == 0)
all0

# proporcion de 'todos uno'
all1 <- mean(rbinom(1000, 8, .5) == 0)
all1

# distribution de 'todos cero'
hist(replicate(1000, mean(rbinom(1000, 8, .5) == 0)), las = 1, main =
'Todos cero', xlab = 'Proporción')
abline(v = all0, col = 2, lwd = 2)

# distribution de 'todos uno'
hist(replicate(1000, mean(rbinom(1000, 8, .5) == 8)), las = 1, main =
'Todos uno', xlab = 'Proporción')
abline(v = all1, col = 2, lwd = 2)

Como veras, la proporcion de 'todos cero' o 'todos uno' varia de manera
considerable, pero en promedio (linea vertical roja) es cercano al valor
que mencionas.

Mas informacion en ?set.seed,  ?rbinom, ?hist t ?replicate.

Saludos cordiales,
Jorge.-



2015-03-25 6:46 GMT+11:00 Our Utopy <ourutopy en gmail.com>:

> Hola de nuevo, ya empiezo a ser pesado ¿no? bueno, no importa porque
> aprendemos todos. Eso, al menos, me parece.
>
> Hoy estuve estudiando en R el tema de la aleatoriedad. Veo que hay
> múltiples posibilidades pero me están chocando mucho. Encuentro que el
> generador de números pseudo aleatorios  es más pseudo de lo que debería.
>
> Me explico, quiero generar 0 y 1 aleatorios. Estoy trabajando con una
> martingala.
>
> Uso x <- runif(1,0,1) y después según cual sea la probabilidad de ganar
> elijo el punto de división del intervalo apropiado.
>
> Bien, los modelos que extraigo no me parecen coherentes. Me gustaban mucho
> más los obtenidos con el aleatorio() de Microsoft Excel. O eso me parece
>
> ¿Por qué? pues porque la probabilidad de sacar 8 "ceros" u 8 "unos"
> seguidas es 0.003 y a mi me está pareciendo que lo hace mucho más a menudo
> de lo debido.
>
> Y ora cosa, ¿sabéis si cuando apago el ordenador me aparecen al iniciarlo
> de nuevo las mismas series aleatorias, en el mismo orden que en el trabajo
> previo?  y de ser así ¿existe alguna función en R que las haga iniciar por
> otro lugar de la serie interna predefinida?
>
> ¿Sabéis si hay test de aleatoriedad para series predefinidos en R? Es que
> el buscador de funciones ?? me busca en los paquetes incluidos pero no
> instalados pero ... ¿hay algo para buscar en el repositorio?
>
> Gracias ya anticipadas. Un saludo
>
>
>
>
>
> --
> Our Utopy
> http://utopicmaths.blogspot.com.es/
> http://financialmeth.blogspot.com.es/
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>

	[[alternative HTML version deleted]]