Differences

This shows you the differences between two versions of the page.

--- r:social_network_analysis [2023/06/12 03:08] – [Hawthorne study] hkimscil
+++ r:social_network_analysis [2024/11/14 09:02] (current) – [Hawthorne study] hkimscil
@@ Line 8: / Line 8: @@
 {{:r:davis.women.club.csv}}
 <code>
-library(tidiverse)
+# install.packages(c("igraph", "tidyverse"))
+library(igraph)
+library(tidyverse)
 sd <- read.csv("http://commres.net/wiki/_media/r/davis.women.club.csv")
 head(sd)
@@ Line 27: / Line 30: @@
 <code>
+V(g)$type <- bipartite_mapping(g)$type
 types <- V(g)$type
 deg <- degree(g)
@@ Line 178: / Line 182: @@
 >
 </code>
+===== stu x class 처럼 분석한 예 =====
 <code>
 actor_matrix <- bipartite_matrix %*% t(bipartite_matrix)
 event_matrix <- t(bipartite_matrix) %*% bipartite_matrix
 diag(actor_matrix) <- 0
 actor_matrix
+actor_matrix_cff_2 <- ifelse(actor_matrix > 2, actor_matrix, 0) # cuttoff 3 below
+actor_matrix_cff_3 <- ifelse(actor_matrix > 3, actor_matrix, 0) # cuttoff 3 below
 actor_g <- graph_from_adjacency_matrix(actor_matrix,
                                        mode = "undirected",
                                        weighted = TRUE)
+actor_g_cff_2 <- graph_from_adjacency_matrix(actor_matrix_cff_2,
+                                             mode = "undirected",
+                                             weighted = TRUE)
+actor_g_cff_3 <- graph_from_adjacency_matrix(actor_matrix_cff_3,
+                                             mode = "undirected",
+                                             weighted = TRUE)
+V(actor_g)$size <- betweenness(actor_g)
+V(actor_g_cff_2)$size <- betweenness(actor_g_cff_2)
+V(actor_g_cff_3)$size <- betweenness(actor_g_cff_3)
+V(actor_g)$label.cex <- betweenness(actor_g) * 0.2
+V(actor_g_cff_2)$label.cex <- betweenness(actor_g_cff_2) * 0.1
+V(actor_g_cff_3)$label.cex <- betweenness(actor_g_cff_3) * 0.4
 actor_g
+actor_g_cff_2
+actor_g_cff_3
 event_g <- graph_from_adjacency_matrix(event_matrix,
@@ Line 195: / Line 220: @@
                                        weighted = TRUE)
 event_g
-# V(actor_g)$size <- degree(actor_g)
-V(actor_g)$label.cex <- actor_g * 0.1
 windowsFonts(d2coding = windowsFont("D2Coding"))
 windowsFonts(lucida = windowsFont("Lucida Console"))
-windowsFonts(courrier = windowsFont("Courrier New"))
 shape <- c("circle", "square")
 fnts <- c("d2coding", "lucida")
-plot(actor_g, layout = layout_with_graphopt,
+plot(actor_g,
      vertex.shape= shape[as.numeric(V(g)$type) + 1],
-     vertex.label.family= fnts[as.numeric(V(g)$type)+1]
+     vertex.label.family= fnts[as.numeric(V(g)$type)+1],
-)
+     edge.color="red", edge.width=3
+)
+plot(actor_g_cff_2,
+     vertex.shape= shape[as.numeric(V(g)$type) + 1],
+     vertex.label.family= fnts[as.numeric(V(g)$type)+1],
+     edge.color="red", edge.width=3
+)
+plot(actor_g_cff_3,
+     vertex.shape= shape[as.numeric(V(g)$type) + 1],
+     vertex.label.family= fnts[as.numeric(V(g)$type)+1],
+     edge.color="red", edge.width=3
+)
 </code>
-{{:r:pasted:20230612-025234.png}}
+[{{:r:pasted:20230612-081851.png|actor_g}}] \\
+[{{:r:pasted:20230612-082040.png|actor_g_cff_2}}] \\
+[{{:r:pasted:20230612-082055.png|actor_g_cff_3}}]
+===== 다른 방법 =====
 <code>