Test de rachas

En este script se comprueba la diferencia entre generar bits al azar y tratar de hacerlo mentalmente. Se comprueba la aleatoriedad o no aleatoriedad de las series generadas mediante un test no paramétrico: el test de rachas. Se concluye que si se trata de generar aleatoriedad la mente se encuentra con dificultades.

#INICIO -------------------------

rm(list=ls(all=TRUE))

#Cargamos un paquete sobre series temporales

library(tseries)

x<-0:1

n<-50

#Generamos bits al azar

dat<-sample(x,n,replace=T)

#Veamos las longitudes de las rachas

A<-rle(dat);A

#Número de rachas

long<-length(A$lengths);long

#Estas son todas las rachas

A$lengths[1:long]

A$values[1:long]

#Ejecutamos el test de rachas

runs.test(as.factor(dat))

#La serie de bits siguiente se ha generado tratando de hacerlo al azar

datt<-c(0,0,1,0,0,1,0,1,0,1,

1,0,0,1,0,1,0,1,0,0,

1,1,0,1,1,0,0,1,0,1,

0,0,1,1,0,1,1,0,0,1,

0,1,0,0,0,1,1,1,0,1)

AA<-rle(datt);AA

longg<-length(AA$lengths);longg

AA$lengths[1:longg]

AA$values[1:longg]

runs.test(as.factor(datt))

#FIN -------------------------