What a Sample Score Represents
A loading describes how a latent factor relates to measured features; a sample score describes where one sample lies along that factor. Taken together, a sample’s scores form coordinates for summarizing its molecular profile. In multi-omics analysis, those coordinates may include factors shared across views and factors that describe variation specific to one view.
A score is not a directly measured molecular quantity. It is an estimated latent coordinate whose meaning comes from the model, the fitted loadings, and the scaling used for the data. A positive score, for example, means movement in the positive direction of the chosen factor coordinate. Which features increase or decrease in that direction is determined by their loadings.
The equation describes the factor contribution to the fitted profile. In a continuous observation model, residuals account for differences between observed and fitted values. For count or binary data, the factor terms may instead form a linear predictor that is mapped to the observation scale by the model’s link function. Thus, a score should be interpreted on the scale of the model component in which it appears, not automatically as a change in observed counts or probabilities.
Meaningful interpretation requires that samples are aligned across views, the fitted factor structure is appropriate, and the data preprocessing and model are adequate for the measurements. Scores and loadings must also be read together: the same score can correspond to different feature patterns in different views because each view has its own loading matrix. The numerical value depends on the fitted coordinate system and its scaling; it is not an absolute biological unit.
Reading Shared and View-Specific Coordinates
A shared score is a coordinate for a direction represented in multiple views. Its estimated value is the same for a sample wherever that shared factor is used, but the resulting feature contributions need not have the same sign or magnitude across views. For example, a positive shared score could contribute positively to some RNA features and negatively to some protein features if the corresponding loadings have opposite signs.
A view-specific score represents a direction included for one view only. A high value on a view-specific RNA score can help explain that sample’s RNA profile, but it is not itself evidence that the same direction is present in protein or metabolite measurements. It is also not a “less important” score by definition: it may capture substantial, reproducible variation in its own view.
Worked Example: Interpreting a Shared Score Across Two Views
Suppose a fitted model has one shared factor. For a hypothetical sample, its shared score is \(z_i=1.5\). The loading for an RNA feature is \(0.80\), and the loading for a protein feature is \(-0.50\). Assume these features were standardized, so fitted contributions are in standard-deviation units.
For the RNA feature, the shared-factor contribution is
For the protein feature, it is
The sample lies in the positive direction of the shared coordinate. Under the fitted representation, that direction contributes positively to the RNA feature and negatively to the protein feature. It would be incorrect to say that the score means “more of both molecules”: the loadings specify how the shared coordinate maps to each feature.
From Coordinates to Fitted Feature Contributions
The fitted profile is useful because it connects compact sample coordinates to the features that were measured. Each feature’s fitted value is the sum of its contributions from shared and view-specific factors. The same decomposition can be used to compare two samples or to ask which latent directions account for a particular fitted feature value.
This expression also gives a direct way to interpret score differences. Holding the fitted loadings fixed, the difference between two samples’ fitted values for a feature is the loading-weighted difference in their scores. This follows by subtracting the fitted-value equation for sample \(a\) from the equation for sample \(b\):
Thus, a score difference alone does not tell how much a specific feature differs. Its contribution depends on that feature’s loading. This is a model-based decomposition of fitted profiles, not a causal explanation of why the samples differ.
Worked Example: Decomposing an RNA Feature
Consider a hypothetical standardized RNA feature with shared loading \(0.60\) and RNA-specific loading \(0.40\). Sample A has shared score \(z_A=0.50\) and RNA-specific score \(u_A=0.20\). Sample B has \(z_B=1.50\) and \(u_B=-0.30\). For this feature, suppose there are no other factors. Calculate the two fitted contributions and their difference.
For sample A, the shared contribution is \(0.60(0.50)=0.30\), and the RNA-specific contribution is \(0.40(0.20)=0.08\). Therefore,
For sample B, the corresponding contributions are \(0.60(1.50)=0.90\) and \(0.40(-0.30)=-0.12\). Therefore,
The fitted difference is \(0.78-0.38=0.40\) standard deviations. Using the score-difference formula gives the same result: \(0.60(1.50-0.50)+0.40(-0.30-0.20)=0.60-0.20=0.40\). The shared direction accounts for a \(0.60\)-unit increase in the fitted feature, partly offset by a \(-0.20\)-unit change in the view-specific contribution.
Worked Example: Identifying What a View-Specific Score Can Explain
Suppose a hypothetical sample has shared score \(z_i=0.8\), RNA-specific score \(u_i^{(RNA)}=1.2\), and protein-specific score \(u_i^{(Protein)}=-0.4\). An RNA feature has shared loading \(0.5\) and RNA-specific loading \(0.3\). A protein feature has shared loading \(-0.2\) and protein-specific loading \(0.7\). All values are on standardized scales.
For the RNA feature, the shared contribution is \(0.5(0.8)=0.40\), and its RNA-specific contribution is \(0.3(1.2)=0.36\), giving a fitted value of \(0.40+0.36=0.76\) standard deviations.
For the protein feature, the shared contribution is \(-0.2(0.8)=-0.16\), and its protein-specific contribution is \(0.7(-0.4)=-0.28\), giving a fitted value of \(-0.16-0.28=-0.44\) standard deviations.
The RNA-specific score contributes to the RNA fitted value, not to the protein fitted value in this model. The protein fitted value instead uses the shared score and the protein-specific score. The shared score is common to both views, but its effects differ because the loadings differ.
Assumptions, Diagnostics, and Alternatives
A score interpretation is only as useful as the representation that produced it. At minimum, the factor model must match samples across views correctly, use an appropriate observation model and preprocessing, and provide an adequate low-dimensional description of the patterns being interpreted. For a linear reconstruction, the score-loading products must correspond to the fitted component on the stated scale. For non-Gaussian observations, interpretation must respect the model’s linear predictor and link rather than treating that predictor as an observed measurement.
Scores are estimated quantities. Their uncertainty may matter, especially for samples with missing measurements, sparse counts, or weak information for a factor. A point estimate should not be mistaken for a perfectly measured biological state. Nor should scores be compared numerically across separately fitted models unless their coordinates have been aligned and their scales made comparable.
Useful diagnostics include inspecting score distributions and outliers, checking whether a sample’s score is driven by a few unusual observations, and comparing fitted profiles with observed data on the appropriate scale. Assess whether samples with similar score patterns have coherent fitted feature contributions across views. If scores change substantially under reasonable changes in preprocessing, model fit, or sample inclusion, interpretations based on small score differences are fragile. A large score is not automatically an error; it is a reason to check the corresponding data and fitted contributions.
Scores also differ from other low-dimensional representations. PCA scores are coordinates in principal directions for a single analyzed matrix (or a specifically constructed joint matrix); they do not inherently separate shared from view-specific structure. NMF represents data through nonnegative components and sample weights, which can be useful for additive patterns but has a different constraint and interpretation. Separate factor models for each view produce coordinates that need not represent common directions. The shared-and-specific factor model makes the distinction explicit, but its score meanings still depend on its fitted assumptions and parameterization.
Common Pitfall: Treating a Score as a Measured Trait
A common mistake is to describe a positive shared score as “high expression across all omics,” or to interpret a view-specific score as if it were a shared biological signal. Both claims skip the loading map. A score indicates a coordinate; a feature’s loading determines whether moving along that coordinate contributes positively or negatively to that feature, and whether the coordinate is shared or view-specific determines where it enters the model.
Another mistake is to infer causation from a score-feature pattern. A factor model summarizes patterns under a statistical representation. Even if a shared factor is associated with features in several views, that alone does not establish a biological mechanism or causal direction. Use score interpretations as descriptions of fitted multivariate structure, and consider other evidence before making biological claims.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- A sample has a shared score of \(2\). For one feature, the shared loading is \(-0.3\). What is the shared-factor contribution to that feature, and what does its sign mean?
- Can an RNA-specific score directly contribute to a protein fitted profile in the shared-and-specific representation defined here? Explain.
- Two samples differ by \(0.5\) on a shared score. A feature’s loading on that factor is \(0.8\). What is the difference in their fitted contributions from that factor?
- Why is it not generally valid to compare score magnitudes from two separately fitted models as if they were on the same scale?
Solutions. (1) The contribution is \((-0.3)(2)=-0.6\) feature-scale units. The negative sign means the positive score contributes in the negative direction for this feature. (2) No. An RNA-specific score enters the RNA profile; the protein profile uses shared scores and protein-specific scores. (3) The fitted contribution difference is \(0.8(0.5)=0.40\) feature-scale units. (4) The fitted coordinate scales and orientations may differ, so the same numerical score need not represent the same amount or direction of latent variation across models.