График 4 доверительных интервала с barplot2 (gplots)

Question

Я хочу сравнить 4 доверительных интервала. У меня есть пример для этого с двумя интервалами, но не могу заставить его работать с четырьмя бок о бок.

Вот рабочий сценарий на 2 интервала. VMP0$LOAD - мои данные.

CIA <- t.test(VMP0$NLoad)$conf.int 
CIB <- t.test(VMP1$NLoad)$conf.int 
#CIC <- t.test(VMP2$NLoad)$conf.int 
#CID <- t.test(VMP3$NLoad)$conf.int 

lower <- c(CIA[1], CIB[1]) 
upper <- c(CIA[2], CIB[2]) 
library(gplots) 
barplot2(c(mean(VMP0$NLoad), mean(VMP1$NLoad)), 
     plot.ci = TRUE, ci.l = lower,  
     ci.u = upper,col = 2:3) 

Answer 1

Как насчет этого?

library(gplots) 

# Just some example data 
load1 <- 1:20 
load2 <- 1:10 
load3 <- 10:20 
load4 <- 21:10 

# uncomment this to make it work for your case 
#load1 <- VMP0$NLoad 
#load2 <- VMP1$NLoad 
#load3 <- VMP2$NLoad 
#load4 <- VMP3$NLoad 

CIA <- t.test(load1)$conf.int 
CIB <- t.test(load2)$conf.int 
CIC <- t.test(load3)$conf.int 
CID <- t.test(load4)$conf.int 

lower <- c(CIA[1], CIB[1], CIC[1], CID[1]) 
upper <- c(CIA[2], CIB[2], CIC[2], CID[2]) 

barplot2(c(mean(load1), mean(load2), mean(load3), 
      mean(load4)), plot.ci = TRUE, 
     ci.l = lower, ci.u = upper,col = 2:3) 

В следующий раз добавьте свои данные на свой вопрос, чтобы их можно было воспроизвести.