Reported Drosophila courtship song rhythms are artifacts of data analysis
© Stern; licensee BioMed Central Ltd. 2014
Received: 13 February 2014
Accepted: 20 May 2014
Published: 26 June 2014
In a series of landmark papers, Kyriacou, Hall, and colleagues reported that the average inter-pulse interval of Drosophila melanogaster male courtship song varies rhythmically (KH cycles), that the period gene controls this rhythm, and that evolution of the period gene determines species differences in the rhythm’s frequency. Several groups failed to recover KH cycles, but this may have resulted from differences in recording chamber size.
Here, using recording chambers of the same dimensions as used by Kyriacou and Hall, I found no compelling evidence for KH cycles at any frequency. By replicating the data analysis procedures employed by Kyriacou and Hall, I found that two factors - data binned into 10-second intervals and short recordings - imposed non-significant periodicity in the frequency range reported for KH cycles. Randomized data showed similar patterns.
All of the results related to KH cycles are likely to be artifacts of binning data from short songs. Reported genotypic differences in KH cycles cannot be explained by this artifact and may have resulted from the use of small sample sizes and/or from the exclusion of samples that did not exhibit song rhythms.
In 1980, Kyriacou and Hall  reported that, on average, the inter-pulse interval cycled with an amplitude of approximately 3 ms and a period of approximately 1 min (0.017 Hz) in D. melanogaster (Figure 1d) and 30 s (0.03 Hz) in D. simulans (Figure 1e). Following tradition, I will call this periodicity KH cycles. In the same paper, they reported that mutations in the period locus that increased, decreased, and ablated periodicity of the circadian rhythm had similar effects on KH cycles. Later papers reported that transformation of the period locus into individuals carrying a null allele of the period gene restores KH cycles to arrhythmic males ; that females of D. melanogaster and D. simulans mate faster when presented with artificial songs that contain the species-specific inter-pulse interval and KH cycles than when presented with heterospecific songs ; that the difference in KH cycles between D. melanogaster and D. simulans is controlled by a locus or loci on the X chromosome (where the period gene resides) ; that hybrid females of D. melanogaster and D. simulans prefer songs with a combination of intermediate inter-pulse interval and KH cycles ; that D. simulans KH cycles can be transferred into D. melanogaster by genetic transformation of the D. simulans period locus into a D. melanogaster strain carrying a period null mutation ; and that this effect of the period locus can be mapped to a 700 bp region of the largest exon of period, which contains a threonine-glycine repeat region that differs in length between D. simulans and D. melanogaster.
This series of papers has been heavily cited and discussed in the secondary literature, in part because this is one of only two examples of a species-specific behavioral difference that has been traced to a genetic difference. (Variation in vasopressin V1A receptor expression likely contributes to differences in affiliative behavior between vole species [11, 12]). Given the importance of this work, it would be of great interest to follow-up on these initial studies and to identify the molecular mode of action of period in courtship song cycles. However, no other genes have been implicated in control of KH cycles and the molecular mechanism by which variation in the period locus influences KH cycles has not been identified. In fact, there have been no papers published since 1991 providing further evidence for the mode of action of the period gene in courtship song rhythms. In particular, there is still no satisfactory explanation for how the period gene, whose expression cycles approximately diurnally and whose protein product is degraded during daylight hours , could contribute to rhythms that cycle on the order of 30 to 60 s. This fact contrasts strongly with the significant progress over the past 20 years revealing the molecular basis for circadian rhythms [14, 15].
There are two possible reasons for the recent absence of published work on this problem. First, the tools available over the past 20 years may have been insufficient to allow more detailed investigation. For example, tools to allow high throughput analysis of courtship song and detailed analysis of neural circuits in Drosophila have only recently become available [3, 16]. Second, it is possible that KH cycles are not real phenomena, but instead are an artifact of the methods of data collection or analysis utilized in earlier papers. This possibility was raised originally in the late 1980s, when two independent groups reported that they were unable to detect KH cycles [17, 18]. (Also see [19–22]). In addition, my colleagues and I were unable to detect significant evidence for KH cycles in five wild-type strains of D. melanogaster. However, one independent group has reported that they can detect KH cycles . As statisticians have noted, it is likely that in most of the early studies too few pulse events were sampled to allow robust estimates of rhythms with periods on the order of 60 s . Nonetheless, there has been no reasonable explanation for the fact that Kyriacou and Hall and at least one other group have consistently found evidence for KH cycles.
In apparent further support for the existence of Kyriacou and Hall rhythms, Ritchie and colleagues showed that female D. melanogaster prefer synthesized courtship song containing KH cycles over song that does not contain this rhythmicity . This may imply that females have the ability to detect KH cycles in courtship song. However, there is an alternative explanation for their observations. Normal D. melanogaster song contains complex, non-random dynamics of the inter-pulse interval on the scale of seconds  (Figure 1b). Ritchie et al. compared synthesized song that contained KH cycles with synthesized song containing either randomly patterned inter-pulse intervals, constant inter-pulse intervals, or silence. Remarkably, females showed very similar preferences for randomly patterned song, constant song, and silence, even though females are known to prefer song over silence . It is possible that the randomly patterned song and constant inter-pulse interval song in this experiment did not accurately mimic wild-type song. Most importantly, this experiment is missing an additional critical control, song with natural patterns of inter-pulse intervals. Without this control, it is premature to conclude that females prefer song containing KH cycles to song without cycles. It is possible that D. melanogaster females prefer some modulation of inter-pulse interval over song with a constant inter-pulse interval — since variation in the inter-pulse interval is a common feature of real song — but this has not yet been demonstrated.
Since KH cycles are a fascinating phenomenon with implications for speciation and evolutionary biology, I had hoped to follow-up on earlier reports to study the molecular and neurobiological mechanisms regulating KH cycles. Despite efforts to match earlier experimental conditions as closely as possible, I was unable to detect KH cycles. As I show, however, non-significant rhythms of the inter-pulse interval in the frequency range reported by Kyriacou and Hall are generated by binning the data of short songs. These results suggest that KH cycles do not exist and that previous reports reflect artifacts of data analysis.
I searched for KH cycles first by attempting to replicate the experimental paradigm used by Kyriacou and Hall as closely as possible. Unlike many earlier attempts, I employed recording chambers with dimensions that matched those used by Kyriacou and Hall. I recorded song from wild-type D. melanogaster and D. simulans flies and mutants of the period locus (null, short and long). I recorded each courting male for approximately 45 min and, to avoid possible artifacts arising from limited sample sizes, I discarded recordings that included fewer than 1,000 pulses.
Testing for rhythmicity in time series of biological data presents multiple challenges [27, 28]. One classical method for detecting rhythmicity is to calculate the Fourier transform of the time series, which essentially transforms the time series into the frequency domain by representing the original wave form as the sum of a series of sinusoids of different frequencies, power, and phase. Power can then be plotted for each frequency in the form of a periodogram and periodicity can be detected as peaks of strong power. However, the Fourier transform requires equally spaced data. With unequally spaced data, or data sampled regularly but with missing values, the Fourier transform does not have well-defined properties and can generate biased periodograms [27, 29]. Furthermore, the statistical significance of values in these periodograms cannot be determined. Everyone who has worked on KH cycles was aware of, and made efforts to ameliorate, these limitations of classical periodogram analysis [23, 30].
To overcome these difficulties, astronomers developed methods based on the least-squares fit of a series of sinusoids of different frequencies to the data [29, 31]. This method, now called the Lomb-Scargle periodogram, does not require regularly sampled time series and, because the power at different frequencies represents the contribution of that frequency to the variance in the signal, it allows estimation of the statistical significance of the power at a particular frequency [29, 31]. The Lomb-Scargle periodogram is now widely recognized as a powerful method for detecting rhythmicity in biological time series [27, 28, 32–34]. In addition, we have shown previously through simulation  that the Lomb-Scargle periodogram is very sensitive to potential KH cycles, with power greater than 80% to detect putative KH cycles when the signal-to-noise ratio is at least one (Additional file 3: Figure S1). Therefore, I used the Lomb-Scargle periodogram to search for rhythmicity in the inter-pulse interval.
KH cycles observed by previous authors might have resulted from the specific analysis techniques or sampling procedures employed in earlier studies. In different papers, authors used different spectral analysis techniques to search for rhythms, so it is unlikely that KH cycles resulted from any specific analytical technique. However, there are two parameters that are common to most published studies: inter-pulse interval data were binned into 10 s intervals and most recordings were approximately 5 min long or shorter (Additional file 1: Table S1). In addition, most studies did not incorporate a robust method of assessing the statistical significance of rhythmicity . Thus, often the periodogram peak with greatest power was taken as the best estimate of rhythmicity. Sampling just the peaks representing the greatest power in the Lomb-Scargle periodograms did not obviously alter the distribution of frequencies from the sample of all peaks with P < 0.01 (Figure 3c). However, 10 s binning of inter-pulse intervals dramatically altered the distribution of strongest signals in the Lomb-Scargle periodograms, forcing many values to be within the frequency range reported for KH cycles (Figure 3d, compare with colored columns), although few of these signals were statistically significant. This truncation of the high frequency range of the spectrogram resulted simply from the fact that the highest frequency that can be estimated in any time series is equal to half the sampling rate, or 0.05 Hz (1/20 s) for 10 s bins .
To more precisely replicate the data used in previous studies, I analyzed the first five minutes of song following the first train of pulse song. I also sampled non-overlapping 5 min clips from across the entire duration of recordings. Both samples of 5 min of song resulted in poor estimation of low frequency rhythms, forcing approximately 65% of the maximal rhythms to be within the range reported for KH cycles (Figure 3d,e). This is similar to the proportion of songs with KH cycles that was reported in the later publications from Kyriacou and Hall (Figure 2a), suggesting that these reported KH rhythms could be fully explained as artifacts of binning short songs. In addition to these manipulations, I also tested many other manipulations of the data described in , including multiple methods for estimating values for empty bins. None of these other manipulations had a large effect on the distributions illustrated in Figure 4 (not shown).
The pattern observed for the wild-type Canton-S strain of D. melanogaster was also observed for a wild-type strain of D. simulans and for the three period mutations (Figure 4a). Binning short songs constrained the maximal rhythmicity close to the range reported by Kyriacou and Hall (Figure 4a). It is important to note that binning the data, but not shortening the songs, reduced the power of this rhythmicity, resulting in much lower statistical significance for the maximal periodogram peaks (Figure 4a). In fact, fewer than two percent of the maximal periodogram peaks observed in binned data reached a threshold for statistical significance of P <0.05. That is, binning did not generate statistically significant rhythms, it simply constrained the maximal peak of periodograms so that they were more likely to fall in the range reported by Kyriacou and Hall. Since most previous studies did not incorporate robust methods of estimating the statistical significance of the maximal peaks, it is likely that their reported maximal peaks were, like the ones I found, not statistically significant.Randomizing inter-pulse interval data produced patterns of maximal periodogram peaks similar to the original data, but the statistical significance of the peaks was much lower on average in the randomized data than in the original data (Figure 4b). This provides further evidence that the statistically significant high frequency rhythms detected in the original data reflect real fluctuations in inter-pulse intervals on the order of seconds.
While it may be impossible to prove that a phenomenon does not exist, the foundations of science rest upon repeatable observations. KH cycles cannot be replicated. Perhaps more importantly, we now have an explanation for why Kyriacou, Hall and colleagues observed courtship song rhythms in a particular frequency range. Periodogram peaks with periods on the order of 60 s are an artifact of binning data from short songs into 10 s intervals. Most of this apparent rhythmicity does not reach any sensible level of statistical significance. KH cycles appear to represent fluctuations of noise, which is supported by their similarity to patterns generated by randomized data.
My observations do not explain why Kyriacou and Hall observed that different genotypes appeared to produce KH cycles with slightly different average periods nor why females of different species appeared to prefer KH cycles with periods matching that of their own species . One possibility is that the observed average differences between genotypes represented random samples from the full distribution of rhythms generated artificially by the analysis techniques, and that by using relatively small sample sizes (N = 8 ± 3.9 (mean ± standard deviation); Additional file 2: Table S2), Kyriacou and Hall fortuitously observed frequencies consistent with their hypotheses. This possibility is emphasized by the extensive overlap seen in the reconstructed ranges of periods reported in earlier studies (Figure 2b). This small-sample-size effect may have been exacerbated by the exclusion of many songs that did not exhibit rhythmicity (Figure 2a, Additional file 2: Table S2), which was common practice in all papers on this subject.
For decades, the work of Kyriacou and Hall has appeared far ahead of its time and their evidence for the role of the period gene in behavioral evolution has stood as a shining example of how genetic evolution can generate a species-specific behavior. Unfortunately, it appears highly likely that this example is not true. This case highlights some of the challenges inherent in identifying the genetic causes of behavioral evolution. Behavior is often variable on multiple time scales and as a result of multiple genetic and environmental influences. Scoring behavior and detecting species-specific differences requires objective methods for quantifying behavior. Given that most behaviors are likely to have evolved as the result of evolution at multiple loci , large samples of individuals must be assayed. Satisfactory methodologies have only recently become available for scoring a handful of fly behaviors [3, 37–40]. There is now at least some hope that these questions can be addressed in a rigorous manner.
Virgin males and females were collected from the Canton-S strain of D. melanogaster and from a wild-type strain of D. simulans collected in Princeton, NJ, USA in 2009. Individual males were isolated in single wells of a 96-well Scienceware™ 96 deep-well plate (Bel-Art, Wayne, NJ, USA) containing approximately 1 mL of standard cornmeal-based fly food covered with an AirPore Tape Sheet (Qiagen, Hilden, Germany). Males were aged for 5 to 10 days in a 25°C incubator with a 12 hour:12 hour light:dark cycle. Females were group housed in standard Drosophila vials and maintained in the same incubator with males for 1 to 3 days prior to pairing with males for recording.
We built 20 × 10 × 3 mm recording chambers from acrylic, which matched the specifications reported in  and which could be fitted in our 32-channel recording system . Data were collected and segmented as reported previously .
Following Kyriacou and Hall, inter-pulse intervals were defined as times between pulse events that were between 15 and 65 ms long. I excluded songs with fewer than 1,000 pulses and time series were initiated at the pulse train that contained at least 10 inter-pulse intervals, as described in . Data were analyzed using the Lomb-Scargle periodogram with the MATLAB script ‘lomb’ written by Saragiotis and available from . Data were binned by calculating the mean inter-pulse interval for 10 s intervals that each contained at least 10 pulses, as described in . Five-minute clips were generated by dividing 45 min recordings into 5 min segments. For the randomization test, all times between pulses were calculated (including the pauses between pulse trains), and these values were shuffled without replacement. This procedure generated time series that resembled real song, but both the initiation time of pulse trains and the inter-pulse intervals were shuffled. All other analyses were performed using custom MATLAB scripts. The raw song data are available at . The segmented inter-pulse intervals for all individuals are available as individual CSV files at . Complete details of the analysis pipeline, including step-by-step instructions to replicate the results reported here, are provided in a file called AnalysisPipeline.rtf together with all required MATLAB scripts in the zipped folder available at .
I thank Elizabeth Kim for building the new recording chambers and for assistance collecting and recording flies; Michael Young for providing flies carrying alleles of the period locus; and Justin Crocker, Yun Ding, Andrew Lemire, Joshua Mast, Ella Preger, Troy Shirangi and Ben Arthur for helpful discussions. I thank Ben Arthur, Jessica Cande, Roian Egnor, Mala Murthy, Troy Shirangi, and three reviewers for helpful comments on the manuscript.
- Ewing AW, Bennet-Clark HC: The courtship songs of Drosophila. Behaviour. 1968, 31: 288-301. 10.1163/156853968X00298.View ArticleGoogle Scholar
- Greenspan RJ, Ferveur JF: Courtship in Drosophila. Annu Rev Genet. 2000, 34: 205-232. 10.1146/annurev.genet.34.1.205.PubMedView ArticleGoogle Scholar
- Arthur BJ, Sunayama-Morita T, Coen P, Murthy M, Stern DL: Multi-channel acoustic recording and automated analysis of Drosophila courtship songs. BMC Biol. 2013, 11: 11-10.1186/1741-7007-11-11.PubMed CentralPubMedView ArticleGoogle Scholar
- Kawanashi M, Watanabe TK: Genetic variation of courtship song of Drosophila melanogaster and D. simulans. Japan J Genet. 1980, 55: 235-240. 10.1266/jjg.55.235.View ArticleGoogle Scholar
- Turner TL, Miller PM: Investigating natural variation in Drosophila courtship song by the evolve and resequence approach. Genetics. 2012, 191: 633-642. 10.1534/genetics.112.139337.PubMed CentralPubMedView ArticleGoogle Scholar
- Kyriacou CP, Hall JC: Circadian rhythm mutations in Drosophila melanogaster affect short-term fluctuations in the male’s courtship song. PNAS. 1980, 77: 6729-6733. 10.1073/pnas.77.11.6729.PubMed CentralPubMedView ArticleGoogle Scholar
- Zehring WA, Wheeler DA, Reddy P, Konopka RJ, Kyriacou CP, Rosbash M, Hall JC: P-element transformation with period locus DNA restores rhythmicity to mutant, arrhythmic Drosophila melanogaster. Cell. 1984, 39: 369-376. 10.1016/0092-8674(84)90015-1.PubMedView ArticleGoogle Scholar
- Kyriacou CPP, Hall JC: The function of courtship song rhythms in Drosophila. Anim Behav. 1982, 30: 794-801. 10.1016/S0003-3472(82)80152-8.View ArticleGoogle Scholar
- Kyriacou CP, Hall JC: Interspecific genetic control of courtship song production and reception in Drosophila. Science. 1986, 232: 494-497. 10.1126/science.3083506.PubMedView ArticleGoogle Scholar
- Wheeler D, Kyriacou CP, Greenacre ML, Yu Q, Rutila JE, Rosbash M, Hall JC: Molecular transfer of a species-specific behavior from Drosophila simulans to Drosophila melanogaster. Science. 1991, 251: 1082-1085. 10.1126/science.1900131.PubMedView ArticleGoogle Scholar
- Young LJ, Nilsen R, Waymire KG, MacGregor GR, Insel TR: Increased affiliative response to vasopressin in mice expressing the V1a receptor from a monogamous vole. Nature. 1999, 400: 766-768. 10.1038/23475.PubMedView ArticleGoogle Scholar
- Lim MM, Wang Z, Olazabal DE, Ren X, Terwilliger EF, Young LJ: Enhanced partner preference in a promiscuous species by manipulating the expression of a single gene. Nature. 2004, 429: 754-757. 10.1038/nature02539.PubMedView ArticleGoogle Scholar
- Young MW, Kay SA: Time zones: a comparative genetics of circadian clocks. Nat Rev Genet. 2001, 2: 702-715. 10.1038/35088576.PubMedView ArticleGoogle Scholar
- Nitabach MN, Taghert PH: Organization of the Drosophila circadian control circuit. Curr Biol. 2008, 18: R84-R93. 10.1016/j.cub.2007.11.061.PubMedView ArticleGoogle Scholar
- Hardin PE: The circadian timekeeping system of Drosophila. Curr Biol. 2005, 15: R714-R722. 10.1016/j.cub.2005.08.019.PubMedView ArticleGoogle Scholar
- Jenett A, Rubin GM, Ngo TT, Shepherd D, Murphy C, Dionne H, Pfeiffer BD, Cavallaro A, Hall D, Jeter J, Iyer N, Fetter D, Hausenfluck JH, Peng H, Trautman ET, Svirskas RR, Myers EW, Iwinski ZR, Aso Y, DePasquale GM, Enos A, Hulamm P, Lam SC, Li HH, Laverty TR, Long F, Qu L, Murphy SD, Rokicki K, Safford T, et al: A GAL4-driver line resource for Drosophila neurobiology. Cell Rep. 2012, 2: 991-1001. 10.1016/j.celrep.2012.09.011.PubMed CentralPubMedView ArticleGoogle Scholar
- Ewing AW: Cycles in the courtship song of male Drosophila melanogaster have not been detected. Anim Behav. 1988, 36: 1091-1097. 10.1016/S0003-3472(88)80068-X.View ArticleGoogle Scholar
- Crossley SA: Failure to confirm rhythms in Drosophila courtship song. Anim Behav. 1988, 36: 1098-1109. 10.1016/S0003-3472(88)80069-1.View ArticleGoogle Scholar
- Ewing AW: Identification of cycles in the courtship song of Drosophila melanogaster - comment. Anim Behav. 1989, 37: 860-861.View ArticleGoogle Scholar
- Kyriacou CP, Hall JC: Failure to detect cycles in Drosophila mating songs - comment. Anim Behav. 1988, 36: 1110-10.1016/S0003-3472(88)80070-8.View ArticleGoogle Scholar
- Crossley SA: On Kyriacou and Hall's defense of courtship song rhythms in Drosophila. Anim Behav. 1989, 37: 861-863.View ArticleGoogle Scholar
- Kyriacou CP, Hall JC: Spectral analysis of Drosophila courtship song rhythms. Anim Behav. 1989, 37: 850-859.View ArticleGoogle Scholar
- Alt S, Ringo J, Talyn B, Bray W, Dowse H: The period gene controls courtship song cycles in Drosophila melanogaster. Anim Behav. 1998, 56: 87-97. 10.1006/anbe.1998.0743.PubMedView ArticleGoogle Scholar
- Logan IG, Rosenberg J: A referee’s comment on the identification of cycles in the courtship song of Drosophila melanogaster. Anim Behav. 1989, 37: 860-View ArticleGoogle Scholar
- Ritchie MG, Halsey EJ, Gleason JM: Drosophila song as a species-specific mating signal and the behavioural importance of Kyriacou & Hall cycles in D. melanogaster song. Anim Behav. 1999, 58: 649-657. 10.1006/anbe.1999.1167.PubMedView ArticleGoogle Scholar
- Bennet-Clark HC, Ewing AW: Pulse interval as a critical parameter in courtship song of Drosophila melanogaster. Anim Behav. 1969, 17: 755-10.1016/S0003-3472(69)80023-0.View ArticleGoogle Scholar
- Ruf T: The Lomb-Scargle periodogram in biological rhythm research: analysis of incomplete and unequally spaced time-series. Biol Rhythm Res. 1999, 30: 178-201. 10.1076/brhm.184.108.40.2062.View ArticleGoogle Scholar
- Van Dongen HP, Olofsen E, VanHartevelt JH, Kruyt EW: A procedure of multiple period searching in unequally spaced time-series with the Lomb-Scargle method. Biol Rhythm Res. 1999, 30: 149-177. 10.1076/brhm.220.127.116.114.PubMedView ArticleGoogle Scholar
- Lomb NR: Least-squares frequency analysis of unequally spaced data. Astrophys Space Sci. 1976, 39: 447-462. 10.1007/BF00648343.View ArticleGoogle Scholar
- Kyriacou CP, van den Berg MJ, Hall JC: Drosophila courtship song cycles in normal and period mutant males revisited. Behav Genet. 1990, 20: 617-644. 10.1007/BF01065875.PubMedView ArticleGoogle Scholar
- Scargle JD: Studies in astronomical time series analysis. II - Statistical aspects of spectral analysis of unevenly spaced data. Astrophys Journal, Part 1. 1982, 263: 835-853.View ArticleGoogle Scholar
- Glynn EF, Chen J, Mushegian AR: Detecting periodic patterns in unevenly spaced gene expression time series using Lomb-Scargle periodograms. Bioinformatics. 2006, 22: 310-316. 10.1093/bioinformatics/bti789.PubMedView ArticleGoogle Scholar
- Bohn A, Hinderlich S, Hütt M-T, Kaiser F, Lüttge U: Identification of rhythmic subsystems in the circadian cycle of crassulacean acid metabolism under thermoperiodic perturbations. Biol Chem. 2003, 384: 721-728.PubMedView ArticleGoogle Scholar
- Deckard A, Anafi RC, Hogenesch JB, Haase SB, Harer J: Design and analysis of large-scale biological rhythm studies: a comparison of algorithms for detecting periodic signals in biological data. Bioinformatics. 2013, 29: 3174-3180. 10.1093/bioinformatics/btt541.PubMed CentralPubMedView ArticleGoogle Scholar
- Nyquist H: Certain topics in telegraph transmission theory. Trans AIEE. 1928, 47: 617-644.Google Scholar
- Robinson GE, Fernald RD, Clayton DF: Genes and social behavior. Science. 2008, 322: 896-900. 10.1126/science.1159277.PubMed CentralPubMedView ArticleGoogle Scholar
- Dankert H, Wang L, Hoopfer ED, Anderson DJ, Perona P: Automated monitoring and analysis of social behavior in Drosophila. Nat Methods. 2009, 6: 297-303. 10.1038/nmeth.1310.PubMed CentralPubMedView ArticleGoogle Scholar
- Kabra M, Robie AA, Rivera-alba M, Branson S, Branson K: JAABA: interactive machine learning for automatic annotation of animal behavior. Nat Methods. 2013, 10: 64-70.PubMedView ArticleGoogle Scholar
- Mendes CS, Bartos I, Akay T, Márka S, Mann RS: Quantification of gait parameters in freely walking wild type and sensory deprived Drosophila melanogaster. Elife. 2013, 2: e00231-10.7554/eLife.00231.PubMed CentralPubMedGoogle Scholar
- Giancardo L, Sona D, Huang H, Sannino S, Managò F, Scheggia D, Papaleo F, Murino V: Automatic visual tracking and social behaviour analysis with multiple mice. PLoS One. 2013, 8: e74557-10.1371/journal.pone.0074557.PubMed CentralPubMedView ArticleGoogle Scholar
- Lomb normalized periodogram. [http://www.mathworks.com/matlabcentral/fileexchange/22215-lomb-normalized-periodogram]
- Raw data from this paper. [http://research.janelia.org/sternlab/rawData.tar]
- Segmented inter-pulse interval data from this paper. [http://research.janelia.org/sternlab/IPIs.zip]
- MATLAB analysis scripts used in this paper. [http://research.janelia.org/sternlab/KHSongAnalysis.zip]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.