Positional dominance in networks

Theoretical Background

A network can be described as a dyadic variable x ∈ 𝒲^𝒟, where 𝒲 is the value range of the network (in the simple case of unweighted networks 𝒲 = {0, 1}) and 𝒟 = 𝒩 × 𝒜 describes the dyadic domain of actors 𝒩 and affiliations 𝒜. If 𝒜 ≠ 𝒩, we obtain a two-mode network and if 𝒜 = 𝒩 a one-mode network consisting of relations among actors.

Definition

Let x ∈ 𝒲^𝒟 be a network and i, j ∈ 𝒩. We say that i is dominated by j under the total homogeneity assumption, denoted by i ≤ j if x_it ≤ x_jt  ∀t ∈ 𝒩. If there exists a permutation π : 𝒩 → 𝒩 such that x_it ≤ x_jπ(t)  ∀t ∈ 𝒩, we say that i is dominated by j under the total heterogeneity assumption, denoted by i ⪯ j.
It holds that i ≤ j ⟹ i ⪯ j but not vice versa.

More about the positional dominance and the positional approach to network analysis can be found in

Brandes, Ulrik. (2016). Network Positions. Methodological Innovations, 9, 2059799116630650. (link)

Positional Dominance in the netrankr Package

library(netrankr)
library(igraph)
library(magrittr)
set.seed(1886) #for reproducibility

The function positional_dominance can be used to check both, dominance under homogeneity and heterogeneity. In accordance with the analytic pipeline of centrality we use the %>% operator.

data("dbces11")
g <- dbces11

#neighborhood inclusion can be expressed with the analytic pipeline
D <- g %>% indirect_relations(type="adjacency") %>% positional_dominance()

More on the indirect_relations() function can be found in this vignette.

The map parameter of positional_dominance allows to distinguish between dominance under total heterogeneity (map=FALSE) and total homogeneity (map=TRUE). In the later case, all relations can be ordered non-decreasingly (or non-increasingly if the relation describes costs, such as distances) and afterwards checked front to back. Dominance under total homogeneity yields a ranking, if the relation is binary (e.g. adjacent or not).

D <- g %>% 
  indirect_relations(type="adjacency") %>% 
  positional_dominance(map=TRUE)

comparable_pairs(D)

## [1] 1

For cost variables like shortest path distances, the benefit parameter is set to FALSE.

D1 <- g %>% 
  indirect_relations(type="dist_sp") %>% 
  positional_dominance(map=FALSE,benefit=FALSE)

From the definition given in the first section, it is clear that there are always at least as many comparable pairs under the total homogeneity assumption as under total heterogeneity.

D1 <- g %>% 
  indirect_relations(type="dist_sp") %>% 
  positional_dominance(map=FALSE,benefit=FALSE)

D2 <- g %>% 
  indirect_relations(type="dist_sp") %>% 
  positional_dominance(map=TRUE,benefit=FALSE)

c("heterogeneity"=comparable_pairs(D1),
  "homogeneity"=comparable_pairs(D2))

## heterogeneity   homogeneity 
##     0.1636364     0.8727273

Additionally, all dominance relations from the heterogeneity assumption are preserved under total homogeneity. (Note: A ⟹ B is equivalent to ¬A ∨ B)

all(D1!=1 | D2==1) 

## [1] TRUE