Comment créer des attributs de sommet basés sur les attributs des voisins des sommets dans igraph?

Question

Par exemple, comment calculer, pour chaque sommet, le pourcentage de liens dirigés vers l'extérieur vers les mâles?

g <- erdos.renyi.game(20, .3, type=c("gnp"), directed = TRUE) 
V(g)$male <- rbinom(20,1,.5) 
V(g)$male[10] <- NA 

Answer 1

Une possible solution (pas nécessaire optimal) est la suivante (ce qui est une ligne unique, je romps il suffit de faire par souci de lisibilité):

unlist(lapply(get.adjlist(g, mode="out"), 
     function (neis) { 
      sum(V(g)[neis]$male, na.rm=T) 
     } 
))/degree(g, mode="out") 

Maintenant, nous allons le diviser en plus petits pièces. Tout d'abord, nous obtenons la liste d'adjacence du graphique en utilisant get.adjlist(g, mode="out"). Cela vous donne une liste de vecteurs, chaque vecteur contenant les voisins extérieurs d'un sommet. Ensuite, nous appliquons une fonction à chaque vecteur de cette liste en utilisant lapply. La fonction appliquée est la suivante:

function (neis) { 
    sum(V(g)[neis]$male, na.rm=T) 
} 

La fonction prend simplement les voisins d'un noeud dans neis et l'utilise pour sélectionner un sous-ensemble de sommets à partir de l'ensemble de vertex définir V(g). Ensuite, l'attribut male est récupéré pour ce sous-ensemble de vertex et les valeurs sont additionnées, en supprimant les valeurs NA à la volée. Essentiellement, cette fonction vous donne le nombre de mâles dans neis.

Maintenant, revenons à notre expression originale, nous avons appliqué cette fonction à la liste d'adjacence du graphe en utilisant lapply, obtenant une liste de nombres, chaque nombre contenant le nombre de voisins mâles d'un sommet donné. Nous convertissons cette liste en un seul vecteur R en utilisant unlist et le divisons élément par les out-degrés des sommets pour obtenir les ratios.