PCA图

作者: 陈宇乔 | 来源:发表于2019-09-29 15:11 被阅读0次
# 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()

相关文章

网友评论

      本文标题:PCA图

      本文链接:https://www.haomeiwen.com/subject/jgmouctx.html