def bayes_fitting(self):
self = pickle.load(
open(
"/%s/%s_self_para_eval.p" %
(self.output_directory, self.output_prefix), 'rb'))
logging.info('Modeling the mean and variance ...')
if self.genes_varmodeling > 0:
self.mrm = MeanVarModel()
#self.mrm.model_mean_disp_by_lm(self.allgenedict)
self.mrm.model_mean_disp_by_glm(self.allgenedict,
self.output_prefix, self.size_f)
else:
mrm = None
pickle.dump(
self,
open(
"/%s/%s_self_fitting.p" %
(self.output_directory, self.output_prefix), 'wb'))
