def METHOD_NAME():
    marginals_ = [Normal(loc=None, scale=2.), Lognormal(s=1., loc=0., scale=np.exp(1))]
    dist_joint_ = JointIndependent(marginals=marginals_)
    data = np.array([[-0.17126121, 0.91793325], [3.99469089, 7.36946747], [2.565957, 3.60736828]])
    mle_fit = dist_joint_.fit(data=data)
    assert np.round(mle_fit['loc_0'], 3) == 2.130