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 に答える