mynet1 <- sienaNet(array(c(tmp3, tmp4), dim=c(32, 32, 2)))
mydata <- sienaDataCreate(mynet1)
myeff<- getEffects(mydata)
mymodel<- sienaModelCreate(fn = maxlikefn, nsub=2, n3=100, maxlike=TRUE)
ans<- siena07(mymodel, data=mydata, effects=myeff, batch=TRUE)
Run the code above in your browser using DataLab