mynet1 <- sienaNet(array(c(tmp3, tmp4), dim=c(32, 32, 2)))
mydata <- sienaDataCreate(mynet1)
myeff<- getEffects(mydata)
myeff[myeff$effectName=='transitive triplets'&
myeff$type=='eval','include']<- TRUE
mymodel<- sienaModelCreate(fn = simstats0c, nsub=2, n3=100)
ans<- siena07(mymodel, data=mydata, effects=myeff, batch=TRUE)
Run the code above in your browser using DataLab