We investigate how the time evolution of different kinetic Ising models depends on the initial conditions of the dynamics. To this end we consider the simultaneous evolution of two identical systems subjected to the same thermal noise. We derive a master equation for the time evolution of a joint probability distribution of the two systems. This equation is then solved within an effective-field approach. By analysing the fixed points of the master equation and their stability we identify regular and chaotic phases.