Scanpy - Phân tích tế bào đơn
Bộ kỹ năng chạy toàn bộ quy trình scRNA-seq bằng Scanpy: từ QC, chuẩn hóa, giảm chiều, phân cụm đến tìm marker và chú giải loại tế bào.
Dữ liệu · Phân tíchNâng cao★ 33,030⑂ 3,248Điểm AI 8/10Cập nhật lần cuối: 9 thg 8, 2026
Làm được gì
Hướng dẫn Claude thực hiện trọn vẹn quy trình phân tích RNA-seq tế bào đơn (scRNA-seq) dựa trên Scanpy.
- Kiểm soát chất lượng: lọc tế bào/gen, kiểm tra tỉ lệ ty thể, phát hiện doublet bằng Scrublet
- Tiền xử lý: chuẩn hóa, log1p, chọn gen biến thiên cao (HVG), scale/regress
- Giảm chiều & phân cụm: PCA → đồ thị lân cận → UMAP/t-SNE, phân cụm Leiden ở nhiều độ phân giải
- Hiệu chỉnh batch: harmony / bbknn / combat
- Marker & chú giải:
rank_genes_groups, xuất CSV theo nhóm, ánh xạ cụm → loại tế bào, chấm điểm chữ ký gen - Phân tích tiếp theo: gộp pseudobulk cho pydeseq2, trích tập con, vẽ hình chất lượng xuất bản
- Tương tác với R: chuyển đổi tệp Seurat/SingleCellExperiment
.rdssang.h5ad
Các script trong scripts/ đều nhận và xuất .h5ad nên có thể nối chuỗi từng bước, hoặc chạy run_pipeline.py một lệnh duy nhất.
Phù hợp với ai
- Nhà nghiên cứu, học viên cao học ngành tin sinh học làm việc thường xuyên với dữ liệu tế bào đơn
- Người cần khám phá nhanh dữ liệu 10x Genomics hoặc tệp h5ad
- Người đang chuyển từ Seurat (R) sang Scanpy (Python)
- Ai cần pipeline tái lập được (config JSON, lưu điểm kiểm tra trung gian)
Ví dụ sử dụng
- Chạy một lệnh: "Xử lý raw.h5ad từ QC đến phân cụm Leiden và tìm marker" →
python scripts/run_pipeline.py raw.h5ad -o processed.h5ad --resolution 0.5 --scrublet - Tích hợp nhiều mẫu: dữ liệu có hiệu ứng batch theo mẫu, dùng
--batch-key sample --batch-method harmonyđể có UMAP hợp nhất - Chuyển từ Seurat: nhận
pbmc.rdstừ đồng nghiệp, làm theoreferences/r_interop.mdđể chuyển sang.h5adrồi phân tích lại trong Scanpy - Hình cho bài báo:
python scripts/plot.py ann.h5ad --kind dotplot --genes CD3D CD14 --groupby cell_type
· · · Hướng dẫn cài đặt · · ·
Cài vào ứng dụng Claude (không cần terminal)
- Tải tệp ZIP bằng nút bên dưới.
- Trong Claude, mở Settings → Capabilities và bật 'Code execution and file creation'. (chỉ một lần)
- Vào Customize → Skills → + → 'Upload a skill' rồi tải tệp ZIP lên.
Cài vào Claude Code
Để Claude làm — dán câu dưới đây vào Claude Code
Hãy cài skill mình tìm thấy trên Claude Skill Mart. Copy thư mục skills/scanpy từ repo GitHub K-Dense-AI/scientific-agent-skills vào ~/.claude/skills/scanpy/ của mình. Sau khi cài xong, cho mình biết skill này làm được gì trong một câu.
Cài bằng lệnh thủ công
git clone https://github.com/K-Dense-AI/scientific-agent-skills.git && mkdir -p ~/.claude/skills && cp -r scientific-agent-skills/skills/scanpy ~/.claude/skills/⚠ Đây là skill do người khác tạo. Hãy kiểm tra repo gốc trước khi cài.
- Kiểm tra môi trường: chạy
python --version, cần Python 3.12 trở lên (scanpy 1.12 không hỗ trợ ≤3.11). - Cài Scanpy:
uv pip install "scanpy[leiden]"hoặcpip install "scanpy[leiden]"— phần[leiden]cài kèmigraphvàleidenalgcần cho phân cụm. - Tải kho mã:
git clone https://github.com/K-Dense-AI/scientific-agent-skills.git - Sao chép kỹ năng:
mkdir -p ~/.claude/skills && cp -r scientific-agent-skills/skills/scanpy ~/.claude/skills/ - Kiểm tra cấu trúc: trong
~/.claude/skills/scanpy/phải cóSKILL.md,scripts/,references/,assets/; đặc biệt giữ nguyênscripts/_common.py. - Khởi động lại Claude Code, mở phiên mới và yêu cầu ví dụ "dùng scanpy để QC tệp h5ad này".
- Chạy thử:
python ~/.claude/skills/scanpy/scripts/inspect_data.py du_lieu.h5adđể xem bản tóm tắt dữ liệu. - (Tùy chọn) Nếu dùng tệp R: cài R cùng zellkonverter/SeuratDisk theo hướng dẫn trong
references/r_interop.md.
Xem mã nguồn trên GitHub ↗Giấy phép: MIT