Si lanza una moneda y obtiene 268 caras y 98 colas, puede calcular la probabilidad de que la moneda sea justa de varias maneras. Una observación simple y heurística probablemente concluiría que tal moneda es injusta. He calculado el valor p en R con:
> coin <- pbinom(98, 366, 0.5)
> coin*2
[1] 2.214369e-19
Este valor es menor que .05, por lo tanto, rechazamos la hipótesis de que es una moneda justa.
Pero, ¿qué sucedería si se dijera que la misma moneda cayó de lado 676 veces durante el juicio? Heurísticamente, probablemente llegarás a la misma conclusión, pero ¿seguirían siendo válidas las típicas pruebas de monedas justas?
Aquí hay un gráfico para ilustrar el problema:
¿Cuáles son los métodos válidos para probar la hipótesis de que existe la misma probabilidad de que ocurra un evento en las áreas sombreadas?
NOTA: hay 629 movimientos positivos (413 negativos) en la ilustración gráfica.
Código R que genera los datos:
require("quantmod")
ticker <- getSymbols("SLV")[,6]
change <- (ticker - lag(ticker, 24)) / lag(ticker, 24)
change <- na.locf(change, na.rm=TRUE)
# some other calculations
dens <- density(change)
plot(dens)
# some formatting stuff
fuente