r - 从 cdplot() 中提取 x 和 y 值

Question

如何使用 cdplot(y~x) 提取任意 x 值下每个因子 y 的精确概率

>

谢谢

Answer 1

按照 ?cdplot 帮助文件中的示例，您可以执行以下操作...

## NASA space shuttle o-ring failures
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)

## CD plot
result <- cdplot(fail ~ temperature)

这是从 cdplot 输出中获取概率的简单方法。

# Getting the probabilities for each group.
lapply(split(temperature, fail), result[[1]])
$no
 [1] 0.8166854 0.8209055 0.8209055 0.8209055 0.8090438 0.7901473 0.7718317 0.7718317 0.7579343
[10] 0.7664731 0.8062898 0.8326761 0.8326761 0.8905854 0.9185472 0.9626185

$yes
[1] 3.656304e-05 6.273653e-03 1.910046e-02 6.007471e-01 7.718317e-01 7.718317e-01 8.062898e-01

请注意，结果是由cdplot不可见地返回的条件密度函数(在失败级别上累积)，因此我们可以拆分温度通过fail并使用lapply将返回的函数应用于这些值。