# Social Network Analysis

## Definition

Social networks are simply networks of social interactions and personal
relationships. Think about our group of friends and how we got to know
them. Maybe we met them while ago from our schooling, or maybe we met
them through a hobby or through our community. In fact, 72% of all
Internet users are active on social media today, including in social
interactions and developing personal relationships. However, to
understand about social networks, we only don’t need to go through
internet or social media, they may come in variety of form in our daily
life.

## Social Network Analysis

Social network analysis (SNA), also known as network science, is a field
of data analytics that uses networks and graph theory to understand
social structures. We can see network around us such as road network,
Learning SNA allow us to explore insight of various data sources.

## SNA Graph

We all are familier about the graph which simply the collection of non
zero vertex and edge. In order to build SNA graphs, two key components
are required actors and relationships. Here, actors represents the vertex
and relationship means the edge between two actors. Let us write SNA
graph in R code. To do this, we should have igraph package
already install in our R or R studio.

library(igraph)
##
## Attaching package: 'igraph'

## The following objects are masked from 'package:stats':
##
##     decompose, spectrum

## The following object is masked from 'package:base':
##
##     union
g <- graph(c(1,2))
plot(g) In figure, we can see the directed graph from node 1 to node 2. From above
graph, it is also clear that by default it gives directed graph. In
above figure we can not clearly see node and edge so let’s increase its
size and give different color to the node.

plot(g,
vertex.color = 'green',
vertex.size = 40,
edge.color ='red',
edge.size = 20) Now we change the node color and node font size. Add more node on graph for
this we follow the following code.

g <- graph(c(1,2,2,3,3,4,4,1))
plot(g,
vertex.color = 'green',
vertex.size =40,
edge.color = 'red',
edge.size = 20) We got the graph with four vertex 1,2,3,4. Here, we also got directed graph.
Let’s make it for this we need to write directed is false.

g <- graph(c(1,2,2,3,3,4,4,1),directed = FALSE)
plot(g,
vertex.color = 'green',
vertex.size =40,
edge.color = 'red',
edge.size = 20) We got our desirable type of graph. Now let’s move forward. We can give the
number of vertex with out writing them for this look following code.

g <- graph(c(1,2,2,3,3,4,4,1),
directed = F, n=7)
plot(g,
vertex.color = 'green',
vertex.size =40,
edge.color = 'red',
edge.size = 20) In above graph we have given seven vertex numbers. Among them we can see three isolated
nodes. The reason to come isolated node is that we did not specify the
edge or relationship between them. Also from graph we can make sense
that this is not directed graph.

Let's see what will happen if we type only g[].

g[]
## 7 x 7 sparse Matrix of class "dgCMatrix"
##
## [1,] . 1 . 1 . . .
## [2,] 1 . 1 . . . .
## [3,] . 1 . 1 . . .
## [4,] 1 . 1 . . . .
## [5,] . . . . . . .
## [6,] . . . . . . .
## [7,] . . . . . . .

This gives us the 7 by 7 adjacency matrix. Adjacency matrix is matrix
where we give 1 if there is edge between two vertex if not then we
give 0. But in above matrix it gave simply . instead of zero.

Let us try to build our graph by keeping text data in place of number.

g1 <-
graph(c("Binu","Binita","Binita","Rita"
,"Rita","Binu","Binu","Rita", "Anju",
"Binita"))
plot(g1,
vertex.color = "green",
vertex.size = 40,
edge.color = "red",
edge.size = 5) If we want to check the features of g1 we simply type g1 and press control and
enter key. We got following output.

g1
## IGRAPH d030509 DN-- 4 5 --
## + attr: name (v/c)
## + edges from d030509 (vertex names):
##  Binu  ->Binita Binita->Rita   Rita  ->Binu   Binu  ->Rita   Anju  ->Binita

It showed that in graph there are 4 nodes 5 edges. And edges are directed
from
Binu ->Binita, Binita->Rita, Rita ->Binu, Binu ->Rita, Anju ->Binita.

## Degree

Degree means numbers of connection to each node.
Let’s check the degree of graph g1. To check degree we can do
degree(g1) or degree(g1, mode='all').

degree(g1) 
##   Binu Binita   Rita   Anju
##      3      3      3      1

Degree of Binu is 3 similarly Anju has degree 1.

degree(g1, mode='all')
##   Binu Binita   Rita   Anju
##      3      3      3      1

## Diameter

Diameter means number of edged inside and outside
of SND. Now, let’s check the diameter of graph.

diameter(g1, directed = F, weights =
NA)
##  2

We got two diameter. i.e Anju <- Binita <- Rita, Anju <- Binita <- Binu.

## Edge density

Edge density means ecount(g1)/(vcount(g1)*(vcount(g1) -1)). We can
calculate it from following code.

edge_density(g1, loops = F)
##  0.4166667

We got the value of edge density 0.4166667.

## Reciprocity

Total edges = 5

Tied edges = 2

Reciprocity = 2/5 = 0.4

reciprocity(g1)
##  0.4

## Closeness

Now getting closeness of graph.

closeness(g1, mode = "all", weights = NA)
##      Binu    Binita      Rita      Anju
## 0.2500000 0.3333333 0.2500000 0.2000000

From above result we can see that Binita is closest to the others three
persons and Anju is farthest from other three persons.

## Betweenness

Let’s calculate between of g1

betweenness(g1, directed = T, weights = NA)
##   Binu Binita   Rita   Anju
##      1      2      2      0

Binita and Rita has 2 inner edge similarly Binu has one inner edge and
Anju has no inner edge.

## Edge Betweenness

For every pair of vertices in a connected graph, there exists at least
one shortest path between the vertices.

edge_betweenness(g1, directed = T, weights = NA)
##  2 4 4 1 3

Here I have loadedTwitterdata from my local machine.

load("F:/MDS R/termDocMatrix.rdata")
m<- as.matrix(termDocMatrix)
termM <- m %*% t(m)
termM[1:10,1:10]
##               Terms
## Terms          analysis applications code computing data examples introduction
##   analysis           23            0    1         0    4        4            2
##   applications        0            9    0         0    8        0            0
##   code                1            0    9         0    1        6            0
##   computing           0            0    0        10    2        0            0
##   data                4            8    1         2   85        5            3
##   examples            4            0    6         0    5       17            2
##   introduction        2            0    0         0    3        2           10
##   mining              4            7    3         1   50        5            3
##   network            12            0    1         0    0        2            2
##   package             2            1    0         2   12        2            0
##               Terms
## Terms          mining network package
##   analysis          4      12       2
##   applications      7       0       1
##   code              3       1       0
##   computing         1       0       2
##   data             50       0      12
##   examples          5       2       2
##   introduction      3       2       0
##   mining           64       1       6
##   network           1      17       1
##   package           6       1      27

Now we have built a term-term adjacency matrix, where the rows and
columns represents terms, and every entry is the number of
co-occurrences of two terms. Next we can build a graph with
graph.adjacency() from package igraph.

g <- graph.adjacency(termM,weighted = T,mode = 'undirected')
g
## IGRAPH d066e18 UNW- 21 151 --
## + attr: name (v/c), weight (e/n)
## + edges from d066e18 (vertex names):
##   analysis    --analysis     analysis    --code
##   analysis    --data         analysis    --examples
##   analysis    --introduction analysis    --mining
##   analysis    --network      analysis    --package
##   analysis    --positions    analysis    --postdoctoral
##  analysis    --r            analysis    --research
##  analysis    --series       analysis    --slides
##  analysis    --social       analysis    --time
## + ... omitted several edges

Here we have built graph on termDocMatrix. In result we can see the edges
between different dodes.

g <- simplify(g)
g
## IGRAPH d06a87e UNW- 21 130 --
## + attr: name (v/c), weight (e/n)
## + edges from d06a87e (vertex names):
##   analysis    --code         analysis    --data
##   analysis    --examples     analysis    --introduction
##   analysis    --mining       analysis    --network
##   analysis    --package      analysis    --positions
##   analysis    --postdoctoral analysis    --r
##  analysis    --research     analysis    --series
##  analysis    --slides       analysis    --social
##  analysis    --time         analysis    --tutorial
## + ... omitted several edges

Function simplify() in igraph handly removes self-loops from a network.
We can see in previous graph there are 151 edges in second graph there
are only 130 edges. Hence there are 21 self loops they are omitted from
the graph.

## Check Degree of graph and nodes of the graph.

V(g)$label <- V(g)$name
V(g)$label ##  "analysis" "applications" "code" "computing" "data" ##  "examples" "introduction" "mining" "network" "package" ##  "parallel" "positions" "postdoctoral" "r" "research" ##  "series" "slides" "social" "time" "tutorial" ##  "users" V(g)$degree <- degree(g)
V(g)$degree ##  17 6 9 9 18 14 12 20 14 13 8 7 8 17 9 11 15 11 11 16 15 We found degree of graph is 20. ## Histogram on the basis of degree hist(V(g)$degree, breaks = 100,col = 'green', main = "Frequency Of Degree",
xlab = " Degree of vertices", ylab = " Frequency") We can see most of nodes have degree 9 and 11. We all khow what is degree of
graph, number of edges that are incident to the node is called degree of
graph.

## Let’s Plot Graph on the Data.

Let's set set.seed(3952). Where set.seed gives same sample with the same seed value.

set.seed(3952)

layout1 <- layout.fruchterman.reingold(g)

plot(g, layout=layout1) A different layout can be generated with the first line of code below. The
second line produces an interactive plot, which allows us to manually
rearrange the layout

plot(g, layout=layout.kamada.kawai) ## Make it better.

V(g)$label.cex <- 2.2 * V(g)$degree / max(V(g)$degree)+ .2 V(g)$label.color <- rgb(0, 0, .2, .8)

V(g)$frame.color <- NA egam <- (log(E(g)$weight)+.4) / max(log(E(g)$weight)+.4) E(g)$color <- rgb(.5, .5, 0, egam)

E(g)$width <- egam # plot the graph in layout1 plot(g, layout=layout1) Here size of words appear according to their degree. From graph we can clearly see that node mining has maximum degree. ## Community detection Communities are seen as groups, clusters, coherent subgroups, or modules in different fields; community detection in a social network is identifying sets of nodes in such a way that the connections of nodes within a set are more than their connection to other network nodes. comm <- cluster_edge_betweenness(g) ## Warning in cluster_edge_betweenness(g): At community.c:461 :Membership vector ## will be selected based on the lowest modularity score. ## Warning in cluster_edge_betweenness(g): At community.c:468 :Modularity ## calculation with weighted edge betweenness community detection might not make ## sense -- modularity treats edge weights as similarities while edge betwenness ## treats them as distances plot(comm,g) There are dense connection within the group within the community the connection is sparse. prop <- cluster_label_prop(g) plot(prop, g) This is different type of algorithms to detect community which is different from previous one. ## Hubs Nodes with most outer edges. We need to find hub score. hs <- hub_score(g,weights = NA)$vector
hs
##     analysis applications         code    computing         data     examples
##    0.9047777    0.3589289    0.5606314    0.5223206    0.9065063    0.8195307
## introduction       mining      network      package     parallel    positions
##    0.7307838    1.0000000    0.7483791    0.7210610    0.4939614    0.3733995
## postdoctoral            r     research       series       slides       social
##    0.4095660    0.9147530    0.4481802    0.6761093    0.8510808    0.6018664
##         time     tutorial        users
##    0.6761093    0.8899001    0.8342594

## Authority

Nodes with most inner edges. We need to find authority score

as <- authority_score(g, weights = NA)\$vector
as
##     analysis applications         code    computing         data     examples
##    0.9047777    0.3589289    0.5606314    0.5223206    0.9065063    0.8195307
## introduction       mining      network      package     parallel    positions
##    0.7307838    1.0000000    0.7483791    0.7210610    0.4939614    0.3733995
## postdoctoral            r     research       series       slides       social
##    0.4095660    0.9147530    0.4481802    0.6761093    0.8510808    0.6018664
##         time     tutorial        users
##    0.6761093    0.8899001    0.8342594

## Hub in Plot

par(mfrow = c(1,2))
plot(g,vertex.size= hs*50, main = "Hubs",
vertex.label = NA,
vertex.colour = rainbow(50)) ## Authority in Plot

plot(g,vertex.size= as*30, main = "Authorities",
vertex.label = NA,
vertex.colour = rainbow(50)) Hub are expected to contain large number of outgoing link. And authority are
expected to contain large number of incoming link from hubs.

## Application of SNA in Real World

Social network analysis can provide information about the reach of
gangs, the impact of gangs, and gang activity. The approach may also
allow we to identify those who may be at risk of gang-association and/or
being exploited by gangs.

Scroll to top