- Research
- Open access
- Published:

# Statistical methods for detecting and comparing periodic data and their application to the nycthemeral rhythm of bodily harm: A population based study

*Journal of Circadian Rhythms*
**volume 8**, Article number: 10 (2010)

## Abstract

### Background

Animals, including humans, exhibit a variety of biological rhythms. This article describes a method for the detection and simultaneous comparison of multiple nycthemeral rhythms.

### Methods

A statistical method for detecting periodic patterns in time-related data via harmonic regression is described. The method is particularly capable of detecting nycthemeral rhythms in medical data. Additionally a method for simultaneously comparing two or more periodic patterns is described, which derives from the analysis of variance (ANOVA). This method statistically confirms or rejects equality of periodic patterns. Mathematical descriptions of the detecting method and the comparing method are displayed.

### Results

Nycthemeral rhythms of incidents of bodily harm in Middle Franconia are analyzed in order to demonstrate both methods. Every day of the week showed a significant nycthemeral rhythm of bodily harm. These seven patterns of the week were compared to each other revealing only two different nycthemeral rhythms, one for Friday and Saturday and one for the other weekdays.

## Background

Analysis of biological activities that fluctuate throughout the day is common in various fields of medicine. Blood pressure and heart rate as well as the occurrence of acute cardiovascular disease are subject to a twenty-four hour rhythm (also referred to as circadian or nycthemeral rhythm) [1, 2]. This rhythm is also present in episodes of dyspnoea in nocturnal asthma [3], intraocular pressure [4, 5], and hormonal pulses [6–8]. Nycthemeral fluctuations in neurotransmitters and hormones have been discussed as influencing human behavior [9–11]. Suicide as well as parasuicide and violence against the person show day-night variation [12–14]. Assaults presenting to trauma centers display a distinct nycthemeral pattern [8–12]. In this study the nycthemeral rhythm of violent crime rates is analyzed to demonstrate a detection method and a comparison method suitable for twenty-four hour time series, but not limited to this sampling period.

Much mathematical effort was invested to detect and model the dependency on the time of day [15–19]. A classification of the data by identifying similarities and distinctions requires statistical methods [20–25].

The cosinor analysis is a common approach [26] that describes data by a single cosine function with fixed frequency plus a constant (single-harmonic model) yielding the three parameters amplitude, phase and mean [27]. Corresponding parameters were compared one by one to compare two or more time series modeled by cosinor analysis [28, 29]. A multivariate technique is applied in this study aiming to compare several periodic patterns simultaneously. Models allowing more than one frequency (multi-harmonic model) show no graphic equivalent for the parameters amplitude and phase. Multi-harmonic models have been used to describe human core-temperature [18], blood pressure and incidence of angina [23] as well as in the nycthemeral distribution of violent crime rates, although the true waveform of nycthemeral rhythms is still a matter of debate. The purpose of this study is to identify the underlying frequencies and to compare the resulting periodic patterns via Fourier transform. This transform is common use in various fields of medicine [16] as well as other scientific areas. The explained variance of individual oscillations is utilized to detect the inherent periodic patterns of the data.

A modification of the analysis of variance (ANOVA) is used to compare two or more time series with periodic patterns. The typical ANOVA tests whether the means of several groups are equal. The scope of ANOVA is extended to periodic patterns by combining it with Fourier analysis. This new test rejects or confirms equality of multiple oscillating time series.

To demonstrate both methods, the oscillations of violent crimes in Middle Franconia, Bavaria/Germany from 2002 to 2005, were analyzed. Nycthemeral rhythms of bodily harm were identified on all seven days of the week. The seven patterns of the week were compared to each other revealing only two different nycthemeral rhythms. We demonstrate that the nycthemeral rhythms on Friday and Saturday are equal and differ significantly from the rhythms of the other weekdays, which are then equal again.

To compare our method with the cosinor method an analysis of the same data is performed and yields no strong evidence of different rhythms.

The simultaneous comparison of a greater number of nycthemeral rhythms is made possible by the use of the mathematical methods described in this study. A need for such procedures derives from the prospect of developing a prediction model for violent crime rates which is of immediate interest for public services such as social facilities, police departments and hospitals.

The section detection method contains a procedure to find the inherent frequencies of the data, the section Fourier Anova describes the comparison method, the results section illustrates both methods by analyzing nycthemeral rhythms of offenses against the person causing bodily harm and in the conclusion limitations, modifications and alternatives to our methods are discussed.

## Methods

### Detection method

A statistical test for finding the frequencies of oscillating data is described. Using harmonic frequencies the data are modeled as a sum of sine and cosine oscillations and a Fourier transform is performed. In our case the Fourier transform equals an ordinary least squares. All frequencies are tested for significance. The ratio of explained variance of a frequency and remaining variance acts as test statistic. Model selection is carried out by a Bonferroni-Holm Method (see [30]).

Fitting harmonic models to nycthemeral rhythms is a common procedure [31–33]. The detection method is ancillary, its output is used as input for the comparison method (see section Fourier ANOVA). From a numerical vantage point linear least squares with orthonormal regressors are applied. From a linear algebra perspective we choose a specific set of vectors forming an orthonormal basis and change basis. Statistical methods are applied to search for single coordinates of the data (relative to the new basis) that are „large“ compared to the other coordinates. The orthonormal basis ensures independent and normal distributed regression coefficients; thus choosing significant frequencies (i.e. model selection) is straightforward. Furthermore the orthonormal regressors are necessary for our extension of ANOVA described in the section Fourier ANOVA.

The model for our data is

with white noise ϵ. Constant terms are omitted. So a time series sampled *n* times with a fixed sampling interval, homoscedasticity and uncorrelated noise and without a linear trend or missing values is assumed. The regressors have the harmonic frequencies

By this choice the regressors cos(2*πf*
_{
j
}
*t*) and sin(2*πf*
_{
j
}
*t*) are an orthogonal basis of *R*
^{n}. Estimating *a* and *b* with ordinary least squares against the normalized regressors yields independent and normal distributed coefficients. To determine significant frequencies we search for large coefficients *a* and *b* by a method similar to a Wald-statistic and by a Bonferroni-Holm procedure [30].

The null hypotheses are : *a*
_{
j
} = *b*
_{
j
} = 0, or: „no significant periodic pattern with frequency *f*
_{
j
} in the data“. To test these hypotheses

is calculated, mimicking a periodogram. The value can be interpreted as the explained variance of frequency *f*
_{
j
}, furthermore *c*
_{
j
} is invariant under time-shift of the data. Then *c* is sorted in descending order a *F* -distributed test statistic is calculated:

which is tested on the corrected significance level . If *T*
_{
j
} does not exceed the critical value for a specific *j*, then all *T*
_{
i
} with *i > j* are not tested anymore. This test yields a set of significant frequencies. A Fourier approximation ℱ_{
F
} of the data is obtained by evaluating equation 1 using only a subset *F* of the harmonic frequencies (e.g. the significant frequencies) and their corresponding amplitudes:

The Fourier approximation filters the periodic components out of the data; it is a denoising procedure. The data is decomposed in a fundamental frequency and its multiple, the harmonics. The Fourier coefficients indicate the strength, i.e. the amplitude of these oscillations. Usually the fundamental frequency has the highest amplitude and the strength decreases for greater harmonics. The influence of the harmonics can reach from only small adjustments of the fundamental oscillations to generating additional maxima, minima or plateaus.

### Comparison method (Fourier ANOVA)

A statistical test for comparing periodic patterns of grouped data is described. The test determines if the rhythm of the groups are equal or not. The mathematical concept of the ANOVA is transferred to periodic patterns by substituting the mean estimators for Fourier approximations. This test compares the periodic patterns in its entirety. The orthogonal regressors mentioned in the section Detection method are necessary for this test.

Suppose data divided in *k* groups with *n* measurements for every group and denote this data as *x*
_{
t,j
} (*t* = 1... *n, j* = 1... *k*). The *F* distributed ANOVA test statistic for equal means in every group is

To compare not the means but the periodic pattern of every group we substitute the mean estimators for the Fourier approximation (see 5):

The frequencies *F* are chosen as described in the section Detection method: the detection method is applied to every group of *x*. Testing with *d* = |*F*| frequencies the degrees of freedom are *df*
_{1} = 2*dk* - 2*d* and *df*
_{2} = *nk* - 2*dk*.

The test uses the same idea as the ANOVA: Calculate the variance within the groups, i.e. the deviation of the data from its Fourier approximation within every group. Furthermore calculate the variance between the groups, i.e. the deviation the Fourier approximation of the single groups and the Fourier approximation of the whole data. If all groups show the same rhythm then the variance between the groups should have roughly the same magnitude as the variance within the groups. Conversely a large variance between the groups argues for an impact of a group on the rhythm.

In the following we will scrutinize the distribution of the test statistic in equation 7: We show that the test statistic *T*
_{
F
} in equation is *F* distributed. Cochran's Theorem, as stated in [34], yields a *χ*
^{2} distribution of the nominator and the denominator of equation 7. To apply this theroem the test statistic needs a matrix representation.

The Fourier approximation in equation 5 has a matrix representation: For *f* ∈ ℝ define the column vectors

with normalization constants *k*
_{
c
}(*f, n*)*, k*
_{
s
}(*f, n*). Then then Fourier approximation can be written as

Let be this transformation matrix of ℱ_{
F
}, then is a symmetric projection, i.e.

Furthermore pile the columns of the data *x* ∈ *R*
^{n,k}one below the other and call this vector *y*
*∈*
*R*
^{nk}. Define the matrices

and

Because *A*
_{1} and *A*
_{2} are symmetric projections the test statistic *T*
_{
F
} in equation 7 can be written as

Now the test statistic has a representation suitable for Cochran's Theorem. All that is left is the orthogonality assumption for the projections *A*
_{2} - *A*
_{1} and . The specific form of the harmonic frequencies is again utilized: The image of *A*
_{1} is spanned by and (*f* ∈ *F* ). The image of *A*
_{2} is spanned by the vectors and filled up with the zero vector **0** = (0...0) ∈ ℝ^{n}:

By definition of the harmonic frequencies (see equation 2) the following equation holds except for normalization factor:

So the image of *A*
_{1} is a subset of the image of *A*
_{2} and it holds:

This equation shows the orthogonality of the projections of Cochran's Theorem.

## Results

Nycthemeral rhythm of violent crime rates are analyzed to demonstrate both the detection and comparison method.

The study included 15881 crimes of violent behavior (without suicides) which were filed at the Police Department of Middle Franconia, Bavaria/Germany between January 1, 2002 and December 31, 2005, and gathered into the EVioS (Erlangener Violence Studies [35]) data base. Bodily harm as defined in § 223 German Criminal Code is more closely examined. We investigate if the seven days of the week show different nycthemeral rhythms of bodily harm. Data handling and calculations were performed by Microsoft Excel^{®}, Matlab^{®} and R. Significance level was set to 0.05.

In the following, the detection method shows the existence of nycthemeral rhythms of bodily harm on all seven days of the week. A comparison of these seven rhythms reveals only two different nycthemeral rhythms, one describing crime rates on Friday and Saturday, the other on Sunday to Thursday. In order to analyze a more homogeneous sample, only crimes committed by male offenders and not occurring on holidays such as New Year's Eve are further surveyed; this sample consists of 11402 cases. The investigated data *x* ∈ ℝ^{24 × 7} are the number of violent acts *x*(*h, d*) at a specific hour *h*
*∈* {1... 24} and „day“ *d* ∈ {1... 7}. We define the first „day“ as the 24 hours starting Sunday at 9:00 a.m. and denote it with *d*1. This definition is adapted to the data: at 9:00 a.m. violent crime rates of all seven days are similar and a renewal of the time series occurs (see Figure 1). Furthermore the second „day“ *d*2 is defined as the 24 hours starting Monday at 9:00 a.m. lasting till Tuesday 9:00 a.m. and so forth.

The histogram in Figure 2 shows the distribution of violent crimes per „day“ with 95% confidence intervals. In particular the number of crimes on *d*6 and *d*7 are distinct. We are interested in the nycthemeral rhythm and not in total numbers; so we normalize the data by dividing the number of crimes at „day“ *d* and hour *h* by the number of crimes on „day“ *d*. The normalized data are called *y* ∈ ℝ^{24 × 7}. So every column of *y* sums up to 1 and thus can be interpeted as relative frequency of crimes.

The assumptions of our model in equation 1 are satisfied by the data *y*: There is no trend or missing values and a constant time between two consecutive samples. *x* consists of count data, so *x*(*h, d*) follows a Poisson distribution and the normalized data *y*(*h, d*) are well approximated by a normal distribution. The sequence *y*(*h*, *d*)_{
h = 1... 24,d = 1... 7 }is assumed to be independent, because sites of crimes are spatially separated or offenders don't even know each other. Homoscedasticity (constant variance of the residuals) and Poisson distributions do not make a good match: For Poisson random variable the mean equals the variance and we assume a oscillating number of crimes. So the residuals will not automatically be homoscedastic and are afterwards tested for „whiteness“ by a Kolmogorov-Smirnov test [36], a Lilliefors test [37] (both for normal distribution), a Breusch-Godfrey test [38, 39] and a Wald Wolvowitz runs test [40] (for absence of autocorrelation, the latter is applied to the signs of the residuals). The data are also tested for stationary cycles by a Canova-Hansen [41] test and a Kwiatkowski-Phillips-Schmidt-Shin test [42]. We also divided the data in 10 disjoint random subsamples to avoid testing hypotheses suggested by the data.

Applying the detection method to the columns of *y* reveals significant nycthemeral rhythm on every „day“. All seven „days“ showed significant periods of length 24 and 12 hours except *d*3 and *d*4, which showed only a significant 24 hour period. So every „day“ shows a nycthemeral rhythm of bodily harm. Note that by analyzing single days of the week, i.e. columns of *y*, which have a length of 24, we restrict our search to the frequency and its integer multiple (see the model in equation 1 and its description). We have two reason for doing so: first we have a priori knowledge: The sun is a zeitgeber for the human biological clock [43], that argues for a 24 hour rhythm. Furthermore the week is the time unit that governs the working life in Germany and separates it in five working days (Monday to Friday) and two weekend days (Saturday and Sunday). Second we get a posteriori knowledge: by applying our detection method to the whole data *y* which revealed no other significant periods, especially no significant period greater than 24 hours and by calculating a periodogram of the data (see Figure 3), which reveals only a day period, a week period and their corresponding harmonics.

Applying the comparison method to *d*1 to *d*7 (frequencies and number of samples *n* = 7·24) generates a p-value smaller than 0.05 (*F* = 18.1639, *df*
_{1} = 24, *df*
_{2} = 140). So there are at least two different periodic patterns in the data. This finding is verified in the 10 randomly-generated subsamples: comparing the period of the subsamples yields p-values within the interval [1.04 · 10^{-10}, 1.3 · 10^{-3}].

Comparing *d*6 and *d*7 (*n* = 2 · 24, *f* as above) yields a p-value of 0.3582 (*F* = 1.1352, *df*
_{1} = 4, *df*
_{2} = 40). So there is no significant difference between *d*6 and *d*7. So Friday and Saturday show the same nycthemeral rhythm of bodily harm. By testing this hypothesis in the 10 subsamples p-values in the interval [0.15, 0.98] are obtained.

We found that nycthemeral rhythm of *d*6 and *d*7 is different from the rhythm of *d*1 to *d*5. For statistical verification, the comparison method was applied to the 93 partitions *P* ⊂ {1...7} of the seven days, that contain at least one element of {1, 2, 3, 4, 5} and at least one element of {6, 7}. So we tested the 93 hypothesis *H*
_{
P
} : „There is no significant difference between the „days“ of partition *P*“. The comparison method yields p-values smaller than 8.5 · 10^{-11}. Bonferroni's inequality yields an upper bound for the p-value of the hypothesis ∪ *H*
_{
P
} („There is at least one of the 93 partitions without significant difference between the nycthemeral rhythms of the „days“ of this partition“): *P* (∪*H*
_{
P
} ) < 93 · 8.5 · 10^{-11}
*<*0.05 and thus reject this hypothesis. We accept the alternative hypothesis: „the „days“ of all 93 partitions have significant different nycthemeral rhythms“.

Comparing only *d*1 to *d*5 (*n* = 5 · 24, *f* as above) yields a p-value of 0.0515 (*F* = 1.7372, *df*
_{1} = 16, *df*
_{2} = 100). Applying this test to the 10 subsamples yields p-values within [0.0457, 0.93], one p-value was lower than 5%. Testing the 26 partitions of {1... 5}, which have at least two elements yields p-values ranged from 0.0047 to 0.9908, none was smaller than Bonferroni-corrected significance level . Altogether we found some significant differences within *d*1 to *d*5, but consider them marginal. So there are only two significantly different nycthemeral rhythms, one describing crime rates on *d*6 and *d*7, the other on *d*1 to *d*5, see Figure 1 for a plot of these two rhythms.

Now the „whiteness“ of the residuals of the fit of *d*1 to *d*5 is tested. Figure 4 shows a quantile-quantile-plot of the residuals against a standard normal distribution, which is almost linear, arguing for normal distributed residuals. A formal test for normal distribution is the Kolmogorov-Smirnov test. Testing the residuals divided by their estimated standard deviation against the standard normal distribution yields *p* = 0.96, *d*
_{
ks
} = 0.0440, *n* = 120.

Autocorrelation of the residuals biases the estimation of the coefficients and is a evidence for a misspecified model. A Breusch-Godfrey test for autocorrelation up to order 23 does also not reject the null hypothesis (*p* = 0.155, = 29.8). So these residuals show no significant autocorrelation.

Stationarity is a property often desired in time series analysis, particular in econometrics [44, 45]. A stationary process fluctuates steadily around a deterministic trend, a nonstationary series is subject to persistent random shocks or can even be transient. If the variables in the regression model are not stationary, then the standard assumptions for asymptotic analysis may not be valid. In other words, the usual F-ratios will not follow a F-distribution, so we cannot validly undertake hypothesis tests about the regression parameters. The Canova-Hansen Test and the Kwiatkowski-Phillips-Schmidt-Shin Test did not reject the null Hypothesis of stationary seasonal cycles. Applying these tests to the residuals of the fit of *d*6 and *d*7 yields the same results (*p* = 0.369, *d*
_{
ks
} = 0.1290, *n* = 48 and *p* = 0.350, = 25.0, no rejection of the null Hypothesis by Canova-Hansen Test and Kwiatkowski-Phillips-Schmidt-Shin Test).

Though our Fourier approximation underestimates the peaked crime rates around midnight the coefficient of determination of the single days is within [0.86, 0.96]. Overall the model is satisfying.

## Conclusion

Two statistical methods that will enlarge the scientists toolbox for analyzing multi-harmonic oscillations were described. As the example demonstrated the methods can be used to detect and compare multi-harmonic patterns in biological rhythm data.

The orthogonality of the sine and cosine vectors is intensively used to calculate the exact distribution of certain test statistics, not just the approximate distribution for large sample sizes. But this orthogonality also limits the set of frequencies in our multi-harmonic model. In this special case our detection method is an extension of the cosinor-method to multi harmonic models. It also includes a model selection process. Our comparison method uses the whole periodic patterns instead of single parameters. This is an enhancement of the commonly used ANOVA with single parameter „mean“. Furthermore the exact distribution of the test statistic is known, not just an approximate or a limiting distribution for large sample sizes. This can in some cases increase the tests power. In addition the method allows a simultaneous comparison of several time series. This allows to test the hypothesis if „at least one time series shows a different rhythm“ without having any a priori knowledge which one could be deviant (this situation can occur if for example the study design or the data does not allow a partition in a control group and a treatment group).

Problems may occur with missing values (no ON-basis), trends in the data (model is not valid) or the choice of the number of samples, when no a priori knowledge of the inherent periods of the data is available. To derive a more robust version of the statistical test use the rank of the residuals instead the residuals analogous to the ANOVA on ranks. Identifying the method's limitations will help improve it and make it more universal, which is one of the reasons for providing a detailed description of the method calculation steps.

Likelihood ratio tests are in common use for model selection or hypothesis testing and could be an alternative to our tests. Least squares estimates of the coefficients coincide with the maximum likelihood estimates, if the residuals are normal distributed and homoscedastic. Our tests confirm, that the residuals have these properties. So there is neither a gain nor a loss in switching to likelihood ratio tests, which are based on maximum likelihood estimates. Furthermore only the limiting distribution of the likelihood ratio test statistic for large sample sizes is known, whereas the exact distribution of our test statistics is specified. The described detection method uses all harmonic frequencies, because potentially all frequencies could be inherent in the data. However this approach can increase the false negative rate of the test, because the corrected significance level becomes too small. So we are using a conservative test. As Albert and Hunsberger [31] point out there is a „wide range of circadian patterns which can be characterized with a few harmonics“ and that they „recommend choosing between one, two, or three harmonics“. We too found only two significant harmonics in our analysis and observed a good coefficient of determination and white noise residuals. So if some frequencies are ruled out by a priori knowledge the detection method can be executed with fewer harmonics to increase the tests power.

We compared our methods with the cosinor method [26], which fits a single cosine wave with a user defined period to the data: coefficient of determination is 0.732 for a 24 hour period and 0.2 for a 12 hour period when fitting Friday and Saturday. Our detection method achieved a coefficient of determination of 0.86. The cosinor method also calculated the amplitude of the 24 hour periods for workdays and weekends: they differed by only 5%. Analyzing the amplitudes of the first harmonic yields overlapping confidence intervals. So the cosinor method gives no strong evidence for different rhythms on workdays and weekends. A significant difference between workdays and weekends is revealed by simultaneously comparing all weekdays as we did in section.

The findings of a 24 hour period on every day could be for example associated with the hormones testosteron and serotonin. Both of them show a nycthemeral rhythm [7, 8] and are linked to violent behavior [46, 47]. The different rhythm on Friday and Saturday could be caused by exogenous factors like increased alcohol consumption [48].

## References

Manfredini R, Boari B, Smolensky MH, Salmi R, la Cecilia O, Maria Malagoni A, Haus E, Manfredini F:

**Circadian variation in stroke onset: identical temporal pattern in ischemic and hemorrhagic events.***Chronobiol Int*2005,**22**(3)**:**417–53.Pierdomenico SD, Bucci A, Costantini F, Lapenna D, Cuccurullo F, Mezzetti A:

**Circadian blood pressure changes and myocardial ischemia in hypertensive patients with coronary artery disease.***J Am Coll Cardiol*1998,**31**(7)**:**1627–34.Martin RJ:

**Nocturnal asthma: circadian rhythms and therapeutic interventions.***Am Rev Respir Dis*1993,**147**(6 Pt 2)**:**S25–8.Renard E, Palombi K, Gronfier C, Pepin JL, Noel C, Chiquet C, Romanet JP:

**Twenty-four hour (Nyctohemeral) rhythm of intraocular pressure and ocular perfusion pressure in normal-tension glaucoma.***Invest Ophthalmol Vis Sci*2009,**51**(2)**:**882–9.Hara T, Tsuru T:

**Increase of peak intraocular pressure during sleep in reproduced diurnal changes by posture.***Arch Ophthalmol*2006,**124**(2)**:**165–8.Buijs RM, Wortel J, Van Heerikhuize JJ, Feenstra MG, Ter Horst GJ, Romijn HJ, Kalsbeek A:

**Anatomical and functional demonstration of a multisynaptic suprachiasmatic nucleus adrenal (cortex) pathway.***Eur J Neurosci*1999,**11**(5)**:**1535–44.Nakayama K:

**Diurnal rhythm in extracellular levels of 5-hydroxyindoleacetic acid in the medial prefrontal cortex of freely moving rats: an in vivo microdialysis study.***Prog Neuropsychopharmacol Biol Psychiatry*2002,**26**(7–8)**:**1383–8.Diver MJ, Imtiaz KE, Ahmad AM, Vora JP, Fraser WD:

**Diurnal rhythms of serum total, free and bioavailable testosterone and of SHBG in middle-aged men compared with those in young men.***Clin Endocrinol (Oxf)*2003,**58**(6)**:**710–7.Placidi GP, Oquendo MA, Malone KM, Huang YY, Ellis SP, Mann JJ:

**Aggressivity, suicide attempts, and depression: relationship to cerebrospinal fluid monoamine metabolite levels.***Biol Psychiatry*2001,**50**(10)**:**783–91.Lidberg L, Tuck JR, Asberg M, Scalia-Tomba GP, Bertilsson L:

**Homicide, suicide and CSF 5-HIAA.***Acta Psychiatr Scand*1985,**71**(3)**:**230–6.Tripodianakis J, Markianos M, Rouvali O, Istikoglou C:

**Gonadal axis hormones in psychiatric male patients after a suicide attempt.***Eur Arch Psychiatry Clin Neurosci*2007,**257**(3)**:**135–9.Preti A, Miotto P:

**Diurnal variations in suicide by age and gender in Italy.***J Affect Disord*2001,**65**(3)**:**253–61.Caracciolo S, Manfredini R, Gallerani M, Tugnoli S:

**Circadian rhythm of parasuicide in relation to violence of method and concomitant mental disorder.***Acta Psychiatr Scand*1996,**93**(4)**:**252–6.Laubichler W, Ruby M:

**Diurnal relations of criminal behavior.***Arch Kriminol*1986,**177**(5–6)**:**176–84.Artis M, Clavel JG, Hoffmann M, Nachane D:

**Harmonic Regression Models: A Comparative Review with Applications.**[http://ideas.repec.org/p/zur/iewwpx/333.html]*IEW - Working Papers iewwp333, Institute for Empirical Research in Economics, University of Zurich*2007.Brown EN, Czeisler CA:

**The Statistical Analysis of Circadian Phase and Amplitude in Constant-Routine Core-Temperature Data.***J Biol Rhythms*1992,**7**(3)**:**177–202.Celec P:

**Analysis of rhythmic variance-ANORVA. A new simple method for detecting rhythms in biological time series.***Biol Res*2004,**37**(4 Suppl A)**:**777–82.Indic P, Brown EN:

**Characterizing the amplitude dynamics of the human core-temperature circadian rhythm using a stochastic-dynamic model.***J Theor Biol*2006,**239**(4)**:**499–506.Teicher MH, Barber NI:

**COSIFIT: an interactive program for simultaneous multioscillator cosinor analysis of time-series data.***Comput Biomed Res*1990,**23**(3)**:**283–95.Klerman EB, Jewett ME:

**Commentary: model building, quantitative testing, and model comparison.***J Biol Rhythms*1999,**14**(6)**:**621–4.Gaffney M, Taylor C, Cusenza E:

**Harmonic regression analysis of the effect of drug treatment on the diurnal rhythm of blood pressure and angina.***Stat Med*1993,**12**(2)**:**129–142.Douglas AS, Strachan DP, Maxwell JD:

**Seasonality of tuberculosis: the reverse of other respiratory diseases in the UK.***Thorax*1996,**51**(9)**:**944–6.Samel A, Wegmann HM, Vejvoda M, Maass H, Gundel A, Schutz M:

**Influence of melatonin treatment on human circadian rhythmicity before and after a simulated 9-hr time shift.***J Biol Rhythms*1991,**6**(3)**:**235–48.Li XM, Levi F:

**Circadian physiology is a toxicity target of the anticancer drug gemcitabine in mice.***J Biol Rhythms*2007,**22**(2)**:**159–66.Goulet G, Mongrain V, Desrosiers C, Paquet J, Dumont M:

**Daily light exposure in morning-type and evening-type individuals.***J Biol Rhythms*2007,**22**(2)**:**151–8.Nelson W, Tong YL, Lee JK, Halberg F:

**Methods for cosinor-rhythmometry.***Chronobiologia*1979,**6**(4)**:**305–23.Zautcke JL, Morris RW, Koenigsberg M, Carmody T, Stein-Spencer L, Erickson TB:

**Assaults from penetrating trauma in the State of Illinois.***Am J Emerg Med*1998,**16**(6)**:**553–6.Tabata M, Takeshima T, Burioka N, Nomura T, Ishizaki K, Mori N, Kowa H, Nakashima K:

**Cosinor analysis of heart rate variability in ambulatory migraineurs.***Headache*2000,**40**(6)**:**457–63.Munakata M, Imai Y, Minami N, Sasaki S, Ichijyo T, Yoshizawa M, Sekino H, Abe K, Yoshinaga K:

**Cosinor analysis of changes in circadian blood pressure rhythm with aging in spontaneously hypertensive rats.***Tohoku J Exp Med*1990,**161:**55–64.Holm S:

**A Simple Sequentially Rejective Multiple Test Procedure.***Scand Stat Theory Appl*1979,**6**(2)**:**65–70.Albert S, Hunsberger PS:

**On analyzing circadian rhythms data using nonlinear mixed models with harmonic terms.***Biometrics*2005,**61:**1115–1120.Greenhouse JB, Kass RE, Tsay RS:

**Fitting nonlinear models with ARMA errors to biological rhythm data.***Stat Med*1987,**6:**167–183.SJ P:

**Harmonic Analysis Applied to Seasonal Variations in Sickness Absence.***J R Stat Soc Ser C Appl Stat*1974,**23**(2)**:**103–120.Hafner R:

*Wahrscheinlichkeitsrechnung und Statistik*. Springer, Wien; 1989.Biermann T, Asemann R, McAuliffe C, Stroebel A, Keller J, Sperling W, Bleich S, Kornhuber J, Reulbach U:

**Relationship between lunar phases and serious crimes of battery: a population-based study.***Compr Psychiatry*2009,**50**(6)**:**573–7.Massey FJ:

**The Kolmogorov-Smirnov Test for Goodness of Fit.***J Am Stat Assoc*1951,**46**(253)**:**68–78.Lilliefors HW:

**On the Kolmogorov-Smirnov Test for Normality with Mean and Variance Unknown.***J Am Stat Assoc*1967,**62**(318)**:**399–402.Breusch T:

**Testing for Autocorrelation in Dynamic Linear Models.***Aust Econ Pap*1979,**17:**334–355.Godfrey L:

**Testing Against General Autoregressive and Moving Average Error Models when the Regressors Include Lagged Dependent Variables.***Econometrica*1978,**46:**1293–1302.Wald A, Wolfowitz J:

**On a test whether two samples are from the same population.***Ann Math Statist*1940,**11:**147–162.Canova F, Hansen BE:

**Are Seasonal Patterns Constant Over Time? A Test for Seasonal Stability.***J Bus Econ Stat*1995,**13**(3)**:**237–52.Kwiatkowski D, Phillips PCB, Schmidt P, Shin Y:

**Testing the null hypothesis of stationarity against the alternative of a unit root: How sure are we that economic time series have a unit root?***J Econom*1992,**54**(1–3)**:**159–178.Roenneberg T, Kumar CJ, Merrow M:

**The human circadian clock entrains to sun time.***Curr Biol*2007,**17**(2)**:**R44–5.Diebold F, Senhadji : [http://econpapers.repec.org/RePEc:wop:pennhp:_054]

*Deterministic vs. Stochastic Trend in U.S. GNP, Yet Again*. Home pages, University of Pennsylvania; 1996.Cochrane J:

**How Big Is The Random Walk in GNP.***J Polit Econ*1988,**96:**893–920.Manuck SB, Flory JD, McCaffery JM, Matthews KA, Mann JJ, Muldoon MF:

**Aggression, impulsivity, and central nervous system serotonergic responsivity in a nonpatient sample.***Neuropsychopharmacology*1998,**19**(4)**:**287–99.Virkkunen M, Rawlings R, Tokola R, Poland RE, Guidotti A, Nemeroff C, Bissette G, Kalogeras K, Karonen SL, Linnoila M:

**CSF biochemistries, glucose metabolism, and diurnal activity rhythms in alcoholic, violent offenders, fire setters, and healthy volunteers.***Arch Gen Psychiatry*1994,**51:**20–7.Sivarajasingam V, Matthews K, Shepherd J:

**Price of beer and violence-related injury in England and Wales.***Injury*2006,**37**(5)**:**388–94.

## Acknowledgements

This work was supported by the Interdisciplinary Center of Clinical Research (IZKF) at the University hospital of the University of Erlangen-Nuremberg. The authors wish to thank Joanne Eysell for proofreading the manuscript.

## Author information

### Authors and Affiliations

### Corresponding author

## Additional information

### Competing interests

The authors declare that they have no competing interests.

### Authors' contributions

AS contributed to the conception and the design of the study, analyzed the data and drafted the manuscript. UR contributed to the conception and the design of the study. TB acquired the data. IK contributed to the analysis. JK contributed to the intellectual content. TG, MB and all other authors read and approved the final version of the article.

## Rights and permissions

## About this article

### Cite this article

Stroebel, A.M., Bergner, M., Reulbach, U. *et al.* Statistical methods for detecting and comparing periodic data and their application to the nycthemeral rhythm of bodily harm: A population based study.
*J Circad Rhythms* **8**, 10 (2010). https://doi.org/10.1186/1740-3391-8-10

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/1740-3391-8-10