r - 无法使用 ggplot2 和 ggpmisc 获得两位小数的 R2 值
问题描述
很抱歉,在这里寻求帮助时,我不知道如何制作可重现的示例。因此,我提供下面的代码和数据;<350 KB 这里(数据:https ://easyupload.io/1r5xuo ),如果你可以看看。
我的问题是:我想要两个 R2 的十进制值(如附图所示),即使我使用了 digits=2,我也无法使用 ggpmisc 的 stat_fit_glance() 函数得到它,但它只显示一个(R2=0.7和 R2=0.9),可能是因为小数点后 1 位的值为零(见图)。
但是,我希望它们为 R2=0.70 和 R2=0.90。
谁能帮我解决这个问题。
library(ggplot2)
library(ggpmisc)
library(dplyr)
library(openxlsx)
library(ggpubr)
# data
data_all <- read.xlsx("./DJFMA_new_files/djfma_hourly_clear_overcast_scatter.xlsx")
colnames(data_all)[1] <- "h"
colnames(data_all)[2] <- "le"
colnames(data_all)[3] <- "ta_ts"
colnames(data_all)[4] <- "q_qRS"
# filter
data_clear <- data_all %>%
filter(type == "Clear-sky")
data_over <- data_all %>%
filter(type == "Overcast")
# Formula
formula1 <- y ~ x
# Plot
ggplot(data=data_all, aes(x=ta_ts, y=h, color=type)) +
geom_hline(yintercept=0, linetype="dashed", color = "grey50", size=.3) +
geom_vline(xintercept=0, linetype="dashed", color = "grey50", size=.3) +
geom_point(alpha=0.3, size=1.5, shape=20) +
labs(x=NULL,
y=expression("H"~ "[W" ~ m^-2~"]")) +
scale_colour_manual("",values = c("Overcast"="#51806a","Clear-sky"="#DD3C51")) +
theme_bw() +
theme(legend.position = c(.25, .85),
legend.background = element_rect(fill = "transparent"),
axis.ticks.length=unit(-0.12, "cm"),
axis.text.y = element_text(margin=margin(5,7,5,5,"pt")),
axis.text.x = element_text(margin=margin(7,5,5,5,"pt"))) +
stat_fit_glance(method = 'lm',
method.args = list(formula = formula1),
#geom = 'text',
label.y = "bottom",
label.x = "right",
aes(label = paste("R^2 ==", round(..r.squared.., digits = 2),
sep = "~")),
parse=TRUE)
解决方案
使用stat_poly_eq()
它可以轻松调整位数。这是使用 CRAN 的当前版本的“ggpmisc”。我最近更新了代码,以便不会删除尾随零。
library(ggpmisc)
library(dplyr)
library(openxlsx)
# data
data_all <- read.xlsx("./R bits and pieces/djfma_hourly_clear_overcast_scatter.xlsx")
colnames(data_all)[1] <- "h"
colnames(data_all)[2] <- "le"
colnames(data_all)[3] <- "ta_ts"
colnames(data_all)[4] <- "q_qRS"
# filter
data_clear <- data_all %>%
filter(type == "Clear-sky")
data_over <- data_all %>%
filter(type == "Overcast")
# Plot
ggplot(data=data_all, aes(x=ta_ts, y=h, color=type)) +
geom_hline(yintercept=0, linetype="dashed", color = "grey50", size=.3) +
geom_vline(xintercept=0, linetype="dashed", color = "grey50", size=.3) +
geom_point(alpha=0.3, size=1.5, shape=20) +
labs(x=NULL,
y=expression("H"~ "[W" ~ m^-2~"]")) +
scale_colour_manual("",values = c("Overcast"="#51806a","Clear-sky"="#DD3C51")) +
theme_bw() +
theme(legend.position = c(.25, .85),
legend.background = element_rect(fill = "transparent"),
axis.ticks.length=unit(-0.12, "cm"),
axis.text.y = element_text(margin=margin(5,7,5,5,"pt")),
axis.text.x = element_text(margin=margin(7,5,5,5,"pt"))) +
stat_poly_eq(label.y = "bottom",
label.x = "right",
rr.digits = 2)
推荐阅读
- flutter - Start slide-in animation on button press flutter
- python - 检测/删除职位发布中的样板文本 [Python]?
- javascript - 获取数据并使用状态挂钩
- java - Selenium 测试返回与 chrome 和 egde 相同的错误 java.net.SocketException
- sql-server - AAD 和 DB 用户的混合身份验证
- vhdl - Quartus 网表优化在状态机中丢失寄存器扇出
- python - 如何根据几列中的指定条件选择数据框中的行?
- mysql - 根据不同列的值更新列的数据
- wordpress - 在 Wordpress 自定义帖子类型管理表单中添加动态更改的下拉列表
- python - 如何快速检查域是否存在?