PCA图

# choose_gene=head(rownames(nrDEG),1000) ## 50 maybe better
# choose_matrix=exprSet[choose_gene,]
choose_matrix<- exprSet
choose_matrix[1:4,1:4]
choose_matrix=t(scale(t(log2(choose_matrix+1)))) 


library(ggfortify)
df=as.data.frame(t(choose_matrix))
df$group=group_list
png('pca.png',res=120)
p=autoplot(prcomp( df[,1:(ncol(df)-1)] ), data=df,colour = 'group')+theme_bw()
print(p)
dev.off()
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容