#--- Prueba de Bondad de Ajuste Ji-Cuadrada ---

#Genera muestra de datos
n<-50
x<-rnorm(n,0,1)
par(mfrow=c(2,1))
hist(x)
boxplot(x,horizontal=T)

#Calcula prueba Ji-Cuadrada con particiones por defecto del comando hist
mu<-mean(x)
sigma<-sd(x)
out<-hist(x,plot=F)
k<-length(out$counts)
delta<-pnorm(out$breaks,mu,sigma)[-1]-pnorm(out$breaks,mu,sigma)[-(k+1)]
qk<-sum((out$counts-delta*n)^2/n/delta)
gl<-k-1-2
1-pchisq(qk,gl)

#Calcula prueba Ji-Cuadrada con particiones fijas dadas por el usuario
inter<-c(-10,-1,0,1,10)
par(mfrow=c(2,1))
hist(x,breaks=inter)
boxplot(x,horizontal=T)
mu<-mean(x)
sigma<-sd(x)
k<-length(inter)-1
out<-hist(x,plot=F,breaks=inter)
delta<-pnorm(inter,mu,sigma)[-1]-pnorm(inter,mu,sigma)[-(k+1)]
qk<-sum((out$counts-delta*n)^2/n/delta)
gl<-k-1-2
1-pchisq(qk,gl)