我使用 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 获得您所寻求的信息，但我们确实知道足够多的信息，能够自行获取。
使用https://stackoverflow.com/a/16819127/1017276 的答案并将此知识与uniroot 结合起来，我们可以求解ecdf 返回的函数，将x 的值返回为y 所需的值。
uniroot(f = function(x) ecdf(data$mm)(v = x) - .996,
interval = range(data$mm, na.rm = TRUE))
uniroot 返回的列表中的root 元素将为您提供所需的值。
【讨论】：