class_docstr: A mixin providing generic (not family-specific) methods for
  variational Bayes mean field fitting.
class_docstr_tok: ["A","mixin","providing","generic","(","not","family-specific",")","methods","for","variational","Bayes","mean","field","fitting","."]
code: def _lp_stats(self, fep_mean, fep_sd, vc_mean, vc_sd):
      tm = np.dot(self.exog, fep_mean)
      tv = np.dot(self.exog**2, fep_sd**2)
      tm += self.exog_vc.dot(vc_mean)
      tv += self.exog_vc2.dot(vc_sd**2)
  
      return tm, tv
  
imports: import numpy
  from scipy.optimize import minimize
  from scipy import sparse
  import statsmodels.base.model
  from statsmodels.iolib import summary2
  from statsmodels.genmod import families
  import pandas
  import warnings
  import patsy
label_desc: Use image node_id 1 for calling the _VariationalBayesMixedGLM obj's underlying member method code with example usage: obj._lp_stats(fep_mean, fep_sd, vc_mean, vc_sd) and returns: tm,  tv