Перекрывающийся доверительный интервал - PullRequest
0 голосов
/ 03 марта 2019

У меня есть три переменные, и у каждой есть доверительный интервал, как следует

variable lowerci beta upperci
a 0.86471152    0.87615515 0.88759878
b 0.966626566 0.970159119 0.973691671
c 0.94946211 0.955502661 0.961543211

Я хочу посмотреть, совпадают ли они (ясно, что нет, но я хочу график).

Как я могу сделать это в R?

Ответы [ 2 ]

0 голосов
/ 03 марта 2019

В дополнение к ответу @Rui Barradas, вот код, чтобы сделать это с базовой графикой:

df <- read.table(text = "variable lowerci beta upperci
a 0.86471152    0.87615515 0.88759878
b 0.966626566 0.970159119 0.973691671
c 0.94946211 0.955502661 0.961543211                  
", header = TRUE, as.is = TRUE)

plot(df[, 3], pch = 19, cex = 2, ylim = range(df[, -1]))
for (i in c(1:nrow(df))) {
    lines(c(i, i), c(df[i, 2], df[i, 4]))
}
0 голосов
/ 03 марта 2019

Если вам нужен только график доверительных интервалов, попробуйте следующее.

library(ggplot2)

ggplot(df1, aes(x = variable, y = beta, ymin = lowerci, ymax = upperci)) +
  geom_pointrange()

enter image description here

Данные.

df1 <- read.table(text = "
variable lowerci beta upperci
a 0.86471152    0.87615515 0.88759878
b 0.966626566 0.970159119 0.973691671
c 0.94946211 0.955502661 0.961543211                  
", header = TRUE)
...