Tracé 4 intervalles de confiance avec barplot2 (gplots)

Question

Je souhaite comparer 4 intervalles de confiance. J'ai un exemple pour le faire avec 2 intervalles, mais je ne peux pas le faire fonctionner avec 4 côte à côte.

Voici le script de travail pour 2 intervalles. Le VMP0$LOAD est mes données.

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

Que pensez-vous de cela?

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) 

La prochaine fois, ajoutez vos données à votre question afin qu'elles puissent être reproduites.