带有双 y 图的 Geom_smooth 函数缺少置信区间

Question

我正在努力制作双 y 轴的图表。 loess 结果没有置信区间，我无法理解原因。

下面我报代码：

library(ggplot2)
library(readxl)

Dati  <- data.frame("r" = c(0.99, 1.42, 2.10, 3.32, 6.09), "Vix" = c(16500, 19200, 22500, 24000, 26000), "OT" = c(23.5, 19, 11, 9, 7), "ref" = c("PU 178", "PU 178", "PU 178", "PU 178", "PU 178"))
attach(Dati)


scaleFactor <- max(Vix) / max(OT)


Graph <- ggplot(Dati, aes(x= r)) +
  geom_point(aes(y= Vix, col=paste0("Vix ", ref)), shape = 1, size = 3.5) +
  geom_smooth(aes(y= Vix, col = paste0("Vix ", ref)),  method="loess", level=0.55, se = TRUE) +  
  geom_point(aes(y= OT * scaleFactor, col=paste0("OT ", ref)), shape = 1, size = 3.5) + 
  geom_smooth(aes(y=OT * scaleFactor, col = paste0("OT ", ref)), method="loess", level=0.55, se = TRUE) +  
  scale_color_manual(values=c('#644196', '#f92410', '#bba6d9',  '#fca49c'),
                     name = "") +

theme(legend.justification = "top") +

  scale_y_continuous(name="Viscosity at 10rpm (mPa s)", sec.axis=sec_axis(~./scaleFactor, name="open time (sec)")) +
  theme(
    axis.title.y.left=element_text(color='#f92410'),
    axis.text.y.left=element_text(color='#f92410'),
    axis.title.y.right=element_text(color='#644196'),
    axis.text.y.right=element_text(color='#644196'),
    legend.position = "none"
  ) + 
  scale_x_continuous(name="ratio A2333/AD5027") 

Graph

结果两行完全没有CI。我认为指定的 level 太大或太小，但也改变它我没有得到 CIs。我觉得5个值太少了，但是我在过去的图表中用5个值没有问题。

有人知道我是否犯了什么错误吗？

下面是我post我得到的图。

做

Answer 1

您的跨度太小（参见），因此估计置信区间的点数太少。例如，如果您这样做：

ggplot(Dati, aes(x= r)) +
geom_point(aes(y= Vix, col=paste0("Vix ",ref)),shape = 1, size = 3.5) +
geom_smooth(aes(y= Vix, col =paste0("Vix ",ref)),  method="loess" ,span=1) +  
geom_point(aes(y= OT * scaleFactor, col=paste0("OT ",ref)), shape = 1, size = 3.5) + 
geom_smooth(aes(y=OT * scaleFactor, col =paste0("OT ",ref) ), method="loess",span=1) +  
scale_color_manual(values=c('#644196', '#f92410', '#bba6d9',  '#fca49c'),
                     name = "") +

theme(legend.justification = "top")

Loess 在这里有点矫枉过正，您可以考虑其他 smooth 并且还可以长时间旋转数据以使其更易于编码：

library(tidyr)
library(dplyr)

Dati %>% mutate(OT = OT*scaleFactor) %>%
pivot_longer(-c(r,ref)) %>%
mutate(name = paste0(name,ref)) %>%
ggplot(aes(x = r,y = value,col = name,fill = name)) +
geom_point(shape = 1, size = 3.5) +
geom_smooth(method="gam",formula = y ~ s(x,k=3),alpha=0.1) +
theme_bw()

或2次多项式：

Dati %>% mutate(OT = OT*scaleFactor) %>%
pivot_longer(-c(r,ref)) %>%
mutate(name = paste0(name,ref)) %>%
ggplot(aes(x = r,y = value,col = name,fill = name)) +
geom_point(shape = 1, size = 3.5) +
geom_smooth(method="lm",formula = y ~ poly(x, 2),alpha=0.1) +
theme_bw()

带有双 y 图的 Geom_smooth 函数缺少置信区间

Missing Confidence Intervals on a Geom_smooth function with double y graph

r

ggplot2

confidence-interval