---
title: 'Consensus MDS: Aggregating Ideal Points'
url: https://www.emergentmind.com/topics/consensus-multidimensional-scaling-comds
type: topic
---

# Consensus MDS: Aggregating Ideal Points

Consensus Multidimensional Scaling (CoMDS) is a distance-based aggregation method that takes in multiple, heterogeneous estimates of political actors’ “ideal points” and returns a single, common embedding that preserves the cross-source structure those estimates share. It is designed for the way applied scholars actually use ideal points in practice: as general indicators of ideological position or extremity, irrespective of the particular data source or behavioral model that generated them. In that sense, CoMDS captures the shared, stable associations of a set of underlying ideal point estimates and can be interpreted as their common spatial representation. In a broader multi-view setting, closely related work formulates the same general problem as multidimensional scaling on multiple input distance matrices, seeking a single consensus embedding from several dissimilarity views [2601.05213] [1605.00286].

## 1. Problem setting and motivation

CoMDS arises from a setting in which researchers have many estimates of candidates’ positions derived from different data sources and models, including roll-call votes, campaign finance, speeches, platforms, social media, endorsements, and surveys. These sources often exhibit weak within-party relationships and strong cross-party separation. Each source mixes the candidate’s general positioning, which is the shared component, with domain-specific idiosyncrasies, which are source-specific components. Different measures are nevertheless used interchangeably in most substantive analyses, even though weak relationships raise questions about the extent to which they capture a shared quantity rather than idiosyncratic, domain-specific factors [2601.05213].

A central motivation is endogeneity. Using an ideal point and an outcome measured from the same domain creates endogeneity concerns through mechanical correlations and measurement co-production. CoMDS addresses this by extracting only the structure shared across sources and delivering a common spatial representation of candidates. In practice, it recovers a consensus embedding that is invariant to each source’s rotation, scale, and shift and tolerant to missingness; it enables relating consensus ideal points to source-domain variables without relying on any single domain’s idiosyncrasies; and it isolates source-specific features via a projection residual, helping clarify what conclusions hinge on the domain rather than the underlying ideology [2601.05213].

The same general problem can be stated outside political measurement. In the multi-view MDS formulation, the inputs are multiple distance or dissimilarity matrices over the same set of objects, and the target is a single consensus embedding shared across all views. This broader framing treats each input distance matrix as one view and asks how to do multidimensional scaling on multiple input distance matrices in a consensus sense [1605.00286].

## 2. Formal model and objective functions

Let $C$ be the set of $n$ candidates, and let $S$ denote the number of sources, indexed by $s=1,\dots,S$. For each source $s$, one may observe either a one-dimensional estimate $\theta_i^{(s)}$ for candidate $i$ or a higher-dimensional embedding $z_i^{(s)} \in \mathbb{R}^{r_s}$. From each source, one constructs a pairwise dissimilarity matrix $D^{(s)}$ with entries $d_{ij}^{(s)} \ge 0$. If $\theta^{(s)}$ is one-dimensional, then
$$
d_{ij}^{(s)} = |\theta_i^{(s)} - \theta_j^{(s)}|.
$$
If $Z^{(s)}$ is multi-dimensional, then
$$
d_{ij}^{(s)} = \|z_i^{(s)} - z_j^{(s)}\|_2
$$
after source-specific standardization [2601.05213].

The CoMDS target is a consensus embedding $Y \in \mathbb{R}^{n \times k}$, with $k=1$ or $2$ in most applications, that preserves these source-specific pairwise dissimilarities up to source-specific diagonal rescalings. Let $v_{ij}^{(s)} \in \{0,1\}$ be a missingness or availability indicator for pair $(i,j)$ in source $s$, and let $w_s \ge 0$ be a source-level weight. The metric CoMDS stress objective is
$$
\min_{Y \in \mathbb{R}^{n \times k}, \{W^{(s)} \in \mathrm{Diag}(k)\}}
\sum_{s=1}^S w_s \sum_{i<j} v_{ij}^{(s)}
\left[
\|W^{(s)}(Y_i - Y_j)\|_2 - d_{ij}^{(s)}
\right]^2.
$$
Here $W^{(s)}$ is a $k \times k$ diagonal matrix allowing each source to rescale the consensus axes, so sources can have different effective units across dimensions [2601.05213].

A non-metric CoMDS variant replaces metric matching with source-specific monotone transformations $f_s$:
$$
\mathrm{Stress}_{\mathrm{nonmetric}}(Y,\{f_s\}) =
\sum_{s=1}^S w_s \sum_{i<j} v_{ij}^{(s)}
\left[
f_s(\|Y_i - Y_j\|_2) - \delta_{ij}^{(s)}
\right]^2,
$$
where $\delta_{ij}^{(s)}$ is a source-provided proximity or dissimilarity and $f_s$ is estimated via isotonic regression mapping consensus distances to source dissimilarities. This variant preserves rank orders rather than metric scales within sources [2601.05213].

Construction of $d_{ij}^{(s)}$ is part of the model definition. For one-dimensional $\theta^{(s)}$, the source is centered and rescaled to unit variance or unit interquartile range before taking absolute differences. For multi-dimensional $Z^{(s)}$, each dimension is standardized within source to unit variance, and Euclidean distances are then used. If one begins from similarities, they may be converted to dissimilarities by $\delta_{ij}^{(s)} = 1 - \mathrm{sim}_{ij}^{(s)}$ or another monotone transformation, followed by non-metric CoMDS [2601.05213].

Identifiability follows standard MDS invariances. MDS-type embeddings are invariant to translation, rotation, and global scale. CoMDS enforces identifiability by centering the columns of $Y$, fixing average squared distance to $1$, and orienting axes for interpretability. For $k=1$, the sign may be set so larger values denote greater conservatism; for $k>1$, Procrustes alignment to a reference or anchoring with known ideologues may be used. Sign flips or rotations across runs with different initializations do not affect distances [2601.05213].

## 3. Estimation, optimization, and computational structure

CoMDS can be estimated by an extension of SMACOF, an iterative majorization algorithm for MDS that guarantees monotone stress decrease and linear convergence. A practical scheme is alternating minimization. In Step A, one updates $Y$ given the current $\{W^{(s)}\}$ using a SMACOF-style update. In Step B, one updates each diagonal $W^{(s)}$ given $Y$ via a weighted nonnegative least squares fit on squared distances. The alternation typically converges in tens of iterations [2601.05213].

For the $Y$-update, define the current fitted distances
$$
\Delta_{ij}^{(s)}(Y) = \|W^{(s)}(Y_i - Y_j)\|_2.
$$
Then construct per-source majorization matrices $B^{(s)}(Y)$ with entries
$$
b_{ij}^{(s)}(Y) =
\begin{cases}
- \dfrac{w_s v_{ij}^{(s)} d_{ij}^{(s)}}{\Delta_{ij}^{(s)}(Y) + \epsilon}
& \text{if } i \ne j \text{ and } \Delta_{ij}^{(s)}(Y) > 0,\\[6pt]
- \sum_{j \ne i} b_{ij}^{(s)}(Y)
& \text{if } i=j.
\end{cases}
$$
Here $\epsilon>0$ is a small constant to avoid division by zero; if $\Delta_{ij}^{(s)}=0$, one sets $b_{ij}^{(s)}=0$. With
$$
B(Y)=\sum_{s=1}^S B^{(s)}(Y), \qquad
V=\mathrm{diag}(V_1,\dots,V_n),
$$
where
$$
V_i = \sum_{s=1}^S w_s \sum_{j \ne i} v_{ij}^{(s)},
$$
the SMACOF update is
$$
Y^{(t+1)} = V^{-1} B(Y^{(t)}) Y^{(t)}.
$$
After the update, one re-centers the columns of $Y$ to mean zero and rescales to fix the global scale; if using $k>1$, one may optionally orthogonalize columns to improve numerical stability [2601.05213].

For the diagonal scaling update, write $W^{(s)}=\mathrm{diag}(w_1^{(s)},\dots,w_k^{(s)})$. The squared fitted distance is
$$
\|W^{(s)}(Y_i - Y_j)\|_2^2
=
\sum_{d=1}^k (w_d^{(s)})^2 (Y_{i,d}-Y_{j,d})^2.
$$
Defining $x_{ij}^{(d)}=(Y_{i,d}-Y_{j,d})^2$ and $\beta_d^{(s)}=(w_d^{(s)})^2$, one fits $\beta^{(s)}$ by weighted nonnegative least squares:
$$
\min_{\beta^{(s)} \in \mathbb{R}_+^k}
\sum_{i<j} v_{ij}^{(s)}
\left[
\sum_{d=1}^k \beta_d^{(s)} x_{ij}^{(d)} - (d_{ij}^{(s)})^2
\right]^2.
$$
Then $w_d^{(s)} = \sqrt{\beta_d^{(s)}}$. This subproblem is convex and can be solved efficiently, including by projected gradient or NNLS [2601.05213].

Initialization and computational details are explicit. One may initialize $Y$ using classical MDS on the average distance matrix $\bar D$ with
$$
\bar d_{ij}
=
\frac{\sum_s w_s v_{ij}^{(s)} d_{ij}^{(s)}}{\sum_s w_s v_{ij}^{(s)}}.
$$
One may initialize $W^{(s)}$ to identity or solve the NNLS once given $Y^0$. A standard stopping rule is relative stress decrease less than $10^{-6}$ or reaching the maximum number of iterations. Complexity per iteration is $O(S n^2 k)$ for $B(Y)$ and $O(S n^2 k)$ for NNLS, often faster in practice. Missingness is handled through $v_{ij}^{(s)}$, so candidates missing from a source are still embedded via other sources [2601.05213].

## 4. Diagnostics, validation, and treatment of endogeneity

CoMDS is accompanied by a suite of diagnostic tools intended to aid practical usage. A basic diagnostic is per-source fit. Define the source-specific sum of squared errors
$$
\mathrm{SSE}_s
=
\sum_{i<j} v_{ij}^{(s)}
\left[
\|W^{(s)}(Y_i-Y_j)\|_2 - d_{ij}^{(s)}
\right]^2
$$
and
$$
\mathrm{SST}_s
=
\sum_{i<j} v_{ij}^{(s)}
\left[
d_{ij}^{(s)} - \mu_s
\right]^2,
$$
with $\mu_s$ the mean of the observed $\{d_{ij}^{(s)}\}$. Then the per-source $R^2$-type fit is
$$
R_s^2 = 1 - \frac{\mathrm{SSE}_s}{\mathrm{SST}_s}.
$$
A second diagnostic is the relative error contribution
$$
\mathrm{RelError}(Z^{(s)}) =
\frac{\mathrm{SSE}_s}{\sum_{r=1}^S \mathrm{SSE}_r},
$$
which lies in $[0,1]$ and sums to $1$ across sources; larger values indicate more source-specific idiosyncrasy and less shared structure [2601.05213].

Leave-one-source-out stability assesses how much the consensus depends on any one source. One estimates CoMDS with all sources to obtain $Y$, then re-estimates after leaving out source $s$ to obtain $Y_{-s}$. For $k=1$, the similarity metric is
$$
\rho(s)=|\mathrm{cor}(Y,Y_{-s})|.
$$
For $k>1$, subspace correlation is defined by
$$
\rho(s)=\frac{1}{k}\sum_{d=1}^k d_i^2,
$$
where the $\{d_i\}$ are singular values of $\mathrm{ortho}(Y)^T \mathrm{ortho}(Y_{-s})$. Values near $1$ indicate stability. Dimension selection proceeds by plotting stress versus $k$ or total $R^2$ versus $k$ and looking for elbows, or by cross-validating on held-out pairs or sources and measuring predictive $\mathrm{SSE}_s$. In political applications, $k=1$ or $k=2$ usually suffices [2601.05213].

Candidate influence and uncertainty can also be examined. Influence may be assessed by recomputing stress after deleting candidate $i$ or by computing the leave-one-candidate-out correlation between full $Y$ and $Y_{-i}$. Uncertainty may be assessed by bootstrapping pairs within sources or bootstrapping sources, then reporting standard errors or confidence bands for $Y_i$. Residual maps based on
$$
r_{ij}^{(s)}=
\|W^{(s)}(Y_i-Y_j)\|_2 - d_{ij}^{(s)}
$$
can reveal systematic residual patterns suggestive of source-specific biases, including agenda or measurement artifacts [2601.05213].

Endogeneity is a principal use-case rather than a secondary consideration. CoMDS mitigates endogeneity by purging source-specific co-production. If one uses campaign-finance-based ideal points and relates them to campaign-finance outcomes, estimates may be mechanically inflated. The consensus embedding is trained to preserve what multiple sources agree on, and source-specific distortions are downweighted through the multi-source stress. Recommended practice includes split-source analyses, in which one re-estimates CoMDS excluding the outcome’s source or assigns that source a small weight, and projection decomposition, in which outcomes are regressed on the shared component only [2601.05213].

Projection decomposition makes the shared-versus-idiosyncratic distinction explicit. After obtaining $Y$, define
$$
P^* = Y(Y^T Y)^{-1}Y^T,
$$
the projection onto the consensus subspace. For source $s$, the shared component of the original coordinates $Z^{(s)}$ is $P^* Z^{(s)}$, and the idiosyncratic residual is
$$
\hat{\nu}^{(s)} = Z^{(s)} - P^* Z^{(s)} = (I-P^*)Z^{(s)}.
$$
By construction, $\hat{\nu}^{(s)}$ is orthogonal to $Y$. This tool quantifies how much a source’s measure deviates from the consensus and can be related to external domain variables to understand source-specific biases [2601.05213].

## 5. Empirical applications in congressional ideal-point estimation

In the U.S. House application for 2016–2024, the method combines NOMINATE, campaign finance scores, and platform positions. The sources are described as NOMINATE, which is roll-call based, one- to two-dimensional, and available for incumbents only; campaign finance scores, including static, dynamic, and DW-DIME variants; and platform positions, based on Wordfish-scaled campaign website issue positions. The analysis includes all candidates with at least two of the three sources, comprising approximately $5{,}389$ candidate-year observations. A one-dimensional consensus embedding displays strong partisan bimodality, and within-party correlations with the consensus are stronger than cross-source correlations alone [2601.05213].

The paper reports three substantive findings. First, on roll-call partisan disloyalty, consensus points reveal that moderates defect more, whereas NOMINATE-only analysis exaggerates magnitudes because of common-source endogeneity. Comparing source-specific measures yields conflicting signs or magnitudes, while the consensus resolves the relationship. Second, on fundraising base, measured as the number of unique donors, the consensus indicates that Democrats’ extremism and Republicans’ moderation have distinct effects; campaign finance scores alone suggest opposite conclusions and much larger magnitudes. Third, on lexical diversity, measured by platform CTTR, the consensus reveals that moderation associates with greater rhetorical sophistication within party; platform-based ideal points alone overstate magnitudes, and other sources sometimes disagree in sign [2601.05213].

These applications define the concrete scope of CoMDS in the political science setting. The method is intended for cases in which ideal points are used as domain-agnostic measures of ideology or extremity across heterogeneous sources, especially when the researcher plans to study outcomes measured in any one domain, including votes, donors, speech, or platforms. A plausible implication is that the method is particularly useful when the substantive target is a stable ideological component rather than a domain-bound behavioral signature, because the projection residual can then be used to study agenda control, rhetoric, donor targeting, or other source-specific phenomena orthogonal to the underlying ideology [2601.05213].

## 6. Relation to multi-view MDS, alternative methods, limitations, and extensions

A broader consensus-MDS formulation is given by Multi-View Multidimensional Scaling (MVMDS), which takes $V=M \ge 2$ distance matrices $D^{(v)}=[\delta_{ij}^{(v)}] \in \mathbb{R}^{n \times n}$ over the same set of objects and learns a single consensus embedding $X \in \mathbb{R}^{n \times p}$ together with view weights $\alpha^{(v)}$. Its objective is
$$
\min_{\alpha^{(v)},X}
\sum_{v=1}^M (\alpha^{(v)})^\gamma
\sum_{i<j} w_{ij}\big(\delta_{ij}^{(v)}-d_{ij}(X)\big)^2
$$
subject to
$$
\sum_{v=1}^M \alpha^{(v)} = 1,
\qquad
0 \le \alpha^{(v)} \le 1,
$$
where $\gamma>1$ is a weight controller. Given $X$, the weights admit the closed-form update
$$
\alpha^{(v)}
=
\frac{(J^{(v)})^{1/(1-\gamma)}}
{\sum_{v'=1}^M (J^{(v')})^{1/(1-\gamma)}},
$$
with
$$
J^{(v)}=\sum_{i<j} w_{ij}\big(\delta_{ij}^{(v)}-d_{ij}(X)\big)^2.
$$
The $X$-step uses a SMACOF-style majorization update $X = V^+BZ$, and with no missing entries the update simplifies to
$$
X =
\frac{1}{N\sum_{v=1}^M (\alpha^{(v)})^\gamma} BZ.
$$
This formulation makes consensus and complementarity explicit: views with lower per-view stress receive higher weight, but $\gamma>1$ prevents trivial collapse to a single view [1605.00286].

The political-science CoMDS and the multi-view MVMDS formulation share a common structural idea: a single embedding is estimated so that multiple dissimilarity structures are jointly respected. CoMDS differs by emphasizing source-specific diagonal rescalings, tolerance to missingness, diagnostic decompositions, and endogeneity mitigation in applications involving heterogeneous ideal point estimators. MVMDS, by contrast, is strictly metric MDS in the paper’s formulation and does not use non-metric monotone transformations. This suggests that CoMDS occupies a more application-specific position within a broader family of consensus MDS methods [2601.05213] [1605.00286].

Alternative methods have distinct assumptions and failure modes.

| Method | Pros | Cons |
|---|---|---|
| Bayesian IRT joint scaling | coherent probabilistic inference | misspecification risk; heavy computation; stringent functional-form assumptions |
| MD2S | explicit decomposition | assumes normality and linearity; sensitive to rotations/scales |
| PCA on stacked measures | simple and fast | sensitive to source rotations/scales; overweights sources with more dimensions |
| Joint or bridge alignments | useful with a gold standard | requires anchors or overlap and a trusted reference |
| JIVE and related consensus methods | shares CoMDS’ spirit | typically assumes linear decompositions |

CoMDS has three stated advantages. It is invariant to source rotations, scales, and shifts because it operates on distances and learns source-specific diagonal scales. It is flexible with missingness and balances across sources rather than over-weighting high-dimensional ones. It does not impose a parametric behavioral model and therefore respects each source’s domain-specific generative assumptions. Its stated limitations are equally clear: it uses only pairwise dissimilarities rather than raw data likelihoods and has no built-in probabilistic uncertainty; the choice of distance and weights matters and can become problematic if a source is very noisy or adversarial; and the algorithm converges to local minima, so initialization matters, even though SMACOF tends to be stable in practice [2601.05213].

Several extensions are proposed. Dynamic CoMDS estimates time-indexed consensus embeddings $Y^{(t)}$ with temporal smoothness penalties, for example
$$
\min \sum_t \mathrm{Stress}_t(Y^{(t)})
+ \lambda \sum_t \|Y^{(t)}-Y^{(t-1)}\|_F^2.
$$
Hierarchical or multi-level CoMDS introduces group-level embeddings and penalties for nested actors. Incorporating uncertainty can be done by weighting pair contributions by inverse variance or adjusting source weights to reflect reliability. Probabilistic formulations treat observed distances as noisy observations of $\|W^{(s)}(Y_i-Y_j)\|$ with additive Gaussian noise and estimate by maximum likelihood or Bayesian posterior sampling while retaining SMACOF-like majorization for MAP. Theoretical guarantees are local rather than global: the SMACOF majorization step guarantees non-increasing stress and linear convergence to a local minimum, and the NNLS updates for $W^{(s)}$ are convex and converge to a global minimum for each fixed $Y$ [2601.05213].

Source: https://www.emergentmind.com/topics/consensus-multidimensional-scaling-comds