require(ergm) Y<-as.matrix(read.table("CaltechY.dat")) m<-ergmMPLE(network(Y)~ edges + sender + receiver, fitmodel = TRUE)