mymodel <- sienaModelCreate(findiff=TRUE, fn=simstats0c)
mynet1 <- sienaNet(array(c(s501, s502, s503), dim=c(50, 50, 3)))
mydata <- sienaDataCreate(mynet1)
myeff <- getEffects(mydata)
print01Report(mydata, myeff)
Run the code above in your browser using DataLab