## Not run:
# data(binned_zmat)
# out <- repfdr(pbz,bz,"replication",
# control = em.control(pi.initial = c(0.48,rep(0.02,26)),
# verbose = TRUE, nr.threads = 1))
# # iterations are printed; run bit slower (1 thread)
# ## End(Not run)
Run the code above in your browser using DataLab