data <- data %>% subset(rowSums(data)>0) ##去掉有0的
data <- t(data)
group=c(rep("normal",24),rep("Drought",15))
data.pca <- prcomp(data, scale. = T)
library(factoextra)
fviz_pca_ind(data.pca, ## repel=T,自动调整文本位置
col.ind=group,
mean.point=F,
addEllipses = T,
legend.title="Groups",repel=T,
ellipse.type="confidence",
ellipse.level=0.9,
palette = c("#CC3333", "#339999","#44cef6","#00e500"))+
theme(panel.border = element_rect(fill=NA,color="black", size=1, linetype="solid"))