An analysis of interactions between different physiological control systems may only possible with correlation functions if the signals have similar spectral distributions. Interactions between such signals can be modelled in state space rather than observation space, i.e we model interactions after first translating the observations into a common domain. Coupled Hidden Markov Models (CHMM) are such state-space models. They form a natural extension to standard Hidden Markov Models. In this paper we perform CHMM parameter estimation under a Bayesian paradigm, using Gibbs sampling, and in a Maximum Likelihood framework, using the expectation maximisation algorithm. The performance differences between the estimators are demonstrated on simulated data as well as biomedical data. We show that the proposed method gives meaningful results when comparing two different signals, such as respiration and EEG.