具有多個重疊錯誤欄的條形圖

我正在嘗試創建具有多個錯誤欄的錯誤。事情是這樣的：http://flyordie.sin.khk.be/r/histogram%20error%20bars.PNG 我有以下數據集： http://flyordie.sin.khk.be/r/output.csv 具有多個重疊錯誤欄的條形圖

我一直在使用GGPLOT2和點陣圖形試過，但沒有發現任何東西，適合我的需要。

我目前用於顯示條形圖代碼是這樣的：

data <- read.csv("c:/output.csv") 
data 

par(las=3) 
barplot(data$PlateId, 
    height=data$HC.Maximum, 
    names.arg=data$PlateId, 
    col="lightblue")

，並顯示出最高errorbars我用這個代碼庫（GGPLOT2）

limits <- aes(ymax = qc$HC.Maximum, ymin = qc$HC.Minimum) 

p <- ggplot(qc, aes(colour=HC.Median,x=PlateId)) 
p + geom_bar(position="dodge")+ geom_errorbar(limits,position="dodge")

但我對如何不知道把它們放在同一個圖形上（如我的例子）

的數據：

qc <- structure(list(row = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = "Row", class = "factor"), 
    ID = 1:14, PlateId = c(35276L, 35279L, 35280L, 35281L, 35282L, 
    35290L, 35291L, 35292L, 35293L, 35294L, 35295L, 35296L, 35297L, 
    35298L), LC.Median = c(439688.495, 509376.055, 475218.99, 
    497368.215, 481801.9, 468603.43, 494713.175, 459047.385, 
    482819.47, 495162.31, 449592.51, 460564.95, 478715.915, 452293.465 
    ), LC.Stdev = c(52290.12229, 49648.49436, 55743.10306, 62002.53552, 
    46908.66149, 52489.615, 48016.94019, 52082.23899, 47934.37133, 
    58977.84845, 45827.62648, 53514.21095, 49638.98286, 139686.144 
    ), LC.Minimum = c(279610.16, 423651.45, 356422.31, 411639.77, 
    397362.84, 345178.07, 406073.72, 352834.86, 339035.77, 369554.11, 
    348688.39, 357341.56, 370463.11, 210367.91), LC.Maximum = c(498195.9, 
    630648.53, 614625.78, 686737.35, 621372.36, 576491.41, 579708.95, 
    580633.28, 580125.9, 622108.73, 530234.87, 563616.65, 614936.33, 
    730272.63), HC.Median = c(507356.465, 553226.525, 447067.77, 
    452223.76, 453439.37, 422491.755, 447438.8, 435034.635, 446148.105, 
    438089.69, 466748.63, 440005.81, 454927.74, 483599.71), HC.Stdev = c(65355.46121, 
    72762.07338, 80118.37641, 43653.99318, 73389.12355, 62590.47601, 
    46421.36678, 62822.88532, 61175.4241, 64418.56174, 63101.2232, 
    68166.51814, 61256.74139, 87354.9441), HC.Minimum = c(381552.05, 
    391124.94, 280614.72, 395454.12, 291433.84, 252579.15, 331661.03, 
    296223.64, 240262.37, 299431.98, 375224.27, 278780.87, 310275.66, 
    213170.04), HC.Maximum = c(626483.6, 635111.41, 555357.3, 
    528822.8, 534172.42, 514927.42, 538385.26, 533024.74, 524973.99, 
    544335.94, 564954.87, 572206.98, 547489.1, 565338.09), zPrime = c(-3.96, 
    -23.73, -7.88, -5.81, -5.32, -5.54, -4.48, -7.98, -6.99, 
    -5.63, -22.54, -33.83, -11.92, -17.44), Sb = c(1.17, 1.03, 
    0.91, 0.91, 0.89, 0.89, 0.9, 0.92, 0.92, 0.89, 1.04, 0.98, 
    0.95, 1.09), Sn = c(1.37, 0.3, -0.83, -0.76, -1.22, -1.01, 
    -1.08, -0.74, -0.86, -0.95, 0.31, -0.2, -0.52, 0.27)), .Names = c("row", 
"ID", "PlateId", "LC.Median", "LC.Stdev", "LC.Minimum", "LC.Maximum", 
"HC.Median", "HC.Stdev", "HC.Minimum", "HC.Maximum", "zPrime", 
"Sb", "Sn"), class = "data.frame", row.names = c(NA, -14L))

來源

2011-04-20 Sir Ksilem

公共服務公告：我用dput（）添加了數據（） – Andrie 2011-04-20 10:56:04

當創建具有多個層的地塊，我的方法，如下所示：

在初始調用定義共同美學ggplot
在每個附加層定義附加美學

請注意，我修改了您的代碼，因爲我無法讓您的示例正常工作：

提供一個明確的binwidth=1到geom_bar刪除警告
取出position=dodge，因爲這是默認和冗餘
供應的明確stat=identity到geom_bar

代碼：

ggplot(qc, aes(x=PlateId)) + 
    geom_bar(aes(y=HC.Median), binwidth=1, stat="identity", fill="cyan") + 
    geom_errorbar(aes(ymin=HC.Minimum, ymax=LC.Minimum), colour="red") + 
    geom_errorbar(aes(ymin=LC.Maximum, ymax=HC.Maximum), colour="purple")

enter image description here

來源

2011-04-20 10:57:37 Andrie

謝謝，這正是我一直在尋找的。 – 2011-04-20 11:41:05

只需添加另一個geom_errorbar。

#Rename limits to limits_hi 
limits_hi <- aes(ymax = qc$HC.Maximum, ymin = qc$HC.Minimum) 

#Define the other error bar 
limits_lo <- aes(ymax = qc$LC.Maximum, ymin = LC.Minimum) 

#I'm not quite sure what you want in the bars; see if this looks right 
p <- ggplot(qc, aes(factor(PlateId), HC.Median)) 

p + 
    geom_bar(position="dodge") + 
    geom_errorbar(limits_hi, position="dodge", colour = "red") + 
    geom_errorbar(limits_lo, position="dodge", colour = "blue") + 
    opts(axis.text.x = theme_text(angle = 30))

來源

2011-04-20 10:39:44

我看到了，我只是誤解了ggplot2的功能。非常感謝您的回答。 – 2011-04-20 11:40:42

具有多個重疊錯誤欄的條形圖

回答

相關問題