Chi Square Analysis using for loop in R

by

A sample of your data would be appreciated, but I think this will work for you. First, create a combination of all columns with combn. Then write a function to use with an apply function to iterate through the combos. I like to use plyr since it is easy to specify what you want for a data structure on the back end. Also note you only need to compute the chi square test once for each combination of columns, which should speed things up quite a bit as well.

library(plyr)

combos <- combn(ncol(Dat),2)

adply(combos, 2, function(x) {
  test <- chisq.test(Dat[, x[1]], Dat[, x[2]])

  out <- data.frame("Row" = colnames(Dat)[x[1]]
                    , "Column" = colnames(Dat[x[2]])
                    , "Chi.Square" = round(test$statistic,3)
                    ,  "df"= test$parameter
                    ,  "p.value" = round(test$p.value, 3)
                    )
  return(out)

})

More Related Contents:

  1. Transform the form of the function's result and get the maximum value
  2. Force the origin to start at 0
  3. Capitalize the first letter of both words in a two word string
  4. How to delete multiple values from a vector?
  5. Frequency count of two column in R
  6. Creating arbitrary panes in ggplot2
  7. Selecting multiple odd or even columns/rows for dataframe
  8. Rolling mean (moving average) by group/id with dplyr
  9. Difference: “Compile PDF” button in RStudio vs. knit() and knit2pdf()
  10. How do I convert a factor into date format?
  11. Deleting reversed duplicates with R
  12. Pass a vector of variable names to arrange() in dplyr
  13. How to select a CRAN mirror in R
  14. percentage on y lab in a faceted ggplot barchart?
  15. agrep: only return best match(es)
  16. Fill missing combinations in a dataframe
  17. Cannot install an R package from Github
  18. Conditional binary join and update by reference using the data.table package
  19. How convert decimal to POSIX time
  20. predict.lm() with an unknown factor level in test data
  21. Delete entries with only one observation in a group
  22. Submit form with no submit button in rvest
  23. R: xtable and dates
  24. Using R to read out excel-colorinfo
  25. Existing function for seeing if a row exists in a data frame?
  26. dynamically add function to r6 class instance
  27. Convert comma separated string to integer in R
  28. Difference between Boolean operators && and & and between || and | in R
  29. Interpolating timeseries
  30. Image not showing in Shiny app R

Leave a Comment