Filter feature table by discriminant score

Wrapper around get_discriminant_features() that returns only features with consensus score ≥ threshold.

get_subset_discriminant(
  tibble_rbims,
  analysis,
  metadata,
  group_col,
  score_min = 1,
  min_presence = 2,
  feature_col = NULL
)

Arguments

tibble_rbims

A tibble returned by consensus_rank() (or rbims_consensus() ...), containing at least the columns feature and score.

analysis

Character: "KEGG", "Pfam", "INTERPRO", "dbCAN", or "MEROPS"

metadata

MAG-level metadata (rownames = MAG identifiers)

group_col

Metadata column to discriminate (e.g. "Depth", "Class", "Phylum")

score_min

Minimum consensus score (default = 2)

min_presence

Minimum number of MAGs in which a feature must appear

feature_col

Column name containing feature IDs

Value

Filtered tibble (same format as input) with the full discriminant object attached as attribute "rbims_disc".