mynet1 <- sienaDependent(array(c(s501, s502, s503), dim=c(50, 50, 3)))
mybeh <- sienaDependent(s50a, type="behavior")
mydata <- sienaDataCreate(mynet1, mybeh)
myeff <- getEffects(mydata)
myeff <- includeGMoMStatistics(myeff, egoX_gmm, interaction1="mybeh")
myeff
Run the code above in your browser using DataLab