performs one VE-step, that is, update of tau and rho
VEstep(VE, mstep, data, modelFamily, directed, fix.iter = 5)
list with variational parameters tau and rho
list with current model parameters and additional auxiliary terms
data vector in the undirected model, data matrix in the directed model
probability distribution for the edges. Possible values:
Gauss
, Gamma
booelan to indicate whether the model is directed or undirected
maximal number of iterations for fixed point equation
updated list VE
with variational parameters tau and rho