def get_global_suff_stats(self, Data, SS, LP, **kwargs):
    ev = LP['edge_variational']
    E, K = ev.shape
    sb_ss0 = np.sum(ev[Data.ind0,:], axis=0)
    sb_ss1 = np.sum(ev[Data.ind1,:], axis=0)
    SS.setField('sb_ss0', sb_ss0, dims=('K'))
    SS.setField('sb_ss1', sb_ss1, dims=('K'))
    return SS

  ######################################################### Global Params
  #########################################################   M-step
Calculate and return sufficient statistics.
dataset/ETHPy150Open daeilkim/refinery/refinery/bnpy/bnpy-dev/bnpy/obsmodel/BernRelObsModel.py/BernRelObsModel.get_global_suff_stats