Estoy tratando de comparar para el mas sencillo (o sea Test de Box-Pierce de aleatoriedad), entre lo que da la función Box.test() y digamos hacerlo a pedal. No me doy cuenta como concluir... va el código R
# Generación de datos
n <- 2000
t1=seq(1,n)#genera un vector
X <- rnorm(n,0,5)# genera n puntos con distribucion N(0,25)
plot(t1,X)
hist(X) # compruebo normalidad
# Primer test de Box-Pierce
Box.test(X)
# da un Pvalue alto siempre, lo que es de esperar
# a pedal
Q <- 0
H <- 10
L <- acf(X)
for(i in 1:H){
Q <- Q + L$acf[i]
}
Q <- n*Q
print(Q)
alfa <- 0.01
qchisq((1-alfa),H)
# No sé como compararlo con lo que da Box.test ???
