记念今天画了一张PCA好图

2014年看一篇论文对PCA很着迷，恁是花了很长时间硬懂了大概。今天终于做出图来了。

-----------------------------------------------

View(X201706D)

rawdata <- X201706D

df <- rawdata[,c(3,4,5,6,8)]

pcadata<-prcomp(df)

score <- pcadata$x #generate scores of observations

newdata<- data.frame(rawdata,score) # merge score with the rawdata

attach(newdata)

library(ggplot2)

group<- Group

mydata<-data.frame(PC1,PC2,group,ID)

ggplot(data=mydata,aes(PC1,PC2,color=factor(group)))+

geom_point()+

geom_text(aes(label=mydata$ID))+ #add labels to plots

guides(colour="none")+

stat_ellipse(lwd=1.5)

--------------------------------------------------

数据不好，但是有图有方法：