def test_GGGM2(verbose=0):
G = GGGM()
sx = 10000
x = nr.randn(sx)
G.init_fdr(x)
G.estimate(x)
assert(G.mixt[1]>0.9)
def test_GGGM1(verbose=0):
G = GGGM()
sx = 10000
x = np.array([float(st.t.rvs(5)) for i in range(sx)])
G.init_fdr(x)
G.estimate(x)
if verbose:
G.parameters()
assert(np.absolute(G.mean)<0.1)