  We demonstrate how to test for conditional independence of two variables with
categorical data using Poisson log-linear models. The size of the conditioning
set of variables can vary from 0 (simple independence) up to many variables. We
also provide a function in R for performing the test. Instead of calculating
all possible tables with for loop we perform the test using the log-linear
models and thus speeding up the process. Time comparison simulation studies are
presented.
