1) split into two functions FDRreg and BayesFDRreg
2) BayesFDRreg has EB and FB modes... the former is fit in C++, the latter in R because of buggy mixture-model sampling
3) In the EB version of BayesFDRreg, EM is used to estimate the mixing density