PCA

比较 stats::procmp() 和 vegan::rda() 两种方法

2018-11-02 Gao
上次提及 PCA 分析的方法有很多种。那不同方法之间的得到的结果会有差异吗? 最近采用 PCA 分析 RNA-seq 样本之间的差异,得到了下面的结果。 […] 生成一个含有 1000 个基因, 27 个样品的数据集. 这 27 个样品来自于 3 个基因型(WT, Mutant1, Mutant2), 3 种处理(CK, Trt1, Trt2), 共分为 \(3 * 3 = 9\) … 阅读全文 →

常见排序分析方法及R语言实现

2018-10-22 gaoch
相信大家在做微生物多样性研究时经常听到PCA分析、PCoA分析,NMDS分析,CCA分析,RDA分析。它们对物种(或基因、功能)的分析具有重要作用,因而频频出现在16S测序及宏基因组测序中。以上分析本质上都属于排序分析(Ordination analysis)。 排序分析(ordination analysis),最早是生态学(ecology)中研究群落(communities)的一大类多元分析手 … 阅读全文 →

排序分析结果作图

2018-10-22 gaoch
𝕝ibrary(vegan) ## Loading required package: permute ## Loading required package: lattice ## This is vegan 2.6-4 data("varespec") pca <- rda(varespec) 首先看一下结果: summary(pca) ## ## Call: ## … 阅读全文 →