1

我想知道 R 是否能够仅打印（或突出显示）输入数据集的帕累托最优点。

例如，在下面的二维图中，您可以观察到一组 50 个点。我希望能够用不同的颜色打印帕累托最优点。在此示例中，考虑将两个维度最小化。

http://i42.tinypic.com/jso7ma.png

有小费吗？

编辑：

根据评论，我使用以下代码实现了预期的结果：

n <- 40
d <- data.frame(
  x = rnorm(n),
  y = rnorm(n)
)
# We want the "extreme" points in the following plot
par(mar=c(1,1,1,1))
plot(d, axes=FALSE, xlab="", ylab="")
for(i in 1:n) {
  polygon( c(-10,d$x[i],d$x[i],-10), c(-10,-10,d$y[i],d$y[i]), 
  col=rgb(.9,.9,.9,.2))
}


d <- d[ order(d$x, decreasing=FALSE), ]
result <- d[1,]
for(i in seq_len(nrow(d))[-1] ) {
  if( d$y[i] < result$y[nrow(result)] ) {
    result <- rbind(result, d[i,])  # inefficient
  } 
}
points(result, cex=2, pch=15)
4

1 回答 1

1

这是使用chull（一组点的凸包）的解决方案。警告：未经测试。

# assume d is your matrix or data frame of points
ch <- d[chull(d), ]
ch[apply(ch, 1, function(p) !any(ch[, 1] < p[1] & ch[, 2] < p[2])), ]

您可以通过替换chull为您自己的凸包算法（例如此处的其中一个）并将条件测试扩展为适合，将其扩展到 >2 维。

于 2013-06-17T16:37:06.827 回答