A family of distances for preference–approvals

A preference–approval on a set of alternatives consists of a weak order on that set and, additionally, a cut-off line that separates acceptable and unacceptable alternatives. In this paper, we propose a new method for defining the distance between preference–approvals taking into account jointly the disagreements in preferences and approvals for each pair of alternatives. The proposed distance is compared to the existing distance functions to deal with clustering problems. Specifically, we prove that our metric improves the estimated clusters in terms of both stability and accuracy.


Introduction
In social choice theory, preference rankings and approvals are two popular ways to collect the preferences of a group of agents on a set of alternatives. Preference rankings order the alternatives from best to worst without distinguishing between acceptable and unacceptable alternatives. That is, if a is ranked above b, we can only infer that a is preferred to b, but we cannot infer anything about their absolute acceptability. In contrast, the approval voting system (Brams & Fishburn, 1978) consists of separating the set of acceptable alternatives from the set of unacceptable alternatives without considering preferences neither over acceptable nor over unacceptable alternatives.
Preference rankings and approval voting are related, but they are basically different types of information and cannot be inferred from each other.
In this paper, we focus on preference-approval structures. They combine preferences over the alternatives, through a weak order, and establish which alternatives are acceptable (Brams, 2008, Chapter 3;Brams & Sanver, 2009;Sanver, 2010). In preference-approval structures, voters can pay attention to which alternatives are acceptable and simultaneously rank-order them. Voters may either rank-order unacceptable alternatives or avoid declaring their preferences about them 1 by (implicitly) showing indifference between these alternatives 2 .
The distance and correlation between two preference rankings are of particular interest within this framework. Kendall's correlation coefficient (Kendall, 1948) is likely the most well-known ranking metric. Kendall's correlation coefficient is a measure of similarity that can be linearly transformed into a measure of dissimilarity (i.e., the Kendall's distance), which counts the number of pairwise disagreements between two rankings. Emond and Mason (2002, p. 20) demonstrated that when indifference between alternatives is allowed (weak orders), Kendall's distance violates the triangle inequality. Moreover, the Kendall's correlation of the all ties ranking with any other weak order is undefined, resulting in a 0/0. The Spearman's distance (Spearman, 1987), another famous ranking metric, between two rank vectors is calculated by taking the square root of the sum of the squared rank differences. Spearman's distance uses rank values as if they were mathematical variates, which leads to anomalous behaviour. Indeed, Spearman distance suffers from sensitivity to irrelevant alternatives (see Emond, 1997, p. 4;Emond & Mason, 2000, p. 16). In short, including additional irrelevant objects in the ranking exercise may alter the maximum agreement solution. Kemeny and Snell (1962) took a different approach to this problem. They defined a set of additional axioms that should be applicable to any distance measure between two weak orders, and introduced a distance that satisfies all these constraints. Besides the classical properties (positivity, symmetry, identity of indiscernibles and triangle inequality), the Kemeny distance is not affected by a relabeling of alternatives (neutrality), and it is consistent in measurement as the number of objects varies (i.e. it is not sensible to irrelevant alternatives). The Kemeny distance is a city block and a geodesic distance in the permutation polytope (see Heiser, 2004). It takes the shortest path between two rankings.
Generally, extensions to ranking measures have mainly focused on the definition of weighted distances (see García-Lapresta & Pérez-Román, 2010, Albano & Plaia 2021. In the last years, there has been a dramatic increase in recent publications about preference-approval structures and the introduction of consensus and distance measures in that setting. Erdamar et al. (2014) introduced a family of distances in the preference-approval setting, and they applied them to measure the consensus in that framework. Kamwa (2019) studied the propensity of the preference-approval voting of electing the Condorcet winner/loser when they exist. Dong et al. (2021) established some axioms implying the existence of a distinct distance function of preference-approval systems. They investigated a preferences aggregation model in the context of group decision-making based on the proposed axiomatic distance function. Kruger and Sanver (2021) investigated the compatibility between ordinal and evaluative approaches to social choice theory under two weak assumptions: respect for unanimity and independence of evaluation of each alternative. They claimed that there is an incompatibility between the two, and described some options whenever the second assumption is relaxed. Long et al. (2021) developed a two-stage consensus reaching method for multi-attribute group decision making problems with preference-approval structures, promoting the efficiency of consensus reaching. Barokas and Sprumont (2021) extended the classing Borda count to rank alternative in preference-approval setting, constructing an axiomatization of a new aggregation procedure called broken Borda rule.
In this paper, we propose a new distance for preference-approvals, following the axiomatic approach of the Kemeny distance. However, while the Kemeny distance can only consider the preference-discordance, our approach takes into account the approval-discordance as well, and use an aggregation function to combine the two types of information for each pair of alternatives.
We show that using, as an aggregation function, the family of weigthed power means (a class of weighted quasiarithmetic means) brings the benefit of many interesting properties. The final aggregated distance will thus be the sum of the pairwise preference-approval discordances. Furthermore, we show that our distance respects the fundamental properties to be defined as a metric and that, under certain assumptions, it has a precise geometric interpretation.
Our proposal can be regarded as the generalization of the Erdamar et al. (2014) distance measure, with the two coinciding for a specific parameter setting. However, we show that the proposed distance family has some advantages over the existing one as it is more versatile and performs better in cluster analysis.
Finally, the proposed metric is used to cluster a set of preference-approvals into homogeneous groups, considering the whole 2-dimensional universe of preference-approvals and a real case study.
The paper is organized as follows. Section 2 is devoted to introduce basic notation, preference-approvals and the codifications used throughout the paper. Section 3 includes our proposal for measuring distances between preference-approvals and some results. Section 4 offers some applications to the clustering task. Finally, Sect. 5 concludes the paper with some remarks.

Preliminaries
Let X = {x 1 , . . . , x n } a finite set of alternatives, with n ≥ 2. A weak order (or complete preorder) on X is a complete and transitive binary relation on X . A linear order on X is an antisymmetric weak order on X . With W (X ) and L(X ) we denote the set of weak and linear orders on X , respectively. Given R ∈ W (X ), with and ∼ we denote the asymmetric and the symmetric parts of R, respectively: Given a set Y , with P(Y ) we denote its power set, i.e., I ∈ P(Y ) ⇔ I ⊆ Y . In turn, with #Y we denote the cardinality of Y .

Preference-approval
Consider that a set of voters V = {v 1 , . . . , v m }, with m ≥ 2, have to express their opinions over X . We assume that each voter ranks the alternatives in X by means of a weak order and, additionally, assesses each alternative as either acceptable or unacceptable by partitioning X into A, the set of acceptable alternatives, and U = X \ A, the set of unacceptable alternatives, where A and U can be empty sets.
We also assume the following consistency condition: given two alternatives x i and x j , if x j is acceptable and x i is ranked above x j , then x i should be acceptable as well.
Definition 1 A preference-approval on X is a pair (R, A) ∈ W (X ) × P(X ) satisfying the following condition: With R(X ) we denote the set of preference-approvals on X .
Remark 1 If (R, A) ∈ R(X ), then the following conditions are satisfied: We now illustrate preference-approval structures through the following example.

Example 1 Let us consider
It means that alternatives in the same row are indifferent, alternatives in upper rows are preferred to those located in lower rows, alternatives above the line are acceptable, i.e., A = {x 1 , x 2 , x 4 , x 6 }, and those below the line are unacceptable, i.e., U = {x 3 , x 5 , x 7 , x 8 }. Table 1 includes the number of possible approvals, linear orders, weak orders and preference-approvals when the number of alternatives is n = 2, 3, . . . , 10.
It is well-known that the total number of approvals (subsets of X ) and linear orders are 2 n and n!, respectively. The number of weak orders is n!(log 2 e) n+1 /2 (see Good, 1980). The formula for calculating the number of preference-approvals has never been defined in the literature. For the first time, the exact number of preference-approvals for n = 2, 3, . . . , 10 alternatives is reported herein in Table 1. The formula to compute the exact number of preference-approvals on a set of n alternatives is where S (r ) n is a Stirling integer (number) of the second kind defined by David and Barton (1962, p. 294), Abramowitz and Stegun (1964, p. 824) and more thoroughly in Fisher and Yates (1953, p. 78), while r denotes the number of distinct positions in a weak order on n alternatives, also known as buckets. For example, considering four alternatives, if two are tied for first place and the other two are tied for third place, we can say that the number of distinct positions, or buckets, is two. Table 2 shows the quotients between preference-approvals and approvals. In turn, Table 3 shows the quotients between preference-approvals and weak orders.
It is clear that the expressivity of voters explodes with preference-approvals.

Codifications
Assigning positions to alternatives in linear orders is trivial because indifferences among distinct alternatives are not allowed. Given R ∈ L(X ), the position of each alternative x i ∈ X in R is defined through the mapping P R : X −→ {1, . . . , n} that assigns 1 to the alternative ranked first, 2 to the alternative ranked second, and so on.
There are different ways of assigning positions to the alternatives in weak orders. One of them is used by García-Lapresta and Pérez-Román (2011) and it is based on Smith (1973), Black (1976) and Cook and Seiford (1982).
The sign function, sgn : R −→ {−1, 0, 1}, is defined as Taking into account May (1952) and Fishburn (2015), we now introduce an index that codifies the order between two alternatives in a weak order R ∈ W (X ): It is worth noting that the index O R (x i , x j ) is also known in the literature as a i j (see Kemeny & Snell, 1962, p. 11) or score matrix (see Emond & Mason, 2002). In this paper, to avoid confusion with the approval-discordance notation of Eq. (8), we chose to use the notation O R (x i , x j ).

The proposal
Given two preference-approvals (R 1 , A 1 ), (R 2 , A 2 ) ∈ R(X ) and two generic alternatives x i , x j ∈ X , we now introduce two indices that measure the discordances between these alternatives with respect to preference and approvals, respectively.
The preference-discordance between x i and x j is defined as Taking into account Eqs. (5), (6) can be defined in an equivalent and simpler way: and therefore, p i j ∈ {0, 0.5, 1}.
The approval-discordance between x i and x j is defined as and again a i j ∈ {0, 0.5, 1}. In both cases, the values of 0, 0.5 and 1 indicate a null, moderate and high discordance, respectively. In order to generate a global measure of discordance between two alternatives, we consider an aggregation function (see Beliakov et al., 2007;Grabisch et al., 2009;Ramík & Vlach, 2012, Sect. 2, among others).

Definition 2 Given an aggregation function
Among the huge variety of aggregation functions, in this proposal we consider a class of weighted quasiarithmetic means 3 : the family of weighted power means, h : where λ ∈ [0, 1] and r > 0.

5.
Comparability: h is increasing in r .
Notice that the inputs of h in Eq. (9) are the pairs of 0, 0.5, 1. In Tables 4 and 5 we show the values of h for these pairs and different values of the parameter r for λ = 0.5, 0.75, respectively.
According to Tables 4 and 5 , the parameter r governs the penalty for each pair of values. Indeed, as r increases, so does the value of h( p i j , a i j ). As a result, taking an excessively large r value results in very similar penalties and reduces the weight of high discordance compared to moderate discordance.
Taking into account Eq. (9) with the aggregation function h in Eq. (10), we now introduce the family of distances on preference-approvals that we analyze in the present paper.
Definition 3 Given λ ∈ [0, 1] and r > 0, the distance associated with λ and r is the

Fig. 1 Preference-approval plane
Remark 4 When r = 2 and λ = 0.5, the geometric interpretation of h( p i j , a i j ) is related to the Euclidean distance. Figure 1 reports the preference-approval plane, that is a Euclidean plane having on the x-axis the preference-discordance, p i j , and on the y-axis the approval-discordance, a i j .
If r = 2 and λ = 0.5, then h( p i j , a i j ) is proportional to the Euclidean distance between ( p i j , a i j ) and the origin, This means that the aggregation function h can be interpreted as a proper distance in the preference-approval plane. As a result, the point of greatest discordance, (1, 1), will be the farthest from the origin of the axes. Conversely, (0, 0) represents the point of greatest agreement. The red segments in Fig. 1 (11)] can be interpreted as the sum of n·(n−1) 2 Euclidean distances in the preference-approval plane. That is, Proposition 1 D r λ is a metric on R(X ) for all λ ∈ (0, 1) and r ≥ 1. That is, for all (R 1 , A 1 ), (R 2 , A 2 ) ∈ R(X ) the following conditions are satisfied 4 : , then D r λ does not verify the identity of indiscernibles, hence it is not a metric.
Proposition 2 demonstrate that our proposal can be considered as a generalization of the preference-approval distance proposed by Erdamar et al. (2014).
Given two preference-approvals , is generated from the preference distance and the approval distance marginally, and eventually aggregate them by a convex combination.
The authors measure the disagreement between preferences by using the Kemeny metric (Kemeny, 1959), d K : Or, equivalently, by considering Eq. (5): In turn, the approval disagreement is measured through the Hamming metric (Hamming, 1950), d H : In order to aggregate d K and d H as a global distance, the two metrics are normalized to the same codomain [0, 1] via dividing by their maximum distances.
The mappings d R : The two normalized distances are eventually aggregated in a final preference-approval distance, d λ : where λ ∈ [0, 1] is a parameter used to control the relative relevance of the two components. Taking into account Eqs. (12) and (13), (14) can be re-written as Note that Proposition 2 is valid for weighted power means. They are the proper weighted quasiarithmetic means that allow us to generalize the distance between preference-approvals introduced by Erdamar et al. (2014).
In Proposition 2, we have shown that D r λ = d λ when r = 1. We now show that is not true if r = 1.

Clustering tasks
This section shows how the proposed distance can be used to study the universe of preferenceapprovals and to determine clusters.
Section 4.1 examines the universe of preference approvals in the case of two alternatives in order to observe how the values of r and λ affect the creation of homogeneous clusters. Afterwards, the influence of the two parameters r and λ when the number of alternatives n varies is investigated. Section 4.2 provides an application on real data, to investigate how the countries of the European Union can be clustered into groups, according to their preference-approvals on nine alternatives concerning social values. The dataset used comes from the Eurobarometer website 5 .
Increasing the value of λ emphasizes the discordance in the preference part, and modifies the relationships between the corresponding preference-approvals. Indeed, when λ = 0.75, there is an increase in the intensity of the distances at the top-right hand side of the graph, which concerns the triples The hierarchical relationship between objects is reported in Fig. 4; the dendrograms show how the hierarchical clustering of the eight preference-approvals changes based on D r λ . Figure 4 shows that the value of λ strongly influences the hierarchical aggregation of preference-approvals.
A similar analysis can be carried out by varying the value of r . In Fig. 5 the distances between the corresponding preference-approvals, for r = 2 and λ = 0.5 are shown.
for all (R 1 , A 1 ), (R 2 , A 2 ) ∈ R(X ). This is due to h being increasing in r . The dendrograms between preference-approvals objects are reported in Fig. 6.

Fig. 6
Hierarchical clustering dendrogram for 2 alternatives, r = 1 (left), r = 2 (right) and λ = 0.5 Table 6 Cophenetic dendrogram correlations for n = 2, r = 1, 1.5, 2, 5, 10 and λ = 0.5 n = 2 1 1 .  Figure 6 shows that an increase in r contributes differently (with respect to an increase in λ) to the change of the hierarchical aggregation structure. In fact, the two dendrograms merge preference-approvals in the same way. What changes is the "height" at which there is the aggregation or, in other words, the distance to be tolerated to aggregate two preferenceapprovals. Note that this happens only for two alternatives.
Tables 6, 7, 8 and 9 show that dendrogram correlations are strictly related to the values of r and n. Overall, the correlations between dendrograms tend to decrease as r increases. This is especially evident when we examine the first column of each table, which reports the correlation between dendrograms obtained with r = 1 and dendrograms obtained with r = 1.5, 2, 5, 10. In terms of the number of alternatives, it should be noted that as n increases, the dendrogram correlations generally decrease with an oscillatory trend.
In other words, Tables 6, 7, 8 and 9 highlight that the parameter r has a considerable influence, not only on the resulting values of the proposed distance D r λ , but also on the cluster structure discovered among the observations of the preference-approvals universe. Specifically, as n increases and the expressiveness of the voters explodes (Table 1), so does the discriminating power of r , allowing different clustering structures to be highlighted. Indeed, the proposed family of distances D r λ is more flexible than the existing one, and it ultimately comes down to a new parameter that can be exploited in various applications, such as maximizing the goodness of a clustering procedure.
To explore further this issue, let us consider a simulation study on the universe of 5 alternatives, which involves three steps: • generate four groups of clustered preference-approvals; • apply a hierarchical clustering algorithm for different values of r .
• compute an external validation index, the Adjusted Rand index (Hubert & Arabie, 1985), to investigate which value of r maximises the similarity between the estimated and the theoretical clusters.
Therefore, we aim to find the value of r that provides more reliable clusters, i.e. clusters that are more consistent with the data-generating process. The number of preference-approvals (on five alternatives) generated within each cluster was determined by randomly drawing four values from a normal distribution N (50, 4) and converting them into integer numbers.
Orderings and approvals were generated individually and merged to produce the final set of preference-approvals. Specifically, orderings within each sub-partition were generated from a Mallows Model (Mallows, 1957), which was one of the earliest probability models suggested for rankings and it is still widely used in theoretical and applied research. It is an exponential model defined by a central permutation σ 0 and a dispersion parameter θ . When θ = 0, σ 0 represents the mode of the distribution, i.e., the permutation with the highest probability of being generated. The probability of any other ranking decays exponentially with increasing distance to the central permutation. The dispersion parameter θ controls the steepness of this decline. The θ values for our simulation studies are {0, 0.5, 1, 1.5, 2}. Assuming that σ is a generic ranking, the probability for this ranking is function of θ , and it is given by: where d is a ranking distance measure and ψ(θ) is a normalization constant. We generated rankings assuming the Kemeny distance d K . The cluster central permutations, σ 0 , used in the analysis are reported in Table 10.
Approvals, within each cluster are generated from four multinomial distributions, with probability vectors, p ik , described in Table 11. Specifically, p ik is the probability to draw i approved alternatives into the k-th cluster.
After deriving clusters, the adjusted Rand index (Hubert & Arabie, 1985) is used to assess their goodness. The adjusted Rand index is a measure of the similarity between two set of clusterings; it is the corrected-for-chance version of the Rand index (Rand, 1971). The correction uses the predicted similarity of all pair-wise comparisons between clusterings described by a random model to generate a baseline. Although the Rand Index can only provide values between 0 and +1 (0 when the two data clusterings do not agree on any pair of points, and 1 when data clusterings are exactly the same), the modified Rand Index can return negative values if the index is lower than the expected similarity of all pair-wise comparisons between clusterings specified by a random model.
The results (Table 12) are obtained by averaging the adjusted Rand index over ten randomly generated datasets for each value of θ . Table 12 shows that, except for the case θ = 1.5, our measure D r λ with r = 1 results in higher average adjusted Rand indices. Thus, r = 1 allows the true clustered structure of data to be found more accurately and provides more accurate clusters.

A real data application
This subsection shows how the proposed metric can be used to perform cluster analysis on real data retrieved from the Eurobarometer website.
Since 1973, Eurobarometer has undertaken a series of public opinion polls on behalf of the European Commission and other European Union (EU) institutions. These polls cover a wide range of topics concerning the EU and its member countries. The data utilized in these analyses are specifically from question Q5 of the poll titled "Defending Democracy, Empowering citizens. Public Opinion at the legislature's midpoint" 7 .
A group of voters, divided by countries, was asked to indicate which of the following values should the European Parliament defend as a matter of priority: • x 1 : Equality between women and men.
• x 2 : The fight against discrimination and for the protection of minorities. • x 3 : Tolerance and respect for diversity in society. • x 4 : Solidarity between EU Member States and between its regions. • x 5 : Solidarity between the EU and poor countries in the world. • x 6 The protection of human rights in the EU and worldwide. • x 7 : Freedom of religion and belief. • x 8 : Freedom of movement. • x 9 : Freedom of speech and thought.
As a result, data are stored in a table (see Table 14) with 27 rows (one row for each EU member country) and 9 columns (each column representing an alternative of X = {x 1 , . . . , x 9 }). The total number of votes cast by the i-th country in favor of the j-th alternative is shown in the table's generic cell i j.
In order to transform the original table into a set of preference-approvals, preferences and approvals need to be derived. For each country, the alternatives are ranked in order of popularity, beginning with the one that received the most votes and ending with the one that received the fewest. Furthermore, in order to generate a vector of approvals, those alternatives that received more votes than the national average were deemed acceptable.
To run the cluster analysis, the distance matrix 27 × 27 was constructed using Eq. (11). All the alternatives seem important in this example, so a distinction between acceptable and unacceptable alternatives should not be interpreted as a distinction between valuable and not valuable, but instead as a distinction between more and less urgent. For this reason, λ = 0.75 was chosen in order to emphasize preference differences more than approvals.
A cluster-wise measure of cluster stability (Hennig, 2007) is used to jointly discover the optimal value of r and the optimal number of clusters k. Stability refers to the property of a meaningful and valid cluster that does not change easily when the data set is perturbed in a non-essential way. That is, when applied to many datasets collected from the same data distribution, a reliable clustering method should produce similar partitions. The cluster stability method (Hennig, 2007) employs three steps: 1. use various strategies to resample new data sets from the original and apply the hierarchical clustering method to each of them; 2. for every given original cluster, find the most similar cluster using the Jaccard coefficient (Jaccard, 1901) in the new data set and record the similarity value; 3. assess the cluster stability of every single cluster by the mean similarity taken over the resampled data sets.
The average cluster-wise stability is shown in Fig. 7 as a function of r (for k = 2, 3, 4 clusters). The procedure suggests that the most stable cluster configuration is k = 2 and r = 2. It is worth noting that, regardless of the value of k, r > 1 always leads to improved cluster stability. Indeed, with two clusters (k = 2) the value of r that maximizes stability is r = 2. Whereas with three or four clusters, the optimal solution is r = 4. In addition, as the number of clusters k increases, the average stability decreases. For several reasons, stability is a particularly relevant cluster validation measure in this example for determining the best value of r . First, it is not possible to use external validation measures in this case as the true clustered structure of the EU countries is not known. At the same time, most internal validation measures employ the distance between observations (D r λ ) to assess the goodness of clusters. However, this may be an issue in our instance since the distance between observations (D r λ ) is influenced by r . Therefore, to determine which value of r yields more accurate clusters, a metric that is independent of r is desirable. Furthermore, cluster stability has been examined both theoretically and practically (Hennig, 2007;Von Luxburg, 2010;Ullmann et al., 2022), and it has been shown to be capable of distinguishing between meaningful stable and spurious clusters. Figures 8 and 9 show the resulting dendrogram and clusters, respectively, obtained with k = 2 and r = 2.
The clustering procedure suggests that the EU countries can be separated into two large groups. Cluster 1 is mainly made up of Western European countries, whereas Cluster 2 of Eastern European countries.
To provide a more in-depth picture of how the EU countries express their views on the nine alternatives proposed, the two preference-approvals that represent the two clusters, that we call representative preference-approvals, are shown in Eq. (18).
To obtain the representative preference-approvals that summarize each cluster, preferences and approvals need to be aggregated. In each cluster, the set of preferences is combined into a unique weak order by deriving the average position for each alternative and ranking them according to it. Note that this aggregation method is equivalent to the the Borda count (Borda, 1781) extended to weak orders (see Smith, 1973;Black, 1976, Cook & Seiford, 1982. In our example, the extended Borda count assigns a score to each alternative, for each country, the number of alternatives ranked below plus half of the number of alternatives that are indifferent to it: Similarly, the set of approvals are combined into a unique approval vector by taking the average approval for each alternative, and then considering those alternatives whose average approval is greater than 0.5 as approved.
Cluster 1 x Cluster 2 It is worth noting that x 6 and x 9 , namely, "The protection of human rights in the EU and worldwide" and "Freedom of speech and thought", respectively, are above the approval line in the two representative preference-approvals, indicating that they can be considered very urgent. Regarding x 1 , that is "Equality between women and men", it is ranked at the top of the representative preference-approval of Cluster 1, while it is just below the approval line in The bold values are the smallest distances between each country and the two clusters the Cluster 2 representative preference-approval. Similarly, x 4 , that is "Solidarity between the EU Member States and between its regions", is ranked fourth (above the approval line) in Cluster 2. Still, it is the first alternative below the approval line in Cluster 1. Furthermore, Cluster 2 prioritizes x 8 , that is "Freedom of movement", which is at the end of the preferenceapproval of Cluster 1. Finally, in both the two representative preference-approvals, x 7 , that is "Freedom of religion and belief", is ranked last. Table 15 reports the D 2 0.75 distances of each country to the representative cluster preference-approvals.
It should be noted that, except for Greece, each country is closer to the preference-approval of its own cluster than the other. Despite being reasonable, this result is not trivial since the technique for obtaining the cluster preference-approval does not involve D r λ . Some countries can be considered central in their clusters as they are very close to the representative preference-approval, e.g. Belgium (0.092), Austria (0.096), Malta (0.096) for Cluster 1, and the Czech Republic (0.036), Lithuania (0.094), Hungary (0.072) and Slovenia (0.105) for Cluster 2. As a rule of thumb, the greater the distance from the own cluster preference-approval, the more the country disagrees with the other countries in its cluster. Finally, it is worth noting that some countries, such as Ireland, Italy and Greece, are located in the middle of the two clusters, as they have a similar distances to the two cluster preferenceapprovals.

Concluding remarks
In social choice theory, preference rankings and approvals are two popular ways to collect the preferences of a group of agents on a set of alternatives. In the preference-approval setting, each agent, in addition to ordering a set of alternatives from best to worst, submits a cut-off line to distinguish between acceptable and unacceptable. Within this framework, in this paper, we propose a new distance for preference-approvals, following the approach of the Kemeny distance.
Given two preference-approvals and two alternatives, we introduce two indices that measure the discordances between these alternatives with respect to preference and approvals, and an aggregation function belonging to the class of weighted power means to define a new distance. This new distance depends on two parameters. The effect of these parameters on the distance is analyzed and described through some heatmaps. The proposed distance can be used to study the universe of preference-approvals and to determine clusters of voters: how the two parameters characterizing the distances affect the clustering process is shown with some dendrograms and by the cophenetic correlations among them. We have shown that the new distance family offers some advantages compared to the existing distance function. Specifically, through a simulation study and the adjusted Rand index, we have proved that D r λ with r = 1 allows the true clustered structure of data to be found more accurately. Similarly, through a cluster-wise stability index, we have shown that D r λ with r = 1 produces more stable clusters on the real data example.
In future work, axiomatizing the new family of distance functions might prove important. Moreover, future research should examine consensus measures based on distances between preference-approvals (see Erdamar et al., 2014), algorithms to determine representative preference-approvals efficiently (see D'Ambrosio, 2017), clustering on alternatives (see González del Pozo et al., 2017), and also reaching consensus processes (see Palomares et al., 2014;Chao et al., 2021, among others).
Acknowledgements The authors are grateful to two anonymous reviewers for their useful comments and suggestions, and also to the Spanish Agencia Estatal de Investigación (project PID2021-122506NB-I00) and the University of Palermo (projects: FFR_D16_PLAIA and FFR_D16_SCIANDRA) for their financial support.
Funding Open access funding provided by Università degli Studi di Palermo within the CRUI-CARE Agreement.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. λ · ( p i j ) r + (1 − λ) · (a i j ) r 1 r + λ · ( p i j ) r + (1 − λ) · (a i j ) r 1 r r .