#install.packages("pheatmap")
rt1=read.table("lssoExp",sep="\t",header=T,row.names=1,check.names=F) #读取文件
rt1=as.data.frame(t(rt1[,3:ncol(rt1)]))
rt=log2(10^rt1+0.001)
outpdf="heatmap.pdf"
library(pheatmap)
Type=read.table("clusterCliGroup.Sig.txt",sep="\t",header=T,row.names=1,check.names=F)
Type=Type[order(Type$risk),]
colnames(rt)=gsub("\\.","-",colnames(rt))
rownames(Type)=gsub("\\.","-",rownames(Type))
rt=rt[,row.names(Type)]
pdf(outpdf,height=7.5,width=10)
pheatmap(rt, annotation=Type,
color = colorRampPalette(c("blue", "white", "red"))(50),
cluster_cols =F,
fontsize=8,
fontsize_row=8,
scale="row",
show_colnames=F,
fontsize_col=3)
dev.off()
临床相关性热图
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- 刘小泽写于19.9.6-第四单元第三讲:多个基因集相关性热图笔记目的:根据生信技能树的单细胞转录组课程探索smar...
- 在进行基因分析中,我们时常会对样本或基因之间的相关性进行分析,虽然R语言中的cor函数可以进行计算,但并没有提供合...