散点拟合作图
数据格式:
代码:
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