Contents

RA plot

AnRA plot, like its cousin, the MA plot, is a re-scaled and (45-degree) rotated version of a simple two-dimensional scatter plot of a versus b where a and b are equal-length vectors of positive measurements. This rescaling and rotation allows for better visibility and emphasis of important outliers points that vary between the two measurement conditions.^[1] Essentially it is a plot of the log ratio [R] vs the average log [A] of each pairing of the elements of a and b. Unlike an MA plot, however, because the RA plot takes non-negative integer counts as input, it must employ work-arounds to include mathematically invisible points (such as points where one or both element(s) of the pair is zero).

If we modify our original a (orb) vector via:

${\displaystyle a={\begin{cases}a+\varepsilon ,&{\text{if }}a=0\\a,&{\text{if }}a>0\end{cases}}}$

where

${\displaystyle 0<\varepsilon <0.5}$

then R and A can be defined as:

${\displaystyle R=\log _{2}(a/b)}$

${\displaystyle A={\frac {1}{2}}\times (\log _{2}a+\log _{2}b)}$

R, like M, is plotted on the y-axis and represents a log (fold change) ratio between a and b. A is plotted on the x-axis and represents the average abundance for a coordinate pair. The RA plot provides a quick overview of the distribution and size of a dataset consisting of non-zero counts.

The acronym prefix "R.A." is sometimes pronounced as the one syllable word "ray" because of the plot's strong resemblance to a geometric ray. This characteristic arrow-like shape derives from two key features: on the right at the vector origin, a long asymptotic tail, and on the left (forming the arrow head) two (often dense) patches of condition-unique points.

Because a large portion of the pairs of a and b contain zeros in one or both conditions, they are impossible to plot as-is on a log scale. Other MA plotting functions artificially include these condition-unique points in the plot by spreading them vertically as a "smear" on the left or horizontally as a "rug" at the very top and bottom of the plot. In an RA plot, by contrast, the uniques are included via addition a small epsilon factor (between .1 and .5) which places them in a more statistically appropriate location in the plot.

MA plot with condition-unique and zero points as a "smear" (via the edgeR Bioconductor package)

Another problem with plotting this (or any) type of count data is overplotting which is solved in the RA plot by jittering the points out away from each other but no so far as to merge with other coordinates. The result of this feature is a patchwork-like appearance to the plot that fades away as the A increases.

An RA plot: many points have identical coordinates and are hidden from each other

The caroline CRAN R package contains the only known implementation of an RA plot. However, the meta-transcriptomics "manta" R package provides a wrapper around this RA plot implementation and is used for assessing fold change in transcription of genes (the points) while simultaneously visualizing each gene's taxonomic distributions as individual pie chart points.^[2]

library(caroline)
a <- rnbinom(n=10000, mu=5, size=2)
b <- rnbinom(n=10000, mu=5, size=2)

raPlot(a, b)

• MA plot
• DNA microarray
• Bland–Altman plot