Why Compositional Measurements Need Relative Coordinates
Many omics measurements are recorded as parts of a whole: for example, the proportions of identified molecules in a sample or feature abundances reported after normalization to a fixed total. In such data, the total is constrained, so an increase in one part's proportion can occur because that part increased, because other parts decreased, or both. A proportion alone does not establish an absolute change in the amount of that feature.
Log-ratio transforms express the information in these data as comparisons between parts. They are useful when the scientifically meaningful quantities are relative abundances, rather than the total amount measured. They do not make every omics data set compositional: the interpretation depends on how the measurements were produced and what the scientific question concerns.
The centered log-ratio (CLR) transform compares each part with the geometric mean of all parts. The additive log-ratio (ALR) transform compares each non-reference part with one chosen reference. Both require strictly positive entries, and both use natural logarithms here. A log-ratio is dimensionless: the units cancel in the ratio.
What the Coordinates Mean
A CLR coordinate is positive when its part is larger than the geometric mean of the parts in that sample, and negative when it is smaller. It is not a log of an absolute abundance. An ALR coordinate is positive when its part exceeds the selected reference and negative when it is smaller. For either transform, a difference of coordinates has a direct ratio interpretation.
Thus the difference between two CLR coordinates is the log of the ratio between those parts; the geometric mean cancels. For ALR coordinates with the same reference \(r\), subtracting two coordinates also cancels that reference: \(\operatorname{alr}_r(x)_i-\operatorname{alr}_r(x)_j=\ln(x_i/x_j)\). If the same feature and reference are compared across samples, the difference in that ALR coordinate is the log of the change in the feature-to-reference ratio.
Both transforms ignore a common multiplicative scale. If every part is multiplied by a positive constant \(c\), then \(g(cx)=c\,g(x)\), so \(\operatorname{clr}(cx)_i=\ln(cx_i/(c\,g(x)))=\operatorname{clr}(x)_i\). Also, \(\operatorname{alr}_r(cx)_i=\ln(cx_i/(cx_r))=\operatorname{alr}_r(x)_i\). This scale invariance is appropriate when the total is not part of the intended analysis, but it also means the transformed coordinates cannot recover that total.
Worked Example: Calculating Centered Log-Ratio Coordinates
A hypothetical sample has three positive part measurements \(x=(2,8,32)\) in a common arbitrary abundance unit. Calculate its CLR coordinates and interpret them.
First calculate the geometric mean:
Now divide each part by 8 and take natural logarithms:
The first part is one-quarter of the geometric mean, so its coordinate is \(\ln(1/4)=-1.3863\). The second equals the geometric mean, and the third is four times the geometric mean. The first-to-third log ratio is \(-1.3863-1.3863=-2.7726=\ln(2/32)\), confirming that a difference of CLR coordinates gives the log ratio between those two parts.
Additive Log-Ratios and the Reference Choice
An ALR representation uses \(D-1\) coordinates because the reference part is compared with itself only implicitly: its ratio to itself is 1, with log equal to 0. Choosing a different reference changes the values and labels of the coordinates. It does not change the underlying pairwise ratios, which can be recovered from coordinate differences or from the original composition.
Worked Example: Changing the ALR Reference
A hypothetical three-part sample is \(x=(12,3,6)\), with parts measured in the same arbitrary unit. Calculate the ALR coordinates using part 2 as reference, then part 1. Interpret why the numerical coordinates differ.
With part 2 as reference, the coordinates are:
These say that part 1 is four times part 2 and part 3 is twice part 2. With part 1 as reference, the coordinates for parts 2 and 3 are:
The coordinate values change because the comparisons changed. The underlying ratio between parts 2 and 3 remains available: using part 2 as reference, \(\ln(x_3/x_2)=0.6931\); using part 1 as reference, subtract the part 2 coordinate from the part 3 coordinate to obtain \(-0.6931-(-1.3863)=0.6931\). Both representations therefore agree that part 3 is twice part 2.
A chosen reference should be positive and measured comparably across samples. It should also be scientifically interpretable as a useful denominator. A low or technically unstable reference can make many ALR coordinates noisy, since every one divides by it. A CLR avoids selecting one reference part, but its denominator depends on all included parts; changing the feature set can therefore change CLR coordinates.
Constraints, Reconstruction, and a Short Derivation
The CLR coordinates always sum to zero. This follows directly from the definition:
The final equality holds because \(g(x)^D=\prod_{i=1}^{D}x_i\). Consequently, CLR coordinates lie in a \(D-1\)-dimensional subspace; they are linearly dependent, even though there are \(D\) coordinates. This is a property of the representation, not evidence of a calculation error. Treating all CLR coordinates as independent can lead to singular covariance matrices in methods that require an invertible covariance matrix.
ALR has exactly \(D-1\) coordinates and does not have this same zero-sum constraint. It is invertible as a representation of relative composition once the reference is specified. Set the reference value to 1 and each other value to the exponential of its ALR coordinate, then divide all values by their sum if proportions summing to 1 are needed. For example, if the ALR coordinates relative to part \(r\) are \(y_i\), set \(q_r=1\), \(q_i=e^{y_i}\), and \(p_i=q_i/\sum_{j=1}^{D}q_j\). Then \(p_i/p_r=e^{y_i}=x_i/x_r\), so the reconstructed composition has the same ratios as \(x\).
CLR can also be converted back to relative proportions: exponentiate each coordinate and close the resulting positive values to sum to 1. The zero-sum property implies that the geometric mean of the exponentiated CLR coordinates is 1; closure supplies the desired total. Neither inverse recovers the original absolute total, because that information was removed by scale invariance.
Worked Example: Scale Invariance Across Two Samples
Two hypothetical samples have measurements \(x=(2,8,32)\) and \(y=(20,80,320)\), in the same units and with the same three parts. The second vector is ten times the first. Compare their CLR coordinates.
For \(x\), the geometric mean is \(g(x)=8\), giving CLR coordinates \((-1.3863,0,1.3863)\), as calculated above. For \(y\):
The coordinates are identical because all three measurements changed by the same factor and their ratios did not change. The transform correctly records no relative change. It cannot distinguish these samples' different totals; if that total is scientifically important, it must be retained and analyzed separately rather than inferred from the log-ratios.
Assumptions, Diagnostics, and Alternatives
Before applying either transform, verify that the coordinates answer the scientific question and that the parts define a comparable composition for each sample. In particular, the selected parts should have consistent meanings across samples, be recorded on a comparable scale, and be strictly positive. The transform is mathematically defined for positive measurements, but mathematical eligibility alone does not establish that relative interpretation is scientifically appropriate.
- Check positivity. A zero or negative entry cannot be used in a logarithmic ratio as defined here. Do not silently replace it with a small positive value; handling zeros and detection limits requires a deliberate method.
- Check the parts and their measurement scale. Inconsistent feature definitions, units that differ across parts, or different sets of included parts can undermine comparisons. A common multiplicative factor applied to every part cancels, but part-specific calibration differences do not generally cancel.
- For ALR, inspect the reference. If the reference is near zero, unstable, or affected by detection limits, ratios involving it can become extreme or highly variable. Assess whether conclusions are sensitive to plausible reference choices.
- For CLR, record the feature set. Because the geometric mean includes every included part, adding or removing features can shift coordinates. Use a consistent feature set when comparing samples and document it.
The main alternatives answer different questions. An ALR is compact and can give a direct, interpretable comparison to a chosen reference, but conclusions expressed as individual coordinates depend on that choice. A CLR treats parts symmetrically and avoids selecting one denominator, but its coordinates are constrained to sum to zero and depend on the included parts. A direct pairwise log ratio, \(\ln(x_i/x_j)\), focuses on a particular comparison without defining a full transformed vector. By contrast, applying a log transform to each raw measurement separately does not itself remove a sample-wide scale: relative interpretation requires comparing logged parts or using an appropriate normalization strategy. The earlier material on normalizing proteomic and metabolomic intensities addresses intensity normalization; a log-ratio transform is not a substitute for deciding whether the measurements should be interpreted compositionally.
Common Mistakes and Practical Use
- Reading a log-ratio as an absolute effect. A higher CLR coordinate means the part is larger relative to the sample's geometric mean, not necessarily that its absolute abundance increased.
- Assuming the total is preserved. Multiplying every part by the same positive factor leaves both transforms unchanged. Keep separate total information when it matters.
- Calling the CLR coordinates independent. They sum to zero, so there is an exact linear dependence among the \(D\) coordinates.
- Ignoring an ALR reference change. Coordinate values from different references are not directly interchangeable. State the reference and interpret every coordinate as a comparison with it.
- Applying logarithms to zeros without a plan. The formulas require positive entries; an undefined value is not repaired by ignoring the issue.
A practical choice is to use CLR when a symmetric representation of the specified set of parts is useful and downstream methods can accommodate its linear dependence. Use ALR when a meaningful, stable reference supports clear comparisons and reference dependence is acceptable. In either case, report the parts included, the reference if applicable, and whether the measurements are treated as relative composition. These details determine what the resulting coordinates mean.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- For \(x=(1,4,16)\), calculate the geometric mean and all three CLR coordinates.
- For \(x=(5,10,20)\), calculate the ALR coordinates using part 3 as the reference. What does the first coordinate mean?
- Why do CLR coordinates sum to zero, and what is the practical implication of that constraint?
- If every part in a sample is multiplied by 7, what happens to its CLR and ALR coordinates? Explain.
- Name one diagnostic concern specific to ALR and one that applies to both CLR and ALR.
Solutions. (1) \(g(x)=(1\cdot4\cdot16)^{1/3}=4\). The coordinates are \(\ln(1/4)=-1.3863\), \(\ln(4/4)=0\), and \(\ln(16/4)=1.3863\). (2) \(\operatorname{alr}_3(x)=(\ln(5/20),\ln(10/20))=(-1.3863,-0.6931)\). The first coordinate says part 1 is one-quarter of part 3. (3) The sum is \(\ln(\prod_i x_i/g(x)^D)=\ln(1)=0\), because \(g(x)^D=\prod_i x_i\). The coordinates are linearly dependent, so a covariance matrix using all \(D\) coordinates is singular. (4) Neither changes: the common factor cancels from every ratio and from the CLR geometric-mean ratio. (5) For ALR, a reference that is unstable or near zero can make all ratios involving it unstable. For both transforms, entries must be positive; zeros make the logarithmic ratios undefined.