在R中，如何找到等于cdplot上某个y值的x值？反之亦然？

Question

我有一个cdplot，我试图找到我的x值，其中分布（或y值）= .5并且找不到可行的方法。另外，我希望在我的x值为0时找到y值，并希望找到该等式，如果它不同。

我无法提供我的代码，因为它依赖于具有大型数据帧的已保存工作区。我将以此为例：

fail <- factor(c(2, 2, 2, 2, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1,1, 2, 1, 1, 1, 1, 1),levels = 1:2, labels = c("no", "yes")) temperature <- c(53, 57, 58, 63, 66, 67, 67, 67, 68, 69, 70, 70,70, 70, 72, 73, 75, 75, 76, 76, 78, 79, 81) cdplot(fail ~ temperature)

所以我不需要快速而肮脏的方法来解决这个具体的例子，我需要一个可以应用到我自己的工作区的代码。

Answer 1

如果捕获cdplot的返回值，则会得到一个可用于查找这些值的函数。

CDP = cdplot(fail ~ temperature
uniroot(function(x) { CDP$no(x) - 0.5}, c(55,80))
> uniroot(function(x) { CDP$no(x) - 0.5}, c(55,80))
$root
[1] 62.34963
$f.root
[1] 3.330669e-16