没有循环。
下面给出了2个分布(y的值)之间的最大差异。
x.e <- ecdf(x)
y.e <- ecdf(y)
mx <- max(x,y)
max(abs(x.e(seq(0,mx,mx/1000)) - y.e(seq(0,mx,mx/1000))))
如何获得达到此最大值的位置(x的值)。
答案 0 :(得分:0)
@docendo discimus,@ danas.zuokas,thx求助。 :)
所以它是这样的:
x.e <- ecdf(x)
y.e <- ecdf(y)
mx <- max(x,y)
steps <- seq(0,mx,mx/1000)
dif <- abs(x.e(steps) - y.e(steps))
y.max <- max(dif)
x.max <- steps[which.max(dif)]