Nothing
"cortest.bartlett" <-
function(R,n=NULL,diag=TRUE) {
#message("Bartlett's test that correlation matrix is an identity matrix")
if (dim(R)[1] != dim(R)[2]) {n <- dim(R)[1]
message("R was not square, finding R from data")
R <- cor(R,use="pairwise")}
p <- dim(R)[2]
if(!is.matrix(R) ) R <- as.matrix(R) #converts data.frames to matrices if needed
if(is.null(n)) {n <- 100
warning("n not specified, 100 used") }
if(diag) diag(R) <- 1 #this will make tests of factor residuals correct
detR <- det(R)
statistic <- -log(detR) *(n -1 - (2*p + 5)/6)
df <- p * (p-1)/2
pval <- pchisq(statistic,df,lower.tail=FALSE)
bartlett <- list(chisq = statistic, p.value =pval, df =df)
return(bartlett) }
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.