笔记—R语言做矩阵散点图
#install.packages("corrgram")
#install.packages("GGally")
install.packages("seriation")
install.packages("wesanderson")
library(corrgram)
library(ggplot2)
library(GGally)
library(wesanderson)#配色包
data <- read.table(file.choose(), header = T,sep = '\t')#以tab制表符分隔
new_data <- data[,-1]#去掉第一例
View(new_data)
# 得到变量之间的相关系数
#data_cor <- cor(new_data)
#cor()无法给出p值
#Hmisc包里的rcorr()函数能够同时给出相关系数以及显著性水平p-value。rcorr(x, type = c(“pearson”,“spearman”))。
#order = TRUE 时,相关矩阵将使用主成分分析法对变量重排序,这将使得二元变量的关系模式更为明显
#。order = FALSE 时,变量按原来的顺序输出。
# 第二种风格的图
# 添加字体
windowsFonts(TNM = windowsFont("Times New Roman"))
pearson <- cor(x= new_data[6:11],y = new_data[1:5],method = "pearson")#计算相关性系数
View(pearson)
ggpairs_theme <- theme_bw()+theme(panel.grid.major = element_blank(), #删除主要网格线
panel.grid.minor=element_blank(), #删除次要网格线
panel.border = element_rect(colour = "grey30", fill = NA, size=0.25),
#使用无填充、0.25磅黑色边框的方格
axis.title.x.bottom=element_text(size=5, face = "bold",family ="TNM",angle = 90),#axis.title 轴的标签
axis.text.x.bottom = element_text(size=6.5,angle = 90, family ="TNM",face = "bold",),#face字体("plain", "italic", "bold",加粗 "bold.italic"),
axis.text.y.left= element_text(size =6.3,family ="TNM",face = "bold",),
strip.background = element_blank())+
theme(strip.text=element_text(size = 7,family="TNM"),strip.text.x = element_text(size = 7,family="TNM",),strip.text.y = element_text(size = 7,family="TNM"))
#theme_bw()+白色背景加黑色网格线主题
#theme_gray()默认主题
p<- ggpairs(new_data,)
#Adiag <- wrap(ggally_densityDiag,size= 1, color ="#FF9999" )#指定颜色
#Alower <- wrap(ggally_points, color = "#99CCFF")
ggplot(as.data.frame(new_data))
p <-ggpairs(new_data,upper = list(continuous=wrap(ggally_cor,digits = 2,
size=3, family="TNM",face = "bold")), )+
ggpairs_theme
p
# lower = list(continuous = Alower),加在ggpairs里指定对角线和下三角颜色
#diag = list(continuous = Adiag),