r - 使用 NADA 库中的 plot(cenros) 命令在同一图表上显示多个概率图

Answer 1

在我看来，“ros”类对象的绘图方法记录不充分。ylim 和 xlim 都设置为参数的范围，这是以防止用户提供替代范围的方式完成的。您需要重新编写绘图方法才能更改此设置。NADA 中绘图方法的帮助页面上没有记录这一事实。

如果您想模仿用于设置 xlim 和 ylin 参数的方法，您需要查看 plot.ros 函数是如何做到的，然后推广到两个模型拟合：

    uncen = x$modeled[!x$censored]
    cen = x$modeled[x$censored]
    pp.uncen.nq = qnorm(x$pp[!x$censored])
    pp.cen.nq = qnorm(x$pp[x$censored])
     ymin = min(c(uncen, cen))
     ymax = max(c(uncen, cen))
     xmin = min(c(pp.uncen.nq, pp.cen.nq))
     xmax = max(c(pp.uncen.nq, pp.cen.nq))

因此，您需要保存模型拟合，然后提取modeled和 'pp' 参数并对其进行处理，直到您可以将某些内容传递给修改后的绘图方法。在尝试了 S4 调用方法但没有成功之后，我决定构建一个可以通过此调用看到的绘图方法的外部化版本会更容易：

showMethods("plot", class="ros", includeDefs=TRUE)