@abner12 wrote:
expdyad<-function(exgrp)
{
trial <- vector(“list”,100)
for ( i in 1:100 )
{
rg <- erdos.renyi.game(n = vcount(exgrp), p.or.m = graph.density(exgrp), directed = TRUE)
trial[[i]] <- dyad.census(rg)
Reduce("+",trial/100)
}
return(data.frame(dyad.census=dyad.census()))
}Dsplit <- split(df1, df1$village)
sapply(Dsplit, function(x) expdyad(graph_from_data_frame(x)))
Posts: 1
Participants: 1