Prueba de esferidad de Barlett

Formulación de Hipótesis:

-H0: Las variables no estan correlacionadas.

-H1: Las variables estan correlacionadas.

t.bartlett <- function(data){

n <- dim(data)[1]
p <- dim(data)[2]
R <- cor(data)

chi.square <- -(n-1-(2*p+5)/6) * log(det(R))

cat("\n","Prueba de Esferecidad de Bartlett","\n","\n")
cat('chi.square = ', round(chi.square,3) ,
', df = ', (p^2-p)/2,
', p-value = ', 1-pchisq(chi.square,(p^2-p)/2))

}