def fit_null_model(self):
null_model = MixedModel(self.pedigrees,
outcome=self.outcome,
fixed_effects=self.fixed_effects)
null_model.add_genetic_effect()
null_model.fit_model()
null_model.maximize(method=self.maximization, verbose=self.verbose)
return null_model
def fit_alternative_model(self):
ibd_model = MixedModel(self.pedigrees,
outcome=self.outcome,
fixed_effects=self.fixed_effects)
ibd_model.add_genetic_effect()
ranef = RandomEffect(self.analysis_individuals,
'IBD',
incidence_matrix='eye',
covariance_matrix=self.ibd_matrix)
ibd_model.add_random_effect(ranef)
ibd_model.fit_model()
ibd_model.maximize(verbose=self.verbose, method=self.maximization)
return ibd_model
