mynet1 <- sienaNet(array(c(s501,s502,s503),dim=c(50,50,3)))
mynet2 <- sienaNet(s50a,type='behavior')
mydata <- sienaDataCreate(mynet1, mynet2)
mydata <- sienaDataConstraint(mydata, mynet2, mynet1, "higher", FALSE)
Run the code above in your browser using DataLab