gets the connected set, then

m2.trace.estimate(sim, model0 = NA, hetero = FALSE)