r语言plsda分析代码附数据

r语言plsda分析代码附数据
木羽长弓.R
Administrator
Sun Feb 12 17:35:58 2017
# 每一列代表一位患者的多个数据，不同颜色代表了不同的分组（当然目前我们进行整理，如果您那边需要我把一个组的患者放在一起的话
# ，我这边再整理一下），我是想通过PLS-DA挖掘下不同组别患者间存在差异的指标有
哪些
library(mixOmics)
## Loading required package: MASS
## Loading required package: lattice
## Loading required package: ggplot2
##
## Loaded mixOmics 6.0.0
##
## Visit /doc/212628404.html, for more details about our methods.
## Any bug reports or comments? Notify us at mixomics at math.univ-toul ouse.fr or https:///doc/212628404.html,/klecao/pac kage-mixomics/issues
##
## Thank you for using mixOmics!
# ## First example
# data(breast.tumors)
# X <- breast.tumors$gene.exp
# Y <- breast.tumors$sample$treatment
#
# plsda.breast <- plsda(X, Y, ncomp = 2)
# palette(c("red", "blue"))
# col.breast <- as.numeric(as.factor(Y))
# plotIndiv(plsda.breast, /doc/212628404.html,s = TRUE, col = col.breast ,ellipse =
TRUE)
# legend('bottomleft', c("After", "Before"), pch = c(16, 16),
# col = unique(col.breast), cex = 1, pt.cex = c(1.2, 1.2),
# title = "Treatment")
# palette("default")
setwd("E:/service/2017/2 11 木羽长弓")
#读取数据
data=read.csv("数据1.csv")
head(data)
## X X101 X103 X106 X107 X124 X132 X141 X142 X144 X151 X157 X158 X1 61
## 1 101 5 3 4 3 4 5 2 8 4 3 4 3 3
## 2 102 0 0 0 0 0 0 0 1 0 0 0 0 0
## 3 103 0 0 0 0 0 0 0 0 0 0 0 0 0
## 4 104 12 12 13 9 17 16 10 19 13 14 13 10 18
## 5 105 0 0 0 0 0 0 0 1 0 2 0 0 0
## 6 106 3 3 1 2 2 2 1 3 2 2 2 3 2
## X164 X166 X168 X172 X173 X175 X180 X730 X731 X732 X741 X748 X753 X 758
## 1 4 7 5 2 2 4 5 4 4 4 2 4 4 4
## 2 0 0 0 0 0 0 0 0 2 0 0 0 0 0
## 3 0 1 0 0 0 0 0 0 0 0 0 0 0 0
## 4 14 16 19 11 12 17 21 22 18 11 17 12 20 12
## 6 2 4 3 1 2 3 3 3 1 2 2 1 1 1
## X778 X788 X.1 X104 X105 X109 X129 X130 X143 X155 X162 X170 X171 X1 74
## 1 3 4 NA 6 4 4 2 3 4 3 4 2 5 2
## 2 0 0 NA 0 0 0 0 0 0 0 0 0 0 0
## 3 0 0 NA 2 2 3 3 2 4 2 1 2 3 2
## 4 15 12 NA 24 22 21 15 19 24 24 26 20 23 14
## 5 0 0 NA 0 0 0 0 1 0 0 0 0 0 0
## 6 3 2 NA 5 2 10 13 10 13 6 6 8 11 9
## X178 X733 X734 X743 X749 X760 X767 X771 X.2 X102 X110 X111 X119 X1 31
## 1 4 4 5 5 4 4 5 2 NA 4 5 3 4
## 2 0 0 0 0 0 0 0 0 NA 3 3 2 2 3
## 3 3 4 3 2 3 1 3 3 NA 1 1 2 3 2
## 4 20 19 20 22 23 18 22 16 NA 4 4 2 2 2
## 5 1 0 0 1 0 0 0 0 NA 2 1 2 2 2
## 6 6 11 8 7 11 8 10 10 NA 2 3 2 2 2
## X135 X145 X152 X159 X160 X167 X169 X179 X740 X744 X746 X751 X756 X 757
## 1 3 2 4 2 4 2 2 3 1 2 2 3 2 2
## 2 3 2 5 2 5 3 2 3 2 2 3 3 3 2
## 3 2 6 3 2 1 1 3 2 4 2 2 3 4 3
## 4 2 0 1 4 1 3 3 1 0 1 2 2 0 3
## 5 1 1 2 1 2 1 2 1 1 1 1 1 2 2
## 6 3 2 3 3 2 1 2 2 1 1 2 2 2 2
## X764 X768 X770 X772 X777 X780 X785 X794 X.3 X108 X134 X154 X747 X7 63
## 1 3 4 4 2 2 3 3 3 NA 0 1 2 3 2
## 2 4 4 5 3 2 2 3 3 NA 0 0 0 0 0
## 4 3 3 2 0 2 2 2 2 NA 2 3 4 2 3
## 5 1 1 1 2 1 2 2 1 NA 7 6 8 8 7
## 6 2 2 4 2 1 1 2 1 NA 4 4 4 2 1
## X773 X776 X787 X.4 X736 X742 X745 X752 X762 X774 X775 X786 X789 X.
5 X738
## 1 3 2 2 NA 5 6 6 4 6 6 3 4 4 N A 1
## 2 0 0 0 NA 3 5 3 4 5 3 2 5 4 N A 2
## 3 4 4 5 NA 2 1 2 2 2 1 1 1 1 N A 4
## 4 3 3 4 NA 7 7 5 6 7 5 7 8 7 N A 6
## 5 6 4 5 NA 7 9 5 7 7 6 8 6 8 N。