def lme_regression(output_dir: str,
table: pd.DataFrame, tree: skbio.TreeNode,
metadata: Metadata, formula: str,
groups: str) -> None:
res = mixedlm(table=table, metadata=metadata._dataframe,
formula=formula, groups=groups)
res.fit()
lme_summary(output_dir, res, tree)
def lme_regression(output_dir: str,
table: pd.DataFrame, tree: skbio.TreeNode,
metadata: Metadata, formula: str,
groups: str) -> None:
if np.any(table.var(axis=0) == 0):
message = ('Detected zero variance balances - '
'double check your table for unobserved features.')
raise UserWarning(message)
res = mixedlm(table=table, metadata=metadata.to_dataframe(),
formula=formula, groups=groups)
res.fit()
lme_summary(output_dir, res, tree)