web-dev-qa-db-ja.com

Rのマルチヒストグラムプロットに重要な凡例を追加する

以下のプロットに重要な凡例を追加するにはどうすればよいですか

enter image description here

右上隅のどこかに、2つの短い水平カラーバーが付いた重要な凡例が必要です。赤いバーは「形成外科がうまくいかなかった」、青いバーは「ドイツ」と表示されます。

次のコードを使用してプロットを作成しました。

bar2 <- read.table("div/ana-mut[...]/barriers-set-2.dat", sep=" ")
bar2val <- c(bar2$V1, bar2$V2)
bar3 <- read.table("div/ana-mut[...]/barriers-set-3.dat", sep=" ")
bar3val <- c(bar3$V1, bar3$V2)
p1 <- hist(subset(bar2val, bar2val < 30), breaks=30)
p2 <- hist(subset(bar3val, bar3val < 30), breaks=30)
plot(p1, col=rgb(1,0,0,8/9), main="Barrier distribution", xlab="Barrier [kcal/mol]", ylab="Mutant count")
plot(p2, col=rgb(0,0,1,8/9), add=T)

ヒントをいただければ幸いです。

rkeyhistogramlegend
9
TMOTTM

Legendコマンドはトリックを行います：

legend("topright", c("Germany", "Plastic"), col=c("blue", "red"), lwd=10)

2つの短い水平カラーバーを取得するには、標準の線を使用しますが、線の太さを増やします。 Rolandが指摘したように、fill引数を使用することもできます。

legend("topright", c("Germany", "Plastic"), fill=c("blue", "red"))

見る ?legend 詳細については。

enter image description here

18
csgillespie

これが代替ソリューションです（コードについては以下を参照してください）

enter image description here

# some semi-random data …
df <- structure(list(Germany = c(1L, 3L, 6L, 1L, 2L), Plastic = c(2L, 
5L, 4L, 2L, 3L)), .Names = c("Germany", "Plastic"), class = "data.frame", row.names = c(NA, 
-5L))

# Expand right side of clipping rect to make room for the legend
par(xpd=T, mar=par()$mar+c(0,0,0,4))

# Graph data (transposing the matrix) using heat colors,  
# put 10% of the space between each bar, and make labels  
# smaller with horizontal y-axis labels
barplot(t(df), main="Barrier distribution", xlab="Barrier [kcal/mol]", ylab="Mutant count", 
   col=c("blue", "red"), space=0.1, cex.axis=0.8, las=1,
   names.arg=c("Mon","Tue","Wed","Thu","Fri"), cex=0.8) 

# Place the legend at (4,9) using blue and red
legend(4, 9, names(df), lwd=4, col=c("blue", "red"))

# Restore default clipping rect
par(mar=c(5, 4, 4, 2) + 0.1)
2
Eric Fail