[{"data":1,"prerenderedAt":-1},["ShallowReactive",2],{"skill:01M6AFJZ340ZXP1AADVXFE3YXD":3,"skill-same:01M6AFJZ340ZXP1AADVXFE3YXD":98},{"skillNo":4,"name":5,"nameZh":6,"description":7,"descriptionZh":6,"repositoryFullName":8,"ownerLogin":9,"ownerAvatarUrl":10,"ownerAccountNo":11,"ownerDisplayName":12,"repositoryNo":13,"starsCount":14,"licenseSpdxId":15,"licenseFileUrl":6,"sourceHtmlUrl":16,"categories":17,"keywords":23,"publishedAt":39,"linkOnly":40,"markdownContent":41,"markdownBody":42,"frontMatter":43,"refName":44,"commitSha":45,"skillDirectory":46,"forksCount":47,"repositoryHtmlUrl":48,"repositoryHomepageUrl":49,"ownerHtmlUrl":50,"ownerBlogUrl":51,"ownerDescription":52,"updatedAt":53,"install":54,"resources":67},"01M6AFJZ340ZXP1AADVXFE3YXD","pydeseq2",null,"Differential gene expression analysis (Python DESeq2). Identify DE genes from bulk RNA-seq counts, Wald tests, FDR correction, volcano\u002FMA plots, for RNA-seq analysis.","davila7\u002Fclaude-code-templates","davila7","https:\u002F\u002Favatars.githubusercontent.com\u002Fu\u002F6216945?v=4","01M6ADEWWFPS191NFM0Y1EX1JH","Daniel Avila","01M6ADJHGM78PHQRD3XZTFDVKK",33481,"MIT","https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates\u002Fblob\u002Fmain\u002Fcli-tool\u002Fcomponents\u002Fskills\u002Fscientific\u002Fpydeseq2\u002FSKILL.md",[18],{"categoryNo":19,"categoryName":20,"taxonomyType":21,"taxonomyTypeDesc":22},"01JBSEEDCAT000000000000023","Skill与智能体开发","SCENARIO","使用场景",[5,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38],"python","analysis","RNA","gene expression","tests","expression","DESeq2","volcano","de","rna-seq","bulk","seq","for","gene","fdr","2026-11-10 10:12:30",false,"---\nname: pydeseq2\ndescription: \"Differential gene expression analysis (Python DESeq2). Identify DE genes from bulk RNA-seq counts, Wald tests, FDR correction, volcano\u002FMA plots, for RNA-seq analysis.\"\n---\n\n# PyDESeq2\n\n## Overview\n\nPyDESeq2 is a Python implementation of DESeq2 for differential expression analysis with bulk RNA-seq data. Design and execute complete workflows from data loading through result interpretation, including single-factor and multi-factor designs, Wald tests with multiple testing correction, optional apeGLM shrinkage, and integration with pandas and AnnData.\n\n## When to Use This Skill\n\nThis skill should be used when:\n- Analyzing bulk RNA-seq count data for differential expression\n- Comparing gene expression between experimental conditions (e.g., treated vs control)\n- Performing multi-factor designs accounting for batch effects or covariates\n- Converting R-based DESeq2 workflows to Python\n- Integrating differential expression analysis into Python-based pipelines\n- Users mention \"DESeq2\", \"differential expression\", \"RNA-seq analysis\", or \"PyDESeq2\"\n\n## Quick Start Workflow\n\nFor users who want to perform a standard differential expression analysis:\n\n```python\nimport pandas as pd\nfrom pydeseq2.dds import DeseqDataSet\nfrom pydeseq2.ds import DeseqStats\n\n# 1. Load data\ncounts_df = pd.read_csv(\"counts.csv\", index_col=0).T  # Transpose to samples × genes\nmetadata = pd.read_csv(\"metadata.csv\", index_col=0)\n\n# 2. Filter low-count genes\ngenes_to_keep = counts_df.columns[counts_df.sum(axis=0) >= 10]\ncounts_df = counts_df[genes_to_keep]\n\n# 3. Initialize and fit DESeq2\ndds = DeseqDataSet(\n    counts=counts_df,\n    metadata=metadata,\n    design=\"~condition\",\n    refit_cooks=True\n)\ndds.deseq2()\n\n# 4. Perform statistical testing\nds = DeseqStats(dds, contrast=[\"condition\", \"treated\", \"control\"])\nds.summary()\n\n# 5. Access results\nresults = ds.results_df\nsignificant = results[results.padj \u003C 0.05]\nprint(f\"Found {len(significant)} significant genes\")\n```\n\n## Core Workflow Steps\n\n### Step 1: Data Preparation\n\n**Input requirements:**\n- **Count matrix:** Samples × genes DataFrame with non-negative integer read counts\n- **Metadata:** Samples × variables DataFrame with experimental factors\n\n**Common data loading patterns:**\n\n```python\n# From CSV (typical format: genes × samples, needs transpose)\ncounts_df = pd.read_csv(\"counts.csv\", index_col=0).T\nmetadata = pd.read_csv(\"metadata.csv\", index_col=0)\n\n# From TSV\ncounts_df = pd.read_csv(\"counts.tsv\", sep=\"\\t\", index_col=0).T\n\n# From AnnData\nimport anndata as ad\nadata = ad.read_h5ad(\"data.h5ad\")\ncounts_df = pd.DataFrame(adata.X, index=adata.obs_names, columns=adata.var_names)\nmetadata = adata.obs\n```\n\n**Data filtering:**\n\n```python\n# Remove low-count genes\ngenes_to_keep = counts_df.columns[counts_df.sum(axis=0) >= 10]\ncounts_df = counts_df[genes_to_keep]\n\n# Remove samples with missing metadata\nsamples_to_keep = ~metadata.condition.isna()\ncounts_df = counts_df.loc[samples_to_keep]\nmetadata = metadata.loc[samples_to_keep]\n```\n\n### Step 2: Design Specification\n\nThe design formula specifies how gene expression is modeled.\n\n**Single-factor designs:**\n```python\ndesign = \"~condition\"  # Simple two-group comparison\n```\n\n**Multi-factor designs:**\n```python\ndesign = \"~batch + condition\"  # Control for batch effects\ndesign = \"~age + condition\"     # Include continuous covariate\ndesign = \"~group + condition + group:condition\"  # Interaction effects\n```\n\n**Design formula guidelines:**\n- Use Wilkinson formula notation (R-style)\n- Put adjustment variables (e.g., batch) before the main variable of interest\n- Ensure variables exist as columns in the metadata DataFrame\n- Use appropriate data types (categorical for discrete variables)\n\n### Step 3: DESeq2 Fitting\n\nInitialize the DeseqDataSet and run the complete pipeline:\n\n```python\nfrom pydeseq2.dds import DeseqDataSet\n\ndds = DeseqDataSet(\n    counts=counts_df,\n    metadata=metadata,\n    design=\"~condition\",\n    refit_cooks=True,  # Refit after removing outliers\n    n_cpus=1           # Parallel processing (adjust as needed)\n)\n\n# Run the complete DESeq2 pipeline\ndds.deseq2()\n```\n\n**What `deseq2()` does:**\n1. Computes size factors (normalization)\n2. Fits genewise dispersions\n3. Fits dispersion trend curve\n4. Computes dispersion priors\n5. Fits MAP dispersions (shrinkage)\n6. Fits log fold changes\n7. Calculates Cook's distances (outlier detection)\n8. Refits if outliers detected (optional)\n\n### Step 4: Statistical Testing\n\nPerform Wald tests to identify differentially expressed genes:\n\n```python\nfrom pydeseq2.ds import DeseqStats\n\nds = DeseqStats(\n    dds,\n    contrast=[\"condition\", \"treated\", \"control\"],  # Test treated vs control\n    alpha=0.05,                # Significance threshold\n    cooks_filter=True,         # Filter outliers\n    independent_filter=True    # Filter low-power tests\n)\n\nds.summary()\n```\n\n**Contrast specification:**\n- Format: `[variable, test_level, reference_level]`\n- Example: `[\"condition\", \"treated\", \"control\"]` tests treated vs control\n- If `None`, uses the last coefficient in the design\n\n**Result DataFrame columns:**\n- `baseMean`: Mean normalized count across samples\n- `log2FoldChange`: Log2 fold change between conditions\n- `lfcSE`: Standard error of LFC\n- `stat`: Wald test statistic\n- `pvalue`: Raw p-value\n- `padj`: Adjusted p-value (FDR-corrected via Benjamini-Hochberg)\n\n### Step 5: Optional LFC Shrinkage\n\nApply shrinkage to reduce noise in fold change estimates:\n\n```python\nds.lfc_shrink()  # Applies apeGLM shrinkage\n```\n\n**When to use LFC shrinkage:**\n- For visualization (volcano plots, heatmaps)\n- For ranking genes by effect size\n- When prioritizing genes for follow-up experiments\n\n**Important:** Shrinkage affects only the log2FoldChange values, not the statistical test results (p-values remain unchanged). Use shrunk values for visualization but report unshrunken p-values for significance.\n\n### Step 6: Result Export\n\nSave results and intermediate objects:\n\n```python\nimport pickle\n\n# Export results as CSV\nds.results_df.to_csv(\"deseq2_results.csv\")\n\n# Save significant genes only\nsignificant = ds.results_df[ds.results_df.padj \u003C 0.05]\nsignificant.to_csv(\"significant_genes.csv\")\n\n# Save DeseqDataSet for later use\nwith open(\"dds_result.pkl\", \"wb\") as f:\n    pickle.dump(dds.to_picklable_anndata(), f)\n```\n\n## Common Analysis Patterns\n\n### Two-Group Comparison\n\nStandard case-control comparison:\n\n```python\ndds = DeseqDataSet(counts=counts_df, metadata=metadata, design=\"~condition\")\ndds.deseq2()\n\nds = DeseqStats(dds, contrast=[\"condition\", \"treated\", \"control\"])\nds.summary()\n\nresults = ds.results_df\nsignificant = results[results.padj \u003C 0.05]\n```\n\n### Multiple Comparisons\n\nTesting multiple treatment groups against control:\n\n```python\ndds = DeseqDataSet(counts=counts_df, metadata=metadata, design=\"~condition\")\ndds.deseq2()\n\ntreatments = [\"treatment_A\", \"treatment_B\", \"treatment_C\"]\nall_results = {}\n\nfor treatment in treatments:\n    ds = DeseqStats(dds, contrast=[\"condition\", treatment, \"control\"])\n    ds.summary()\n    all_results[treatment] = ds.results_df\n\n    sig_count = len(ds.results_df[ds.results_df.padj \u003C 0.05])\n    print(f\"{treatment}: {sig_count} significant genes\")\n```\n\n### Accounting for Batch Effects\n\nControl for technical variation:\n\n```python\n# Include batch in design\ndds = DeseqDataSet(counts=counts_df, metadata=metadata, design=\"~batch + condition\")\ndds.deseq2()\n\n# Test condition while controlling for batch\nds = DeseqStats(dds, contrast=[\"condition\", \"treated\", \"control\"])\nds.summary()\n```\n\n### Continuous Covariates\n\nInclude continuous variables like age or dosage:\n\n```python\n# Ensure continuous variable is numeric\nmetadata[\"age\"] = pd.to_numeric(metadata[\"age\"])\n\ndds = DeseqDataSet(counts=counts_df, metadata=metadata, design=\"~age + condition\")\ndds.deseq2()\n\nds = DeseqStats(dds, contrast=[\"condition\", \"treated\", \"control\"])\nds.summary()\n```\n\n## Using the Analysis Script\n\nThis skill includes a complete command-line script for standard analyses:\n\n```bash\n# Basic usage\npython scripts\u002Frun_deseq2_analysis.py \\\n  --counts counts.csv \\\n  --metadata metadata.csv \\\n  --design \"~condition\" \\\n  --contrast condition treated control \\\n  --output results\u002F\n\n# With additional options\npython scripts\u002Frun_deseq2_analysis.py \\\n  --counts counts.csv \\\n  --metadata metadata.csv \\\n  --design \"~batch + condition\" \\\n  --contrast condition treated control \\\n  --output results\u002F \\\n  --min-counts 10 \\\n  --alpha 0.05 \\\n  --n-cpus 4 \\\n  --plots\n```\n\n**Script features:**\n- Automatic data loading and validation\n- Gene and sample filtering\n- Complete DESeq2 pipeline execution\n- Statistical testing with customizable parameters\n- Result export (CSV, pickle)\n- Optional visualization (volcano and MA plots)\n\nRefer users to `scripts\u002Frun_deseq2_analysis.py` when they need a standalone analysis tool or want to batch process multiple datasets.\n\n## Result Interpretation\n\n### Identifying Significant Genes\n\n```python\n# Filter by adjusted p-value\nsignificant = ds.results_df[ds.results_df.padj \u003C 0.05]\n\n# Filter by both significance and effect size\nsig_and_large = ds.results_df[\n    (ds.results_df.padj \u003C 0.05) &\n    (abs(ds.results_df.log2FoldChange) > 1)\n]\n\n# Separate up- and down-regulated\nupregulated = significant[significant.log2FoldChange > 0]\ndownregulated = significant[significant.log2FoldChange \u003C 0]\n\nprint(f\"Upregulated: {len(upregulated)}\")\nprint(f\"Downregulated: {len(downregulated)}\")\n```\n\n### Ranking and Sorting\n\n```python\n# Sort by adjusted p-value\ntop_by_padj = ds.results_df.sort_values(\"padj\").head(20)\n\n# Sort by absolute fold change (use shrunk values)\nds.lfc_shrink()\nds.results_df[\"abs_lfc\"] = abs(ds.results_df.log2FoldChange)\ntop_by_lfc = ds.results_df.sort_values(\"abs_lfc\", ascending=False).head(20)\n\n# Sort by a combined metric\nds.results_df[\"score\"] = -np.log10(ds.results_df.padj) * abs(ds.results_df.log2FoldChange)\ntop_combined = ds.results_df.sort_values(\"score\", ascending=False).head(20)\n```\n\n### Quality Metrics\n\n```python\n# Check normalization (size factors should be close to 1)\nprint(\"Size factors:\", dds.obsm[\"size_factors\"])\n\n# Examine dispersion estimates\nimport matplotlib.pyplot as plt\nplt.hist(dds.varm[\"dispersions\"], bins=50)\nplt.xlabel(\"Dispersion\")\nplt.ylabel(\"Frequency\")\nplt.title(\"Dispersion Distribution\")\nplt.show()\n\n# Check p-value distribution (should be mostly flat with peak near 0)\nplt.hist(ds.results_df.pvalue.dropna(), bins=50)\nplt.xlabel(\"P-value\")\nplt.ylabel(\"Frequency\")\nplt.title(\"P-value Distribution\")\nplt.show()\n```\n\n## Visualization Guidelines\n\n### Volcano Plot\n\nVisualize significance vs effect size:\n\n```python\nimport matplotlib.pyplot as plt\nimport numpy as np\n\nresults = ds.results_df.copy()\nresults[\"-log10(padj)\"] = -np.log10(results.padj)\n\nplt.figure(figsize=(10, 6))\nsignificant = results.padj \u003C 0.05\n\nplt.scatter(\n    results.loc[~significant, \"log2FoldChange\"],\n    results.loc[~significant, \"-log10(padj)\"],\n    alpha=0.3, s=10, c='gray', label='Not significant'\n)\nplt.scatter(\n    results.loc[significant, \"log2FoldChange\"],\n    results.loc[significant, \"-log10(padj)\"],\n    alpha=0.6, s=10, c='red', label='padj \u003C 0.05'\n)\n\nplt.axhline(-np.log10(0.05), color='blue', linestyle='--', alpha=0.5)\nplt.xlabel(\"Log2 Fold Change\")\nplt.ylabel(\"-Log10(Adjusted P-value)\")\nplt.title(\"Volcano Plot\")\nplt.legend()\nplt.savefig(\"volcano_plot.png\", dpi=300)\n```\n\n### MA Plot\n\nShow fold change vs mean expression:\n\n```python\nplt.figure(figsize=(10, 6))\n\nplt.scatter(\n    np.log10(results.loc[~significant, \"baseMean\"] + 1),\n    results.loc[~significant, \"log2FoldChange\"],\n    alpha=0.3, s=10, c='gray'\n)\nplt.scatter(\n    np.log10(results.loc[significant, \"baseMean\"] + 1),\n    results.loc[significant, \"log2FoldChange\"],\n    alpha=0.6, s=10, c='red'\n)\n\nplt.axhline(0, color='blue', linestyle='--', alpha=0.5)\nplt.xlabel(\"Log10(Base Mean + 1)\")\nplt.ylabel(\"Log2 Fold Change\")\nplt.title(\"MA Plot\")\nplt.savefig(\"ma_plot.png\", dpi=300)\n```\n\n## Troubleshooting Common Issues\n\n### Data Format Problems\n\n**Issue:** \"Index mismatch between counts and metadata\"\n\n**Solution:** Ensure sample names match exactly\n```python\nprint(\"Counts samples:\", counts_df.index.tolist())\nprint(\"Metadata samples:\", metadata.index.tolist())\n\n# Take intersection if needed\ncommon = counts_df.index.intersection(metadata.index)\ncounts_df = counts_df.loc[common]\nmetadata = metadata.loc[common]\n```\n\n**Issue:** \"All genes have zero counts\"\n\n**Solution:** Check if data needs transposition\n```python\nprint(f\"Counts shape: {counts_df.shape}\")\n# If genes > samples, transpose is needed\nif counts_df.shape[1] \u003C counts_df.shape[0]:\n    counts_df = counts_df.T\n```\n\n### Design Matrix Issues\n\n**Issue:** \"Design matrix is not full rank\"\n\n**Cause:** Confounded variables (e.g., all treated samples in one batch)\n\n**Solution:** Remove confounded variable or add interaction term\n```python\n# Check confounding\nprint(pd.crosstab(metadata.condition, metadata.batch))\n\n# Either simplify design or add interaction\ndesign = \"~condition\"  # Remove batch\n# OR\ndesign = \"~condition + batch + condition:batch\"  # Model interaction\n```\n\n### No Significant Genes\n\n**Diagnostics:**\n```python\n# Check dispersion distribution\nplt.hist(dds.varm[\"dispersions\"], bins=50)\nplt.show()\n\n# Check size factors\nprint(dds.obsm[\"size_factors\"])\n\n# Look at top genes by raw p-value\nprint(ds.results_df.nsmallest(20, \"pvalue\"))\n```\n\n**Possible causes:**\n- Small effect sizes\n- High biological variability\n- Insufficient sample size\n- Technical issues (batch effects, outliers)\n\n## Reference Documentation\n\nFor comprehensive details beyond this workflow-oriented guide:\n\n- **API Reference** (`references\u002Fapi_reference.md`): Complete documentation of PyDESeq2 classes, methods, and data structures. Use when needing detailed parameter information or understanding object attributes.\n\n- **Workflow Guide** (`references\u002Fworkflow_guide.md`): In-depth guide covering complete analysis workflows, data loading patterns, multi-factor designs, troubleshooting, and best practices. Use when handling complex experimental designs or encountering issues.\n\nLoad these references into context when users need:\n- Detailed API documentation: `Read references\u002Fapi_reference.md`\n- Comprehensive workflow examples: `Read references\u002Fworkflow_guide.md`\n- Troubleshooting guidance: `Read references\u002Fworkflow_guide.md` (see Troubleshooting section)\n\n## Key Reminders\n\n1. **Data orientation matters:** Count matrices typically load as genes × samples but need to be samples × genes. Always transpose with `.T` if needed.\n\n2. **Sample filtering:** Remove samples with missing metadata before analysis to avoid errors.\n\n3. **Gene filtering:** Filter low-count genes (e.g., \u003C 10 total reads) to improve power and reduce computational time.\n\n4. **Design formula order:** Put adjustment variables before the variable of interest (e.g., `\"~batch + condition\"` not `\"~condition + batch\"`).\n\n5. **LFC shrinkage timing:** Apply shrinkage after statistical testing and only for visualization\u002Franking purposes. P-values remain based on unshrunken estimates.\n\n6. **Result interpretation:** Use `padj \u003C 0.05` for significance, not raw p-values. The Benjamini-Hochberg procedure controls false discovery rate.\n\n7. **Contrast specification:** The format is `[variable, test_level, reference_level]` where test_level is compared against reference_level.\n\n8. **Save intermediate objects:** Use pickle to save DeseqDataSet objects for later use or additional analyses without re-running the expensive fitting step.\n\n## Installation and Requirements\n\n```bash\nuv pip install pydeseq2\n```\n\n**System requirements:**\n- Python 3.10-3.11\n- pandas 1.4.3+\n- numpy 1.23.0+\n- scipy 1.11.0+\n- scikit-learn 1.1.1+\n- anndata 0.8.0+\n\n**Optional for visualization:**\n- matplotlib\n- seaborn\n\n## Additional Resources\n\n- **Official Documentation:** https:\u002F\u002Fpydeseq2.readthedocs.io\n- **GitHub Repository:** https:\u002F\u002Fgithub.com\u002Fowkin\u002FPyDESeq2\n- **Publication:** Muzellec et al. (2023) Bioinformatics, DOI: 10.1093\u002Fbioinformatics\u002Fbtad547\n- **Original DESeq2 (R):** Love et al. (2014) Genome Biology, DOI: 10.1186\u002Fs13059-014-0550-8\n","# PyDESeq2\n\n## Overview\n\nPyDESeq2 is a Python implementation of DESeq2 for differential expression analysis with bulk RNA-seq data. Design and execute complete workflows from data loading through result interpretation, including single-factor and multi-factor designs, Wald tests with multiple testing correction, optional apeGLM shrinkage, and integration with pandas and AnnData.\n\n## When to Use This Skill\n\nThis skill should be used when:\n- Analyzing bulk RNA-seq count data for differential expression\n- Comparing gene expression between experimental conditions (e.g., treated vs control)\n- Performing multi-factor designs accounting for batch effects or covariates\n- Converting R-based DESeq2 workflows to Python\n- Integrating differential expression analysis into Python-based pipelines\n- Users mention \"DESeq2\", \"differential expression\", \"RNA-seq analysis\", or \"PyDESeq2\"\n\n## Quick Start Workflow\n\nFor users who want to perform a standard differential expression analysis:\n\n```python\nimport pandas as pd\nfrom pydeseq2.dds import DeseqDataSet\nfrom pydeseq2.ds import DeseqStats\n\n# 1. Load data\ncounts_df = pd.read_csv(\"counts.csv\", index_col=0).T  # Transpose to samples × genes\nmetadata = pd.read_csv(\"metadata.csv\", index_col=0)\n\n# 2. Filter low-count genes\ngenes_to_keep = counts_df.columns[counts_df.sum(axis=0) >= 10]\ncounts_df = counts_df[genes_to_keep]\n\n# 3. Initialize and fit DESeq2\ndds = DeseqDataSet(\n    counts=counts_df,\n    metadata=metadata,\n    design=\"~condition\",\n    refit_cooks=True\n)\ndds.deseq2()\n\n# 4. Perform statistical testing\nds = DeseqStats(dds, contrast=[\"condition\", \"treated\", \"control\"])\nds.summary()\n\n# 5. Access results\nresults = ds.results_df\nsignificant = results[results.padj \u003C 0.05]\nprint(f\"Found {len(significant)} significant genes\")\n```\n\n## Core Workflow Steps\n\n### Step 1: Data Preparation\n\n**Input requirements:**\n- **Count matrix:** Samples × genes DataFrame with non-negative integer read counts\n- **Metadata:** Samples × variables DataFrame with experimental factors\n\n**Common data loading patterns:**\n\n```python\n# From CSV (typical format: genes × samples, needs transpose)\ncounts_df = pd.read_csv(\"counts.csv\", index_col=0).T\nmetadata = pd.read_csv(\"metadata.csv\", index_col=0)\n\n# From TSV\ncounts_df = pd.read_csv(\"counts.tsv\", sep=\"\\t\", index_col=0).T\n\n# From AnnData\nimport anndata as ad\nadata = ad.read_h5ad(\"data.h5ad\")\ncounts_df = pd.DataFrame(adata.X, index=adata.obs_names, columns=adata.var_names)\nmetadata = adata.obs\n```\n\n**Data filtering:**\n\n```python\n# Remove low-count genes\ngenes_to_keep = counts_df.columns[counts_df.sum(axis=0) >= 10]\ncounts_df = counts_df[genes_to_keep]\n\n# Remove samples with missing metadata\nsamples_to_keep = ~metadata.condition.isna()\ncounts_df = counts_df.loc[samples_to_keep]\nmetadata = metadata.loc[samples_to_keep]\n```\n\n### Step 2: Design Specification\n\nThe design formula specifies how gene expression is modeled.\n\n**Single-factor designs:**\n```python\ndesign = \"~condition\"  # Simple two-group comparison\n```\n\n**Multi-factor designs:**\n```python\ndesign = \"~batch + condition\"  # Control for batch effects\ndesign = \"~age + condition\"     # Include continuous covariate\ndesign = \"~group + condition + group:condition\"  # Interaction effects\n```\n\n**Design formula guidelines:**\n- Use Wilkinson formula notation (R-style)\n- Put adjustment variables (e.g., batch) before the main variable of interest\n- Ensure variables exist as columns in the metadata DataFrame\n- Use appropriate data types (categorical for discrete variables)\n\n### Step 3: DESeq2 Fitting\n\nInitialize the DeseqDataSet and run the complete pipeline:\n\n```python\nfrom pydeseq2.dds import DeseqDataSet\n\ndds = DeseqDataSet(\n    counts=counts_df,\n    metadata=metadata,\n    design=\"~condition\",\n    refit_cooks=True,  # Refit after removing outliers\n    n_cpus=1           # Parallel processing (adjust as needed)\n)\n\n# Run the complete DESeq2 pipeline\ndds.deseq2()\n```\n\n**What `deseq2()` does:**\n1. Computes size factors (normalization)\n2. Fits genewise dispersions\n3. Fits dispersion trend curve\n4. Computes dispersion priors\n5. Fits MAP dispersions (shrinkage)\n6. Fits log fold changes\n7. Calculates Cook's distances (outlier detection)\n8. Refits if outliers detected (optional)\n\n### Step 4: Statistical Testing\n\nPerform Wald tests to identify differentially expressed genes:\n\n```python\nfrom pydeseq2.ds import DeseqStats\n\nds = DeseqStats(\n    dds,\n    contrast=[\"condition\", \"treated\", \"control\"],  # Test treated vs control\n    alpha=0.05,                # Significance threshold\n    cooks_filter=True,         # Filter outliers\n    independent_filter=True    # Filter low-power tests\n)\n\nds.summary()\n```\n\n**Contrast specification:**\n- Format: `[variable, test_level, reference_level]`\n- Example: `[\"condition\", \"treated\", \"control\"]` tests treated vs control\n- If `None`, uses the last coefficient in the design\n\n**Result DataFrame columns:**\n- `baseMean`: Mean normalized count across samples\n- `log2FoldChange`: Log2 fold change between conditions\n- `lfcSE`: Standard error of LFC\n- `stat`: Wald test statistic\n- `pvalue`: Raw p-value\n- `padj`: Adjusted p-value (FDR-corrected via Benjamini-Hochberg)\n\n### Step 5: Optional LFC Shrinkage\n\nApply shrinkage to reduce noise in fold change estimates:\n\n```python\nds.lfc_shrink()  # Applies apeGLM shrinkage\n```\n\n**When to use LFC shrinkage:**\n- For visualization (volcano plots, heatmaps)\n- For ranking genes by effect size\n- When prioritizing genes for follow-up experiments\n\n**Important:** Shrinkage affects only the log2FoldChange values, not the statistical test results (p-values remain unchanged). Use shrunk values for visualization but report unshrunken p-values for significance.\n\n### Step 6: Result Export\n\nSave results and intermediate objects:\n\n```python\nimport pickle\n\n# Export results as CSV\nds.results_df.to_csv(\"deseq2_results.csv\")\n\n# Save significant genes only\nsignificant = ds.results_df[ds.results_df.padj \u003C 0.05]\nsignificant.to_csv(\"significant_genes.csv\")\n\n# Save DeseqDataSet for later use\nwith open(\"dds_result.pkl\", \"wb\") as f:\n    pickle.dump(dds.to_picklable_anndata(), f)\n```\n\n## Common Analysis Patterns\n\n### Two-Group Comparison\n\nStandard case-control comparison:\n\n```python\ndds = DeseqDataSet(counts=counts_df, metadata=metadata, design=\"~condition\")\ndds.deseq2()\n\nds = DeseqStats(dds, contrast=[\"condition\", \"treated\", \"control\"])\nds.summary()\n\nresults = ds.results_df\nsignificant = results[results.padj \u003C 0.05]\n```\n\n### Multiple Comparisons\n\nTesting multiple treatment groups against control:\n\n```python\ndds = DeseqDataSet(counts=counts_df, metadata=metadata, design=\"~condition\")\ndds.deseq2()\n\ntreatments = [\"treatment_A\", \"treatment_B\", \"treatment_C\"]\nall_results = {}\n\nfor treatment in treatments:\n    ds = DeseqStats(dds, contrast=[\"condition\", treatment, \"control\"])\n    ds.summary()\n    all_results[treatment] = ds.results_df\n\n    sig_count = len(ds.results_df[ds.results_df.padj \u003C 0.05])\n    print(f\"{treatment}: {sig_count} significant genes\")\n```\n\n### Accounting for Batch Effects\n\nControl for technical variation:\n\n```python\n# Include batch in design\ndds = DeseqDataSet(counts=counts_df, metadata=metadata, design=\"~batch + condition\")\ndds.deseq2()\n\n# Test condition while controlling for batch\nds = DeseqStats(dds, contrast=[\"condition\", \"treated\", \"control\"])\nds.summary()\n```\n\n### Continuous Covariates\n\nInclude continuous variables like age or dosage:\n\n```python\n# Ensure continuous variable is numeric\nmetadata[\"age\"] = pd.to_numeric(metadata[\"age\"])\n\ndds = DeseqDataSet(counts=counts_df, metadata=metadata, design=\"~age + condition\")\ndds.deseq2()\n\nds = DeseqStats(dds, contrast=[\"condition\", \"treated\", \"control\"])\nds.summary()\n```\n\n## Using the Analysis Script\n\nThis skill includes a complete command-line script for standard analyses:\n\n```bash\n# Basic usage\npython scripts\u002Frun_deseq2_analysis.py \\\n  --counts counts.csv \\\n  --metadata metadata.csv \\\n  --design \"~condition\" \\\n  --contrast condition treated control \\\n  --output results\u002F\n\n# With additional options\npython scripts\u002Frun_deseq2_analysis.py \\\n  --counts counts.csv \\\n  --metadata metadata.csv \\\n  --design \"~batch + condition\" \\\n  --contrast condition treated control \\\n  --output results\u002F \\\n  --min-counts 10 \\\n  --alpha 0.05 \\\n  --n-cpus 4 \\\n  --plots\n```\n\n**Script features:**\n- Automatic data loading and validation\n- Gene and sample filtering\n- Complete DESeq2 pipeline execution\n- Statistical testing with customizable parameters\n- Result export (CSV, pickle)\n- Optional visualization (volcano and MA plots)\n\nRefer users to `scripts\u002Frun_deseq2_analysis.py` when they need a standalone analysis tool or want to batch process multiple datasets.\n\n## Result Interpretation\n\n### Identifying Significant Genes\n\n```python\n# Filter by adjusted p-value\nsignificant = ds.results_df[ds.results_df.padj \u003C 0.05]\n\n# Filter by both significance and effect size\nsig_and_large = ds.results_df[\n    (ds.results_df.padj \u003C 0.05) &\n    (abs(ds.results_df.log2FoldChange) > 1)\n]\n\n# Separate up- and down-regulated\nupregulated = significant[significant.log2FoldChange > 0]\ndownregulated = significant[significant.log2FoldChange \u003C 0]\n\nprint(f\"Upregulated: {len(upregulated)}\")\nprint(f\"Downregulated: {len(downregulated)}\")\n```\n\n### Ranking and Sorting\n\n```python\n# Sort by adjusted p-value\ntop_by_padj = ds.results_df.sort_values(\"padj\").head(20)\n\n# Sort by absolute fold change (use shrunk values)\nds.lfc_shrink()\nds.results_df[\"abs_lfc\"] = abs(ds.results_df.log2FoldChange)\ntop_by_lfc = ds.results_df.sort_values(\"abs_lfc\", ascending=False).head(20)\n\n# Sort by a combined metric\nds.results_df[\"score\"] = -np.log10(ds.results_df.padj) * abs(ds.results_df.log2FoldChange)\ntop_combined = ds.results_df.sort_values(\"score\", ascending=False).head(20)\n```\n\n### Quality Metrics\n\n```python\n# Check normalization (size factors should be close to 1)\nprint(\"Size factors:\", dds.obsm[\"size_factors\"])\n\n# Examine dispersion estimates\nimport matplotlib.pyplot as plt\nplt.hist(dds.varm[\"dispersions\"], bins=50)\nplt.xlabel(\"Dispersion\")\nplt.ylabel(\"Frequency\")\nplt.title(\"Dispersion Distribution\")\nplt.show()\n\n# Check p-value distribution (should be mostly flat with peak near 0)\nplt.hist(ds.results_df.pvalue.dropna(), bins=50)\nplt.xlabel(\"P-value\")\nplt.ylabel(\"Frequency\")\nplt.title(\"P-value Distribution\")\nplt.show()\n```\n\n## Visualization Guidelines\n\n### Volcano Plot\n\nVisualize significance vs effect size:\n\n```python\nimport matplotlib.pyplot as plt\nimport numpy as np\n\nresults = ds.results_df.copy()\nresults[\"-log10(padj)\"] = -np.log10(results.padj)\n\nplt.figure(figsize=(10, 6))\nsignificant = results.padj \u003C 0.05\n\nplt.scatter(\n    results.loc[~significant, \"log2FoldChange\"],\n    results.loc[~significant, \"-log10(padj)\"],\n    alpha=0.3, s=10, c='gray', label='Not significant'\n)\nplt.scatter(\n    results.loc[significant, \"log2FoldChange\"],\n    results.loc[significant, \"-log10(padj)\"],\n    alpha=0.6, s=10, c='red', label='padj \u003C 0.05'\n)\n\nplt.axhline(-np.log10(0.05), color='blue', linestyle='--', alpha=0.5)\nplt.xlabel(\"Log2 Fold Change\")\nplt.ylabel(\"-Log10(Adjusted P-value)\")\nplt.title(\"Volcano Plot\")\nplt.legend()\nplt.savefig(\"volcano_plot.png\", dpi=300)\n```\n\n### MA Plot\n\nShow fold change vs mean expression:\n\n```python\nplt.figure(figsize=(10, 6))\n\nplt.scatter(\n    np.log10(results.loc[~significant, \"baseMean\"] + 1),\n    results.loc[~significant, \"log2FoldChange\"],\n    alpha=0.3, s=10, c='gray'\n)\nplt.scatter(\n    np.log10(results.loc[significant, \"baseMean\"] + 1),\n    results.loc[significant, \"log2FoldChange\"],\n    alpha=0.6, s=10, c='red'\n)\n\nplt.axhline(0, color='blue', linestyle='--', alpha=0.5)\nplt.xlabel(\"Log10(Base Mean + 1)\")\nplt.ylabel(\"Log2 Fold Change\")\nplt.title(\"MA Plot\")\nplt.savefig(\"ma_plot.png\", dpi=300)\n```\n\n## Troubleshooting Common Issues\n\n### Data Format Problems\n\n**Issue:** \"Index mismatch between counts and metadata\"\n\n**Solution:** Ensure sample names match exactly\n```python\nprint(\"Counts samples:\", counts_df.index.tolist())\nprint(\"Metadata samples:\", metadata.index.tolist())\n\n# Take intersection if needed\ncommon = counts_df.index.intersection(metadata.index)\ncounts_df = counts_df.loc[common]\nmetadata = metadata.loc[common]\n```\n\n**Issue:** \"All genes have zero counts\"\n\n**Solution:** Check if data needs transposition\n```python\nprint(f\"Counts shape: {counts_df.shape}\")\n# If genes > samples, transpose is needed\nif counts_df.shape[1] \u003C counts_df.shape[0]:\n    counts_df = counts_df.T\n```\n\n### Design Matrix Issues\n\n**Issue:** \"Design matrix is not full rank\"\n\n**Cause:** Confounded variables (e.g., all treated samples in one batch)\n\n**Solution:** Remove confounded variable or add interaction term\n```python\n# Check confounding\nprint(pd.crosstab(metadata.condition, metadata.batch))\n\n# Either simplify design or add interaction\ndesign = \"~condition\"  # Remove batch\n# OR\ndesign = \"~condition + batch + condition:batch\"  # Model interaction\n```\n\n### No Significant Genes\n\n**Diagnostics:**\n```python\n# Check dispersion distribution\nplt.hist(dds.varm[\"dispersions\"], bins=50)\nplt.show()\n\n# Check size factors\nprint(dds.obsm[\"size_factors\"])\n\n# Look at top genes by raw p-value\nprint(ds.results_df.nsmallest(20, \"pvalue\"))\n```\n\n**Possible causes:**\n- Small effect sizes\n- High biological variability\n- Insufficient sample size\n- Technical issues (batch effects, outliers)\n\n## Reference Documentation\n\nFor comprehensive details beyond this workflow-oriented guide:\n\n- **API Reference** (`references\u002Fapi_reference.md`): Complete documentation of PyDESeq2 classes, methods, and data structures. Use when needing detailed parameter information or understanding object attributes.\n\n- **Workflow Guide** (`references\u002Fworkflow_guide.md`): In-depth guide covering complete analysis workflows, data loading patterns, multi-factor designs, troubleshooting, and best practices. Use when handling complex experimental designs or encountering issues.\n\nLoad these references into context when users need:\n- Detailed API documentation: `Read references\u002Fapi_reference.md`\n- Comprehensive workflow examples: `Read references\u002Fworkflow_guide.md`\n- Troubleshooting guidance: `Read references\u002Fworkflow_guide.md` (see Troubleshooting section)\n\n## Key Reminders\n\n1. **Data orientation matters:** Count matrices typically load as genes × samples but need to be samples × genes. Always transpose with `.T` if needed.\n\n2. **Sample filtering:** Remove samples with missing metadata before analysis to avoid errors.\n\n3. **Gene filtering:** Filter low-count genes (e.g., \u003C 10 total reads) to improve power and reduce computational time.\n\n4. **Design formula order:** Put adjustment variables before the variable of interest (e.g., `\"~batch + condition\"` not `\"~condition + batch\"`).\n\n5. **LFC shrinkage timing:** Apply shrinkage after statistical testing and only for visualization\u002Franking purposes. P-values remain based on unshrunken estimates.\n\n6. **Result interpretation:** Use `padj \u003C 0.05` for significance, not raw p-values. The Benjamini-Hochberg procedure controls false discovery rate.\n\n7. **Contrast specification:** The format is `[variable, test_level, reference_level]` where test_level is compared against reference_level.\n\n8. **Save intermediate objects:** Use pickle to save DeseqDataSet objects for later use or additional analyses without re-running the expensive fitting step.\n\n## Installation and Requirements\n\n```bash\nuv pip install pydeseq2\n```\n\n**System requirements:**\n- Python 3.10-3.11\n- pandas 1.4.3+\n- numpy 1.23.0+\n- scipy 1.11.0+\n- scikit-learn 1.1.1+\n- anndata 0.8.0+\n\n**Optional for visualization:**\n- matplotlib\n- seaborn\n\n## Additional Resources\n\n- **Official Documentation:** https:\u002F\u002Fpydeseq2.readthedocs.io\n- **GitHub Repository:** https:\u002F\u002Fgithub.com\u002Fowkin\u002FPyDESeq2\n- **Publication:** Muzellec et al. (2023) Bioinformatics, DOI: 10.1093\u002Fbioinformatics\u002Fbtad547\n- **Original DESeq2 (R):** Love et al. (2014) Genome Biology, DOI: 10.1186\u002Fs13059-014-0550-8\n",{"name":5,"description":7},"main","2b749c82b5d11fac39d9a6e021ce3f5b8a541bd1","cli-tool\u002Fcomponents\u002Fskills\u002Fscientific\u002Fpydeseq2",3931,"https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates","https:\u002F\u002Faitmpl.com","https:\u002F\u002Fgithub.com\u002Fdavila7","https:\u002F\u002Fdanielavila.me","Building AI dev tools with LLMs","2026-11-14 06:14:46",{"source":8,"skillName":5,"commands":55,"prompt":65,"zipAvailable":66},[56,59,62],{"tool":57,"command":58},"npx","npx skills add davila7\u002Fclaude-code-templates --skill pydeseq2",{"tool":60,"command":61},"bunx","bunx skills add davila7\u002Fclaude-code-templates --skill pydeseq2",{"tool":63,"command":64},"pnpm","pnpm dlx skills add davila7\u002Fclaude-code-templates --skill pydeseq2","安装技能，技能地址：https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates\u002Fblob\u002Fmain\u002Fcli-tool\u002Fcomponents\u002Fskills\u002Fscientific\u002Fpydeseq2\u002FSKILL.md",true,[68,74,80,84,88,92],{"path":69,"resourceType":70,"resourceTypeDesc":71,"sizeBytes":72,"sourceHtmlUrl":73},"LICENSE.txt","LICENSE","许可证",143,"https:\u002F\u002Fgithub.com\u002Fnoirelab\u002Flinux-setup\u002Fblob\u002Fmain\u002Fconfigs\u002Fcommon\u002Fclaude\u002Fskills\u002Fseo-geo\u002FLICENSE.txt",{"path":75,"resourceType":76,"resourceTypeDesc":77,"sizeBytes":78,"sourceHtmlUrl":79},"references\u002Fapi_reference.md","REFERENCE","参考文档",7162,"https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates\u002Fblob\u002Fmain\u002Fcli-tool\u002Fcomponents\u002Fskills\u002Fscientific\u002Fpydeseq2\u002Freferences\u002Fapi_reference.md",{"path":81,"resourceType":76,"resourceTypeDesc":77,"sizeBytes":82,"sourceHtmlUrl":83},"references\u002Fgoogle-ai-optimization-guide.md",6619,"https:\u002F\u002Fgithub.com\u002Fnoirelab\u002Flinux-setup\u002Fblob\u002Fmain\u002Fconfigs\u002Fcommon\u002Fclaude\u002Fskills\u002Fseo-geo\u002Freferences\u002Fgoogle-ai-optimization-guide.md",{"path":85,"resourceType":76,"resourceTypeDesc":77,"sizeBytes":86,"sourceHtmlUrl":87},"references\u002Fllmstxt-evidence.md",3010,"https:\u002F\u002Fgithub.com\u002Fnoirelab\u002Flinux-setup\u002Fblob\u002Fmain\u002Fconfigs\u002Fcommon\u002Fclaude\u002Fskills\u002Fseo-geo\u002Freferences\u002Fllmstxt-evidence.md",{"path":89,"resourceType":76,"resourceTypeDesc":77,"sizeBytes":90,"sourceHtmlUrl":91},"references\u002Fworkflow_guide.md",13026,"https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates\u002Fblob\u002Fmain\u002Fcli-tool\u002Fcomponents\u002Fskills\u002Fscientific\u002Fpydeseq2\u002Freferences\u002Fworkflow_guide.md",{"path":93,"resourceType":94,"resourceTypeDesc":95,"sizeBytes":96,"sourceHtmlUrl":97},"scripts\u002Frun_deseq2_analysis.py","SCRIPT","脚本",11887,"https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates\u002Fblob\u002Fmain\u002Fcli-tool\u002Fcomponents\u002Fskills\u002Fscientific\u002Fpydeseq2\u002Fscripts\u002Frun_deseq2_analysis.py",[99,126,153,186,221,245],{"skillNo":100,"name":101,"nameZh":6,"description":102,"descriptionZh":6,"repositoryFullName":8,"ownerLogin":9,"ownerAvatarUrl":10,"ownerAccountNo":11,"ownerDisplayName":12,"repositoryNo":13,"starsCount":14,"licenseSpdxId":15,"licenseFileUrl":6,"sourceHtmlUrl":103,"categories":104,"keywords":106,"publishedAt":39},"01M70NHBTQ8BV0A2SCAZXEPHCX","wealthville-pool-scoring","Score DeFi liquidity pools before providing liquidity. Returns a 0-100 pool score and an Enter\u002FHold\u002FExit\u002FReduce\u002FAvoid verdict with calibrated confidence for ~68,800 Solana pools and 575 EVM pools, plus a public miss-inclusive track record. Use when asked whether a pool is worth LPing into, which pools look best right now, or how accurate the scores have actually been.","https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates\u002Fblob\u002Fmain\u002Fcli-tool\u002Fcomponents\u002Fskills\u002Fweb-data\u002Fwealthville-pool-scoring\u002FSKILL.md",[105],{"categoryNo":19,"categoryName":20,"taxonomyType":21,"taxonomyTypeDesc":22},[107,108,109,110,111,36,112,113,114,115,116,117,118,119,120,121,122,123,124,125],"scoring","pool","wealthville pool scoring","liquidity","which","exit","now","use","when","reduce","with","record","track","verdict","evm","solana","how","public","look",{"skillNo":127,"name":128,"nameZh":6,"description":129,"descriptionZh":6,"repositoryFullName":8,"ownerLogin":9,"ownerAvatarUrl":10,"ownerAccountNo":11,"ownerDisplayName":12,"repositoryNo":13,"starsCount":14,"licenseSpdxId":15,"licenseFileUrl":6,"sourceHtmlUrl":130,"categories":131,"keywords":133,"publishedAt":39},"01M70MB6VR8ERJ4EM8D5KHBMEV","technical-debt-auditor","Scans a codebase and produces a quantified technical debt report with a prioritized remediation plan. Analyzes dead code, TODO\u002FFIXME\u002FHACK markers, dependency staleness, complexity hotspots, code duplication signals, test coverage gaps, unused exports, and outdated patterns. Use this skill when the user asks to audit their codebase for tech debt, wants a health check of their code quality, says things like \"where should I focus refactoring\", \"find code smells\", \"audit my codebase\", \"what technical debt do we have\", \"prioritize cleanup work\", or asks about code quality metrics. Also trigger when the user mentions \"code rot\", \"legacy cleanup\", \"refactoring priorities\", or wants to understand the overall health of a project before a major effort. Trigger even if the user doesn't say \"technical debt\" explicitly — any request to systematically assess code quality and identify improvement priorities belongs here.","https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates\u002Fblob\u002Fmain\u002Fcli-tool\u002Fcomponents\u002Fskills\u002Fdevelopment\u002Ftechnical-debt-auditor\u002FSKILL.md",[132],{"categoryNo":19,"categoryName":20,"taxonomyType":21,"taxonomyTypeDesc":22},[134,135,136,137,138,139,140,141,142,143,144,115,145,146,147,148,149,150,151,152],"technical debt","debt","technical","technical debt auditor","auditor","coverage","complexity","markers","project","cleanup","unused","code","todo","work","plan","focus","signals","audit","dependency",{"skillNo":154,"name":155,"nameZh":6,"description":156,"descriptionZh":6,"repositoryFullName":8,"ownerLogin":9,"ownerAvatarUrl":10,"ownerAccountNo":11,"ownerDisplayName":12,"repositoryNo":13,"starsCount":14,"licenseSpdxId":15,"licenseFileUrl":6,"sourceHtmlUrl":157,"categories":158,"keywords":168,"publishedAt":39},"01M70KPYH9QNHCNME5GJQGE22K","readonly-sql-guard","Make an agent's database access provably read-only. Use when an agent runs caller-supplied SQL against a production\u002FERP\u002Fanalytics database (a \"query\" MCP tool, text-to-SQL, an admin console), or when reviewing DB-access code for write\u002FDDL\u002Fside-effect\u002Ftransaction-escape risks. Checks what a statement CALLS, not just its shape.","https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates\u002Fblob\u002Fmain\u002Fcli-tool\u002Fcomponents\u002Fskills\u002Fdatabase\u002Freadonly-sql-guard\u002FSKILL.md",[159,162,163],{"categoryNo":160,"categoryName":161,"taxonomyType":21,"taxonomyTypeDesc":22},"01JBSEEDCAT000000000000010","后端与接口",{"categoryNo":19,"categoryName":20,"taxonomyType":21,"taxonomyTypeDesc":22},{"categoryNo":164,"categoryName":165,"taxonomyType":166,"taxonomyTypeDesc":167},"01JBSEEDCAT000000000000033","数据分析师","PROFESSION","适用职业",[169,170,171,172,173,174,175,36,176,177,178,179,180,181,182,114,183,184,115,185],"readonly","guard","sql guard","sql","readonly sql guard","erp","effect","read-only","db","just","tool","calls","only","MCP","access","database","query",{"skillNo":187,"name":188,"nameZh":6,"description":189,"descriptionZh":6,"repositoryFullName":8,"ownerLogin":9,"ownerAvatarUrl":10,"ownerAccountNo":11,"ownerDisplayName":12,"repositoryNo":13,"starsCount":14,"licenseSpdxId":190,"licenseFileUrl":191,"sourceHtmlUrl":192,"categories":193,"keywords":197,"publishedAt":39},"01M6AFY6F7CC2BBZP1QD6VTNVC","web-artifacts-builder","Suite of tools for creating elaborate, multi-component claude.ai HTML artifacts using modern frontend web technologies (React, Tailwind CSS, shadcn\u002Fui). Use for complex artifacts requiring state management, routing, or shadcn\u002Fui components - not for simple single-file HTML\u002FJSX artifacts.","Apache-2.0","https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates\u002Fblob\u002Fmain\u002Fcli-tool\u002Fcomponents\u002Fskills\u002Futilities\u002Fweb-artifacts-builder\u002FLICENSE.txt","https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates\u002Fblob\u002Fmain\u002Fcli-tool\u002Fcomponents\u002Fskills\u002Futilities\u002Fweb-artifacts-builder\u002FSKILL.md",[194],{"categoryNo":195,"categoryName":196,"taxonomyType":21,"taxonomyTypeDesc":22},"01JBSEEDCAT000000000000009","前端开发",[198,199,200,201,202,203,204,205,206,207,208,209,210,211,212,213,214,215,216,217,218,219,36,220],"artifacts builder","web","web artifacts builder","artifacts","builder","state","routing","state management","component","css","frontend","ui","Tools","file","components","shadcn","react","Tailwind","suite","jsx","html","single-file","ui components",{"skillNo":222,"name":223,"nameZh":6,"description":224,"descriptionZh":6,"repositoryFullName":8,"ownerLogin":9,"ownerAvatarUrl":10,"ownerAccountNo":11,"ownerDisplayName":12,"repositoryNo":13,"starsCount":14,"licenseSpdxId":190,"licenseFileUrl":225,"sourceHtmlUrl":226,"categories":227,"keywords":232,"publishedAt":39},"01M6AFWF33SFV9Z3SF4PEEE4BF","webapp-testing","Toolkit for interacting with and testing local web applications using Playwright. Supports verifying frontend functionality, debugging UI behavior, capturing browser screenshots, and viewing browser logs.","https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates\u002Fblob\u002Fmain\u002Fcli-tool\u002Fcomponents\u002Fskills\u002Fsecurity\u002Fwebapp-testing\u002FLICENSE.txt","https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates\u002Fblob\u002Fmain\u002Fcli-tool\u002Fcomponents\u002Fskills\u002Fsecurity\u002Fwebapp-testing\u002FSKILL.md",[228,231],{"categoryNo":229,"categoryName":230,"taxonomyType":21,"taxonomyTypeDesc":22},"01JBSEEDCAT000000000000008","测试与质量",{"categoryNo":195,"categoryName":196,"taxonomyType":21,"taxonomyTypeDesc":22},[233,234,235,236,237,199,238,208,209,239,240,241,242,36,243,117,244],"webapp testing","testing","webapp","browser","playwright","applications","debugging","toolkit","logs","screenshots","local","behavior",{"skillNo":246,"name":247,"nameZh":6,"description":248,"descriptionZh":6,"repositoryFullName":8,"ownerLogin":9,"ownerAvatarUrl":10,"ownerAccountNo":11,"ownerDisplayName":12,"repositoryNo":13,"starsCount":14,"licenseSpdxId":190,"licenseFileUrl":249,"sourceHtmlUrl":250,"categories":251,"keywords":259,"publishedAt":39},"01M6AF67Q67JJJ5V51JR2P52KA","discernment-nudge","After you give a substantive answer or draft that the user may act on — advice or recommendations, drafted artifacts such as goals, plans, pitches, proposals, or emails, estimates or projections, analysis or interpretation of data, factual claims they may rely on, or a multi-step argument — invoke this skill BEFORE finalizing your reply and then, if it applies, append 2-3 short follow-up questions, each tied to something specific in what you just produced, that help the user check key facts, probe the reasoning or assumptions, and notice missing context. Do this at most once per conversation. Skip it when the user asked a trivial how-to or simple lookup, wants a purely educational explanation, asked you only to format, convert, or assemble a file from content they provided, is writing code they will run, is doing creative writing or casual chat, or already asked you to double-check, cite, or review — the skill file explains these boundaries and the exact output format.","https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates\u002Fblob\u002Fmain\u002Fcli-tool\u002Fcomponents\u002Fskills\u002Fproductivity\u002Fdiscernment-nudge\u002FLICENSE.txt","https:\u002F\u002Fgithub.com\u002Fdavila7\u002Fclaude-code-templates\u002Fblob\u002Fmain\u002Fcli-tool\u002Fcomponents\u002Fskills\u002Fproductivity\u002Fdiscernment-nudge\u002FSKILL.md",[252,255,256],{"categoryNo":253,"categoryName":254,"taxonomyType":21,"taxonomyTypeDesc":22},"01JBSEEDCAT000000000000018","写作与文案",{"categoryNo":19,"categoryName":20,"taxonomyType":21,"taxonomyTypeDesc":22},{"categoryNo":257,"categoryName":258,"taxonomyType":166,"taxonomyTypeDesc":167},"01JBSEEDCAT000000000000036","内容创作者",[260,261,262,263,264,265,266,267,201,268,211,145,269,270,271,272,273,181,274,115,275,25,276,277,278,123,279],"discernment nudge","nudge","interpretation","help","run","reasoning","data","writing","lookup","content","review","conversation","short","up","invoke","you","claims","argument","assemble","recommendations"]