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)
--------------------------------------------------
数据不好,但是有图有方法: