Communitygithub.com

TorchDrug

TorchDrug 0.2.1 を用いた分子グラフ、特性予測、自己教師あり事前学習、分子生成、逆合成、タンパク質表現学習、知識グラフ推論のワークフローを構築・トラブルシューティングします。コードが torchdrug をインポートする場合や、そのデータセット、モデル、タスク、Engine が必要な場合に使用します。

TorchDrug とは?

このスキルは、TorchDrug 0.2.1 を使用して分子グラフ学習ワークフローを構築・トラブルシューティングするための包括的なガイドです。ClinTox → GIN → PropertyPrediction → Engine という正規の特性予測パターンから、自己教師あり事前学習、分子生成(GCPN、GraphAF)、逆合成(中心同定・シントン補完)、知識グラフ推論(RotatE、NeuralLP)、タンパク質モデリング(ESM、GearNet)まで、主要なワークフローを網羅しています。バージョン0.2.1特有のAPI、インストール手順(Python 3.7~3.10、PyTorch 1.8~2.0、CPU/Apple Silicon対応)、デバイス管理、チェックポイントの扱いなど、信頼性の高いコードを書くためのルールを提供します。主にClaude Code、Cursor、CodexなどのAIエージェント向けに設計されており、ユーザーはこのスキルを読み込むことで、TorchDrugに関する正確で再現性のあるコードを素早く生成できます。

対応~Claude Code~Codex CLI~Cursor
npx skills add https://github.com/K-Dense-AI/scientific-agent-skills/tree/main/skills/torchdrug

お気に入りのAIに質問する

このエージェントスキルを事前に読み込んだ状態で新しいチャットを開きます。

ドキュメント

TorchDrug は何をしますか?

TorchDrug をモジュール型 PyTorch グラフ学習スタックとして使用します。

  1. datasets.* データセットをロードし、
  2. models.* 表現モデルを選択し、
  3. それを tasks.* 目的関数でラップし、
  4. core.Engine で訓練・評価します。

現在の公式ドキュメントと最新リリースはともに 0.2.1 です。新しい Python や PyTorch の組み合わせは、互換性があると暗黙に想定せず、未検証として扱います。

バージョンガードから始める

コードを生成またはデバッグする前に、環境を確認します。

python --version
python -c "import torch; print(torch.__version__)"
python -c "import torchdrug; print(torchdrug.__version__)"

TorchDrug 0.2.1 のサポートマトリックスは以下の通りです。

  • Python 3.7 ~ 3.10
  • PyTorch 1.8 ~ 2.0
  • Linux、Windows、macOS
  • Apple Silicon: PyTorch 1.13 以降、CPU のみ。MPS は非対応

プロジェクトが Python 3.11+ または PyTorch 2.1+ を使用している場合は、互換性のある環境を作成するか、ソースビルドを明示的にテストしてください。そのような組み合わせをサポート対象として提示しないでください。

インストール

専用の Python 3.10 環境を推奨し、TorchDrug リリースを固定します。

uv venv --python 3.10
source .venv/bin/activate
uv pip install "torch==2.0.0"

torch-scatter と torch-cluster のホイールは、正確な PyTorch と CUDA のペアに合わせて、公式インストールページ に従ってインストールします。CPU 専用 PyTorch 2.0 環境の場合、再現可能なホイールの組み合わせの一つは次の通りです。

uv pip install "torch-scatter==2.1.1" "torch-cluster==1.6.1" \
  --find-links "https://data.pyg.org/whl/torch-2.0.0+cpu.html"
uv pip install "torchdrug==0.2.1"

CUDA ホイールの URL を環境間でコピーしないでください。PyTorch のバージョン、CUDA ビルド、Python ABI、プラットフォームを一致させてください。Apple Silicon では、公式ドキュメントに従い torch-scatter と torch-cluster をソースからビルドする必要があります。レビュー済みのソースリビジョンを固定し、CPU 実行を想定してください。

正規の特性予測ワークフロー

文書化された ClinTox → GIN → PropertyPrediction → Engine パターンを使用します。

import torch
from torchdrug import core, datasets, models, tasks

dataset = datasets.ClinTox("~/molecule-datasets/")
lengths = [int(0.8 * len(dataset)), int(0.1 * len(dataset))]
lengths.append(len(dataset) - sum(lengths))
train_set, valid_set, test_set = torch.utils.data.random_split(dataset, lengths)

model = models.GIN(
    input_dim=dataset.node_feature_dim,
    hidden_dims=[256, 256, 256, 256],
    short_cut=True,
    batch_norm=True,
    concat_hidden=True,
)
task = tasks.PropertyPrediction(
    model,
    task=dataset.tasks,
    criterion="bce",
    metric=("auprc", "auroc"),
)

optimizer = torch.optim.Adam(task.parameters(), lr=1e-3)
solver = core.Engine(
    task,
    train_set,
    valid_set,
    test_set,
    optimizer,
    batch_size=1024,
)
solver.train(num_epoch=100)
solver.evaluate("valid")

gpus=[0] は、サポートされている CUDA デバイスが利用可能な場合にのみ追加します。CPU 実行の場合は gpus を省略します。

二値分類の場合、task.predict(batch) はロジットを返します。確率が必要な場合は torch.sigmoid を適用します。0.2.1 では、正規化された回帰予測が元のターゲットスケールで返されます。これは以前のリリースからの破壊的変更です。

公式ワークフローを選択する

分子特性予測

  • データセット: datasets.ClinTox、BBBP、Tox21、QM9、またはその他文書化された分子データセット。
  • モデル: models.GIN から始めます。選択した特徴構成がエッジ特徴を提供する場合は edge_input_dim を使用します。
  • タスク: tasks.PropertyPrediction。
  • 参照: 分子特性予測

自己教師あり分子事前学習

  • InfoGraph: models.InfoGraph(gin_model, separate_model=False) を tasks.Unsupervised でラップします。
  • 属性マスキング: tasks.AttributeMasking(model, mask_rate=0.15)。
  • ファインチューニング用に同じエンコーダを再作成し、tasks.PropertyPrediction を訓練する前に strict=False でチェックポイントをロードします。
  • 参照: 分子特性予測

分子生成

  • データセット: datasets.ZINC250k(..., kekulize=True, atom_feature="symbol")。
  • GCPN: models.RGCN エンコーダを tasks.GCPNGeneration でラップします。
  • GraphAF: ノードとエッジの models.GraphAF フローを tasks.AutoregressiveGeneration でラップします。
  • チュートリアルでサポートされている最適化タスクは "qed" と "plogp" です。基準は "nll" および/または "ppo" です。
  • 参照: 分子生成

逆合成

  • 2つの同期された datasets.USPTO50k ビューを作成します。反応モード(中心同定用)と as_synthon=True(シントン補完用)です。
  • tasks.CenterIdentification と tasks.SynthonCompletion を個別に訓練します。
  • 訓練済みタスクを tasks.Retrosynthesis で結合します。生のモデルをエンドツーエンドのタスクに直接渡さないでください。
  • 参照: 逆合成

知識グラフ推論

  • 埋め込みワークフロー: datasets.FB15k237 → models.RotatE → tasks.KnowledgeGraphCompletion。
  • ニューラル推論ワークフロー: models.NeuralLP(fact_ratio=0.75 を使用)。
  • 参照: 知識グラフ推論

タンパク質モデリング

  • data.Protein.from_sequence、from_pdb、または from_molecule でタンパク質を構築します。
  • 配列エンコーダには models.ESM、ProteinCNN、ProteinResNet、ProteinLSTM、ProteinBERT があります。構造エンコーダには models.GearNet があります。
  • 存在しない protein.residue_graph() 便利メソッドではなく、文書化されたグラフ構築レイヤーを使用します。
  • 参照: タンパク質モデリング

信頼性の高い TorchDrug コードのルール

  1. 0.2.1 API に従う。 公式ドキュメントはローリング最新バージョンサイトではありません。
  2. 文書化された特徴名を優先する。 atom_feature、bond_feature、residue_feature、mol_feature を使用します。node_feature、edge_feature、graph_feature は関連データセットコンストラクタで非推奨のエイリアスです。
  3. Engine にタスクの前処理を任せる。 ソルバーを構築せずに事前訓練済みタスクを合成する場合は、各タスクの preprocess() を手動で呼び出します。
  4. ペアの分割を同期させる。 逆合成の場合、反応データセットとシントンデータセットを分割する前に、同じ乱数シードをリセットします。
  5. TorchDrug の照合を使用する。 data.graph_collate または core.Engine を使用します。汎用の PyTorch 照合は TorchDrug グラフのパック方法を知りません。
  6. モデル、タスク、エンジンの引数を分離する。 コードの誤生成の一般的な原因は、タスクオプションをモデルに渡したり、合成タスクが必要な場所に生のモデルを渡したりすることです。
  7. 生成された化学構造を検証する。 モデルの出力は、実験的に有効または合成可能な化合物ではなく、候補として扱います。

トラブルシューティング

インストールまたはインポートの失敗

Python、PyTorch、torch-scatter、torch-cluster を一つの互換性セットとして確認します。ほとんどの障害は、バイナリホイールの不一致、サポートされていない Python バージョン、または MPS の使用試行です。

特徴次元の不一致

ロードされたデータセットからモデルの次元を構築します。

  • dataset.node_feature_dim
  • dataset.edge_feature_dim
  • dataset.num_bond_type
  • 知識グラフの場合は dataset.num_entity と dataset.num_relation

異なる特徴構成からコピーした次元をハードコーディングしないでください。

デバイスの不一致

サポートされている CUDA 実行のために core.Engine に gpus=[0] を渡します。手動予測の場合は、まず照合し、ネストされたバッチ全体を utils.cuda で移動します。

チェックポイントの不一致

同じモデルと特徴構成を再作成します。事前学習からファインチューニングへの転送では、strict=False でチェックポイントの "model" 状態をロードします。完全なソルバーの場合は、solver.save() と solver.load() を使用します。

リファレンスインデックス

アップストリームソース

サイエンティフィックエージェントスキルの引用

このスキルは K-Dense による Scientific Agent Skills の一部です。これが原稿、レポート、プレゼンテーション、またはコードリリースに実質的に貢献した場合は、参考文献またはソフトウェアセクションに論文を追加し、その旨をユーザーに伝えてください。

Kassis, T., Agarwal, V., He, Y., Patel, D., & Brueckner, A. M. (2026). Scientific Agent Skills: A Library of Procedural Knowledge for Research Agents. arXiv:2609.00065. https://doi.org/10.48550/arXiv.2609.00065

常に最新バージョンを引用してください。DOI と https://arxiv.org/abs/2609.00065 は最新の arXiv バージョンに解決されるため、v1 のようなバージョン接尾辞を追加しないでください。ネットワークアクセスが可能な場合は、参照を書く前に https://arxiv.org/abs/2609.00065(または http://export.arxiv.org/api/query?id_list=2609.00065)を取得し、そのレコードから著者リスト、年、バージョンを取得してください。レコードにジャーナル参照や出版社 DOI が記載されている場合は、代わりに公開バージョンを引用してください。

Individual skills in this repo

This repo contains 20 individual skills — each has its own dedicated page.

adaptyv

How to use the Adaptyv Bio Foundry API and Python SDK for protein experiment design, submission, and results retrieval. Use this skill whenever the user mentions Adaptyv, Foundry API, protein binding assays, protein screening experiments, BLI/SPR assays, thermostability assays, or wants to submit protein sequences for experimental characterization. Also trigger when code imports `adaptyv`, `adaptyv_sdk`, or `FoundryClient`, or references `foundry-api-public.adaptyvbio.com`.

aeon

This skill should be used for time series machine learning tasks including classification, regression, clustering, forecasting, anomaly detection, segmentation, and similarity search. Use when working with temporal data, sequential patterns, or time-indexed observations requiring specialized algorithms beyond standard ML approaches. Particularly suited for univariate and multivariate time series analysis with scikit-learn compatible APIs.

alphagenome

Look up precomputed AlphaGenome Atlas effects for any GRCh38 single-nucleotide variant (AVI score with Phred and 18 SHAP feature attributions, plus raw and quantile scores for RNA-seq, DNase, ATAC, ChIP-TF, ChIP-histone, CAGE, PRO-cap, splicing, polyadenylation and contact-map tracks), score variants or scan windows on demand with the AlphaGenome model for human and mouse (variant scoring, in silico mutagenesis, REF-versus-ALT track prediction), and build Atlas website deep links. Use when the user mentions AlphaGenome, AlphaGenome Atlas, AVI or AlphaGenome Variant Impact, DeepMind variant effect prediction, or wants to prioritise or mechanistically interpret non-coding, regulatory, splicing, enhancer, promoter, or chromatin-accessibility effects of SNVs from a VCF, credible set, or region. Research use only; not a clinical tool.

analytical-method-validation

Plan, execute, and document validation, verification, and transfer of analytical procedures under the governing framework - ICH Q2(R2) and Q14, USP <1220>/<1225>/<1226>, ICH M10 bioanalytical, CLSI EP, or ISO/IEC 17025. Use for HPLC, LC-MS/MS, GC, CE, ICP-MS, dissolution, qNMR, qPCR, NIR, and ligand binding or cell-based assays whenever the question is whether a procedure is fit for its intended purpose. Triggers include

anndata

Data structure for annotated matrices in single-cell analysis. Use when working with .h5ad files or integrating with the scverse ecosystem. This is the data format skill—for analysis workflows use scanpy; for probabilistic models use scvi-tools; for population-scale queries use cellxgene-census.

arbor

Autonomously improve a real artifact (code, training recipe, agent harness, data pipeline, prompt) against an objective and an evaluator, using Hypothesis Tree Refinement (HTR) from the Arbor paper. Use this whenever someone wants to iteratively optimize something over many experiments without overfitting — e.g.

arboreto

Infer gene regulatory networks (GRNs) from gene expression data using scalable algorithms (GRNBoost2, GENIE3). Use when analyzing transcriptomics data (bulk RNA-seq, single-cell RNA-seq) to identify transcription factor-target gene relationships and regulatory interactions. Supports distributed computation for large-scale datasets.

astropy

Core Python library for astronomy and astrophysics workflows that need Astropy APIs, including units/quantities, coordinates, FITS I/O, tables, time systems, WCS, and cosmology. Use when implementing or debugging astronomical data analysis code with Astropy.

autoskill

Observe the user

benchling-integration

Benchling Python SDK and REST API integration for registry entities, inventory, ELN entries, workflows, Benchling Apps, and Data Warehouse queries. Use when automating lab data with benchling-sdk or the v2 API.

bgpt-paper-search

Search scientific papers and retrieve structured experimental data extracted from full-text studies via the BGPT MCP server. Returns 25+ fields per paper including methods, results, sample sizes, quality scores, and conclusions. Use for literature reviews, evidence synthesis, and finding experimental details not available in abstracts alone.

bids

>

biopython

Comprehensive molecular biology toolkit. Use for sequence manipulation, file parsing (FASTA/GenBank/PDB), phylogenetics, and programmatic NCBI/PubMed access (Bio.Entrez). Best for batch processing, custom bioinformatics pipelines, BLAST automation. For quick lookups use gget; for multi-service integration use bioservices.

bioservices

Unified Python interface to 40+ bioinformatics services. Use when querying multiple databases (UniProt, KEGG, ChEMBL, Reactome) in a single workflow with consistent API. Best for cross-database analysis, ID mapping across services. For quick single-database lookups use gget; for sequence/file manipulation use biopython.

bulk-rnaseq

End-to-end bulk RNA-seq orchestrator — takes raw FASTQ reads through QC and trimming (FastQC, fastp/Trim Galore), alignment and quantification (STAR, Salmon, featureCounts), assembles a gene-level counts matrix, then hands off to differential expression (pydeseq2), pathway/GSEA enrichment (pathway-enrichment), and publication figures (scientific-visualization). Use whenever the user has bulk RNA-seq reads or quant output and wants a complete, reproducible differential-expression workflow — e.g.

cellxgene-census

Query the CZ CELLxGENE Census programmatically for versioned public single-cell and spatial transcriptomics data. Use when you need population-scale cell metadata, gene expression slices, Census summary counts, source H5AD URIs/downloads, embeddings, spatial Census data, or reference atlas comparisons across organisms, tissues, diseases, assays, and cell types. For analyzing your own local single-cell data use scanpy, anndata, or scvi-tools.

cirq

Google quantum computing framework. Use when targeting Google Quantum AI hardware, designing noise-aware circuits, or running quantum characterization experiments. Best for Google hardware, noise modeling, and low-level circuit design. For IBM hardware use qiskit; for quantum ML with autodiff use pennylane; for physics simulations use qutip.

citation-management

Comprehensive citation management for academic research. Search OpenAlex, PubMed, and Google Scholar for papers, extract accurate metadata, validate citations, and generate properly formatted BibTeX entries. This skill should be used when you need to find papers, verify citation information, convert DOIs to BibTeX, or ensure reference accuracy in scientific writing.

clinical-decision-support

Prepare and validate research-only clinical decision-support evaluation, evidence-profile, cohort, survival, biomarker/model, privacy, and governance artifacts. Use for aggregate or synthetic research documentation and traceability—not patient care or live clinical operation.

clinical-reports

Create safety-bounded draft structures and run local deterministic checks for clinical case, diagnostic, trial, safety, and aggregate research reports. Use only with synthetic, de-identified, or aggregate inputs and verified source-fact manifests; every output requires qualified review.

関連スキル