Question

我想比较4个置信区间。我有一个例子，它有2个间隔，但不能使它与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)

下次将您的数据添加到您的问题中，以便可以重现。

用barplot2（gplots）绘制4个置信区间

1 个答案: