r - 在绘图中注释/标记区域

Question

我有以下图表：

该图（红色字母除外）由以下代码生成：

boundary12 <- c(0.01, 0.005, 0.002, 0.001)
boundary23 <- c(0.1, 0.05, 0.02, 0.01)
boundary34 <- c(0.5, 0.6, 0.7, 0.55)

plot(boundary12, type="lines", ylim=c(0,1), xlab="Pollution level", ylab="Probability of event")
lines(boundary23)
lines(boundary34)

问题：

1）如果我生成很多图表，有没有办法半自动地把红色字母放进去？

2) 有没有办法给区域着色，例如区域 1 是蓝色，区域 2 是黄色，区域 3 是绿色？

标签： rplotpolygon

Answer 1

可以通过计算多边形质心来自动放置注释。

plot(boundary12, type="n", ylim=c(0,1), xlab="Pollution level",
  ylab="Probability of event")

poly12 <- data.frame(x=c(1:4, 4:1), y=c(boundary12, rev(boundary23)))
poly23 <- data.frame(x=c(1:4, 4:1), y=c(boundary23, rev(boundary34)))
poly34 <- data.frame(x=c(1:4, 4:1), y=c(boundary34, rep(1, 4)))

polygon(poly12, border="#00000066", col="skyblue")
polygon(poly23, border="#00000066", col="gold")
polygon(poly34, border="#00000066", col="darkolivegreen3")

library(maps)
centroid <- maps:::centroid.polygon

p12 <- centroid(poly12)
p23 <- centroid(poly23)
p34 <- centroid(poly34)

textcoor <- rbind(p12, p23, p34)
textcoor[1, 1] <- textcoor[1 ,1] * 0.8 # to smooch nr1 a bit to the left

text(textcoor, labels=1:3, col="darkblue", cex=1.2)