function(x, catg)
{
	n <- length(x)
	k <- length(catg)
	mu <- mean(x)
	sig <- std(x)
	ax <- x
	ax <- code(x, catg, ax)
	frey <- as.matrix(table(ax), nrow = 
		length(catg))
	exp <- c(1:k)
	exp[1] <- n * pnorm(catg[1], mu, sig
		)
	for(i in 2:k) {
		exp[i] <- n * (pnorm(catg[i],
			mu, sig) - pnorm(
			catg[i - 1], mu, sig
			))
	}
	exp <- as.matrix(exp, nrow = k)
	chi <- sum(frey^2/exp) - n
	pval <- 1 - pchisq(chi, k)
	res <- list(chi, pval)
	names(res) <- c("CHISQ", "P-value")
	res
}
