[英]Confidence intervals for ggplot
如果要在同一圖上繪制兩組數據,是否有任何方法可以僅獲取其中一個數據集而不獲取另一個數據集的置信區間? 例如,回歸僅顯示對數據組A有意義,而對B則不顯着,但是您仍然希望在同一圖中直觀地描繪A和B的數據,而置信區間僅在重要組A周圍。
您可以有選擇地選擇要傳遞給回歸繪圖儀的數據。
考慮以下示例:
set.seed(10)
#Make sample data
df <- data.frame(
group=rep(c("A","B"), each=10),
X = rep(1:10, 2))
df$Y <- 2*df$X + runif(20, -20, 20) #Create y values with lots of noise
#Reduce the noise for group A
df[df$group == "A", "Y"] <- 2*df[df$group == "A", "X"] + rnorm(10)
#Compare regression p-values
coef(summary(lm(Y ~ X, data=df[df$group == "A", ])))[, 4] #p < 0.05 for group A
# (Intercept) X
#1.577943e-01 5.411004e-09
coef(summary(lm(Y ~ X, data=df[df$group == "B", ])))[, 4] #p > 0.05 for group B
#(Intercept) X
# 0.7338232 0.1309030
#Graph all points, coloring by group. Add a regression line for group A only.
ggplot(df, aes(x=X, y=Y, colour= group)) + theme_bw() +
geom_point(size=2.5) +
geom_smooth(data = df[df$group == "A",], method="lm")
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.