def ols_regression(output_dir: str, table: pd.DataFrame, tree: skbio.TreeNode, metadata: Metadata, formula: str) -> None: res = ols(table=table, metadata=metadata._dataframe, formula=formula) res.fit() ols_summary(output_dir, res, tree)
def ols_regression(output_dir: str, table: pd.DataFrame, tree: skbio.TreeNode, metadata: Metadata, formula: 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 = ols(table=table, metadata=metadata.to_dataframe(), formula=formula) res.fit() ols_summary(output_dir, res, tree)