0

我想比较 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)
4

1 回答 1

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)

下次将您的数据添加到您的问题中，以便可以复制它。

于 2015-03-17T15:59:37.197 回答