跳到正文
cccskills
免费在 GitHub 公开

bioconductor-pcamethods

Provides Bayesian PCA, Probabilistic PCA, Nipals PCA, Inverse Non-Linear PCA and the conventional SVD PCA. A cluster based method for missing value estimation is included for comparison. BPCA, PPCA and NipalsPCA may be used to perform PCA o

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

实测评分 未实测

查看安装方法

包含的文件(1)

  • SKILL.md5.8 KB

SKILL.md(原文)

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

pcaMethods

Dependencies & Environment

Package-intrinsic requirements from the Bioconductor landing page — reproduce in any R environment.

  • Version: 2.4.0 · Bioconductor: 3.23 · R: ≥ 4.6
  • Depends: Biobase
  • Imports: BiocGenerics, Rcpp, MASS
  • System requirements: Rcpp
  • Install: BiocManager::install("pcaMethods")

When to Use

  • Handling Outliers: When performing PCA on metabolite or microarray data corrupted with extreme values using a robust singular value decomposition (method="robustPca").
  • Missing Value Imputation: When estimating missing values (NA) in incomplete datasets using Probabilistic PCA (method="ppca") or Nipals (method="nipals").
  • Algorithm Comparison: When comparing multiple PCA formulations (e.g., "svd", "ppca", "bpca", "svdImpute", "nipals", "nlpca") on the same ExpressionSet or matrix using a unified pca() interface.

When NOT to Use

  • Sparse Single-Cell Data: For single-cell RNA-seq data with high dropout rates, use scry or GLM-PCA instead because they model the count distribution directly.
  • Complete Datasets: For standard, complete matrices without missing values or outliers, use base R's prcomp or irlba because they are faster and computationally lighter.

Data Requirements

  • Input Format: A numeric matrix, data.frame, or an ExpressionSet (from the Biobase package).
  • Missing Values: Missing data should be represented as NA.
  • Normalization: Data is typically mean-centered and scaled (e.g., using the prep() function with scale="none" and center=TRUE) prior to running PCA.

Key Parameters

  • method ("svd"): The PCA algorithm to use (e.g., "svd", "robustPca", "ppca", "bpca", "svdImpute", "nipals", "nlpca").
  • nPcs (2): Number of principal components to calculate.
  • center (TRUE): Logical indicating whether to mean-center the variables.
  • scale ("none"): Scaling method to apply in the prep() function.
  • maxSteps (1000): Maximum number of iterations for iterative algorithms like NLPCA.
  • fold (10): Number of folds for cross-validation when using the Q2() function.
  • evalPcs (1:5): The principal components to evaluate when estimating error with kEstimate().
  • em ("nrmsep"): The error metric to use during cross-validation in kEstimate().

Best Practices

  • Pre-processing: Always pre-process data using the prep() function to handle centering and scaling before running pca().
  • Determine Components: Use cross-validation via the Q2() function to evaluate the optimal number of principal components to retain.
  • Outlier Management: When dealing with extreme outliers, set them to NA and use method="ppca", or use method="robustPca" directly on the complete data.
  • Error Estimation: Estimate imputation error using kEstimate() with cross-validation to ensure the chosen method and number of PCs are appropriate.

Common Pitfalls

  • Failing to Center Before Outlier Injection: Creating artificial outliers for testing on un-centered data, which shifts original means and prevents objective comparison. Fix: Use scale(..., center=TRUE, scale=FALSE) before injecting outliers.
  • PPCA Instability: PPCA failing to converge or misestimating values in rare cases. Fix: Use a while loop checking if sum(abs(Q2)) > 1 to re-run the Q2() function if PPCA becomes unstable.
  • Non-Numeric Input: Passing a matrix containing character columns or factors to pca(), causing the algorithms to fail. Fix: Ensure the input matrix or ExpressionSet contains only numeric values.

Alternatives

  • missMDA: For PCA and imputation on incomplete datasets using regularized iterative PCA.
  • impute: For K-Nearest Neighbors imputation (impute.knn) of missing microarray data.
  • FactoMineR: For advanced exploratory multivariate data analysis on complete datasets.
  • irlba: For fast, truncated SVD and PCA on extremely large, sparse matrices.

Citations

  • Liu, L., Hawkins, D.M., Ghosh, S. and Young, S.S. Robust singular value decomposition analysis of microarray data. PNAS, 2003;100:13167–13172.
  • Hawkins, D.M., Liu, L. and Young, S.S. Robust Singular Value Decomposition. National Institute of Statistical Sciences, 2001, Tech Report 122.

References

<!-- biomate-cta -->

Run this on BioMate

This skill is the knowledge layer — when, why, and how to use pcamethods. To run this analysis on your own data with managed compute, automated QC, and reproducible outputs, use BioMate — free to start.

▶ Open pcamethods on BioMate →

评价

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

评论

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

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

提交后会立即公开。请勿填写邮箱、电话等个人信息。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 的全部技能

报告此技能的问题