¿No prestó atención? Señale con el dedo aquí -
1. Introducción
En la literatura, verá algunos diagramas de correlación de triángulos rotados , busque en Internet y los compartirá con usted.
2 dibujos
library(ggplot2)
library(grid)
# 加载内置数据集
data('mtcars')
# 计算相关性系数
corda <- data.frame(cor(mtcars))
# 上三角操作
corda[upper.tri(corda)] <- NA
# # 计算相关性系数
# corda <- data.frame(cor(mtcars))
#
# # 下三角操作
# corda[lower.tri(corda)] <- NA
# 加载R包
library(reshape2)
library(tidyverse)
# 增加行名列
corda$y <- rownames(corda)
# 宽数据转长数据
da <- melt(data = corda) %>% na.omit()
# 因子化
da$variable <- factor(da$variable,levels = unique(da$variable))
da$y <- factor(da$y,levels = unique(da$y))
# squre
p <- ggplot(da) +
# 矩形图层
geom_tile(aes(x = variable,y = y),fill = 'white',
show.legend = F,
color = 'grey80',size = 1) +
# 点图层
geom_point(aes(x = variable,y = y,fill = value,size = value),
show.legend = F,
shape = 21,color = 'black') +
theme_minimal() +
# 主题调整
theme(panel.grid = element_blank(),
aspect.ratio = 1,
axis.text = element_text(color = 'black',size = 20),
axis.text.x = element_text(angle = 45,hjust = 0),
axis.text.y = element_text(angle = 45,hjust = 1)) +
# x轴标签位置
scale_x_discrete(position = 'top') +
# 点颜色
scale_fill_gradientn(colors = colorRampPalette(c("#019267", "white", "red"))(10)) +
# 点大小范围
scale_size(range = c(7,14)) +
# 图例设置
guides(size = 'none',fill = guide_colorbar(title = 'Corr',
barwidth = 1.5,barheight = 15,
frame.colour = 'black',
ticks.colour = "black")) +
xlab('') + ylab('')
print(p, vp = viewport(width = unit(0.5, "npc"),
height = unit(0.5, "npc"), angle = -45))
rectángulo:
p <- ggplot(da) +
# 矩形图层
geom_tile(aes(x = variable,y = y,fill = value),
show.legend = F,
color = 'grey80',size = 1) +
# # 点图层
# geom_point(aes(x = variable,y = y,fill = value,size = value),
# show.legend = F,
# shape = 21,color = 'black') +
theme_minimal() +
# 主题调整
theme(panel.grid = element_blank(),
aspect.ratio = 1,
axis.text = element_text(color = 'black',size = 16),
axis.text.x = element_text(angle = 45,hjust = 0),
axis.text.y = element_text(angle = 45,hjust = 1)) +
# x轴标签位置
scale_x_discrete(position = 'top') +
# 点颜色
scale_fill_gradientn(colors = colorRampPalette(c("#398AB9", "white", "red"))(10)) +
# 点大小范围
scale_size(range = c(7,14)) +
xlab('') + ylab('')
print(p, vp = viewport(width = unit(0.8, "npc"),
height = unit(0.8, "npc"), angle = -45))
3 terminando
Se ve bien.
Bienvenido al grupo de intercambio Shengxin. Agrégueme a WeChat y también lo llevaré al chat grupal de WeChat . El código de datos se cargó al grupo QQ, bienvenido a unirse y descargar.老俊俊生信交流群
Código QR del grupo :
Lao Junjun WeChat:
Planeta del conocimiento:
Entonces, ¿estudiaste hoy?
Eso es todo por el intercambio de hoy, ¡estén atentos para el próximo!
Finalmente, bienvenidos a todos a compartir y reenviar , ¡sus Me gusta son mi aliento y afirmación !
Si crees que te ayuda mucho, ¡disfruta de una taza de agua feliz!
revisión anterior
◀ Rsamtools procesamiento por lotes de archivos bam
◀ GSEApy para análisis y visualización de enriquecimiento
◀ pysam lee archivos bam para preparar datos de control de calidad de Ribo-seq
◀ rueda de sankey para dibujar un hermoso diagrama de sankey
◀ ggplot dibuja diagrama de violín + diagrama de caja + prueba estadística
◀ Análisis ascendente de datos de Ribo-seq
◀ Reflexionando sobre el problema de la extracción de secuencias
◀ Aprenda a hacer gráficos con Cell Reports-CLIP-seq visualización de datos
◀ enriquecimiento de m6A en el centro del pico
◀ Distribución del metagen m6A Coordenadas corregidas
◀...