Skip to main content
Glama

diversity_pca

Performs PCA on genotype dosage data to reveal population structure, outputting sample coordinates, variance explained, and outlier flags. Optionally labels groups from metadata for plot coloring.

Instructions

Principal component analysis of population structure.

Runs PCA on the alt-allele dosage matrix (monomorphic markers dropped, missing mean-imputed, Patterson scaling). Writes pca_coords.csv (per-sample PC coordinates) and reports variance explained plus any PC1/PC2 outlier samples (beyond outlier_sd SD). Pass metadata_tsv + group_column to add a group column (population label per sample) for colouring the PC plot. For large sets pass method="allelematrix" + max_markers to avoid a full VCF export.

Input Schema

TableJSON Schema
NameRequiredDescriptionDefault
methodNoGenotype source: 'vcf' (full export, cached) or 'allelematrix' (paged, server-side subset).vcf
regionNoRestrict analysis to a genomic window: 'chrom' or 'chrom:start-end' (1-based).
id_columnNoColumn in the metadata TSV holding the individual/accession id (default 'individual').individual
outlier_sdNoFlag points more than this many standard deviations from the mean.
output_dirNoDirectory for the output CSV(s) (default ./gigwa_results/<module>/).
max_markersNoCap analysis to the first N markers in canonical Gigwa search order; omit to use all.
group_columnNoColumn in the metadata TSV holding the group/population label.
metadata_tsvNoPath to a metadata TSV (import_metadata format) used to define groups.
n_componentsNoNumber of principal components to compute.
variant_set_db_idYesBrAPI variantSetDbId identifying the run (MODULE§project§run) -- copy the exact string from list_variant_sets / list_content, never assemble one by hand: the middle segment is a numeric project index, not the project's name, and a wrong guess fails with an opaque HTTP 500 rather than a clear error.

Output Schema

TableJSON Schema
NameRequiredDescriptionDefault
resultYes

Schema Changelog

Changes observed during successful MCP inspections. Dates show when Glama detected each change.

  1. Changed2 schema fields changedv1.9.1
    • changedInput schema / properties / max_markers / description
      Previous value: -"Cap the number of markers analysed (evenly-spaced subsample); omit to use all."New value: +"Cap analysis to the first N markers in canonical Gigwa search order; omit to use all."
    • changedInput schema / properties / variant_set_db_id / description
      Previous value: -"BrAPI variantSetDbId identifying the run (MODULE§project§run); from list_variant_sets / list_content."New value: +"BrAPI variantSetDbId identifying the run (MODULE§project§run) -- copy the exact string from list_variant_sets / list_content, never assemble one by hand: the middle segment is a numeric project index, not the project's name, and a wrong guess fails with an opaque HTTP 500 rather than a clear error."
  2. First observedv1.4.16

TDQS

A4.7/5.0
Behavior5/5

Does the description disclose side effects, auth requirements, rate limits, or destructive behavior?

No annotations are present, so the description carries the full burdens. It discloses preprocessing decisions (monomorphic markers dropped, missing mean-impured, Patterson scaling), output artifacts (pca_coords.csv, variance explained, outlier flags), and the trade-off between VCF export and allelematrix streaming. This is unusually transparent.

Agents need to know what a tool does to the world before calling it. Descriptions should go beyond structured annotations to explain consequences.

Conciseness5/5

Is the description appropriately sized, front-loaded, and free of redundancy?

Four compact sentences: purpose, algorithm and outputs, optional grouping, and large-set guidance. Every sentence earns its place and there is no restatement of schema fields.

Shorter descriptions cost fewer tokens and are easier for agents to parse. Every sentence should earn its place.

Completeness5/5

Given the tool's complexity, does the description cover enough for an agent to succeed on first attempt?

With 10 parameters and no annotations, the description covers the core algorithm, outputs, optional metadata integration, and scale guidance. The rest of parameter meaning is fully supplied by the 100%-coverage input schema, and return values are covered by the output schema.

Complex tools with many parameters or behaviors need more documentation. Simple tools need less. This dimension scales expectations accordingly.

Parameters4/5

Does the description clarify parameter syntax, constraints, interactions, or defaults beyond what the schema provides?

Schema coverage is 100%, so the baseline is 3. The description adds coupling semantics: metadata_tsv + group_column together produce a group column for plotting, and method='allelematrix' + max_markers together avoid full VCF export. It does not describe region, id_column, or output_dir, but the schema already documents those completely.

Input schemas describe structure but not intent. Descriptions should explain non-obvious parameter relationships and valid value ranges.

Purpose5/5

Does the description clearly state what the tool does and how it differs from similar tools?

The description opens with 'Principal component analysis of population structure' and immediately clarifies the mechanism: 'Runs PCA on the alt-allele dosage matrix'. This names a specific algorithm and resource, and the PCA framing distinguishes it from other diversity/structure tools like diversity_structure or diversity_kinship.

Agents choose between tools based on descriptions. A clear purpose with a specific verb and resource helps agents select the right tool.

Usage Guidelines4/5

Does the description explain when to use this tool, when not to, or what alternatives exist?

It gives actionable conditional guidance: add metadata_tsv + group_column for group labelling/colouring, and use method='allelematrix' + max_markers for large sets to avoid full VCF export. No explicit when-not/alternative sibling routing is provided, but the practical usage context is clear.

Agents often have multiple tools that could apply. Explicit usage guidance like "use X instead of Y when Z" prevents misuse.

Latest Blog Posts

MCP directory API

We provide all the information about MCP servers via our MCP API.

curl -X GET 'https://glama.ai/api/mcp/v1/servers/gkanogiannis/Gigwa-MCP'

If you have feedback or need assistance with the MCP directory API, please join our Discord server