Scripted example of using chisq to mine for relations between nominal variables

chisq_mining.r

testvars = c(6,7,9,10,11,12,13,14,16, 17,18,19,20,21,23,24,25,26,384,375,376,386,385,387,388)

resultlist = c()
for (i in testvars) {
  xsq = chisq.test(big.dataset[,i], big.dataset$DV_3lvls)$statistic
  varname = names(big.dataset)[i]
  tab = xtabs(~DV_3lvls + big.dataset[,i], data=big.dataset)
  resultlist = rbind(resultlist, list(chisq=xsq, testvar=varname, xtab=tab))
}