xyplot với khoảng tin cậy (ô và ô cửa sổ) được thêm vào trong r

Question

Tôi muốn thêm khoảng tin cậy 99% và 95% trong XYplot.

Sau đây là dữ liệu:

X <- 1:20 
Y <- c(0.5, 1.4, 2.3, 3.4, 4.5, 
     3.3, 3.0, 2.1, 1.5, 0, 
     0, 3.4, 4.5, 6.7, 5.3, 2.8, 
     0.5, 3.4, 3.5, 3.7) 
mydata <- data.frame (X, Y) 

Tôi muốn xác định giá trị lớn nhất của Y và giá trị X tương ứng là vị trí của trung bình trong hộp và cốt truyện râu ria. Bất cứ khi nào giá trị của Y giảm 1 điểm (trái hoặc phải) là khoảng tin cậy 99% (nằm bên trong hộp) và bất cứ khi nào Y giảm xuống còn 2 (cả trái và phải), các vị trí tương ứng trong x sẽ được biểu diễn bởi râu ria.

âm mưu mong muốn:

Giải thích.

tương ứng x giá trị max (Y) = 6,7 tương ứng giá trị x để hộp bên trái = 6,7 - 1, hộp đúng = 6,7-1 tương ứng giá trị x để râu ria trái = 6,7-2, Whisker đúng = 6,7 - 2

Answer 1

Nội dung sau đây sẽ giúp bạn bắt đầu. Nó sử dụng hình chữ nhật chứ không phải là một bwplot, và nó cũng không làm nội suy.

Tạo dữ liệu:

library(ggplot2) 

dat <- data.frame(
    x = 1:20, 
    y = c(0.5, 1.4, 2.3, 3.4, 4.5, 3.3, 3.0, 2.1, 1.5, 0, 0, 3.4, 4.5, 6.7, 5.3, 2.8, 0.5, 3.4, 3.5, 3.7) 
) 

Viết một hàm trả về 5 điểm mong muốn:

getRange <- function(x, a=1, b=2){ 
    maxy <- max(x) 
    xMax <- which.max(x) 
    x2 <- max(which(x[1:xMax] <= (maxy-a))) 
    x1 <- max(which(x[1:x2] <= (maxy-b))) 
    x3 <- xMax + min(which(x[-(1:xMax)] < (maxy+a))) 
    x4 <- x3 + min(which(x[-(1:x3)] < (maxy+b))) 
    data.frame(x1=x1, x2=x2, max=xMax, x3=x3, x4=x4) 
} 

Lấy giá trị phạm vi và cốt truyện:

rr <- getRange(dat$y, 1, 3) 

ggplot(dat, aes(x, y)) + geom_line() + geom_point() + 
    geom_rect(data=rr, aes(xmin=x2, xmax=x3, NULL, NULL), 
       ymin=-Inf, ymax=Inf, fill="blue", alpha=0.25) + 
    geom_rect(data=rr, aes(xmin=x1, xmax=x4, NULL, NULL), 
       ymin=-Inf, ymax=Inf, fill="blue", alpha=0.25)