#Program to create network data in R from SAS 
# load R packages 
library(foreign)
library(statnet)
library(network)
#Read in the edge-level data from external file
s_edge <- read.xport('c:/jwm/conferences/columbia/workshop/s_edge.xpt')
#create a network object from the Node names of edgelist
s_friends <- network(s_edge[,c("N1" ,"N2" )])
summary(s_friends)
#Now create the node attributes 
s_nodes <- read.xport('c:/jwm/conferences/columbia/workshop/s_nodes.xpt')
s_friends %v% "FEMALE" <- s_nodes[,7 ]
s_friends %v% "NID" <- s_nodes[,5 ]
s_friends %v% "S3" <- s_nodes[,1 ]
s_friends %v% "S9" <- s_nodes[,2 ]
s_friends %v% "S48" <- s_nodes[,3 ]
s_friends %v% "S64" <- s_nodes[,4 ]
s_friends %v% "WHITE" <- s_nodes[,6 ]
list.vertex.attributes(s_friends)
#Now assign edge attributes 
set.edge.attribute(s_friends,"n1" ,s_edge[,1 ])
set.edge.attribute(s_friends,"n2" ,s_edge[,2 ])
list.edge.attributes(s_friends)
#Now assign dyad values 
s_dycov <- read.xport('c:/jwm/conferences/columbia/workshop/s_dycov.xpt')
s_clubs <- network(s_dycov[,c("N1" ,"N2" )])
summary(s_clubs)
set.edge.attribute(s_clubs,"n1" ,s_dycov[,1 ])
set.edge.attribute(s_clubs,"n2" ,s_dycov[,2 ])
set.edge.attribute(s_clubs,"ovlpec" ,s_dycov[,3 ])