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()

你可能感兴趣的:(PCA图)