跳到正文
cccskills
免费在 GitHub 公开

bioconductor-phipdata

PhIPData defines an S4 class for phage-immunoprecipitation sequencing (PhIP-seq) experiments. Buliding upon the RangedSummarizedExperiment class, PhIPData enables users to coordinate metadata with experimental data in analyses. Additionally, PhIPData provides specialized methods to subset and identify beads-only samples, subset objects using virus aliases, and use existing peptide libraries to populate object parameters.

中文简介正在准备中。当前显示原文说明。

查看安装方法

包含的文件(1)

  • SKILL.md5.3 KB

SKILL.md(原文)

安装前可以查看交给 AI 代理的指令内容。

PhIPData

Workflows

Standard Workflow

PhIPData defines an S4 class for phage-immunoprecipitation sequencing (PhIP-seq) experiments. Buliding upon the RangedSummarizedExperiment class, PhIPData enables users to coordinate metadata with experimental data in analyses. Additionally, PhIPData provides specialized methods to subset and identify beads-only samples, subset objects using virus aliases, and use existing peptide libraries to populate object parameters.

library(PhIPData)

# 1. Define dimensions and simulate data matrices
n_samples <- 5L
n_peps <- 10L
counts_dat <- matrix(sample(1:1e6, n_samples*n_peps, replace = TRUE), nrow = n_peps)
logfc_dat <- matrix(rnorm(n_samples*n_peps, mean = 0, sd = 10), nrow = n_peps)
prob_dat <- matrix(rbeta(n_samples*n_peps, shape1 = 1, shape2 = 1), nrow = n_peps)

# 2. Define metadata
peptide_meta <- data.frame(pep_id = paste0("pep_", 1:n_peps))
sample_meta <- data.frame(
  sample_name = paste0("sample", 1:n_samples), 
  group = c("ctrl", "beads", "trt", "trt", "trt")
)

# 3. Construct PhIPData object and subset beads-only samples
phip_obj <- PhIPData(counts_dat, logfc_dat, prob_dat, peptide_meta, sample_meta)
beads_only <- subsetBeads(phip_obj)

Note: Inputs are parallel matrices of counts, log fold-changes, probabilities, and metadata data frames; output is a coordinated PhIPData object.

When to Use

  • Storing PhIP-seq experimental results: Manage raw read counts, log2 estimated fold-changes, and enrichment probabilities in a single coordinated container.
  • Identifying and subsetting control samples: Quickly isolate beads-only control samples using the specialized subsetBeads() wrapper.
  • Managing peptide libraries and virus aliases: Query and subset specific viral species using setAlias(), getAlias(), and standard subset() operations.

When NOT to Use

  • For general RNA-seq differential expression: Use DESeq2 or edgeR because PhIPData is specifically tailored for PhIP-seq antibody enrichment assays.
  • For single-cell transcriptomics: Use SingleCellExperiment because PhIPData is designed for bulk phage-display sequencing experiments.
  • For mass spectrometry-based proteomics: Use QFeatures or MSnbase because PhIPData is designed for high-throughput sequencing of phage libraries.

Data Requirements

  • Assay Matrices: Three parallel matrices of identical dimensions: counts (non-negative numeric values), logfc (log2 estimated fold-changes), and prob (enrichment probabilities/p-values).
  • Peptide Metadata: A DataFrame or data.frame containing peptide identifiers, and optionally start/end positions (pos_start, pos_end).
  • Sample Metadata: A DataFrame or data.frame containing sample names and group designations (e.g., "beads", "ctrl", "trt").

Key Parameters

  • counts (matrix): Raw read counts or pseudocounts for each peptide and sample.
  • logfc (matrix): Log2 estimated fold-changes compared to negative controls.
  • prob (matrix): Probabilities or p-values associated with enriched antibody response.
  • withDimnames (TRUE): In librarySize() and propReads(), controls whether output vectors/matrices retain row and column names.

Best Practices

  • Ensure that the sample metadata contains a column named group with a "beads" level to enable automatic identification of beads-only control samples via subsetBeads().
  • Use makeLibrary() to save peptide annotations as reusable templates, avoiding redundant re-importing of large peptide metadata files.
  • Coerce PhIPData objects to DGEList using as(phip_obj, "DGEList") when transitioning to edgeR for differential enrichment analysis.

Common Pitfalls

  • Missing peptide start/end positions: If pos_start and pos_end are missing from the peptide metadata, the constructor replaces them with 0 and prints a warning.
  • Missing required assays: Attempting to construct a PhIPData object without counts, logfc, or prob will result in empty matrices being initialized for the missing assays.
  • Case-sensitivity in aliases: Keys in the alias database (e.g., "hiv" vs "HIV") are case-sensitive; ensure consistent casing when calling getAlias().

Alternatives

  • SummarizedExperiment: For general-purpose coordinated representation of rectangular biological assay data.
  • RangedSummarizedExperiment: For genomic-coordinate-linked assay data.
  • edgeR (DGEList): For standard differential count analysis.

Citations

  • Chen A, Scharpf R, Ruczinski I (2026). PhIPData: A Container for PhIP-Seq Experiments. R package version 1.12.0.

References

评价

还没有评价。欢迎分享你的使用感受。

评论

无需登录即可发表。欢迎分享使用感受、提问,或说明在哪些环境下无法运行。

还没有评论。来写第一条吧。

提交后会立即公开。请勿填写邮箱、电话等个人信息。URL 最多 2 个。

0 / 1000 字

此表单受 reCAPTCHA 保护,适用 Google 的隐私权政策和服务条款。

同一仓库的技能

The package contains functions for exploratory oligonucleotide array analysis. The dependence on tkWidgets only concerns few convenience functions. 'affy' is fully functional without it.

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

This tool uses the label-tree function from HyPhy to annotate a phylogenetic tree. It allows users to select a subset of leaves using either a regular expression or a list of sequence names, and then apply a specified label to these selecte

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

This package provides a client for the Bioconductor AnnotationHub web resource. The AnnotationHub web resource provides a central location where genomic files (e.g., VCF, bed, wig) and other resources from standard locations (e.g., UCSC, En

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

apeglm provides Bayesian shrinkage estimators for effect sizes for a variety of GLM models, using approximation of the posterior for individual coefficients.

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

ASURAT is a software for single-cell data analysis. Using ASURAT, one can simultaneously perform unsupervised clustering and biological interpretation in terms of cell type, disease, biological process, and signaling pathway activity. Inputting a single-cell RNA-seq data and knowledge-based databases, such as Cell Ontology, Gene Ontology, KEGG, etc., ASURAT transforms gene expression tables into original multivariate tables, termed sign-by-sample matrices (SSMs).

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

AUCell allows to identify cells with active gene sets (e.g. signatures, gene modules...) in single-cell RNA-seq data. AUCell uses the "Area Under the Curve" (AUC) to calculate whether a critical subset of the input gene set is enriched with

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

查看 bioMate-AI 的全部技能

报告此技能的问题