Minimum example to remove edges connecting nodes of interest

remove_edge.R

set.seed(1234)
g <- sample_gnp(n=10, 0.15) #make_ring(10)
V(g)$name <- letters[1:10]
plot(g)
rm_e <- combinations(n=3, r=2, v=c('b', 'i', 'c'))
rm_e_do <- apply(rm_e, 1, function(r) are.connected(g, r[1], r[2]))
rm_e <- rm_e[rm_e_do, ] %>% apply(., 1, paste, collapse='|')
par(mfrow=c(1, 2))
plot(g)
plot(g - edge(rm_e))
dev.off()