#Contrast the Simmelian ties in the Coleman friendship network with the "raw" ties
data(coleman)
fall<-coleman[1,,] #Fall ties
spring<-coleman[2,,] #Spring ties
sim.fall<-simmelian(coleman[1,,]) #Fall Simmelian ties
sim.spring<-simmelian(coleman[2,,]) #Spring Simmelian ties
par(mfrow=c(2,2))
gplot(fall,main="Nominations in Fall")
gplot(spring,main="Nominations in Spring")
gplot(sim.fall,main="Simmelian Ties in Fall")
gplot(sim.spring,main="Simmelian Ties in Spring")
#Which ties shall survive?
table(fall=gvectorize(fall),spring=gvectorize(spring)) #Fall vs. spring
table(sim.fall=gvectorize(sim.fall),spring=gvectorize(spring))
sum(fall&spring)/sum(fall) #About 58% of ties survive, overall...
sum(sim.fall&spring)/sum(sim.fall) #...but 74% of Simmelian ties survive!
sum(sim.fall&sim.spring)/sum(sim.fall) #(About 44% stay Simmelian.)
sum(sim.fall&sim.spring)/sum(sim.spring) #39% of spring Simmelian ties were so in fall
sum(fall&sim.spring)/sum(sim.spring) #and 67% had at least some tie in fall
Run the code above in your browser using DataLab