r - “lines()”表现得像“polygon()”R？

Question

我试图在红色回归线周围画两条黑线。但是该lines()命令绘制的东西更像是一条polygon()而不是一条简单的线（见下图代码）。

我想知道是否有一个解决方法可以简单地在回归线周围画两条线（即不确定区间），或者我错过了什么？

library(rstanarm) 
data(kidiq)
d <- kidiq  

fit <- stan_glm(kid_score ~ mom_iq,
           data = d,   
           prior = normal(0, 2.5),  
           prior_intercept = normal(0, 10),  
           prior_aux = cauchy(0, 100)) 

plot(kid_score ~ mom_iq, data = d, type = "n")
abline(fit, col = 2)

pred_lin <- posterior_linpred(fit)

loop <- length(d$mom_iq)
I <- matrix(NA, loop, 2)
for(i in 1:loop){
I[i,] = quantile(pred_lin[,i], c(.025, .975))
}
lines(d$mom_iq, I[,1], lty = 2)
lines(d$mom_iq, I[,2])

Answer 1

尝试订购的data.frame，如：

a <- cbind(d$mom_iq, I[,1])
a <- a[order(a[,1]),]
lines(a)

所以你也可以写：

lines(sort(d$mom_iq), I[,2][order(d$mom_iq)])

或者简单地说：

apply(I, 2, function(x) lines(sort(d$mom_iq), x[order(d$mom_iq)]))