///|
/// ANCOMBC module demo - Microbiome differential abundance with bias correction.
///
/// Run: moon run examples/ancombc_demo
fn main {
  println("=== ANCOMBC Module Demo ===")
  println("Analysis of Composition of Microbiomes with Bias Correction")
  println("Reference: Lin & Peddada (2020) Nat Commun 11:5144")
  println("")

  // ------------------------------------------------------------------
  // 1. Load sample data
  // ------------------------------------------------------------------
  println("=== 1. Load Sample Data ===")
  let data = @bio.ancombc_sample_data()
  println("Samples: " + data.n_samples().to_string())
  println("Features: " + data.n_features().to_string())
  let groups = data.groups()
  println("Groups: " + groups.length().to_string())
  for g in groups {
    println("  - " + g)
  }
  println("")
  let lib_sizes = data.library_sizes()
  println("Library sizes:")
  for i in 0..