Dummy example. Only degree-corrected LBM is available yet.

help("DcLbm")
simu <- rlbm(50, 10, rep(1 / 5, 5), rep(1 / 10, 10), matrix(runif(50), 5, 10))
model = DcLbm(alpha = 1)
sol = greed(X = simu$x, model=model, alg = Hybrid())