R 如何获取ks图的数据

我用Gumbel拟合绘制了天气数据的Kolmogorov-Smirnov测试图

ks.gumbel(data$mm, 63.2899, 30.06080, alternative = "two.sided", plot = TRUE)

我想得到绘制在y=0.996的曲线的x值

---

有什么想法吗？谢谢

查看

ks.gumbel

的源代码，绘图由

plot（ecdf（x））

生成，但函数不返回绘图中的值。因此，直接从ks.gumbel获得您想要的信息是不可能的，但我们知道的足够多，可以自己获得

使用at的答案，并将此知识与

uniroot

相结合，我们可以解决

ecdf

返回的函数，以将

x

的值返回到所需的

y

值

uniroot(f = function(x) ecdf(data$mm)(v = x) - .996,
        interval = range(data$mm, na.rm = TRUE))

---

由

uniroot

返回的列表的

root

元素将为您提供所需的值。

请指定使用的软件包。我使用了reliaR软件包