A statistical system is disclosed to analyze multivariate response of circadian rhythms in a crossover design, in which response variables are continuously monitored to evaluate the therapeutic effect of a regimen on circadian rhythms such as blood pressure and blood sugar. The methods determine the alteration of not only amplitude but also correlation of multiple circadian rhythms under the influence of a regimen.
