Raleigh plots, or Rayleigh plots (also called circlegrams and closely related to circular histograms,[1] phasor diagrams,[2] and wind roses[3]), are statistical graphics that serve as graphical representations for a Raleigh test that map a mean vector to a circular plot. Raleigh plots have many applications in the field of chronobiology, such as in studying butterfly migration patterns or protein and gene expression, and in other fields such as geology, cognitive psychology, and physics.
Raleigh plots was first introduced by Lord Rayleigh. The concept of Raleigh plots evolved from Raleigh tests, also introduced by Lord Rayleigh in 1880. The Rayleigh test is a popular statistical test used to measure the concentration of data points around a circle, identifying any unimodal bias in the distribution.[4] Rayleigh plots emerged from this analysis as a means to illustrate the nature of the distribution.
In a Raleigh plot, each individual is assigned a unit vector with a corresponding angle. These unit vectors are averaged together in a Raleigh plot into the mean vector. The length of the mean vector is determined by r (or R), the mean resultant length. r is a measure of concentration, ranging in value between 0 and 1. If the individual angles of the unit vectors are tightly clustered, then the r value will be closer to 1 while if they are widely scattered, then r will be closer to zero.[5] The mean vector is positioned in the center of a circle. Dashes along the circumference of this circle denote desired values. Examples include angles from magnetic north (zero degrees) going clockwise (e.g., 90 degrees from magnetic north, or eastward); times of day, which can also be described in zeitgeber time and circadian time; and phase. Dots on the circumference are usually used to indicate individual unit vectors and their respective angle in regard to the values being measured. Raleigh plots can also use more than one mean vector, particularly if one wants to display the mean vector for different tested groups in the study or to compare mean vectors between groups.
The example to the left is a Raleigh plot that has a high r value. Blue and yellow dots indicate individuals from different groups being tested, and the position of the dots indicate in which angle from magnetic north each tested individual is traveling. Due to the high overall concentration of individuals going at an angle between zero degrees and ninety degrees, the mean vector is much longer. Compare the figure to the left with the figure below:
In this second example (to the right), the Raleigh plot has a low r value. Both yellow and blue dots are spread along the circumference of the circle, indicating many individuals are traveling at different angles. The largest cluster of individuals, a group traveling between 180 degrees and 270 degrees, causes the mean vector to be pointed at an angle in that direction. Notably, due to the variability in the direction within this group, the mean vector is much smaller.
Raleigh plots have been used in chronobiology studies on the biological clocks behind monarch butterfly migration patterns. They are particularly relevant for studying sun compass orientation in migrating butterflies.
In butterfly migration studies, the Raleigh plot maps the orientation of the butterfly when allowed to fly, where the circumference is marked as a compass, with north (N) at the top position. Given the plotted data points, a mean r vector is drawn to indicate the mean orientation of the butterflies in a particular condition.[6] [7]
In his studies on the neurobiology of butterfly migration, Steven M. Reppert observes the oriented flight behavior of monarch butterflies. Reppert explains how Raleigh plots are used to handle butterfly orientation data and as tools for the data analysis.
In a 2012 study by Reppert and colleagues on the sufficiency of an antenna for proper time compensation and sun compass orientation in the monarch butterfly, Raleigh plots were used to present the mean flight orientation of butterflies subjected to different study conditions. Along the edge of the circle, degrees 0º (magnetic north) through 360º are shown, the orientation of each butterfly is marked with a dot, and a mean vector is drawn to represent the mean flight orientation recorded.[8] In a 2016 annual review on the neurobiology of monarch butterfly migration, Reppert, Guerra, and Merlin also use Raleigh plots to present butterfly orientation data. The plots were used in the study of the time-compensated sun compass in monarch butterflies.[9] A 2018 review by Reppert and de Roode on the mechanisms of monarch butterfly migration also used Raleigh plots, or circlegrams, to represent butterfly orientation data. Each dot indicates the orientation in which a butterfly individual flew continuously for 5 minutes or longer, and the vector points in the mean direction (some degrees from north) with a magnitude proportional to the mean orientation.[10]
Raleigh plots can be used to visualize circadian rhythms in protein or gene expression, and how their phases are affected by other variables or induced conditions.
Jennifer Mohawk, then a postdoc at the University of Texas Southwestern Medical Center, used multiple Raleigh plots to illustrate PER2::LUC expression in her 2013 paper "Methamphetamine and Dopamine Receptor D1 Regulate Entrainment of Murine Circadian Oscillators."[11] Specifically, Mohawk investigated how injections of methamphetamine and D1 antagonist SCH-23390 would shift the peak time of PER2 expression in the liver, lung, pituitary gland, and salivary gland. In these plots, the Raleigh plot can be interpreted as a 24 hour clock with CT 0 at the top of the circle and CT 12 at the bottom of the circle. Each arrow represents the average peak phase of PER2::LUC expression of each group. The strength of the phase clustering is symbolized by the length of the arrow, meaning stronger clustering or closer data points resulting in longer arrows. The individual data points are plotted on the outside of the circle and their unique color and shape resemble the different groups of conditions. At ZT 7 is a pink box that shows the timing of the methamphetamine injection. Mohawk and collaborators compared the angle of the vector, or the mean phase, between the different groups in order to determine if methamphetamine injections induced statistically significant phase changes of PER2::LUC expression within the different glands and organs.
Similarly, Tsedey Mekbib, then a PhD student at the Morehouse School of Medicine, utilized Raleigh plots to depict how the knockdown of SIAH2 impacted the rhythmic expression of all other genes in both males and females. After profiling the entire transcriptome via RNA sequence on the liver at frequent time intervals, the expression of peak timing for all rhythmically expressed genes was plotted on Raleigh plots for each group. These Raleigh plots contain vectors that represent the average peak phase. However, instead of differing the length of the vectors to illustrate the variability in the data points, Mekbib and the other collaborators added a +/- 95% confidence interval that is represented by red range along the circle. In addition to the typical Raleigh plot, the left half of the circle is shaded darker to better visualize the night phase occurring between CT 12 and CT 24.[12]
Raleigh plots can be used to visualize the circadian rhythms in protein or gene expression in the suprachiasmatic nucleus (SCN) specifically, and how they might influence other peripheral tissues.
Raleigh plots were used by Elizabeth Maywood, an English researcher at the MRC Laboratory of Molecular Biology, to visualize how pacemaking activity and synchrony between host SCN cells lacking vasoactive intestinal peptide (VIP) can be restored with a wild-type SCN graft. These plots show vectors that represent the phase of the host SCN cells, measured by PER2::LUC expression. Each plot has a value representing the mean vector length with time points where the cell phases are closer in sync having a value closer to 1. Maywood and collaborators showed that VIP-null host SCN cells synchrony deteriorated over time based on the mean vector length of the Raleigh plots decreasing, and concluded that paracrine signaling from an introduced wild-type SCN graft is sufficient to restore the synchrony between SCN cell pacemakers based on the mean vector length increasing after the graft was introduced.[13]
Mariko Izumo, while working at the University of Texas Southwestern Medical Center, used Raleigh plots to assess the effect of knocking out BMAL1 in the SCN on the circadian expression of PER2, measured using PER2::LUC expression, in the SCN and in peripheral tissues. Izumo and others found that knocking out BMAL1 led to the desynchronization and dampening of PER2 expression in peripheral tissues with Raleigh plots showing different mean phases of the rhythm of PER2 expression. They also show that light/dark cycles and feeding can restore synchrony in peripheral tissues with Raleigh plots showing the same mean phases of the rhythm of PER2 expression.[14]
Similarly, Yongli Shan, also while working at the University of Texas Southwestern Medical Center, used Raleigh plots to show that BMAL1 knockout AVP SCN neurons and VIP SCN neurons show a loss in rhythm in expression of PER2. BMAL1 knockout neurons have data points around the circular diagram and a small mean vector r, while wild-type neurons have a larger mean vector r and data points closer together. Additionally, Shan and others show that intercellular connections to the rest of the SCN was sufficient in restoring rhythmicity in BMAL1 knockout VIP neurons but not AVP neurons with an increase in the mean vector r in VIP neurons but not in AVP neurons.[15]
Raleigh plots or variations on Raleigh plots are used in fields beyond chronobiology. While Raleigh plots visualize a mean vector for data, the variational plots that are closely related to Raleigh plots may visualize histogram data in spokes on a circular chart.
In geology, circular histogram plots or rose diagrams can be used to characterize tectonic plate movements. For example, they may be used to visualize the frequency and direction of fault line motion.[16] In meteorology and climate studies, wind roses are used to present data on the direction, duration, and speed of winds that occurred at a given location. For example, in wind roses released by the Midwestern Regional Climate Center, the length of a spoke in a particular direction, representing a histogram bin, is proportional to the duration of time for which the wind was blowing in that direction, with different colors to show wind speed categories.
Additionally, Raleigh plots can be used in cognitive psychology. Joëlle Provasi, professor at École pratique des hautes études, used them to explain the response of children with or without a lesion in their cerebellum due to surgically removed tumor to a rhythmic stimulus. Provasi and others show a Raleigh plot depicting responses that are close to the stimulus as data points around the top at 0 degrees with a mean vector value close to 1, and a Raleigh plot depicting responses that are irregular with dots spread around the circular plot and a mean vector value closer to 0.[17]