1 Introduction
Political scientists are often interested in making inferences about the ideology of legislators based on ideal-point models. For example, scholars have sought to make inferences about party influence (Snyder and Groseclose Reference Snyder and Groseclose2000; McCarty, Poole, and Rosenthal Reference McCarty, Poole and Rosenthal2001), the dimensionality of voting (e.g., Rosenthal Reference Rosenthal1992; Roberts, Smith, and Haptonstahl Reference Roberts, Smith and Haptonstahl2016), pivot and cartel theories of lawmaking (e.g., Chiou and Rothenberg Reference Chiou and Rothenberg2003; Covington and Bargen Reference Covington and Bargen2004), the identity of the most liberal or conservative legislators (Clinton, Jackman, and Rivers Reference Clinton, Jackman and Rivers2004a), and ideological change over time (e.g., Epstein et al. Reference Epstein, Martin, Quinn and Segal2007; Poole Reference Poole2007).
The bulk of these studies have used parametric ideal-point models that specify particular functional forms for the utility functions and error distributions. It is well known that the choice of parametric assumptions can affect the estimates and inferences derived from data. This has been demonstrated in the ideal-point context (e.g., Ho and Quinn Reference Ho and Quinn2010; Krehbiel and Peskowitz Reference Krehbiel and Peskowitz2015) and has led studies to different conclusions on questions such as the dimensionality of congressional voting (see Poole and Rosenthal Reference Poole and Rosenthal1991; Heckman and Snyder Reference Heckman and Snyder1997). Yet scholars rarely have theoretical or empirical grounds for choosing a particular functional form.
The only alternative approach has been to use Optimal Classification (OC), a model-free estimation method introduced by Poole (Reference Poole2000). Because OC “is not a statistical model” (Poole Reference Poole2000, p. 211), statistical inference is not possible.Footnote 1 OC is also not a consistent estimator of the rank order of the ideal points under quadratic utility.Footnote 2 Thus, for many purposes, OC cannot be used in place of a parametric model.
The main contribution of this paper is to introduce a method for ideal-point inference and estimation under much weaker assumptions than previous approaches. It requires no parametric assumptions, making it the first approach to nonparametric ideal-point inference.Footnote 3 Its estimates of the rank order of ideal points are also consistent under the same weak, nonparametric assumptions.
This method allows for formal hypothesis tests of whether legislators have the same ordering on different groups of bills without the validity of the tests depending on knowledge of the parametric form, which is almost always unavailable. Because there are no assumptions made about the distribution of bill characteristics in the two groups, these inferences remain valid even if there are also differences in agenda control or other selection effects affecting the bills being voted on in each group. For example, this method can be applied to whether a legislator’s voting behavior changed after a hypothesized point in time, such as a party switch or change in public opinion. Similarly, it can be applied to testing hypotheses regarding dimensionality by testing whether a common ideological dimension can capture voting on different groups of bills.
The model consists of two key assumptions: unidimensionality and concave utility.Footnote 4 Aside from unidimensionality, the model’s assumptions are far weaker than in most ideal-point models. Not only is the model not restricted to quadratic utility or errors drawn from particular parametric families, but the utility functions need not be symmetric and the error distributions can vary across votes provided the errors follow the same distribution across legislators for each vote. The utility functions can be any concave function as long as the function is shared by all legislators, aside from differences in their ideal points.Footnote 5
The rest of the paper is organized as follows. The next section introduces the nonparametric ideal-point model used throughout this paper. Sections 3 and 4 then present methods for estimation and inference, respectively, under this model. Section 5 explores the properties of these methods through Monte Carlo simulations. Section 6 presents two applications of the model to Supreme Court data. The last section gives a concluding discussion of the methods and their possible extensions.
2 Model
Consider a legislature of
$K$
legislators voting on
$N$
bills. Assume legislators have unidimensional spatial preferences and vote sincerely according to the stochastic utility function given by
where
$v_{ij}$
is the vote of legislator
$i$
on bill
$j$
,
$v_{ij}=1$
indicates a vote of “yea,”
$v_{ij}=0$
indicates a vote of “nay,”
$y_{j1}$
is the policy outcome of the bill,
$y_{j0}$
is the status quo,
$x_{i}$
is the ideal point of legislator,
$\unicode[STIX]{x1D701}_{ij1}$
and
$\unicode[STIX]{x1D701}_{ij0}$
are stochastic errors, and
$G_{j}$
is the joint distribution of
$\unicode[STIX]{x1D701}_{ij1}$
and
$\unicode[STIX]{x1D701}_{ij0}$
for bill
$j$
. Assume that
$f$
is a strictly concave function with a maximum at zero.Footnote
6
Further, assume that the errors,
$\unicode[STIX]{x1D701}_{ij1}$
and
$\unicode[STIX]{x1D701}_{ij0}$
, are independent across bills and legislators and identically distributed for each bill. Finally, assume legislators vote “yea” when indifferent.Footnote
7
Note that they only differ in their ideal points,
$x_{i}$
, and not in
$f$
, which is assumed to be the same for all legislators. Note also that the distribution
$G_{j}$
of the stochastic utility components,
$\unicode[STIX]{x1D701}_{ij1}$
and
$\unicode[STIX]{x1D701}_{ij0}$
, is assumed to be the same for all legislators for any given bill but can vary across bills.
The probability of legislator
$i$
voting “yea” on bill
$j$
is given by
Let
$H_{j}$
be the cumulative distribution function of
$\unicode[STIX]{x1D701}_{ij0}-\unicode[STIX]{x1D701}_{ij1}~\left(H_{j}\left(z\right)=\text{Pr}\left(\unicode[STIX]{x1D701}_{ij0}-\unicode[STIX]{x1D701}_{ij1}\leqslant z\right)\right)$
. Thus,
The results in this paper rest on the implication of this model that
$\text{Pr}\left(v_{ij}=1\right)$
is monotonic in
$x_{i}$
. It is also worth noting that, if
$y_{j0}\neq y_{j1}$
and
$H$
is strictly increasing, then
$\text{Pr}\left(v_{ij}=1\right)$
is strictly monotonic in
$x_{i}$
. These two conditions require only that the policy outcome of the bill differs from the status quo and that the error distribution has support over the entire real line.
The following theorem establishes these results:
Theorem 1. If
$f$
is strictly concave,
$H_{j}\left(f\left(z-y_{j1}\right)-f\left(z-y_{j0}\right)\right)$
is monotonic in
$z$
. Further, if
$y_{j0}\neq y_{j1}$
and
$H$
is strictly increasing, then
$H_{j}\left(f\left(z-y_{j1}\right)-f\left(z-y_{j0}\right)\right)$
is strictly monotonic in
$z$
.
Proof. Consider the three cases
$y_{j0}=y_{j1}$
,
$y_{j0}<y_{j1}$
, and
$y_{j0}>y_{j1}$
.
- Case 1:
$y_{j0}=y_{j1}$
. In this case,
(6)
$$\begin{eqnarray}H_{j}\left(f\left(z-y_{j1}\right)-f\left(z-y_{j0}\right)\right)=H_{j}\left(f\left(z-y_{j0}\right)-f\left(z-y_{j0}\right)\right)=0.\end{eqnarray}$$
Since this constant in
$z$
, it is also monotonic in
$z$
.- Case 2:
$y_{j0}<y_{j1}$
. Take any
$x_{i}<x_{k}$
. Since
$f$
is strictly concave,
$tf\left(z_{1}\right)+\left(1-t\right)f\left(z_{2}\right)<f\left(tz_{1}+\left(1-t\right)z_{2}\right)$
for all
$z_{1}\neq z_{2}$
and
$t\in \left(0,1\right)$
. Let
$t^{\prime }=\frac{y_{j1}-y_{j0}}{x_{k}-x_{i}+y_{j1}-y_{j0}}$
. Since
$y_{j1}-y_{j0}>0$
and
$x_{k}-x_{i}>0$
,
$t^{\prime }\in \left(0,1\right)$
. Thus, (7)
$$\begin{eqnarray}\displaystyle & & \displaystyle \displaystyle t^{\prime }f\left(x_{i}-y_{j1}\right)+\left(1-t^{\prime }\right)f\left(x_{k}-y_{j0}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad <f\left(t^{\prime }\left(x_{i}-y_{j1}\right)+\left(1-t^{\prime }\right)\left(x_{k}-y_{j0}\right)\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\left(\frac{y_{j1}-y_{j0}}{x_{k}-x_{i}+y_{j1}-y_{j0}}\right)\left(x_{i}-y_{j1}\right)+\left(1-\frac{y_{j1}-y_{j0}}{x_{k}-x_{i}+y_{j1}-y_{j0}}\right)\left(x_{k}-y_{j0}\right)\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\left(\frac{y_{j1}-y_{j0}}{x_{k}-x_{i}+y_{j1}-y_{j0}}\right)\left(x_{i}-y_{j1}\right)+\left(\frac{x_{k}-x_{i}}{x_{k}-x_{i}+y_{j1}-y_{j0}}\right)\left(x_{k}-y_{j0}\right)\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\frac{\left(y_{j1}-y_{j0}\right)\left(x_{i}-y_{j1}\right)+\left(x_{k}-x_{i}\right)\left(x_{k}-y_{j0}\right)}{x_{k}-x_{i}+y_{j1}-y_{j0}}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\frac{\left(y_{j1}x_{i}-y_{j0}x_{i}-y_{j1}y_{j1}+y_{j0}y_{j1}\right)+\left(x_{k}x_{k}-x_{i}x_{k}-x_{k}y_{j0}+x_{i}y_{j0}\right)}{x_{k}-x_{i}+y_{j1}-y_{j0}}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\frac{y_{j1}x_{i}-y_{j1}y_{j1}+y_{j0}y_{j1}+x_{k}x_{k}-x_{i}x_{k}-x_{k}y_{j0}}{x_{k}-x_{i}+y_{j1}-y_{j0}}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\frac{y_{j1}x_{i}-y_{j1}y_{j1}+y_{j0}y_{j1}+x_{k}x_{k}-x_{i}x_{k}-x_{k}y_{j0}+x_{k}y_{j1}-x_{k}y_{j1}}{x_{k}-x_{i}+y_{j1}-y_{j0}}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\frac{y_{j1}\left(-x_{k}+x_{i}-y_{j1}+y_{j0}\right)+x_{k}\left(x_{k}-x_{i}+y_{j1}-y_{j0}\right)}{x_{k}-x_{i}+y_{j1}-y_{j0}}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\frac{x_{k}\left(x_{k}-x_{i}+y_{j1}-y_{j0}\right)-y_{j1}\left(x_{k}-x_{i}+y_{j1}-y_{j0}\right)}{x_{k}-x_{i}+y_{j1}-y_{j0}}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\frac{\left(x_{k}-y_{j1}\right)\left(x_{k}-x_{i}+y_{j1}-y_{j0}\right)}{x_{k}-x_{i}+y_{j1}-y_{j0}}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(x_{k}-y_{j1}\right).\end{eqnarray}$$
Similarly,
(8)
$$\begin{eqnarray}\displaystyle & & \displaystyle \displaystyle \left(1-t^{\prime }\right)f\left(x_{i}-y_{j1}\right)+t^{\prime }f\left(x_{k}-y_{j0}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad <f\left(\left(1-t^{\prime }\right)\left(x_{i}-y_{j1}\right)+t^{\prime }\left(x_{k}-y_{j0}\right)\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\left(1-\frac{y_{j0}-y_{j1}}{x_{i}-x_{k}+y_{j0}-y_{j1}}\right)\left(x_{i}-y_{j1}\right)+\left(\frac{y_{j0}-y_{j1}}{x_{i}-x_{k}+y_{j0}-y_{j1}}\right)\left(x_{k}-y_{j0}\right)\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\left(\frac{x_{i}-x_{k}}{x_{i}-x_{k}+y_{j0}-y_{j1}}\right)\left(x_{i}-y_{j1}\right)+\left(\frac{y_{j0}-y_{j1}}{x_{i}-x_{k}+y_{j0}-y_{j1}}\right)\left(x_{k}-y_{j0}\right)\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\left(x_{i}-x_{k}\right)\left(x_{i}-y_{j1}\right)+\frac{\left(y_{j0}-y_{j1}\right)\left(x_{k}-y_{j0}\right)}{x_{i}-x_{k}+y_{j0}-y_{j1}}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\left(x_{i}x_{i}-x_{k}x_{i}-x_{i}y_{j1}+x_{k}y_{j1}\right)+\frac{\left(y_{j0}x_{k}-y_{j1}x_{k}-y_{j0}y_{j0}+y_{j1}y_{j0}\right)}{x_{i}-x_{k}+y_{j0}-y_{j1}}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\frac{y_{j0}x_{k}-y_{j0}y_{j0}+y_{j1}y_{j0}+x_{i}x_{i}-x_{k}x_{i}-x_{i}y_{j1}}{x_{i}-x_{k}+y_{j0}-y_{j1}}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\frac{y_{j0}x_{k}-y_{j0}y_{j0}+y_{j1}y_{j0}+x_{i}x_{i}-x_{k}x_{i}-x_{i}y_{j1}+x_{i}y_{j0}-x_{i}y_{j0}}{x_{i}-x_{k}+y_{j0}-y_{j1}}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\frac{y_{j0}\left(-x_{i}+x_{k}-y_{j0}+y_{j1}\right)+x_{i}\left(x_{i}-x_{k}+y_{j0}-y_{j1}\right)}{x_{i}-x_{k}+y_{j0}-y_{j1}}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\frac{x_{i}\left(x_{i}-x_{k}+y_{j0}-y_{j1}\right)-y_{j0}\left(x_{i}-x_{k}+y_{j0}-y_{j1}\right)}{x_{i}-x_{k}+y_{j0}-y_{j1}}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(\frac{\left(x_{i}-y_{j0}\right)\left(x_{i}-x_{k}+y_{j0}-y_{j1}\right)}{x_{i}-x_{k}+y_{j0}-y_{j1}}\right)\nonumber\\ \displaystyle & & \displaystyle \displaystyle \quad =f\left(x_{i}-y_{j0}\right).\end{eqnarray}$$
Adding these two results together yields
(9)
$$\begin{eqnarray}f\left(x_{i}-y_{j1}\right)+f\left(x_{k}-y_{j0}\right)<f\left(x_{k}-y_{j1}\right)+f\left(x_{i}-y_{j0}\right).\end{eqnarray}$$
So,
(10)
$$\begin{eqnarray}f\left(x_{i}-y_{j1}\right)-f\left(x_{i}-y_{j0}\right)<f\left(x_{k}-y_{j1}\right)-f\left(x_{k}-y_{j0}\right).\end{eqnarray}$$
Since
$H$
is nondecreasing, (11)
$$\begin{eqnarray}H_{j}\left(f\left(x_{i}-y_{j1}\right)-f\left(x_{i}-y_{j0}\right)\right)\leqslant H_{j}\left(f\left(x_{k}-y_{j1}\right)-f\left(x_{k}-y_{j0}\right)\right)\end{eqnarray}$$
with strict inequality holding if
$H$
is strictly increasing. So, by definition of monotonicity,
$H_{j}\left(f\left(z-y_{j1}\right)-f\left(z-y_{j0}\right)\right)$
is monotonic in
$z$
with strict monotonicity if
$H$
is strictly increasing.- Case 3:
$y_{j0}<y_{j1}$
. Take any
$x_{i}<x_{k}$
. The same argument as in Case 2 with
$y_{j1}$
and
$y_{j0}$
swapped establishes that (12)
$$\begin{eqnarray}H_{j}\left(f\left(x_{i}-y_{j0}\right)-f\left(x_{i}-y_{j1}\right)\right)\leqslant H_{j}\left(f\left(x_{k}-y_{j0}\right)-f\left(x_{k}-y_{j1}\right)\right).\end{eqnarray}$$
Thus, negating both sides of the inequality yields
(13)
$$\begin{eqnarray}H_{j}\left(f\left(x_{i}-y_{j1}\right)-f\left(x_{i}-y_{j0}\right)\right)\geqslant H_{j}\left(f\left(x_{k}-y_{j1}\right)-f\left(x_{k}-y_{j0}\right)\right)\end{eqnarray}$$
with strict inequality holding if
$H$
is strictly increasing. So, by definition of monotonicity,
$H_{j}\left(f\left(z-y_{j1}\right)-f\left(z-y_{j0}\right)\right)$
is monotonic in
$z$
with strict monotonicity if
$H$
is strictly increasing.
Therefore, in all cases,
$H_{j}\left(f\left(z-y_{j1}\right)-f\left(z-y_{j0}\right)\right)$
is monotonic in
$z$
. If
$y_{j0}\neq y_{j1}$
and
$H$
is strictly increasing, then
$H_{j}\left(f\left(z-y_{j1}\right)-f\left(z-y_{j0}\right)\right)$
is strictly monotonic in
$z$
.◻
This model bares some similarity to a nonparametric item-response theory model developed by Mokken (Reference Mokken1971): the monotone homogeneity model. The crucial difference is that the above model assumes that the probability of a certain vote choice is merely a monotonic function of a legislator’s ideal point and can, thus, be either nondecreasing or nonincreasing. In contrast, the monotone homogeneity model assumes a priori knowledge of which answer is correct and, therefore, nondecreasing in a student’s ability. Unfortunately, such a simple approach is not applicable to the more general problem of nonparametric ideal-point estimation, where researchers cannot be presumed to know whether a “yea” vote is the more “conservative” or “liberal” vote choice.Footnote 8 Mokken’s (Reference Mokken1971) approach also requires an absence of missing values, further limiting its relevance to analyzing roll-call data.
3 Estimation
Under the model described in the previous section, knowing which of two legislators is to the right of the other does not tell us which is more likely to vote “yea” on a given bill. This is in sharp contrast to the monotone heterogeneity model, where knowing which student has higher ability directly implies which is more likely to answer a question correctly.
Most approaches to ideal-point estimation have addressed the need to know the polarity of each bill and other bill-specific parameters by simultaneously estimating these parameters.Footnote 9 Introducing these parameters can complicate both estimation and inference. For example, maximum likelihood estimators may be inconsistent when such parameters are included (Neyman and Scott Reference Neyman and Scott1948). Further, even aside from these complications, estimating the relationship between ideal points and vote probabilities is impractical here because the parameter space is infinite.
In this paper, I adopt a different approach that avoids estimating bill-specific parameters. The key insight for both estimation and inference is that, while we cannot learn much from the votes of a pair of legislators alone, the votes of two pairs can be informative. For any two pairs of legislators, if we condition on each pair containing one “yea” vote and one “nay” vote, then the conditional probability that the more conservative of the legislators in each pair voted together and the more liberal legislators in each pair voted together is greater than the probability that the more liberal of the legislators of each pair voted with the more conservative of the other.
More formally, if
$x_{i}<x_{k}$
and
$x_{\ell }<x_{m}$
, then either
$\text{Pr}\left(v_{kj}=1\right)\,-\,\text{Pr}\left(v_{ij}=1\right)\geqslant 0$
and
$\text{Pr}\left(v_{mj}=1\right)-\text{Pr}\left(v_{\ell j}=1\right)\geqslant 0$
or
$\text{Pr}\left(v_{kj}=1\right)-\text{Pr}\left(v_{ij}=1\right)\leqslant 0$
and
$\text{Pr}\left(v_{mj}=1\right)-\text{Pr}\left(v_{\ell j}=1\right)\leqslant 0$
. For notational convenience, let
$p_{i}=\text{Pr}\left(v_{ij}=1\right)$
and
$q_{i}=\text{Pr}\left(v_{ij}=0\right)=1-\text{Pr}\left(v_{ij}=1\right)$
. In either case,
So,
By independence of
$v_{ij}$
,
$v_{kj}$
,
$v_{\ell j}$
, and
$v_{mj}$
, this implies
Thus, if we consider only the votes where
$v_{ij}\neq v_{kj}$
and
$v_{\ell j}\neq v_{mj}$
, the expected count of votes on which
$v_{kj}=v_{mj}$
is greater than or equal to the expected count of votes on which
$v_{kj}=v_{mj}$
. Let
$$\begin{eqnarray}v_{im\mid k\ell }^{j}=\left\{\begin{array}{@{}ll@{}}1\quad & \text{if }v_{kj}=v_{mj}\wedge v_{ij}\neq v_{kj}\wedge v_{\ell j}\neq v_{mj}\\ 0\quad & \text{otherwise}\end{array}\right.\end{eqnarray}$$
so that
$v_{im\mid k\ell }^{j}$
indicates that legislators
$i$
and
$m$
voted together and legislators
$k$
and
$\ell$
voted together but each pair voted differently from the other on bill
$j$
. Under weak regularity conditions,Footnote
10
the law of large numbers implies that
Thus, if we consider only bills for which legislators
$i$
and
$k$
voted differently and legislators
$\ell$
and
$m$
voted differently, the fraction of these in which legislators
$k$
and
$m$
voted together should be greater than
$\frac{1}{2}$
.
Let
$\prec$
be the binary relation giving the true ordering of the ideal points where
$i\prec k$
indicates that legislator
$i$
is to the left of legislator
$k$
(i.e.,
$x_{i}<x_{k}$
). Let
$\hat{\prec }$
be the binary relation giving the estimated ordering relation where
$i\,\hat{\prec }\,k$
indicates that legislator
$i$
is estimated to be to the left of legislator
$k$
. If
$\hat{\prec }$
is a consistent estimator of
$\prec$
, then
must imply
$i\hat{\prec }k\Leftrightarrow \ell \hat{\prec }m$
for sufficiently large
$N$
and all
$i\neq k\neq \ell \neq m$
. Thus, a consistent estimator will minimize
$\sum _{j=1}^{N}v_{km\mid i\ell }^{j}-v_{im\mid k\ell }^{j}$
for sufficiently large
$N$
and all
$i\neq k\neq \ell \neq m$
such that
$i\hat{\prec }k\wedge \ell \hat{\prec }m$
. Therefore, it also minimizes the sum of this quantity over all possible groups of legislators satisfying
$i\hat{\prec }k\wedge \ell \hat{\prec }m$
,
where
$\mathbb{1}\left(x\right)$
is the indicator function (
$\mathbb{1}\left(x\right)=1$
if
$x$
is true and
$0$
otherwise) for sufficiently large
$N$
.
This suggests choosing as an estimator the relation
$\hat{\prec }$
that minimizes Formula (28). Identifiability also requires constraining the order of two legislators. Given this constraint and the necessary regularity condition (see footnote 25), this estimator is unique, at least for sufficiently large
$N$
, and therefore consistent.
Computing this estimator is difficult. There is no apparent alternative to an exhaustive search, but this requires considering
$\frac{1}{2}K!$
possible orderings for
$K$
legislators and summing over
$2\binom{K}{4}$
choices of legislators satisfying
$i\hat{\prec }k\wedge \ell \hat{\prec }m$
. Thus, computing Formula (28) for
$\frac{1}{2}K!$
possible orderings requires computing the inner sum of Formula (28)
$K!\binom{K}{4}$
times.Footnote
11
This is feasible for a sufficiently small number of legislators, but cases where
$n>10$
seem impractical.Footnote
12
Thus, for larger numbers of legislators, a different approach is necessary. However, if we choose two legislators,
$\ell$
and
$m$
, and consider only the votes in which they disagree, then the remaining
$K-2$
can be ordered based on the fraction of the votes in which each voted with
$m$
instead of
$\ell$
. As was the case when we considered two pairs of legislators, the expectation of the fraction of the votes in which each voted with
$m$
is increasing in their ideal points if
$\ell \prec m$
and decreasing if
$m\prec \ell$
.
If we repeat this process for each possible pair of legislators, we then only need to combine these orderings of each of the
$\binom{K}{2}$
sets of legislators into a single ordering for all legislators. Doing so can be done by creating a
$\binom{K}{2}\times K$
matrix with each row giving a numeric ranking of a different set of
$K-2$
with the two held-out legislators given the mean rank, then mean-centering each row, and finally performing a singular value decomposition (SVD). If all rankings on each subset of
$\binom{K}{2}$
legislators are correct, as will be true asymptotically, the order of the values in the first right-singular vector (i.e., the one corresponding to the largest singular value) should agree with the true ordering of the ideal points up to a reversal. Further, if it is too burdensome to compute rankings for all
$\binom{K}{2}$
groups of legislators, this can be replaced by a smaller number of groups of
$K-2$
chosen at random.
Although this estimator is also consistent under the same regularity conditions, it can only be based on roll-call items with no missing votes. In contrast, the alignment of two pairs of legislators can be computed using only the votes of those four legislators. Thus, the original estimator proposed in this paper is able to use any roll-call items in which at least four legislators voted. Since roll-call data often contain few or no items without any missing votes, the SVD-based approach as outlined above is not feasible but can be modified by computing the fraction of times one legislator agrees with another with missing values imputed by the mean agreement on nonmissing votes. This estimator need not be consistent but is likely to return similar results to the other, consistent estimators provided the fraction of missing votes is not too large.
4 Inference
Suppose we have two groups of roll-call items and wish to test whether a single ideological dimension could explain voting on both groups—or, equivalently, we wish to test whether the legislators’ ideal points maintained the same ordering on both items. The same approach used for estimation of choosing two pairs of legislators and focusing on those votes where each pair disagrees also provides a natural approach for performing inference. Indeed, the pairs-of-pairs approach is, in some ways, a more natural fit to inference than estimation and will provide a means to test the equal-ordering hypothesis without needing to estimate an ordering or to impose identifiability constraints.
Consider a two-sample extension of the model in Section 2 in which
$K$
legislators vote on two separate groups of bills, the first consisting of
$N_{A}$
items and the second of
$N_{B}$
. Let
$A$
be the set of indices of the first group, such that bill
$j$
is in group A if and only if
$j\in A$
and let
$B$
be the set of indices for group B, where
$A$
and
$B$
are disjoint. Assume the order of the legislators’ ideal points is given by the binary relation
$\prec _{A}$
for the first group of
$N_{A}$
bills and
$\prec _{B}$
for the second.Footnote
13
Finally, although most of the results in this section can be derived without this assumption, it will simplify matters to assume that the bill-specific parameters (
$y_{j1}$
,
$y_{j0}$
, and
$G_{j}$
) are jointly independently and identically distributed within each group, with distinct distributions for the bill parameters for each group and with no restrictions placed on the forms of these distributions. Thus, the bill parameters for one group may be drawn from a quite different distribution than the bill parameters in the other.
I aim to construct a test of the null hypothesis that the legislators have the same ordering on both groups—that is, that
$i\prec _{A}k\Leftrightarrow i\prec _{B}k$
for all
$i$
and
$k$
.
Let
indicate the marginal probability that
$\text{Pr}\left(v_{kj}=v_{mj}\mid v_{ij}\neq v_{kj}\wedge v_{\ell j}\neq v_{mj}\right)$
marginalized over the distribution of the bill parameters for group A, with
$\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(B\right)}$
defined analogously for group B. Thus, given the number of votes in which
$v_{ij}\neq v_{kj}\wedge v_{\ell j}\neq v_{mj}$
, the conditional distribution for the number of votes in which
$v_{kj}$
also equals
$v_{mj}$
(i.e., the sum of
$v_{km\mid i\ell }^{j}$
over
$A$
) is the sum of
$\sum _{j\in A}v_{km\mid i\ell }^{j}+v_{im\mid k\ell }^{j}$
trials each with success probability
$\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(A\right)}$
, where, as in the previous section,
$v_{km\mid i\ell }^{j}=\mathbb{1}\left(v_{kj}=v_{mj}\wedge v_{ij}\neq v_{kj}\wedge v_{\ell j}\neq v_{mj}\right)$
. Thus, the distribution is given by
Similarly,
If
$i\prec _{A}k$
and
$\ell \prec _{A}m$
. Under the null hypothesis,
$i\prec _{A}k\Leftrightarrow i\prec _{B}k$
, so
$i\prec _{B}k$
and
$\ell \prec _{B}m$
. From the previous section, this implies
$\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(A\right)}\geqslant \frac{1}{2}$
and
$\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(B\right)}\geqslant \frac{1}{2}$
. The same holds if
$k\prec _{A}i$
and
$m\prec _{A}\ell$
. If, instead,
$k\prec _{A}i$
and
$\ell \prec _{A}m$
or if
$i\prec _{A}k$
and
$m\prec _{A}\ell$
, then
$\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(A\right)}\leqslant \frac{1}{2}$
and
$\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(B\right)}\leqslant \frac{1}{2}$
. Thus, either
$\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(A\right)}\geqslant \frac{1}{2}$
and
$\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(B\right)}\geqslant \frac{1}{2}$
or
$\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(A\right)}\leqslant \frac{1}{2}$
and
$\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(B\right)}\leqslant \frac{1}{2}$
. We, thus, form a p-value for a test of the hypothesis that
$\left(\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(A\right)}\geqslant \frac{1}{2}\wedge \unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(B\right)}\geqslant \frac{1}{2}\right)\vee \left(\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(A\right)}\leqslant \frac{1}{2}\wedge \unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(B\right)}\leqslant \frac{1}{2}\right)$
. Call this value
$\unicode[STIX]{x1D70C}_{ik\ell m}$
.
We want to combine the values of
$\unicode[STIX]{x1D70C}_{ik\ell m}$
across all possible choices of the four indices
$ik\ell m$
. When combining p-values for disjoint groups of four legislators, the p-values are independent. Thus, the p-values can be aggregated via Fisher’s method. This gives
$-2\sum \log \left(\unicode[STIX]{x1D70C}_{ik\ell m}\right)$
, which, under the null hypothesis, will be stochastically dominated by a chi-squared distribution with degrees of freedom equal to twice the number of p-values being aggregated.
If p-values are computed for two groups of legislators that overlap, the p-values are not independent. However, the expectation of the sum of these p-values is unaffected by the lack of independence. Thus,
$E\left[\sum \left(-2\log \left(\unicode[STIX]{x1D70C}_{ik\ell m}\right)-1\right)\right]\leqslant 0$
under the null hypothesis. We can then compute the sample distribution of
$\sum \left(-2\log \left(\unicode[STIX]{x1D70C}_{ik\ell m}\right)-1\right)$
using a stratified nonparametric bootstrap with stratification of the two groups, A and B.Footnote
14
From this, we can test the hypothesis that
$E\left[\sum \left(-2\log \left(\unicode[STIX]{x1D70C}_{ik\ell m}\right)-1\right)\right]\leqslant 0$
, which then also serves as a test of the hypothesis that the orderings of the two groups are the same.
This test as constructed has relatively low power in cases where
$\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(A\right)},\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(B\right)}\ll \frac{1}{2}$
or
$\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(A\right)},\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(B\right)}\gg \frac{1}{2}$
. Under the null hypotheses, all p-values have at least 50% probability of being greater than
$\frac{1}{2}$
. However, in cases where
$\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(A\right)},\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(B\right)}\ll \frac{1}{2}$
or
$\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(A\right)},\unicode[STIX]{x1D70B}_{km\mid i\ell }^{\left(B\right)}\gg \frac{1}{2}$
, p-values may be particularly conservative and have at least 50% probability of being greater than
$\frac{1}{2}$
. To improve its power, we can only include p-values in the sum for groups of legislators where
$\sum _{j\in A}v_{km\mid i\ell }^{j}>\sum _{j\in A}v_{im\mid k\ell }^{j}$
and
$\sum _{j\in B}v_{km\mid i\ell }^{j}<\sum _{j\in B}v_{im\mid k\ell }^{j}$
or where
$\sum _{j\in A}v_{km\mid i\ell }^{j}<\sum _{j\in A}v_{im\mid k\ell }^{j}$
and
$\sum _{j\in B}v_{km\mid i\ell }^{j}>\sum _{j\in B}v_{im\mid k\ell }^{j}$
. These are the groups of legislators where the null hypothesis seems most suspect.
Of course, since we are selecting only a subset of p-values, their distribution will no longer be uniform under the null hypothesis as the p-value is not independent of its inclusion. Thus, without adjustment, we can no longer expect these p-values to be uniformly distributed under the null hypothesis. However, if we replace these p-values with p-values computed conditional on their inclusion—that is, based on the conditional distribution of a test statistic conditional on
$\sum _{j\in A}v_{km\mid i\ell }^{j}>\sum _{j\in A}v_{im\mid k\ell }^{j}$
and
$\sum _{j\in B}v_{km\mid i\ell }^{j}<\sum _{j\in B}v_{im\mid k\ell }^{j}$
or
$\sum _{j\in A}v_{km\mid i\ell }^{j}<\sum _{j\in A}v_{im\mid k\ell }^{j}$
and
$\sum _{j\in B}v_{km\mid i\ell }^{j}>\sum _{j\in B}v_{im\mid k\ell }^{j}$
—then the included p-values will again be uniformly distributed (or stochastically dominated by a uniform distribution) under the null hypothesis.
In some cases, we may wish to test whether a particular legislator changed location relative to any other legislator regardless of whether other legislators changed order. Doing so requires identifying two other legislators who can safely be assumed to have the same ordering in both groups. With this, we can compute the test statistic as above, but only varying
$k$
with
$i$
equal to the legislator whose movement is to be tested and
$\ell$
and
$m$
equal to the legislators identified as nonmoving. The test statistic is, thus,
$\sum _{k}\left(-2\log \left(\unicode[STIX]{x1D70C}_{ik\ell m}\right)-1\right)$
with nonpositive expectation under the null hypothesis. Its distribution can again be estimated via the stratified bootstrap.
5 Monte Carlo Studies
I conducted a series of Monte Carlo simulations to test the properties of the estimator and hypothesis test developed in the previous sections. For all simulations, I use a legislature of size
$K=9$
with legislator
$i$
having ideal point
$x_{i}=i$
. In each simulation, I generate probabilities of voting “yea” according to one of two possible distributions.
The first possible distribution for vote probabilities, which I will term the “uniform bill distribution,” the probability of voting in the conservative direction is
$\unicode[STIX]{x1D6F7}\left(2\left(x_{i}-\unicode[STIX]{x1D6FC}_{j}\right)\right)$
where
$\unicode[STIX]{x1D6F7}$
is the cumulative distribution function of a standard normal distribution, each
$\unicode[STIX]{x1D6FC}_{j}$
is drawn independently from a uniform distribution over the ideal points (
$x_{i}$
), the conservative direction is a “yea” vote with probability
$\frac{1}{2}$
and is independent of any other random variables in the model. This distribution is uniform in the sense that the point at which the probability of voting “yea” is exactly
$\frac{1}{2}$
occurs at the ideal point of each legislator with equal probability.
The second possible distribution, which I will term the “lopsided bill distribution,” the probability of legislator
$i$
voting in the conservative direction is
$1-\left(\frac{1}{2}\right)^{i}$
. Since
$x_{i}=i$
, this probability is strictly increasing in the ideal point of a legislator. As in the previous case, the probability that the conservative direction is a “yea” vote is also randomly determined with
$\frac{1}{2}$
probability of each. In this scenario, the bills being considered are both located toward one end of the spectrum, making the votes of the legislators at the other end more predictable. OC may be less likely to perform well in this scenario because the cutting points—that is, the points separating voters more likely to vote “yea” from those more likely to vote “nay”—occur only at one end of the spectrum.
5.1 Estimation
Table 1 shows the mean estimates and mean squared error based on the uniform bill distribution with 100 votes. Both the new estimator and the OC estimator recover the true rank of the ideal points quite well, although the new estimator has slightly lower mean squared error.
Estimate means for the uniform bill distribution.

Note:
$N=100$
votes. Based on 10,000 Monte Carlo simulations. “New” indicates the estimator discussed in Section 3.
Estimate means for lopsided bill distributions.

Note:
$N=100$
votes. Based on 10,000 Monte Carlo simulations. “New” indicates the estimator discussed in Section 3.
Table 2 shows the analogous results based on the lopsided bill distribution. Here, the legislators further to the right are all likely to vote together. Thus, it is difficult to estimate their relative positions. Both estimators do show higher mean squared error. However, while the new estimator and OC estimators performed relatively similarly with the uniform bill distribution, the new estimator performs significantly better here, with a mean squared error over 80% smaller.
Mean squared error by number of votes

Finally, Table 3 shows mean squared error based on data with different numbers of votes. Again, using the uniform bill distribution, both models perform roughly similarly. However, for data generated using the lopsided bill distribution, the new estimator performs noticeably better. Further, whereas the mean squared error for the new estimator decreases as more votes are added, the mean squared error of OC appears to increase slightly as the number of votes increases.
5.2 Inference
To test the inference, I simulate data for two groups of bills, with vote probabilities for each generated from the uniform bill distribution. For the first simulation, the legislators maintain the same ordering on both sets of bills. Thus, the null hypothesis is true, so the test should not reject at the 5% significance level in more than 5% of simulations. Note that, given the weak assumptions underlying the test, it is likely to be conservative and reject in significantly less than 5% of simulations.
In the next simulation, I choose two neighboring legislators at random and swap their order, which provides as minimal a violation of the null hypothesis as possible. Only one pair of legislators is in a different order in one group of votes as in the other. All other pairs maintain the same ordering. Thus, the null hypothesis is false and can be correctly rejected, but, given that the change is relatively small, then it may be difficult to detect.
Finally, I simulate data in which the ordering of the legislators for the first group of bills is unrelated to the ordering for the second. Here, each possible ordering of the legislators is equally likely, with the ordering chosen for the first group of bills independent of the ordering chosen for the second.
Histograms of p-values for test of same ordering from Monte Carlo simulations.

The results mirror expectations. Figure 1 shows histograms of p-values for each of the three scenarios discussed above and for each of 100, 500, and 1000 votes per group. The p-values in the first column, where the null hypothesis is true, tend to be greater than would be expected under a uniform distribution and the probability of rejection under the null hypothesis is below the nominal level of the test. This indicates the test is more conservative than possible, which is unsurprising.
In the second column, the null hypothesis is false but the orderings differ as little as possible. Here, the power of the test is quite low with 100 votes. However, as the number over voter increases, so does the power of the test. By 1000 votes, the probability of rejecting the null at the 5% significance level is 86.9%.
In the third column, the null hypothesis is false and a legislator’s ideology on one dimension is independent of her ideology on the other. Here, even with only 100 votes per group, the probability of rejecting the null is estimated to be 97.4%. With at least 500 votes per group, the power is over 99.9%.
6 Supreme Court
6.1 Did Justice Blackmun move to the left?
As a second example, consider the question of whether Justice Harry Blackmun’s ideology changed relative to those of the other justices during the 1970s. Blackmun is often hypothesized to have became more liberal during this period. While the tests developed in this paper cannot test this hypothesis in an absolute test, we can test whether Blackmun moved relative to the other justices. For the purposes of this example, I test specifically whether Blackmun’s position relative to the other justices during the Nixon and Ford administrations was the same as his position during the Carter and Reagan administrations.
First, we can test whether a single ordering can explain voting both before 1977 and after. Using the test developed in Section 4, we can strongly reject this hypothesis (
$p=0.0005$
). Thus, the data are incompatible with the justices having the same ordering in both time periods.
But this does not reveal whether Blackmun in particular changed. Choosing two justices to anchor the space permits a test of whether Blackmun in particular changed relative to any of the remaining justices. I use Justice Marshall and Chief Justice Burger. Note that this does not require assuming that Marshall and Burger do not move relative to the other justices. It only requires assuming that Burger is more conservative than Marshall in both time periods. Again, this test strongly rejects a null hypothesis that Blackmun remained on the same side of each of the remaining justices (
$p<0.0001$
).
Of course, this does not establish whether it was Blackmun’s ideology that changed or whether one or more of the other justices moved relative to Blackmun. However, we can test the hypothesis that the remaining justices maintained the same ordering across both time periods. With Blackmun excluded, the p-value becomes
$p=0.71$
. Thus, we see no evidence that the ordering of the remaining justices changed during this time period.
6.2 Is the Supreme Court multidimensional?
First, consider the question of whether the Supreme Court is multidimensional. Using data from the Supreme Court Database (Spaeth et al. Reference Spaeth, Epstein, Andrew D., Segal, Ruger and Benesh2015), I will use the hypothesis test developed in Section 4 to test whether a common dimension can explain voting on pairs of issue areas. I aggregate cases into four groups based on the Supreme Court Database’s coding of issue area: (1) criminal procedure, which includes only cases coded as “criminal procedure,” (2) civil liberties, which includes cases coded as “civil rights,” “First Amendment,” “due process,” and “privacy,” (3) economics, which includes cases coded as “economic activity,” “unions,” and “federal taxation,” and (4) other, which includes all other cases. Because the criminal procedure issue area is sometimes grouped together with the civil liberties issue areas and these issue areas are the ones “conventionally associated with the standard left–right splits on the Court” (Bailey Reference Bailey2013), I test whether a single ordering could explain voting on criminal procedure and civil liberties cases, treated as a single group, and on economics for each decade beginning with the 1950s.Footnote 15 I then similarly test whether a single ordering could explain voting on criminal procedure and economics, on civil liberties and economics, and on criminal procedure and civil liberties.
Tests of the hypothesis that the ideal points of Supreme Court justices are in the same order on cases on two different issue areas.

Note: Data from the Supreme Court Database (Spaeth et al. Reference Spaeth, Epstein, Andrew D., Segal, Ruger and Benesh2015). Civil liberties combines “civil rights,” “First Amendment,” “due process,” and “privacy” issue areas. Economics combines “economic activity,” “unions,” and “federal taxation” issue areas.
Table 4 shows the p-value computed for each comparison in each decade. Several comparisons are significant during the 1960s. Tests for the comparison between economics and criminal procedure, between economics and civil liberties, and between economics and both criminal procedure and civil liberties reject the null hypothesis that the ideal points have a common ordering during the 1960s at the 5% significance level. The test comparing criminal procedure and economics also nearly rejects for the 1970s (
$p=0.066$
).
The p-values reported in the table make no adjustment for the fact that tests are performed for six different decades. However, even when the p-values for each comparison are adjusted using the Holm procedure so that the family wise error rate—that is, the collective probability of falsely rejecting the null across all six decades—for each comparison across all six decades is less than 5%, the null hypothesis can be rejected for the comparisons between economics and civil liberties (adjusted p-value of
$0.039$
) and between economics and both criminal procedure and civil liberties (adjusted p-value of
$0.009$
).Footnote
16
Thus, the extent to which the Court is multidimensional—that is, the extent to which different orderings are needed to explain Supreme Court voting on different issue areas—may depend on the time period in question. For the 1960s, we can reject the hypothesis that the justices had the same ordering on cases involving economics as on criminal procedure and civil liberties. On the other hand, there is no evidence that multiple dimensions are needed for other decades. Note that, while the number of cases the Court considered was slightly lower during the 1990s and 2000s, this alone does not appear to explain the difference, both because there were roughly equal numbers of cases considered for each issue area during the 1960s and 1980s and because the significance of the results persist for the 1960s even when using shorter, five-year periods of time.
7 Conclusion
The nonparametric methods developed in this paper require only that voting is unidimensional, that the shape of the deterministic component of utility is concave and identical for all voters aside from differences in their ideal points, and that the errors are independent and identically distributed across legislators. The underlying model makes no parametric assumptions about the shape of the utility functions or error distribution.
Under these modeling assumptions, the estimator is consistent. In Monte Carlo simulations, it outperformed OC in terms of mean squared error by small margins in one case and larger margins in another. In some scenarios where OC is inconsistent, this is a significant improvement. Nonetheless, in most scenarios, both estimators are likely to return relatively similar results.
The method for nonparametric inference is more novel. Given two groups of bills, this method permits researchers to test the hypothesis that legislators’ ideal points have the same ordering on one group as on the other. This is not quite the same as testing whether the legislators have the same ideal points, as the ideal points could differ while still maintaining the same order. Nonetheless, a change in the order of the ideal points does imply a change in the ideal points themselves.
There are also several reasons why testing whether the orderings are the same can be preferable. Recovering cardinal information about ideal points requires the sort of stronger assumptions provided by parametric models, but researchers are rarely confident in the accuracy of these assumptions. Previous work has suggested using only ordinal information from ideal-point models because cardinal information can be very sensitive to parametric assumptions (Ho and Quinn Reference Ho and Quinn2010). Further, changes in the ordering may be substantively meaningful while at least small changes in the ideal points may not be.
The primary limiting assumption of the model in this paper is of unidimensionality. The remaining assumptions are significantly weaker than those adopted in other approaches to ideal-point estimation or inference. Most existing research on ideal-point estimation, such as Poole and Rosenthal (Reference Poole and Rosenthal1991), Heckman and Snyder (Reference Heckman and Snyder1997), and Clinton, Jackman, and Rivers (Reference Clinton, Jackman and Rivers2004b), requires assuming a particular parametric form for the utility functions and error distributions. The only other approaches, such as Poole (Reference Poole2000), do not use a stochastic model, removing the possibility of statistical inference. As far as I am aware, this paper is the first to provide an approach to provably consistent estimation and valid statistical inference for the order of ideal points without requiring knowledge of the parametric form of the utility functions or error distributions.
While the model cannot easily be extended to include a multidimensional Euclidean policy space, as this would render rank order a meaningful concept, other approaches to dealing with a single set of bills with multiple dimensions are possible. I suggest following the direction of Lauderdale and Clark (Reference Lauderdale and Clark2012) and Bonica (Reference Bonica2011) as well as nonparametric item-response theory by assuming legislators have unidimensional preferences for any bill but for which their ordering may change from one bill to another. Thus, one might be able to divide a set of bills into subsets that involve different dimensions, with legislators having different dimensions on each and the determination of which bills belong to the same subset perhaps informed by the sort of hypothesis test of equal ordering discussed in this paper. I leave this extension as a subject for future research.
Supplementary material
For supplementary material accompanying this paper, please visit https://doi.org/10.1017/pan.2017.38.






