JEvalMstep function

Evaluation of criterion J

Evaluation of the criterion J to verify the convergence of the VEM algorithm

JEvalMstep(VE, mstep, data, directed, sparse, method = "hist")

Arguments

• VE: Results of the previous VE for iterative computation

• mstep: Results of the previous mstep for iterative computation

  • mstep$sum_rhotau : N_Q vector (not needed in the function)
  • mstep$sum_rhotau_obs : N_Q vector
  • mstep$logintensities.ql : N_Q x Dmax matrix
  • m.step$beta : N_Q vector

• data: Data same of mainVEM

• directed: Boolean for directed (TRUE) or undirected (FALSE) case

• sparse: Boolean for sparse (TRUE) or not sparse (FALSE) case

• method: List of string. Can be "hist" for histogram method or "kernel" for kernel method