散点拟合作图

数据格式:

代码:

library(openxlsx)
library(ggplot2)
library(ggpmisc)

dat = read.xlsx("./散点拟合数据.xlsx",
                sheet = 1,)

data_x_name = names(dat)[2]   # 获取数据列名
data_y_name = names(dat)[3]   # 获取数据列名

names(dat)[1] = 'SampleID'    # 重命名列名
names(dat)[2] = 'value_x'       # 重命名列名
names(dat)[3] = 'value_y'   # 重命名列名

# 画图
p = ggplot(data = dat,mapping=aes(x= value_x,
                                  y=value_y,))+
    
    geom_point(size = 2,
               alpha=1,
               color = 'black',
               shape = 16,
               stroke = 1,
               )+

    geom_smooth(method='lm',formula = y~x,colour='red',alpha=0.1)+
    stat_poly_eq(aes(label = paste(..eq.label..,..rr.label..,..p.value.label.., sep = "~~~")), 
                 formula = y ~ x, 
                 parse = TRUE, 
                 label.x = 0.1,  # 标签的位置坐标
                 label.y = 1,   # 标签的位置坐标
                 color = 'black', 
                 size = 5,
                 fontface = 'bold',
                 family = "serif",  # 设置字体新罗马字体
                 )+
    labs(x = data_x_name,  # 设置x轴标签
         y = data_y_name,  # 设置y轴标签
         ) +  # 设置图例
    theme_classic()  # 经典主题

p

数据:

sampleID diversity Multifunctionality group
TP1 0.939099 0.612299 1
TP2 0.875949 0.53845 1
TP3 0.824453 0.411662 1
TP4 0.601341 0.266782 2
TP5 0.700021 0.283706 2
TP6 0.733696 0.287109 2
本文最后更新于2025-01-23如有失效,请留言
THE END