AI Agents in Statistical Genetics Coding: A Practical Guide
2026-05-04
Statistical genetics is a coding-intensive field. On any given day you might be writing R scripts for quality control, building GWAS pipelines in Python, submitting batch jobs to an HPC cluster, or debugging statistical models — across projects spanning GWAS meta-analysis, whole-genome sequencing, rare variant association testing, and functional genomics.
Over the past year, AI coding agents have fundamentally changed how I work. This post covers the tools I use, the concepts that make them effective, and practical patterns for statistical genetics specifically, including how to work safely with sensitive genomic data.
▶ View as slides — 20-slide deck for lab meetings
Institutional compliance
Before listing the tools: if you work at MGB or the Broad Institute, GitHub Copilot is currently the only AI coding tool approved under both institutions’ data governance policies, covered by the Microsoft enterprise agreement. Claude Code, Codex, and other third-party tools operate under personal subscriptions and have not been cleared under either institutional agreement — including Broad.
This does not mean those tools are unsafe by design — it means they have not gone through institutional review. Until that changes, use Copilot for any work in these institutional environments.
If you work at another institution, check with your IT or compliance team before using any AI tool in a research context. The core question to ask: Is this tool covered by an enterprise agreement that includes a BAA (Business Associate Agreement) or equivalent data governance terms?
The AI toolkit
VS Code: the hub
All my AI tools live inside VS Code, which I use both locally and via remote SSH to cluster head nodes. This makes the IDE the natural integration point.
GitHub Copilot
GitHub Copilot is available through institutional and enterprise GitHub licenses — check with your institution’s IT department. If you’re at a research university or major institute, there’s a good chance it’s already covered. This makes it the lowest-friction entry point.
Copilot in VS Code has come a long way. The Copilot Chat panel now does most of what you’d get from a dedicated AI CLI: it can read files, make edits, run terminal commands, and work through multi-step tasks in your project. Recent additions make it especially useful:
- Autopilot mode runs tasks end-to-end without asking for approval on every step. VS Code added configurable approval levels — including an “approve all” option — so you can let the agent run uninterrupted, similar to using a CLI tool.
- Auto model lets the agent self-select the best model for each task. The pool includes Claude models, Codex, and Gemini. This is genuinely useful: simpler tasks route to faster, cheaper models, while complex reasoning gets a capable one. Copilot claims this delivers around 10% cost savings, though part of that may simply be cheaper models handling easier tasks.
- Local models are also supported — Qwen and custom endpoints if you have them. Open-source models have caught up substantially for routine tasks: boilerplate generation, simple refactoring, format conversions. They may be a few months behind on complex reasoning, but “a few months behind the best” is often fine for everyday research coding.
One important billing note: starting June 1, 2026, GitHub Copilot is moving to usage-based billing, replacing premium requests with GitHub AI Credits. It is worth checking with your institution’s IT or finance team whether the enterprise agreement covers this or whether individual costs will change.
Claude Code
Claude Code is a terminal-first AI coding agent from Anthropic. I use it as a personal subscription alongside Copilot. It operates in your shell, can read and write files across your project, run commands, and maintain context about your entire codebase. Its strength is in extended, multi-step tasks where having the full project context matters.
Claude Code introduces a few concepts worth knowing even if you use other agents, because the ideas transfer:
- Memory: persistent context that carries across sessions — your project structure, conventions, HPC environment details, which packages you prefer. You don’t re-explain your setup every time.
- Skills: reusable, domain-specific capabilities you can invoke by name. Think of them as encoded expertise for recurring task types.
- Hooks: automated actions triggered by events — for example, running a linter after every file edit, or checking that a generated SLURM script has valid memory and partition settings before submission.
- Multi-agents: spawning parallel agents to explore different parts of a codebase or documentation simultaneously, with results feeding back into the main task.
Other tools
I use Codex for quick one-off generation and explaining unfamiliar code, and Gemini when I need a large context window — processing long documentation, comparing multiple pipeline scripts, or working with lengthy summary statistics files.
Setting up your agent for success: the context file
The single most impactful thing you can do before using any AI agent on a research project is give it explicit context about your environment. For Claude Code, this is a CLAUDE.md file at the project root. Copilot and other tools have equivalents (.github/copilot-instructions.md for Copilot, AGENTS.md for some others).
A good context file tells the agent:
- What the project is about scientifically
- Where things live (data directories, code directories, output paths)
- What environment it’s running in (HPC scheduler, node types, container runtime)
- What tools and packages you use and prefer
- What coding conventions to follow
- What to avoid
Here’s an example from one of my active projects:
## Scientific Focus
- Prioritize statistical and biological correctness over surface-level style
- When reviewing analysis scripts, lead with the science before describing mechanics
## Project Structure
- Pipeline: limma differential expression → gene-level effect sizes →
burden correlation → plots
- Key directories: ref_data/ (reference), results/ (per-analysis outputs),
data/ (intermediate), results/
## HPC Environment (UGER / Grid Engine)
- Submit jobs with qsub; interactive sessions with qrsh
- GPU jobs require three flags: -l gpu=1, -l os=RedHat8, -hard
- Container runtime: podman (rootless, Docker CLI-compatible)
- GPU batch job template:
#!/usr/bin/env bash
#$ -l gpu=1
#$ -l os=RedHat8
#$ -hard
#$ -l h_rt=HH:MM:SS
#$ -l h_vmem=<N>G
podman run --device nvidia.com/gpu=all --rm \
--security-opt=label=disable <image> <command>
## Data Conventions
- Use ENSG IDs throughout (not gene symbols)
## Coding Standards
- R is primary; use data.table and apply family over loops
- Scripts must be self-contained with explicit seeds
- Keep commit messages ≤ 5 words
This kind of file turns generic AI suggestions into project-aware ones. The agent knows to write data.table::fread not read.csv, to target your specific scheduler not a generic one, to use podman not docker, and to reference ENSG IDs in any gene-level code. It takes twenty minutes to write and saves hours of corrections.
The workflow: Explore, Plan, Implement, Commit
Regardless of which AI tool you use, this four-phase pattern produces better results than jumping straight to “write me a script.”
Explore
Before writing code, let the agent understand what already exists. Point it at your project and ask it to find existing implementations, read relevant documentation, or understand how a previous analysis was structured. This prevents reinventing the wheel and ensures new code is consistent with the rest of the project.
# Example prompt to an agent
"I need to run a rare variant burden test on TOPMed data.
Search this project for any existing STAAR or SKAT scripts.
Check the renv lockfile for which R packages are installed.
Look at how previous analyses handled variant annotation."
Plan
Before writing code, outline the approach. What are the inputs and outputs? Which existing functions can be reused? What edge cases matter — missing phenotype data, multi-allelic variants, related individuals? For statistical genetics work, this step is especially important because a subtle mistake in variant filtering or covariate adjustment can invalidate an entire analysis.
Most AI agents support a plan-first mode. Use it. Review the plan before any code is written.
Implement
Write code with the agent’s assistance. At this stage it has context from exploration and planning, so suggestions are project-aware. Tips:
- Be specific about statistical methods. “Run a burden test” is too vague. “Run STAAR-O with MAF < 0.01, adjusting for age, sex, and the first 10 PCs, looping over all annotated genes on chr22” is actionable.
- Let the agent handle boilerplate: SLURM/UGER headers, argument parsing, logging, file-path handling.
- Stay involved in the statistics. The agent can implement any method you name; choosing the right method requires your domain expertise.
Commit
Review the diff, run a sanity check, and commit. A good agent writes a concise commit message, checks that no hardcoded local paths leaked into the script, and verifies the .gitignore excludes data files.
Approval settings: keeping the agent out of dangerous territory
AI agents are most useful when they can act autonomously — reading files, running code, making edits — without you approving every keystroke. But unconstrained autonomy is also how things go wrong: an agent that can delete files, overwrite data, or push to a remote repository without asking can cause damage that is hard or impossible to undo. Approval settings let you define exactly where the agent acts freely and where it stops and asks.
Copilot agent approval settings in VS Code
The core setting is chat.permissions.default, which takes one of three values:
default(Default Approvals): the agent pauses and asks before running terminal commands or modifying files. Safest, and the right starting point for new tasks.autopilot: auto-approves all tool calls and runs until task completion without pausing. Use this only for well-scoped tasks after you have reviewed the plan.autoApprove: bypasses approvals entirely. Avoid this in research environments — it disables the safety net with no granular control.
Beyond the top-level mode, VS Code exposes per-category controls that give you more precision:
// settings.json
{
// auto-approve only safe, read-only commands; rm and curl still prompt
"chat.tools.terminal.autoApprove": ["ls", "cat", "head", "grep", "Rscript", "plink2"],
// only allow edits inside src/ and scripts/ without asking
"chat.tools.edits.autoApprove": ["src/**", "scripts/**"],
// cap how many steps an agent can take before stopping (default: 25)
"chat.agent.maxRequests": 25,
// block the agent from making outbound calls to external domains
"chat.agent.deniedNetworkDomains": ["*"],
"chat.agent.allowedNetworkDomains": ["api.anthropic.com", "api.openai.com"]
}
maxRequests is an underused setting — it acts as a runaway-prevention mechanism. If an agent gets into a loop or starts doing more than you expected, it will stop at 25 steps and hand control back to you. For exploratory tasks raise it; for well-defined ones leave it at the default.
The network domain filtering (chat.agent.deniedNetworkDomains / chat.agent.allowedNetworkDomains) is particularly relevant for genomics work. Denying all domains except the model APIs means the agent cannot exfiltrate data to an unexpected endpoint, even accidentally.
VS Code also offers a sandbox mode (currently preview on macOS and Linux) that runs agent shell commands in an isolated environment. Enable it with chat.agent.sandbox.enabled: "on" for full filesystem isolation, or "allowNetwork" for network access only. Worth enabling when running unfamiliar agent sessions.
For separate control over the planning vs. implementation phases, two settings let you assign different models to each:
{
"chat.planAgent.defaultModel": "claude-opus-4-7", // careful reasoning for planning
"github.copilot.chat.implementAgent.model": "claude-sonnet-4-6" // faster for writing code
}
Claude Code permissions
Claude Code has a permission system in the project’s .claude/settings.json:
{
"permissions": {
"allow": [
"Bash(Rscript:*)",
"Bash(plink2:*)",
"Bash(regenie:*)",
"Read(**)",
"Write(src/**)",
"Write(scripts/**)"
],
"deny": [
"Bash(rm -rf *)",
"Write(/protected/**)"
]
}
}
allow entries let the agent use specific commands or write to specific paths without prompting. deny entries are hard blocks — the agent will not execute them even if asked. A well-written settings file for a genomics project allows running R, PLINK, and REGENIE, allows reading anywhere in the project tree, allows writing to src/ and scripts/, and hard-denies any writes to the /protected/ data directory.
Plan mode: the pre-flight check
Both Copilot and Claude Code have a plan-before-code mode. In Copilot, select the Plan agent from the agent picker. In Claude Code, enable plan mode to have the agent describe what it intends to do and wait for your sign-off before writing anything.
For complex tasks — “refactor this entire QC pipeline,” “add chromosome-stratified output to all downstream scripts” — always review the plan first. Look specifically for:
- Files it intends to modify that you did not expect
- Shell commands it plans to run (especially anything involving deletion, overwriting, or external network calls)
- Assumptions about data formats or file locations that might be wrong
Approving the plan does not mean blindly trusting the implementation. Review the diff after implementation, before committing.
What should always require your approval
Regardless of settings, these actions should always stop and ask:
- Any deletion —
rm, overwriting an existing file, clearing a directory - Git operations that affect remote state —
push,force-push, branch deletion - Writes to protected data directories — any path containing your actual genomic or phenotypic data
- External API calls — anything that sends data outside your environment
- HPC job submission —
qsub,sbatch, or equivalent: a misfired array job is expensive to cancel and can hold cluster resources for hours
The goal is not to make the agent slow. It is to ensure that irreversible or high-impact actions always have a human in the loop.
Working with sensitive genomic data
This deserves its own section because the stakes are higher.
The core rule: do not let AI agents read your data files directly. Genetic and health data are often subject to data use agreements (DUAs), IRB restrictions, and institutional policies. Even if the agent runs locally, the risk of data being embedded in conversation logs, sent to a remote API, or cached somewhere unexpected is real.
The practical workflow instead:
-
You inspect the data yourself. Check file headers and directory structure manually:
head -3 pheno.txt ls -lh /protected/data/topmed/ zcat variants.vcf.gz | head -20 -
Paste the metadata — not the data — into the agent. Give it the column names, the file format, the variable descriptions, the directory tree. This is usually enough for it to write correct code.
My phenotype file has these columns (tab-separated, header row): FID, IID, AGE, SEX, BMI, LDL, HDL, TG There are ~8,000 rows. File path: /protected/topmed/pheno.txt My genotype data is PLINK2 format (pgen/psam/pvar) at: /protected/topmed/geno/chr{1..22}.* Write a REGENIE step 1 command that uses this phenotype file and includes age, sex, and BMI as covariates. -
The agent writes code; you run it. Review the generated script, then execute it in your own environment. The data never leaves your cluster.
-
For exploratory analysis on results (summary statistics, QQ plots, effect size distributions): these are typically non-identifiable and lower risk. Even so, check your DUA — some restrict any sharing of intermediate results.
This pattern lets you get full AI assistance on the hardest part (writing correct statistical code) while keeping identifiable data entirely within your controlled environment.
Common workflows where agents shine
GWAS QC pipeline: a sequence of PLINK2 commands for sample and variant QC — missingness, Hardy-Weinberg equilibrium, sex discrepancy, relatedness. Provide the filter criteria; the agent produces a complete, commented shell script.
Rare variant aggregation tests: R scripts using STAAR or SKAT that read annotation files, define variant sets, specify null models, and loop over genomic regions. The boilerplate-to-logic ratio is high — ideal for agent assistance.
Visualization: Manhattan plots, QQ plots, regional association plots, forest plots. These are ggplot2-heavy and benefit from inline suggestions for aesthetics and from agents that can iterate on layout.
HPC job submission: SLURM or UGER scripts that array over chromosomes or genomic regions. The agent handles job array syntax, resource flags, output file naming — things that are easy to get wrong and expensive to debug on a busy cluster. The context file (CLAUDE.md or equivalent) that encodes your scheduler details pays off most here.
Advanced usage
Models vs. agents: understanding what is actually running
This distinction matters practically, not just conceptually.
An LLM (Large Language Model) is the underlying neural network that reads text and generates responses: Claude Opus, Claude Sonnet, Claude Haiku, GPT-4o, Gemini 1.5 Pro, Qwen, and so on. It is stateless — it has no memory between calls, no ability to run code, and no awareness of your file system. By itself, it is a very powerful text-in / text-out function.
An agent is a system built on top of one or more LLMs that adds the infrastructure to actually do things:
- Tools: the ability to read and write files, run shell commands, call APIs, search the web
- Memory / context management: persisting project state across turns, loading relevant files into the prompt
- Action loop: the cycle of generating a plan, executing a tool call, observing the result, and planning the next step
- Orchestration: routing sub-tasks to different models or spawning parallel workers
Claude Code, GitHub Copilot, and OpenClaw are all agents. They each have their own orchestration logic, tool sets, and UX — but they can all be powered by many of the same underlying LLMs. Claude Sonnet running inside Claude Code and Claude Sonnet running inside Copilot is the same model, producing identical raw capabilities, in two different agent wrappers. The agent wrapper determines what the model can act on, not what it can reason about.
Why does this matter for research?
- You can swap the underlying model in most agents without rewriting your prompts or workflow. The context file (CLAUDE.md) works regardless of which model is in the loop.
- The “best” agent for a task is not always the one with the most powerful model — it is the one whose tool set and orchestration fit the task. A lightweight model inside an agent with direct HPC access beats a frontier model that can only suggest commands.
- Cost and speed are model properties, not agent properties. You pay for tokens processed by the model, not by the agent wrapper. Switching to a smaller model within the same agent cuts cost without changing your workflow.
Switching models mid-workflow with cc-switch
In Claude Code, /model (or --model at the CLI) lets you swap the underlying LLM without leaving your session. You can use this deliberately as a cost-and-speed strategy:
# Start a session with a powerful model for exploration and planning
claude --model claude-opus-4-7
# Inside the session, switch to a faster model for implementation
/model claude-sonnet-4-6
# Or switch to the fastest model for pure boilerplate generation
/model claude-haiku-4-5-20251001
A practical pattern for a GWAS pipeline:
- Exploration and planning (Opus or Sonnet): “Here is my project structure. Here are the column names from my phenotype file. Design a REGENIE step 1 and step 2 pipeline that handles related individuals using the KING kinship matrix.” Let the model reason carefully about the design.
- Implementation (Sonnet or Haiku): “Now write the SLURM array script for step 2, one job per chromosome.” This is mostly template-filling — a smaller model is fast and sufficient.
- Review (back to Sonnet): “Read the generated script and check for statistical correctness: are the covariate flags right, is the phenotype file path correct, does the output naming make sense?”
The same principle applies in VS Code Copilot’s Auto model mode, which does this automatically based on task complexity. Explicitly controlling the switch yourself gives you more precision; Auto mode gives you convenience.
OpenClaw: a local-first alternative
OpenClaw is an open-source personal AI assistant designed to run on your own devices rather than in the cloud. Where Claude Code operates primarily in your terminal and VS Code, OpenClaw routes AI interactions through messaging channels you already use — WhatsApp, Telegram, Slack, Discord — and adds voice capabilities on macOS and iOS.
For researchers, the most relevant properties are:
- Local-first: the Gateway (control plane) runs on your own machine. Conversations route through your device rather than a hosted service, which reduces concerns about where data is processed.
- Multi-channel: you can ask a question via Slack while your cluster job is running and get the response in the same thread. Useful for asynchronous workflows on HPC where you are not sitting in front of a terminal.
- Model-agnostic: connects to multiple providers. You can point it at Claude, GPT-4o, a locally-hosted Qwen instance, or a custom endpoint.
- Skills registry: similar in concept to Claude Code’s skills — reusable agent capabilities you can invoke by name.
The tradeoff relative to Claude Code or Copilot is that OpenClaw is less focused on code editing and file manipulation, and more on being an always-available assistant across your communication channels. Think of it as complementary: Claude Code for heads-down coding sessions; OpenClaw for quick questions, job monitoring, or retrieval tasks when you are away from your IDE.
It is also a useful reference implementation if you want to build a custom agent that routes to your own HPC environment — the architecture (Gateway + channel adapters + model backends) maps well onto research infrastructure.
MCP: connecting agents to external tools
The Model Context Protocol (MCP) is an open standard that lets AI agents connect to external data sources and tools through a common interface. Instead of pasting PubMed abstracts into the chat, you install an MCP server and the agent can search PubMed directly. Instead of copy-pasting gnomAD allele frequencies, the agent queries gnomAD itself.
The architecture is straightforward: an MCP server sits between the agent and the external resource, translating the agent’s requests into API calls and returning structured results. Once a server is registered in your agent’s config, you use it like any other tool — you describe what you need in plain language, the agent decides which server to call, and the result comes back into the conversation.
For statistical genetics, a handful of MCP servers are immediately practical:
- PubMed / NCBI: search literature, retrieve abstracts and metadata without leaving the coding session. Useful for verifying a method’s original paper or checking whether a gene has known lipid associations.
- bioRxiv: access preprint metadata and abstracts directly. Particularly useful for tracking recent GWAS results before peer review.
- ChEMBL: query compound bioactivity, targets, and drug mechanisms. Relevant if you work near the drug-gene interface (lipid-lowering targets, Mendelian randomization for therapeutic targets).
- ClinicalTrials.gov: search registered trials by condition or intervention. Useful for translational context around GWAS findings.
- gnomAD (via custom or community servers): retrieve population-level allele frequencies, variant annotations, and constraint metrics without switching browser tabs.
The practical value accumulates in workflow continuity. A typical task without MCP: you’re writing a lncRNA annotation script, pause, open a browser, search gnomAD, copy a frequency, paste it back, lose your train of thought. With a gnomAD MCP server registered, you stay in the agent session and say: “check the gnomAD allele frequency for rs12345 in EUR and AFR populations.” The agent calls the server, the answer appears inline, and you keep coding.
Setting up MCP in Claude Code is done via the claude mcp add command or by editing your project’s .mcp.json. For VS Code Copilot, MCP support has been added in recent versions — check the Copilot settings panel under “MCP Servers.”
The MCP registry lists community-maintained servers across many domains. If nothing exists for a niche resource (say, a specific population biobank), the MCP specification is simple enough to implement a lightweight server yourself in an afternoon.
Cost awareness
With Copilot, your usage is tracked in premium requests — a unit that counts how many “expensive” model interactions you consume from your monthly allowance. The key insight is that different models cost different numbers of premium requests per single interaction:
| Tier | Premium requests used | Examples |
|---|---|---|
| Standard | 1 per interaction | GPT-4.5, Claude Sonnet 4.6 — everyday coding tasks |
| Premium | 3 per interaction | Claude Opus 4.6 — stronger reasoning, planning |
| High-tier / Agentic | 15 per interaction | Claude Opus 4.7, agentic modes — complex multi-step work |
So the “3×” and “15×” labels are literal: one Opus 4.7 interaction consumes 15 premium requests in a single shot. If your monthly plan includes 300 premium requests, that’s 300 standard interactions, 100 premium interactions, or just 20 high-tier interactions before you hit the cap. One long agentic pipeline run with Opus 4.7 can exhaust a significant portion of your monthly allowance.
Practical strategies:
- Default to Standard for most coding tasks: Sonnet-class models handle GWAS QC scripting, plot generation, and boilerplate just as well at 1 premium request per interaction.
- Reserve Premium (3 requests) for planning and architecture: designing a new pipeline, debugging a subtle statistical error, or reviewing a complex script — tasks where stronger reasoning actually pays off.
- Use High-tier (15 requests) deliberately: for genuinely complex multi-step agentic work only. Don’t leave an Opus 4.7 session running overnight on routine tasks.
- Keep context lean: every turn in a session draws from your allowance. Sharing a full VCF header when you only need column names costs premium requests on every subsequent exchange in that conversation.
- Prefer targeted questions: “why is my REGENIE step 2 producing NA p-values for chromosome 6?” uses fewer premium requests than “review my entire pipeline” — and tends to get a better answer.
- Check your usage: GitHub Settings → Billing → Copilot. Worth checking after your first few agentic sessions so you can calibrate.
From June 1, 2026, GitHub Copilot moves to usage-based billing (AI Credits replace premium requests for most plans). Check with your institution’s IT team on how the enterprise agreement handles this transition before running extended high-tier sessions.
Tips and caveats
What AI agents are reliably good at:
- Syntactically correct code in R, Python, and bash
- Translating a described statistical method into working code
- Generating boilerplate (job scripts, argument parsers, logging)
- Refactoring messy scripts into clean, modular functions
- Explaining error messages and suggesting fixes
What they are not good at (yet):
- Choosing the right statistical method for your data. They can implement any method you name; the choice requires your expertise.
- Knowing your cluster’s specifics without being told. The context file is how you bridge this.
- Validating scientific results. A Manhattan plot can be generated correctly but still represent a batch artifact — that judgment is yours.
How to verify AI-generated statistical code:
- Read every line. A wrong sign, a dropped covariate, or an off-by-one in variant filtering can silently produce wrong results.
- Test on a small subset (e.g., chromosome 22) before submitting a 22-chromosome array job.
- Compare against known results or a previous version of the pipeline where possible.
- Check that test statistics, p-values, and effect sizes are in plausible ranges.
Don’t get locked into one provider. This field moves faster than any other area of software I have watched. The model that is clearly best today may be third-best in six months, and a newcomer you haven’t heard of may be leading on code benchmarks by the time you read this. The practical implication: structure your workflow around the interface (VS Code, a context file, the Explore→Plan→Implement→Commit loop) rather than around any specific model. Copilot’s Auto mode helps here — it routes to whichever model performs best for a given task, across Anthropic, OpenAI, Google, and others, without you having to manually track rankings. Outside of Copilot, it is worth spending an afternoon every few months with a new model on a real task you know well, so you have a concrete sense of where the landscape actually stands.
The bottom line: AI agents do not replace statistical genetics expertise — they amplify it. They handle the mechanical parts of coding so you can focus on the science.