99999久久久久久亚洲,欧美人与禽猛交狂配,高清日韩av在线影院,一个人在线高清免费观看,啦啦啦在线视频免费观看www

熱線電話:13121318867

登錄
2019-03-05 閱讀量: 1452
R語言中如何畫熱力圖

基于ggplot的熱力圖繪制

library (ggplot2)
library (reshape2)
setwd("D:/example/") #設(shè)定數(shù)據(jù)的存放路徑
data1<-read.csv('a4.csv',header = TRUE) #讀取文件
data2<- melt (data1, id="Month") #轉(zhuǎn)化成長數(shù)據(jù)集
data2$variable= factor(data2$variable, levels=c('EL10_TMP','EL14_TMP','EL18_TMP','EL10_PRE','EL14_PRE','EL18_PRE','EL10_PDSI','EL14_PDSI','EL18_PDSI')) #設(shè)置順序
data2$Month= factor(data2$Month, levels=c('P.Jun','P.Jul','P.Aug','P.Sep','P.Oct','P.Nov','P.Dec','C.Jan','C.Feb','C.Mar','C.Apr','C.May','C.Jun','C.Jul','C.Aug','C.Sep','C.Oct')) ## 設(shè)置順序
p1<-ggplot(data2,aes(x=Month,y=variable,fill=value)) #熱力圖繪制
p2<-p1+geom_raster()+ scale_fill_gradient2(low="red", high="darkgreen", mid="white") #填充不同的顏色

基于pheatmap包的繪制

library(pheatmap)
library(RColorBrewer)
setwd("D:/example/")
data1<-read.csv('a4.csv',header = TRUE)
pheatmap(data1)
p3<-pheatmap(data1,
color=brewer.pal(11,"PRGn"),
cellwidth = 20,
cellheight = 10,
cluster_rows = FALSE,
cluster_cols = FALSE,
labels_row=c('P.Jun','P.Jul','P.Aug','P.Sep','P.Oct','P.Nov','P.Dec','C.Jan','C.Feb','C.Mar','C.Apr','C.May','C.Jun','C.Jul','C.Aug','C.Sep','C.Oct'),
fontsize = 12,
#display_numbers = TRUE,
number_color = "black",
)

7.1701
4
關(guān)注作者
收藏
評論(0)

發(fā)表評論

暫無數(shù)據(jù)
推薦帖子