library(igraph)
Last compiled on oktober, 2025
In class
# install.packages('formatR')
numbers <- sample(x = (0:1), size = 16, replace = T)
net1 <- matrix(data = numbers, nrow = 4, ncol = 4)
diag(net1) <- NA
net1
#> [,1] [,2] [,3] [,4]
#> [1,] NA 0 1 0
#> [2,] 1 NA 1 0
#> [3,] 1 1 NA 1
#> [4,] 1 0 1 NA
mean(rowSums(net1, na.rm = T))
#> [1] 2
density <- (sum(rowSums(net1, na.rm = T)))/16
density
#> [1] 0.5
# reciprocity
net1 <- (t(net1) + net1)
num_twos <- sum(net1 == 2, na.rm = TRUE)
num_ones_twos <- sum(net1 %in% c(1, 2), na.rm = TRUE)
reciprocity <- num_twos/num_ones_twos
reciprocity
#> [1] 0.6
set.seed(123643)
net_ex <- matrix(sample(0:1, 16, replace = T), nrow = 4, ncol = 4)
diag(net_ex) <- 0
net_ex_un <- net_ex + t(net_ex)
net_ex_un[net_ex_un == 2] <- 1
diag(net_ex_un) <- 0
net_ex_un
#> [,1] [,2] [,3] [,4]
#> [1,] 0 1 1 1
#> [2,] 1 0 1 0
#> [3,] 1 1 0 1
#> [4,] 1 0 1 0
netG <- graph_from_adjacency_matrix(net_ex_un)
transitivity(netG)
#> [1] 0.75
plot <- plot(netG)

net_ex <- (t(net_ex) + net_ex)
num_twos <- sum(net_ex == 2, na.rm = TRUE)
num_ones_twos <- sum(net_ex %in% c(1, 2), na.rm = TRUE)
reciprocity <- num_twos/num_ones_twos
reciprocity
#> [1] 0.6
# ?igraph
netG_d <- graph_from_adjacency_matrix(net_ex)
dyad_census(netG_d)
#> $mut
#> [1] 5
#>
#> $asym
#> [1] 0
#>
#> $null
#> [1] 1
triad_census(netG_d)
#> [1] 0 0 0 0 0 0 0 0 0 0 2 0 0 0 0 2
plot(netG_d)

Homework
Introduction
The principle of homophily refers to the larger likelihood of
observing a positive relationship between two people who are similar. In
the world of academic research, this principle also applies. For
example, we see that researchers with greater physical proximity or the
same gender are more likely to collaborate (Horta et al., 2022). Given
the different positions researchers can hold in universities, this
raises the question to what extent this same goes for researchers of
similar positions.There are several reasons why this might occur.
Researcher of similar posiitons might prefer to work with other
researchers with a similar level (selection), or working together with a
reserahcer of a higher level might influence a researcher to promote
(influence).
If this homophily structure exists, this raises the question what the
consequences are of this structure of collaberations on papers. Other
researchers might prefer to cite papers by researchers of a higher
position. Furthermore, working with researchers of a higher position
might influence these researchers in such a way that their future papers
will be cited more often.
#RQs
(Descriptive) RQ1: To what degree is there position homophily in
publication collaberations between social science researchers in the
Netherlands?
(Explanatory) RQ2: How does the position composition of a
researcher’s egonet (paper collaberations / department influence a
researcher’s citations?
(Explanatory) RQ3: How does the position composition of a paper
influence the number of citations of the paper?
Data needed
Node characteristics: position, citations, affiliation.
Relational attributes: who collaberated with whom, paper
citations.
LS0tDQp0aXRsZTogIldlZWsgMiINCiNiaWJsaW9ncmFwaHk6IHJlZmVyZW5jZXMuYmliDQpkYXRlOiAiMjAyNS0wOS0xMiINCi0tLQ0KDQpgYGB7ciwgZ2xvYmFsc2V0dGluZ3MsIGVjaG89RkFMU0UsIHdhcm5pbmc9RkFMU0UsIHJlc3VsdHM9J2hpZGUnfQ0KbGlicmFyeShrbml0cikNCg0Ka25pdHI6Om9wdHNfY2h1bmskc2V0KGVjaG8gPSBUUlVFKQ0Kb3B0c19jaHVuayRzZXQodGlkeS5vcHRzPWxpc3Qod2lkdGguY3V0b2ZmPTEwMCksdGlkeT1UUlVFLCB3YXJuaW5nID0gRkFMU0UsIG1lc3NhZ2UgPSBGQUxTRSxjb21tZW50ID0gIiM+IiwgY2FjaGU9VFJVRSwgY2xhc3Muc291cmNlPWMoInRlc3QiKSwgY2xhc3Mub3V0cHV0PWMoInRlc3QyIikpDQpvcHRpb25zKHdpZHRoID0gMTAwKQ0KcmdsOjpzZXR1cEtuaXRyKCkNCg0KDQoNCmNvbG9yaXplIDwtIGZ1bmN0aW9uKHgsIGNvbG9yKSB7c3ByaW50ZigiPHNwYW4gc3R5bGU9J2NvbG9yOiAlczsnPiVzPC9zcGFuPiIsIGNvbG9yLCB4KSB9DQoNCmBgYA0KDQpgYGB7ciBrbGlwcHksIGVjaG89RkFMU0UsIGluY2x1ZGU9VFJVRX0NCmtsaXBweTo6a2xpcHB5KHBvc2l0aW9uID0gYygndG9wJywgJ3JpZ2h0JykpDQoja2xpcHB5OjprbGlwcHkoY29sb3IgPSAnZGFya3JlZCcpDQoja2xpcHB5OjprbGlwcHkodG9vbHRpcF9tZXNzYWdlID0gJ0NsaWNrIHRvIGNvcHknLCB0b29sdGlwX3N1Y2Nlc3MgPSAnRG9uZScpDQpgYGANCg0KYGBge3J9DQpsaWJyYXJ5KGlncmFwaCkNCmBgYA0KDQpMYXN0IGNvbXBpbGVkIG9uIGByIGZvcm1hdChTeXMudGltZSgpLCAnJUIsICVZJylgDQoNCjxicj4NCg0KIyBJbiBjbGFzcw0KDQpgYGB7cn0NCiNpbnN0YWxsLnBhY2thZ2VzKCJmb3JtYXRSIikNCg0KbnVtYmVyczwtIHNhbXBsZSh4PSgwOjEpLCBzaXplID0gMTYsIHJlcGxhY2UgPSBUKQ0KbmV0MSA8LSBtYXRyaXgoZGF0YT1udW1iZXJzLCBucm93PTQsIG5jb2w9NCkNCmRpYWcobmV0MSkgPC0gTkENCm5ldDENCg0KbWVhbihyb3dTdW1zKG5ldDEsIG5hLnJtPVQpKQ0KDQpkZW5zaXR5IDwtIChzdW0ocm93U3VtcyhuZXQxLCBuYS5ybT1UKSkpLzE2DQoNCmRlbnNpdHkNCg0KDQojcmVjaXByb2NpdHkNCm5ldDEgPC0gKHQobmV0MSkgKyBuZXQxKQ0KDQpudW1fdHdvcyA8LSBzdW0obmV0MSA9PSAyLCBuYS5ybSA9IFRSVUUpDQoNCm51bV9vbmVzX3R3b3MgPC0gc3VtKG5ldDEgJWluJSBjKDEsIDIpLCBuYS5ybSA9IFRSVUUpDQoNCnJlY2lwcm9jaXR5IDwtIG51bV90d29zIC8gbnVtX29uZXNfdHdvcw0KcmVjaXByb2NpdHkNCg0KDQoNCnNldC5zZWVkKDEyMzY0MykNCm5ldF9leCA8LSBtYXRyaXgoc2FtcGxlKDA6MSwgMTYsIHJlcGxhY2UgPSBUKSwgbnJvdyA9IDQsIG5jb2w9NCkNCmRpYWcobmV0X2V4KSA8LSAwDQpuZXRfZXhfdW4gPC0gbmV0X2V4ICsgdChuZXRfZXgpDQpuZXRfZXhfdW5bbmV0X2V4X3VuID09IDJdIDwtIDENCmRpYWcobmV0X2V4X3VuKSA8LSAwDQpuZXRfZXhfdW4NCg0KbmV0RyA8LSBncmFwaF9mcm9tX2FkamFjZW5jeV9tYXRyaXgobmV0X2V4X3VuKQ0KdHJhbnNpdGl2aXR5IChuZXRHKQ0KcGxvdCA8LSBwbG90KG5ldEcpDQoNCg0KbmV0X2V4IDwtICh0KG5ldF9leCkgKyBuZXRfZXgpDQoNCg0KbnVtX3R3b3MgPC0gc3VtKG5ldF9leCA9PSAyLCBuYS5ybSA9IFRSVUUpDQoNCm51bV9vbmVzX3R3b3MgPC0gc3VtKG5ldF9leCAlaW4lIGMoMSwgMiksIG5hLnJtID0gVFJVRSkNCg0KcmVjaXByb2NpdHkgPC0gbnVtX3R3b3MgLyBudW1fb25lc190d29zDQpyZWNpcHJvY2l0eQ0KDQpgYGANCg0KYGBge3J9DQojP2lncmFwaA0KDQpuZXRHX2QgPC0gZ3JhcGhfZnJvbV9hZGphY2VuY3lfbWF0cml4KG5ldF9leCkNCg0KZHlhZF9jZW5zdXMobmV0R19kKQ0KdHJpYWRfY2Vuc3VzKG5ldEdfZCkNCg0KcGxvdChuZXRHX2QpDQpgYGANCg0KIyBIb21ld29yaw0KDQojIEludHJvZHVjdGlvbg0KDQpUaGUgcHJpbmNpcGxlIG9mIGhvbW9waGlseSByZWZlcnMgdG8gdGhlIGxhcmdlciBsaWtlbGlob29kIG9mIG9ic2VydmluZyBhIHBvc2l0aXZlIHJlbGF0aW9uc2hpcCBiZXR3ZWVuIHR3byBwZW9wbGUgd2hvIGFyZSBzaW1pbGFyLiBJbiB0aGUgd29ybGQgb2YgYWNhZGVtaWMgcmVzZWFyY2gsIHRoaXMgcHJpbmNpcGxlIGFsc28gYXBwbGllcy4gRm9yIGV4YW1wbGUsIHdlIHNlZSB0aGF0IHJlc2VhcmNoZXJzIHdpdGggZ3JlYXRlciBwaHlzaWNhbCBwcm94aW1pdHkgb3IgdGhlIHNhbWUgZ2VuZGVyIGFyZSBtb3JlIGxpa2VseSB0byBjb2xsYWJvcmF0ZSAoSG9ydGEgZXQgYWwuLCAyMDIyKS4gR2l2ZW4gdGhlIGRpZmZlcmVudCBwb3NpdGlvbnMgcmVzZWFyY2hlcnMgY2FuIGhvbGQgaW4gdW5pdmVyc2l0aWVzLCB0aGlzIHJhaXNlcyB0aGUgcXVlc3Rpb24gdG8gd2hhdCBleHRlbnQgdGhpcyBzYW1lIGdvZXMgZm9yIHJlc2VhcmNoZXJzIG9mIHNpbWlsYXIgcG9zaXRpb25zLlRoZXJlIGFyZSBzZXZlcmFsIHJlYXNvbnMgd2h5IHRoaXMgbWlnaHQgb2NjdXIuIFJlc2VhcmNoZXIgb2Ygc2ltaWxhciBwb3NpaXRvbnMgbWlnaHQgcHJlZmVyIHRvIHdvcmsgd2l0aCBvdGhlciByZXNlYXJjaGVycyB3aXRoIGEgc2ltaWxhciBsZXZlbCAoc2VsZWN0aW9uKSwgb3Igd29ya2luZyB0b2dldGhlciB3aXRoIGEgcmVzZXJhaGNlciBvZiBhIGhpZ2hlciBsZXZlbCBtaWdodCBpbmZsdWVuY2UgYSByZXNlYXJjaGVyIHRvIHByb21vdGUgKGluZmx1ZW5jZSkuDQoNCklmIHRoaXMgaG9tb3BoaWx5IHN0cnVjdHVyZSBleGlzdHMsIHRoaXMgcmFpc2VzIHRoZSBxdWVzdGlvbiB3aGF0IHRoZSBjb25zZXF1ZW5jZXMgYXJlIG9mIHRoaXMgc3RydWN0dXJlIG9mIGNvbGxhYmVyYXRpb25zIG9uIHBhcGVycy4gT3RoZXIgcmVzZWFyY2hlcnMgbWlnaHQgcHJlZmVyIHRvIGNpdGUgcGFwZXJzIGJ5IHJlc2VhcmNoZXJzIG9mIGEgaGlnaGVyIHBvc2l0aW9uLiBGdXJ0aGVybW9yZSwgd29ya2luZyB3aXRoIHJlc2VhcmNoZXJzIG9mIGEgaGlnaGVyIHBvc2l0aW9uIG1pZ2h0IGluZmx1ZW5jZSB0aGVzZSByZXNlYXJjaGVycyBpbiBzdWNoIGEgd2F5IHRoYXQgdGhlaXIgZnV0dXJlIHBhcGVycyB3aWxsIGJlIGNpdGVkIG1vcmUgb2Z0ZW4uDQoNCiNSUXMNCg0KKERlc2NyaXB0aXZlKSBSUTE6IFRvIHdoYXQgZGVncmVlIGlzIHRoZXJlIHBvc2l0aW9uIGhvbW9waGlseSBpbiBwdWJsaWNhdGlvbiBjb2xsYWJlcmF0aW9ucyBiZXR3ZWVuIHNvY2lhbCBzY2llbmNlIHJlc2VhcmNoZXJzIGluIHRoZSBOZXRoZXJsYW5kcz8NCg0KKEV4cGxhbmF0b3J5KSBSUTI6IEhvdyBkb2VzIHRoZSBwb3NpdGlvbiBjb21wb3NpdGlvbiBvZiBhIHJlc2VhcmNoZXIncyBlZ29uZXQgKHBhcGVyIGNvbGxhYmVyYXRpb25zIC8gZGVwYXJ0bWVudCBpbmZsdWVuY2UgYSByZXNlYXJjaGVyJ3MgY2l0YXRpb25zPw0KDQooRXhwbGFuYXRvcnkpIFJRMzogSG93IGRvZXMgdGhlIHBvc2l0aW9uIGNvbXBvc2l0aW9uIG9mIGEgcGFwZXIgaW5mbHVlbmNlIHRoZSBudW1iZXIgb2YgY2l0YXRpb25zIG9mIHRoZSBwYXBlcj8NCg0KIyBEYXRhIG5lZWRlZA0KDQotICAgTm9kZSBjaGFyYWN0ZXJpc3RpY3M6IHBvc2l0aW9uLCBjaXRhdGlvbnMsIGFmZmlsaWF0aW9uLg0KDQotICAgUmVsYXRpb25hbCBhdHRyaWJ1dGVzOiB3aG8gY29sbGFiZXJhdGVkIHdpdGggd2hvbSwgcGFwZXIgY2l0YXRpb25zLg0KDQojIFJlZmVyZW5jZXMNCg0KSG9ydGEsIEguLCBGZW5nLCBTLiAmIFNhbnRvcywgSi5NLiBIb21vcGhpbHkgaW4gaGlnaGVyIGVkdWNhdGlvbiByZXNlYXJjaDogYSBwZXJzcGVjdGl2ZSBiYXNlZCBvbiBjby1hdXRob3JzaGlwcy4gU2NpZW50b21ldHJpY3MgMTI3LCA1MjPigJM1NDMgKDIwMjIpLiA8aHR0cHM6Ly9kb2kub3JnLzEwLjEwMDcvczExMTkyLTAyMS0wNDIyNy16Pg0KDQpgYGB7ciBzZXR1cCwgaW5jbHVkZT1GQUxTRX0NCmtuaXRyOjpvcHRzX2NodW5rJHNldChlY2hvID0gVFJVRSkNCmBgYA0K