從一維數據繪製置信區間？

Question

如果我的數據是逗號分隔值（持續時間）的一長排，則R可以通過

summary(dat) 

問題給我很多信息

能R還使一個陰謀有置信區間的數據？

Answer 1

作爲第一款谷歌比賽說明，你可以嘗試做如下：

n<-50 
x<-sample(40:70,n,rep=T) 
y<-.7*x+rnorm(n,sd=5) 
plot(x,y,xlim=c(20,90),ylim=c(0,80)) 
mylm<-lm(y~x) 
abline(mylm,col="red") 
newx<-seq(20,90) 
prd<-predict(mylm,newdata=data.frame(x=newx),interval = c("confidence"), 
level = 0.90,type="response") 
lines(newx,prd[,2],col="red",lty=2) 
lines(newx,prd[,3],col="red",lty=2) 

當然這是可能性之一 另一個例子是：

x <- rnorm(15) 
y <- x + rnorm(15) 
new <- data.frame(x = seq(-3, 3, 0.5)) 

pred.w.clim <- predict(lm(y ~ x), new, interval="confidence") 


# Just create a blank plot region with axes first. We'll add to this 
plot(range(new$x), range(pred.w.clim), type = "n", ann = FALSE) 


# For convenience 
CI.U <- pred.w.clim[, "upr"] 
CI.L <- pred.w.clim[, "lwr"] 

# Create a 'loop' around the x values. Add values to 'close' the loop 
X.Vec <- c(new$x, tail(new$x, 1), rev(new$x), new$x[1]) 

# Same for y values 
Y.Vec <- c(CI.L, tail(CI.U, 1), rev(CI.U), CI.L[1]) 

# Use polygon() to create the enclosed shading area 
# We are 'tracing' around the perimeter as created above 
polygon(X.Vec, Y.Vec, col = "grey", border = NA) 


# Use matlines() to plot the fitted line and CI's 
# Add after the polygon above so the lines are visible 
matlines(new$x, pred.w.clim, lty = c(1, 2, 2), type = "l", col = 
c("black", "red", "red")) 

Example 1

Example 2