因為常被問到，也好久沒來「上班」了，發一個教學文。 問題：怎麼畫出像下面這種圖？也就是二因子的柱狀圖。 http://imgur.com/riRkMGA 比較麻煩的部份是： 1. 同時有線條填色與顏色填色 2. 放在外部的圖示 3. 有 error bar 以下請見該圖的原始碼。 # dat 為輸入的二因子資料 set.seed(1234) dat <- data.frame( y = rnorm(80, 5, 2), x1 = gl(2, 40, labels = c("A","B")), x2 = gl(4, 10, labels = c("1","2","3","4")) ) # 利用 tapply 求出各組之 mean 和 sd y.mean <- with(dat, tapply(y, paste(x1, x2, sep="."), mean)) y.sd <- with(dat, tapply(y, paste(x1, x2, sep="."), sd)) # 為了方便之後畫圖，把 y.mean 和 y.sd 轉成矩陣 # 矩陣各列表示 X2，各欄表示 X1 Y <- matrix(y.mean, 4) E <- matrix(y.sd, 4) colnames(Y) <- colnames(E) <- c("A", "B") rownames(Y) <- rownames(E) <- c("1", "2", "3", "4") # 接下來就開始畫圖 # 過程小小複雜，若有不了解可以逐步進行並觀看圖怎麼變化 # 先開一個畫布 # 因為我想把圖示畫在右側外部，所以 mar 和 xpd 要自行設定 dev.new(width = 5, height = 5) par(mar=c(4, 4, 1, 10) + 0.1, xpd = T) # 接下來是畫上柱狀圖主體 # 因為 X2 == "4" 柱體要畫成同時有「斜線填色」和「顏色填色」 # 所以 barplot 要畫二次 # 以下是第一次，只畫主要文字、框線和顏色填色 # 我習慣先自行算好 ylim # 注意 barplot 的回傳值被我定義為變數 bp（之後有用） bp <- barplot( Y, beside = T, xaxt = "n", ylim = c(0,10), ylab = "Title for y-axis (dat$y)", xlab = "Title for x-axis (dat$x1)", col = c(0,3,0,"yellow") ) # 接下來是疊上 X2 == "4" 柱體的「斜線填色」 # 上一步已經畫過的就不要再畫一次 # 只要畫上 X2 == "4" 柱體的「斜線填色」即可 # 注意到 add = T 就是疊圖的作用 barplot( Y, beside = T, xaxt = "n", ylim = c(0,10), yaxt = "n", ylab = "", xlab = "", col = c(0,1,1,1), density = c(0,0,25,15), angle = c(NULL,NULL,45,135), border = c(0,0,0,0), add = T ) # 這裡只是示範怎麼用 axis() axis(1, colMeans(bp), c("A", "B"), tick = F) # 這裡是畫「只向上方延伸的 error bar」 # 這時候 bp 的值（也就是每個柱的橫軸位置）就很有用了 arrows(bp, Y, bp, Y+E, length = 0.1, angle = 90) # 當然也可以再補畫上「向下方延伸的 error bar」 # 但這個例子不太適當就不做了 # 有需要的話請再加入 # arrows(bp, Y, bp, Y-E, length = 0.1, angle = 90) # 接下來是畫圖示 # 這裡示範如何把圖示畫在右側外部 # 需要與前面 par() 配合才行 # 也因為 X2 == "4" 柱體的同時有「斜線填色」和「顏色填色」 # 所以 legend() 一樣要畫兩次 # 第一次負責文字、線框、顏色填色 legend( "right", legend = c("1", "2", "3", "4"), fill = c(0,3,0,"yellow"), title = "Row name (dat$x2)", text.width = 5.5, box.lwd = 1, col = c(1,1,1,1), inset=c(-0.9,0) ) # 第二次負責 X2 == "4" 的斜線填色就好 # 重覆的東西都不要再畫 # 注意兩次畫 legend 要完全重疊才會成功 # 因此我定了 text.width 之類看似多餘的參數 legend( "right", legend = c("", "", "", ""), fill = c(0,0,1,1), border = c(0,0,0,0), density = c(0,0,25,15), angle = c(NULL,NULL,45,135), title = "", text.width = 5.5, box.lwd = 0, col = c(0,0,0,0), inset=c(-0.9,0) ) # 畫個框框 box() [關鍵字]: 柱狀圖 線條填色 外部圖示 --

※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 125.230.76.169 ※ 文章網址: https://webptt.com/m.aspx?n=bbs/R_Language/M.1450200152.A.E08.html ※ 編輯: andrew43 (125.230.76.169), 12/16/2015 01:25:06 ※ 編輯: andrew43 (125.230.76.169), 12/16/2015 01:29:21