微生物研究R语言做图

R语言|PLS-DA分析绘图示例

2021-06-15  本文已影响0人  维凡生物

偏最小二乘法判别分析(PLS_DA)是一种用于判别分析的多变量统计分析方法,一种根据观察或测量到的若干变量值,来判断研究对象如何分类的常用统计分析方法。对不同处理样本(如观测样本、对照样本)的特性分别进行训练,产生训练集,并检验训练集的可信度。以下是PLS_DA分析绘图的一个小示例。

(本文仅作绘图示例,不涉及分析说明)

install.packages("BiocManager") 
## install mixOmics 
BiocManager::install('mixOmics')
library(mixOmics)

data1=read.table("test.txt", header=T, row.names=1,sep="\t")
group=read.table("group.txt",sep="\t",header=T,row.names=1)
X=t(data1[, rownames(group)])
Y=group$group
##PLS-DA分析,选取2个主成分进行分析
plsda.datatm <-plsda(X, Y, ncomp = 3)
plotIndiv(plsda.datatm, ind.names = FALSE, legend=TRUE,ellipse = TRUE, title="sPLS-DA - final result")
background <- background.predict(plsda.datatm, comp.predicted=2,dist = "max.dist") 
image.png
#plotVar(plsda.datatm) 
plotIndiv(plsda.datatm, comp = 1:2, 
          ind.names = FALSE, title = "Maximum distance",
          legend = TRUE,  background = background,ellipse = TRUE)
#auc.plsda <- auroc(plsda.datatm)
image.png

软文同步于“作图帮”公众号,感兴趣的小伙伴可以关注一下哦~

上一篇下一篇

猜你喜欢

热点阅读