小编use*_*291的帖子

我想在R中绘制以下函数

f(w) = 1/(1-5*e^(-iw))

其中i是-1的平方根.R可以在绘图中处理复数吗？

所以我正在尝试计算帕累托前线(http://en.wikipedia.org/wiki/Pareto_efficiency)R并且能够做到,但是,我无法有效地做到这一点.特别是随着点对的数量增加,计算显着减慢.

所以一般来说,我想做的是检查所有非支配(或支配)对.现在,我这样做的方法是找到所有这样的点对,使得x _i > X 和y _i > Y,其中(x _i,y _i)是单对,X和Y表示所有点x和y.现在,这部分工作非常快并且易于实现,但是,还有可能多个x值可能相同,但它们将具有不同的y值,因此在这种情况下我希望能够识别x值具有最低y值(对于具有相同y值但x值不同的点,反之亦然).

为了说明这一点,这里是来自维基百科的图片:

所以基本上我希望能够识别出红线上的所有点.

这是我的代码可以正常工作,但对于大型数据集效率非常低:

#Example Data that actually runs quickly
x = runif(10000)
y = runif(10000)

pareto = 1:length(x)

for(i in 1:length(x)){
    cond1 = y[i]!=min(y[which(x==x[i])])
    cond2 = x[i]!=min(x[which(y==y[i])])
    for(n in 1:length(x)){
        if((x[i]>x[n]  &  y[i]>y[n]) | (x[i]==x[n] & cond1) | …