2 A Statistical Model for Roll Call Analysis with Non-Separable Preferences

In this section, we describe our statistical model for roll call votes with non-separable preferences. The theoretical foundation of ideal point estimation are spatial voting models that define actors’ utility functions with stochastic disturbances and functional forms for estimation (see, e.g., Clinton et al. Reference Clinton, Jackman and Rivers2004). Statistical models of roll call analysis assume that legislators’ preferences across dimensions are additively separable, so that an actor’s utility is given by the (weighted) sum of deviations along dimensions. Non-separability also allows for the interaction of deviations across dimensions and is part of the canonical weighted Euclidean distance model of spatial voting (Davis et al. Reference Davis, Hinich and Ordeshook1970; Hinich and Munger Reference Hinich and Munger1997).

Budget negotiations are the classical example where non-separability matters in multidimensional decision-making (Hinich and Munger Reference Hinich and Munger1997). Because a single budget limits total spending, preferences over individual issues are rendered non-separable: more spending in one area must be offset by less spending in another area. Overarching budgetary constraints induce quid pro quo trade-offs across issues.

It is likely that non-separability also matters in other settings of multidimensional decision-making (Milyo Reference Milyo2000). In models of roll call votes, proposals can have implications on multiple dimensions representing conflict in distinct policy areas (e.g., as a consequence of issue bundles across areas; Kadane Reference Kadane1972; Kardasheva Reference Kardasheva2013). Therefore, legislators generally take deviations from their unconditional ideal point across dimensions into account when casting their votes. The question is whether these deviations enter actors’ utility functions solely in an additively separable manner or whether joint cross-dimensional interactions of deviations (as implied by non-separability) also figure in their utility functions.

What could induce such a relationship across dimensions? In an electoral choice example, Stoetzer and Zittlau (Reference Stoetzer and Zittlau2020) argue that distinct issue dimensions are related at a broader ideological level in voters’ decisions (see also Stoetzer and Zittlau Reference Stoetzer and Zittlau2015). Therefore, voters may favor a party representing a specific combination of deviations from their unconditional ideal point to an equally distant party representing a different combination of deviations. Similarly, Finke (Reference Finke2009) argues that preferences over distinct areas are non-separable in political negotiations because the final outcome is an aggregated result of competition between negotiators over each area. Negotiators favor certain aggregate directions of an outcome to others as they give and take across areas, resembling the quid pro quo of budgetary negotiations. Comparable rationales could apply to decision-making in legislatures when actors consider the aggregate direction implied by the combination of deviations.

Hence, when proposals have implications on dimensions and these are related at an aggregate level in actors’ considerations, it is theoretically plausible to consider non-separability in their decision-making processes. To incorporate non-separability in ideal point estimation, we assume that the utility of actors follows a weighted quadratic Euclidean distance model (Davis et al. Reference Davis, Hinich and Ordeshook1970; Hinich and Munger Reference Hinich and Munger1997). The spatial utility function for legislator i with ideal point $\boldsymbol {\theta }_{i}$ regarding a proposal j with outcome $\mathbf {p}_{j}$ in a D-dimensional space (both $\boldsymbol {\theta }_{i}, \mathbf {p}_{j} \in \mathbb {R}^{D}$ ) is assumed to be

(1) $$ \begin{align} v_{i}\left(\mathbf{p}_{j}\right) = -[\mathbf{p}_{j}-\boldsymbol{\theta}_{i}]^{\prime}\mathbf{A}[\mathbf{p}_{j}-\boldsymbol{\theta}_{i}]. \end{align} $$

$\mathbf {A}$ is a $D \times D$ weighting matrix that is assumed to be symmetric positive definite. The diagonal entries of $\mathbf {A}$ are dimensional salience weights that reflect the relevance of distances on a particular dimension in a legislator’s utility function. The off-diagonal entries of $\mathbf {A}$ contain non-separability parameters and reflect whether and how distances on one dimension are related to distances on any other dimension in the utility function. If an off-diagonal entry equals zero, preferences across the corresponding pair of dimensions are separable, otherwise they are non-separable.

The derivation of a statistical model for roll call votes from this spatial utility function assumes that legislators compare the utility of a Yea-vote to the utility of a Nay-vote and vote accordingly (Clinton et al. Reference Clinton, Jackman and Rivers2004). We illustrate the implications of non-separability with a two-dimensional example that shows the unconditional ideal point of a legislator $\boldsymbol {\theta }_{i}$ and her indifference contours in Figure 1. We denote $\mathbf {Y}_{j}$ as the Yea position and $\mathbf {N}_{j}$ as the Nay position of a bill in $\mathbb {R}^{D}$ . In both panels, $\mathbf {Y}_{j}$ and $\mathbf {N}_{j}$ are equally distant from $\boldsymbol {\theta }_{i}$ . While $\mathbf {Y}_{j}$ represents deviations in the same direction across dimensions, $\mathbf {N}_{j}$ combines deviations in opposing directions.

Figure 1 Illustration in two dimensions for a proposal with $\mathbf {Y}_{j} = [1,1]$ and $\mathbf {N}_{j} = [1,-1]$ .

In the left panel of Figure 1, the legislator’s preferences are separable and she is indifferent between $\mathbf {Y}_{j}$ and $\mathbf {N}_{j}$ because her utility equally decreases in all directions away from her ideal point. In the right panel, her preferences are non-separable and the off-diagonal elements of $\mathbf {A}$ are positive. Therefore, her utility decreases faster when deviations in the same direction across dimensions are combined than when deviations in opposing directions are combined in a substitutionary manner. $\mathbf {N}_{j}$ is favored over $\mathbf {Y}_{j}$ because of the aggregate direction implied by the combination of deviations across dimensions. Finally, if non-separability were negative, the contours would be flipped across either axis. Then, $\mathbf {Y}_{j}$ would be preferred to $\mathbf {N}_{j}$ as it combines deviations in equal directions across dimensions in a complementary manner.

While this illustration clarifies the mathematical working of non-separability in utility functions, its substantial interpretation depends on the meaning and rotation of the underlying dimensions. For example, the two dimensions may represent conflict on economic and social issues with lower values indicating more left-wing outcomes and higher values indicating more right-wing outcomes on both dimensions. Then, a substitutionary relationship would imply that actors favor outcomes combining more left-wing deviations on one issue dimension with more right-wing deviations on the other dimension, while a complementary relationship favors deviations in the same direction across dimensions. If either the meaning or the rotation of a dimension were altered, the same non-separability parameter would represent an altered combination of favored deviations. Crucially, non-separability can only be meaningfully interpreted relative to the latent dimensions.

Stochastic aspects influence legislators’ decisions beside the systematic component given by the utilities $v_{i}\left (\mathbf {Y}_{j}\right )$ and $v_{i}\left (\mathbf {N}_{j}\right )$ . Following Clinton et al. (Reference Clinton, Jackman and Rivers2004), we assume that errors for the Yea utility $\eta _{ij}$ and the Nay utility $\nu _{ij}$ are jointly normal distributed with $E[\eta _{ij}] = E[\nu _{ij}] = 0$ . We further assume that $Var[\eta _{ij} - \nu _{ij}] = 1$ . Given that the overall utility is $u_{i}\left (\mathbf {Y}_{j}\right ) = v_{i}\left (\mathbf {Y}_{j}\right ) + \eta _{ij}$ and $u_{i}\left (\mathbf {N}_{j}\right ) = v_{i}\left (\mathbf {N}_{j}\right ) + \nu _{ij}$ , legislator i votes for a proposal $y_{ij} =1$ when $u_{i}(Y_{j})> u_{i}(N_{j})$ and opposes a proposal $y_{ij} =0$ otherwise. This results in a probabilistic choice model in which $\Phi (.)$ is the cumulative density function of the standard normal distribution:

(2) $$ \begin{align}P(y_{ij}=1) &= P\left(v_{i}\left(\mathbf{Y}_{j}\right) + \eta_{ij}> v_{i}\left(\mathbf{N}_{j}\right) + \nu_{ij}\right)\nonumber\\[1pt] &= P\left(v_{i}(\mathbf{Y}_{j}) - v_{i}(\mathbf{N}_{j}) > \nu_{ij} - \eta_{ij}\right)\nonumber\\[1pt] &= \Phi\left(v_{i}(\mathbf{Y}_{j}) - v_{i}(\mathbf{N}_{j})\right)\\[1pt] &= \Phi\left(-[\mathbf{Y}_{j}-\boldsymbol{\theta}_{i}]^{\prime}\mathbf{A}[\mathbf{Y}_{j}-\boldsymbol{\theta}_{i}] + [\mathbf{N}_{j}-\boldsymbol{\theta}_{i}]^{\prime}\mathbf{A}[\mathbf{N}_{j}-\boldsymbol{\theta}_{i}]\right)\nonumber\\[1pt] & = \Phi\left( 2 (\mathbf{Y}_{j} - \mathbf{N}_{j}) ^{\prime}\mathbf{A} \boldsymbol{\theta}_{i} -\mathbf{Y}_{j}^{\prime}\mathbf{A} \mathbf{Y}_{j} +\mathbf{N}_{j}^{\prime}\mathbf{A} \mathbf{N}_{j}\right) .\nonumber \end{align} $$

We can define a multidimensional item response theory (IRT) model based on this specification. The first term of the equation is a function of the Yea and Nay positions of the proposal, the weight matrix, and the legislator’s ideal points. Collecting the proposal-specific terms gives a vector of discrimination parameters  $\boldsymbol {\beta }_{j} = 2 (\mathbf {Y}_{j} - \mathbf {N}_{j})$ . The second part of the equation depends on the proposal-specific parameters and the weight matrix, which we denote as the difficulty parameter  $\alpha _{j} = -\mathbf {Y}_{j}^{\prime }\mathbf {A}\mathbf {Y}_{j} +\mathbf {N}_{j}^{\prime }\mathbf {A}\mathbf {N}_{j}$ , resulting in an IRT model:

(3) $$ \begin{align} P(y_{ij}=1) = \Phi\left(\boldsymbol{\beta}_{j}^{\prime}\mathbf{A}\boldsymbol{\theta}_{i} + \alpha_{j} \right). \end{align} $$

In this formulation as an IRT model, non-separability can influence the systematic component via the product of a legislator’s ideal point on one dimension and a proposal’s discrimination parameter on another dimension. Preferences are additively separable if the off-diagonal elements of the weight matrix are zero, so that the cross-dimensional product equals zero too. Then, the systematic component is determined by the (weighted) sum of products between ideal points and discrimination parameters of each dimension.Footnote ³ If off-diagonal elements are not equal to zero, then this analysis is altered: cross-dimensional interactions representing non-separable preferences can influence the systematic component determining legislators’ decisions.

2.1 Identification and Estimation

Multidimensional latent variable models such as the IRT model described above are not identified without further constraints regarding the issues of location, scale, and rotation (Bafumi et al. Reference Bafumi, Gelman, Park and Kaplan2005; Jackman Reference Jackman2001; Rivers Reference Rivers2003). The locational issue arises because a constant can be added and subtracted from either term of the model without altering the result: $\boldsymbol {\beta }_{j}^{\prime }\mathbf {A}\boldsymbol {\theta }_{i} + \alpha _{j} = (\boldsymbol {\beta }_{j}^{\prime }\mathbf {A}\boldsymbol {\theta }_{i} - a) + (\alpha _{j} + a)$ . The scalar issue addresses that any element of the interaction term can be multiplied by a constant if another element is divided by the same constant: $\boldsymbol {\beta }_{j}^{\prime }\mathbf {A}\boldsymbol {\theta }_{i} = (\boldsymbol {\beta }_{j}/b)^{\prime }\mathbf {A}(b\boldsymbol {\theta }_{i})$ . Finally, $\mathbf {A}$ can be decomposed and multiplied with both sides of the interaction, leading to rotational invariance: $\boldsymbol {\beta }_{j}^{\prime }\mathbf {A}\boldsymbol {\theta }_{i} = \boldsymbol {\beta }_{j}^{\prime }\mathbf {R}^{\prime }\mathbf {R}\boldsymbol {\theta }_{i} = (\mathbf {R}\boldsymbol {\beta }_{j})^{\prime }(\mathbf {R}\boldsymbol {\theta }_{i}) = \boldsymbol {\tilde {\beta }}_{j}^{\prime }\boldsymbol {\tilde {\theta }}_{i}$ .Footnote ⁴ Therefore, we employ standard constraints on the means and variances of the ideal points and the discrimination parameters to address the issues of location and scale, while we use a hierarchical modeling approach for the ideal points to address the rotational issue (similar to Bafumi et al. Reference Bafumi, Gelman, Park and Kaplan2005).Footnote ⁵

This identification approach comes with a set of advantages. Importantly, it disentangles distinct quantities of interest: variation in non-separability and dimensional salience in legislative decision-making processes is captured by the weight matrix, whereas variation in the substantive meaning of conflict along dimensions is captured by the modeling of the ideal points. The confirmatory modeling approach regarding the latent dimensions is based on prior knowledge and allows us to assess additional quantities of interest regarding them. The estimated coefficients relating covariates to the ideal points tell us how strongly each covariate influences the position of legislators along each dimension. Comparing the predicted values of $\widehat {\boldsymbol {\theta }}$ based on the systematic component of their parametrization to the estimated values of $\boldsymbol {\theta }$ , we can assess how inductive the covariates are overall for the spatial dispersion of legislators via the $R^{2}$ statistic (Gelman et al. Reference Gelman, Goodrich, Gabry and Vehtari2019).

We estimate the model using Bayesian inference. Given the assumption of independence of the observed roll call votes for bills and legislators conditional on parameters, the item response model has the following likelihood:

(4) $$ \begin{align} L(\boldsymbol{\beta},\alpha,\mathbf{A}, \boldsymbol{\theta} \mid \mathbf{Y}) = \prod_{i=1}^{N} \prod_{j=1}^{J} P(y_{ij}=1)^{y_{ij}} (1-P(y_{ij}=1))^{(1-y_{ij})}, \end{align} $$

where $\boldsymbol {\beta }$ is a $J\times D$ matrix that collects the discrimination parameters $\boldsymbol {\beta }_{j}$ of J proposals in rows, $\alpha $ is a vector capturing proposals’ difficulty parameters, $\mathbf {A}$ is the weight matrix, $\boldsymbol {\theta }$ is an $I\times D$ matrix of ideal points $\boldsymbol {\theta }_{i}$ of I legislators in rows, and $\mathbf {Y}$ is an $I \times J$ matrix of binary recorded votes, with $y_{ij}$ at position $(i,j)$ .

We specify priors for the item parameters $\boldsymbol {\beta }$ and $\alpha $ , the weight matrix $\mathbf {A}$ , and the ideal points $\boldsymbol {\theta }$ . Item parameters are assigned standard normal distributions as prior distributions, so that $\beta _{jd} \sim \mathcal {N}(0,1)$ and $\alpha _{j} \sim \mathcal {N}(0,1)$ . We decompose $\mathbf {A}$ in a manner similar to a variance-covariance matrix composed of a correlation matrix and scaling parameters ( $A=\tau \Omega \tau$ ; see Barnard, McCulloch, and Meng Reference Barnard, McCulloch and Meng2000).Footnote ⁶ This decomposition assures that $\mathbf {A}$ is symmetric positive definite. The correlation matrix component $\boldsymbol {\Omega }$ of the weight matrix $\mathbf {A}$ is drawn from a prior distribution defined in Lewandowski, Kurowicka, and Joe (Reference Lewandowski, Kurowicka and Joe2009) with $\nu = 1$ and the scale parameters $\tau $ from a truncated normal distribution $\tau \sim \mathcal {N}_{+}(0,1)$ . Finally, the ideal point of legislator i on dimension $d \in (1,\ldots ,D)$ , $\theta _{id}$ is drawn from a normal distribution with legislator-specific mean $\theta _{id} \sim \mathcal {N}(\widehat {\theta }_{id}, 1)$ . The parametrization of $\widehat {\theta }_{id}$ is application-specific, and we discuss it in more detail below. We approximate the posterior distribution of our model using Markov Chain Monte Carlo simulations, relying on the No-U-Turn Sampler as implemented in Stan (for code, see Section 1 of the Supplementary Material; Carpenter et al. Reference Carpenter2017), and assess convergence across chains via the $\widehat {R}$ statistic (Gelman and Rubin Reference Gelman and Rubin1992).

3 Monte Carlo Simulations

In this section, we compare the estimates from a separable and non-separable model using Monte Carlo simulations.Footnote ⁷ We focus on two parameters: (i) the dimensional salience weights in estimated matrices $\mathbf {A}$ and (ii) estimated ideal points $\boldsymbol {\theta }$ . Both are of central interest to applied researchers, as they speak to the relevance and substance of dimensions, and the spatial positioning of legislators. More details of the data-generating process and additional results (including results across scenarios, further parameters, and model comparisons) are provided in Section 2 of the Supplementary Material.

Throughout, we specify that 100 legislators consider deviations from their ideal points along two dimensions and vote on 250 roll call votes. We simulate a total of 1,000 legislatures for which we vary the extent of non-separability, the relative salience of the two dimensions in the weight matrix $\mathbf {A}$ , and the correlation between ideal points $\boldsymbol {\theta }$ to accomodate a realistic range of potential scenarios. In line with our confirmatory identification approach, ideal points $\theta _{.d}$ are modeled via $\widehat {\theta }_{.d} = \gamma _{d} z_{d}$ , where $\gamma _{d} = 1$ and $z_{d} \in \{-1, 0, 1\}$ . $z_{d}$ resembles a group membership index, specifying which group of legislators should be positioned toward which end of the latent scale. All other parameters ( $\alpha $ , $\boldsymbol {\beta }$ ) are sampled from standard normal distributions.

A separable specification overestimates the salience of a primary dimension in the presence of non-separability, thereby distorting estimates of dimensional salience in legislators’ roll call votes. Figure 2 displays the deviation in estimated dimensional salience in weight matrices $\mathbf {A}$ across varying levels of non-separability. If preferences are separable and non-separability equals 0, both specifications return equal estimates. However, the separable specification—represented by the dashed line—systematically overestimates the salience of the primary dimension and underestimates that of the secondary dimension at higher absolute levels of non-separability. The non-separable specification—the solid line—does not exhibit this pattern. The non-separable model estimates the true non-separability value accurately (see Figure 2 in the Supplementary Material).

Figure 2 Deviation in estimated dimensional salience in weight matrices $\mathbf {A}$ .

Figure 3 Estimated ideal points $\boldsymbol {\theta }$ .

The estimated ideal points $\boldsymbol {\theta }$ from a separable specification less accurately reflect the original latent parameters if non-separability matters. Figure 3a displays the correlation between true and estimated ideal points across dimensions. Overall, the correlation decreases as non-separability increases. However, this pattern is more pronounced for the separable specification than for the non-separable specification. The separable specification estimates an increasingly distorted representation of legislators’ preferences at higher absolute values of non-separability compared to the non-separable specification.

These distortions—in dimensional salience and ideal points—are attributable to the way by which the separable specification accommodates non-separability: the estimated ideal points become increasingly related across dimensions. This is visible in Figure 3b, which displays the correlation of ideal points across dimensions. Across simulations, the expected correlation equals zero. After estimation, ideal points from the separable specification are correlated in a manner consistent with the direction and size of non-separability. Similar to the other quantities focused on in our discussion, the non-separable specification does not exhibit similar patterns.

These results concentrate on distortions relative to latent parameters that are known in simulations, but unknown otherwise. Therefore, we also consider model fit relative to the data. Out-of-sample fit approximated by the leave-one-out cross-validation (LOO) information criteria shows that the non-separable specification generally performs better (see Figure 7 in the Supplementary Material; Vehtari, Gelman, and Gabry Reference Vehtari, Gelman and Gabry2017), while in-sample predictions show little difference between the two models (see Figure 8 in the Supplementary Material). This is in line with the observation that the separable specification estimates alternative sets of latent parameters to explain the same variance in observed roll call votes. At increasing levels of non-separability, the effective dimensionality of the latent space decreases: if two dimensions are perfect substitutes or complements, indifference contours collapse to a line representing a mixture of the original latent dimensions. In such a scenario, estimates from the separable model resemble an increasingly unidimensional scenario (via deviations in dimensional salience) in which the dominant dimension also reflects variation on the secondary dimension (via correlated ideal points) as it cannot accomodate such dimensionality reduction via non-separability. These patterns are more accentuated when both dimensions are similarly salient. The non-separable specification more accurately portrays legislators’ decision-making processes in such scenarios as characterized by distinct, but related latent dimensions.

4 Applications

We focus on two applications, the U.S. Senate (37 Congresses between 1945 and 2019; Lewis et al. Reference Lewis, Poole, Rosenthal, Boche, Rudkin and Sonnet2020) and the EP (five sessions between 1979 and 2004; Hix et al. Reference Hix, Noury and Roland2006). We compare the results from the non-separable and separable model for each Congress and EP session separately, focusing on the same parameters as in the Monte Carlo discussion: (i) non-separability and dimensional salience as elements of the weight matrix $\mathbf {A}$ and (ii) legislators’ ideal points $\boldsymbol {\theta }$ . More details and results are provided in Sections 3 and 4 in the Supplementary Material for the Senate and the EP, respectively.

In line with prior research, we assume that roll call votes in both settings are determined by two latent dimensions whose salience and substantive meaning may vary over time. In the U.S. Senate, these dimensions reflect interparty conflict between Democrats and Republicans as well as intraparty conflict within parties (e.g., Clausen and Cheney Reference Clausen and Cheney1970; Nye Reference Nye1991; Poole and Rosenthal Reference Poole and Rosenthal1991, Reference Poole and Rosenthal2001, Reference Poole and Rosenthal2007). While the former conflict revolves around economic issues, the substance of the latter varies more over time: in earlier Congresses, intraparty conflict was more strongly associated with racial/social issues. But as party lines became the increasingly important determinant of legislative behavior in the U.S. Senate overall, the substantive meaning of this dimension has come to vary more as it accommodates residual variation in votes not due to partisanship (Poole and Rosenthal Reference Poole and Rosenthal1991, Reference Poole and Rosenthal2001, Reference Poole and Rosenthal2007). The meaning of the two dimensions in the EP is more consistent over time: while one dimension can be characterized as a left-right divide, a second dimension reflects conflict regarding EU integration (Brack Reference Brack2013; Hix et al. Reference Hix, Noury and Roland2006).

Figure 4 U.S. Senate—estimated weight matrices $\mathbf {A}$ .

We parametrize the distribution of legislators’ ideal points along dimensions accordingly. In the United States, Democrats are pit against Republicans to identify the interparty dimension, while Senators from racially liberal coastal states are juxtaposed with those from racially conservative states of the former Confederacy to identify the intraparty dimension (lower values on dimensions represent more liberal positions; Acharya, Blackwell, and Sen Reference Acharya, Blackwell and Sen2016; Key Jr. Reference Key1949). In the EP, members of left-wing and of right-wing party groups stand at opposing ends of the left–right dimension, while members of EU-favorable and EU-sceptic party groups occupy the poles of the EU dimension (lower values represent more left-wing or EU-favorable positions).

Additionally, we leverage that many legislators are incumbents to aid in the identification of the dimensions over time. In the U.S. Senate, around 85% of Senators in one Congress will also be present in the next Congress, while this is the case for around 45% of the members of parliament in a session of the EP. We exploit this by adding the residual of a legislator’s ideal point from her group-level mean in the last legislative session as a individual-level predictor of her ideal point in the current session on the same dimension (akin to a random walk in dynamic models; e.g., Lo Reference Lo2018).Footnote ⁸

4.1 U.S. Senate

We begin our discussion with the U.S. Senate. Elements of the estimated weight matrices ( $\mathbf {A}$ ) of each Congress after 1945 are shown in Figure 4. The estimated non-separability parameters from the non-separable specification are shown in Figure 4a (mean and 95% credible interval). The resulting estimates indicate a substantial degree of non-separability across dimensions in most Congresses. Cross-dimensional considerations prominently figure in legislative behavior in the U.S. Senate, and the underlying general dimensions are not independent in Senators’ roll call votes.

The nonzero non-separability parameters have implications for the estimates of dimensional salience in the decision-making process of Senators as shown in Figure 4b. In line with the Monte Carlo simulations, the separable specification (indicated by the dashed line) overestimates the salience of a dominant dimension and underestimates the weight of the secondary dimension compared to the non-separable model. This pattern is especially visible in the period between the 89th Congress (1965–1967) and the 98th Congress (1983–1985). This period is characterized both by comparatively high levels of non-separability and similarly salient dimensions. Whereas the non-separable specification estimates roughly equal salience weights for both dimensions during this period (indicated by the solid line), the separable specification attributes more weight to the intraparty dimension than to the interparty dimension. Since then, the salience of the interparty dimension has increased and that of the intraparty dimension has decreased in both specifications (in line with, e.g., Poole and Rosenthal Reference Poole and Rosenthal2007).

The influence of non-separability is also visible in the correlation of ideal points across dimensions in Figure 5. Similar to the findings of the Monte Carlo simulations, the separable specification accommodates non-separability by estimating correlated ideal points. From the 86th Congress (1959–1961) to the 102nd Congress (1991–1993), the separable specification estimates a much higher degree of positive correlation across dimensions than the non-separable specification. Accordingly, the correlation between the estimated ideal points across models on the secondary intraparty dimension drops to values of around .6 to .8 in this period (see Figure 10 in the Supplementary Material).

Figure 5 U.S. Senate estimated ideal points $\boldsymbol {\theta }$ —correlation of ideal points across dimensions.

A clarification in the substantive meaning of the intraparty dimension by explicitly incorporating non-separability is visible if we compare the ideal points of Senators across specifications. In Figure 6, the top 10 Senators in terms of their spatial shifts to the extremes of this dimension are shown for the period between the 89th Congress and the 99th Congress, just after the Civil Rights Era.Footnote ⁹ Substantially, the intraparty dimension is more strongly determined by racial/social issues in the non-separable specification: while liberal members of the Republican Party such as Jacob Javits, John Chaffee, and Charles Percy shift downward toward a more liberal position on this dimension, conservative Democrats such as J. James Exon, James Fulbright, and James Sasser shift upward toward a more conservative position. We also find that Senators are more strongly polarized according to whether they stem from racially/socially conservative or liberal states on this dimension in the non-separable rather than the separable specification in this period (see Figures 11 and 12 in the Supplementary Material).

Figure 6 U.S. Senate—shifts in estimated ideal points $\boldsymbol {\theta }$ on the intraparty dimension. $\circ $ = separable specification, $\bullet $ = non-separable specification. The rank of each Senator (out of 137) is shown beside each ideal point.

Taken together, the positive non-separability estimate in the period just after the Civil Rights Era indicates that Senators favored outcomes combining deviations in opposite directions across dimensions (rather than the same direction). This means that proposals bundling more liberal outcomes on one dimension with more conservative outcomes on the other dimension were preferred to those combining more liberal or conservative outcomes on both (because lower values are associated with more liberal preferences and higher values are associated with more conservative preferences on both dimensions in this time). The non-separable specification reflects this quid pro quo between economic and racial/social issues in legislators’ decisions, highlights the presence of distinct, but related dimensions, and clarifies the meaning of the intraparty dimension in this period.

4.2 European Parliament

Non-separability also figures in legislative behavior in the EP. Figure 7 displays the recovered non-separability parameter in the five sessions of the EP. In four sessions, the estimated parameter is positive and the 95% credible intervals do not include zero. Similar to the post–Civil Rights Era in the U.S. Senate, the two dimensions act as substitutes: legislators favor proposals representing deviations in opposite directions on the left/right dimension and the EU integration dimension. Because of the polarity of these dimensions, this means that proposals combining more right-wing outcomes with more EU integration (or more left-wing outcomes with less integration) are preferred to those representing deviations in the same direction. The degree of non-separability has increased in the fourth and fifth sessions, potentially indicating that the association between dimensions has increased in EP roll call voting.

Figure 7 EP—estimated non-separability parameter across sessions.

The comparatively low absolute level of non-separability limits its impact on other estimated parameters. The two sets of estimates are very similar across specifications. This concerns estimated salience parameters (Figure 21 in the Supplementary Material), the correlation between estimated ideal points across specifications (Figure 22 in the Supplementary Material), the estimated coefficients of covariates structuring the ideal points (Figure 26 in the Supplementary Material), as well as the extent to which the ideal points are structured by the systematic component of their parametrization (Figure 25 in the Supplementary Material).