Wie füge ich dem folgenden Plot eine Schlüssellegende hinzu?

Irgendwo in der oberen rechten Ecke möchte ich eine Schlüssellegende mit zwei kurzen horizontalen Farbbalken haben, in denen die rote "Plastische Chirurgie schief gelaufen" und die blaue "Deutschland" sagen soll.

Ich habe den folgenden Code verwendet, um den Plot zu erstellen:

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)

Alle Hinweise wäre sehr dankbar.