如何将图例添加到多个条形图？

Question

我的代码：

#!/usr/bin/R 

par(mfrow=c(4,1)) 

fcm <-c(13.0,12.5,11.8) 
gk <-c(10.9 , 10.5 , 10.2) 
gg <-c(12.0 , 11.0 , 10.8) 
data1 <- rbind(fcm,gk,gg) 
colnames(data1) <- c(5,6,7) 

fcm <-c(1.29,1.25,1.22) 
gk <-c(1.19,1.05,1.00) 
gg <-c(1.10,1.02,1.07) 
data2 <- rbind(fcm,gk,gg) 
colnames(data2) <- c(5,6,7) 

fcm <-c(10.5,9.8,8.8) 
gk <-c(11.0,8.4,10.0) 
gg <-c(2.0,1.8,1.6) 
data3 <- rbind(fcm,gk,gg) 
colnames(data3) <- c(5,6,7) 

fcm <-c(0.33,0.29,0.31) 
gk <-c(0.49,0.40,0.35) 
data4 <- rbind(fcm,gk) 
colnames(data4) <- c(5,6,7) 

barplot(as.matrix(data1),ylim=c(0,20),main="P wave", 
     xlab="number of clusters", ylab="traveltime rms(ms)", 
     col=c("red", "black", "green"), beside=TRUE) 
barplot(as.matrix(data2),ylim=c(0,2),main="MT", 
     xlab="number of clusters", ylab="MT functions", 
     col=c("red", "black", "green"), beside=TRUE) 
barplot(as.matrix(data3),ylim=c(0,20),main="XBI", 
     xlab="number of clusters", ylab="index value", 
     col=c("red", "black", "green"), beside=TRUE) 
barplot(as.matrix(data4),ylim=c(0,0.6),main="NCE", 
     xlab="number of clusters", ylab="index value", 
     col=c("red", "black"), beside=TRUE) 

legend(legend = c("fcm","gk","gg"), fill = c("red", "black", "green")) 

dev.copy(png,"dp.png",width=9,height=6,units="in",res=200) 
dev.off() 

我已经尝试了传说中的右上方，但它重叠image.Anyone知道如何添加传说，位置也不是那么进口只是为了visible.I正在提交这对一个科学杂志我需要一个不错的布局。

Answer 1

您可以使用layout代替par(mfrow...)，然后在需要的位置添加图例。我的每个主要地块占据两行，而传说中的每一行都在底部。

本示例将图例水平放置在底部。不幸的是，我无法上传图片，因为我在防火墙后面。

layout(matrix(c(1, 1, 2, 2, 3, 3, 4, 4, 5), ncol=1)) 

# main plots 
par(mai=rep(0.5, 4)) 
fcm <-c(13.0,12.5,11.8) 
gk <-c(10.9 , 10.5 , 10.2) 
gg <-c(12.0 , 11.0 , 10.8) 
data1 <- rbind(fcm,gk,gg) 
colnames(data1) <- c(5,6,7) 

fcm <-c(1.29,1.25,1.22) 
gk <-c(1.19,1.05,1.00) 
gg <-c(1.10,1.02,1.07) 
data2 <- rbind(fcm,gk,gg) 
colnames(data2) <- c(5,6,7) 

fcm <-c(10.5,9.8,8.8) 
gk <-c(11.0,8.4,10.0) 
gg <-c(2.0,1.8,1.6) 
data3 <- rbind(fcm,gk,gg) 
colnames(data3) <- c(5,6,7) 

fcm <-c(0.33,0.29,0.31) 
gk <-c(0.49,0.40,0.35) 
data4 <- rbind(fcm,gk) 
colnames(data4) <- c(5,6,7) 

barplot(as.matrix(data1),ylim=c(0,20),main="P wave", 
     xlab="number of clusters", ylab="traveltime rms(ms)", 
     col=c("red", "black", "green"), beside=TRUE) 
barplot(as.matrix(data2),ylim=c(0,2),main="MT", 
     xlab="number of clusters", ylab="MT functions", 
     col=c("red", "black", "green"), beside=TRUE) 
barplot(as.matrix(data3),ylim=c(0,20),main="XBI", 
     xlab="number of clusters", ylab="index value", 
     col=c("red", "black", "green"), beside=TRUE) 
barplot(as.matrix(data4),ylim=c(0,0.6),main="NCE", 
     xlab="number of clusters", ylab="index value", 
     col=c("red", "black"), beside=TRUE) 

# legend for blank plot 
par(mai=c(0,0,0,0)) 
plot.new() 
legend(legend = c("fcm","gk","gg"), fill = c("red", "black", "green"), 
    "center", horiz=TRUE)