如何在可靠性图中创建箱柜

Question

我使用

创建了一个带有二项式响应变量的logistic / logit模型

model <- glm(response~predictor1+predictor2+...)

然后我使用predict函数创建一个新的数据框

outcome <-data.frame(predict(model,newdata=IndependentDataSet,type="response"),as.numeric(as.character(Independent$ResponseVariable)))

names(outcome) <- c("Pr","Obs")

我可以使用以下功能之一

plot(verify(data$obs,data$pr),CI=TRUE) 
attribute(verify(data$obs,data$pr)) 

创建一个看起来像这样的图

或

reliability.plot(verify(data$obs,data$pr))

来自

library(verification)

创建可靠性图表。我想知道如何根据具体值分离垃圾箱。例如，我正在评估的模型是基于19％（0.19）的气候学，我希望在（1/3）*气候学，气候学和气候学上升（2/3）对于进行中的垃圾箱。我怎么能这样做？

此外，我看到这些垃圾箱表示为圆圈，其大小与该垃圾箱中数据的百分比成正比。有谁知道如何在R中制作更美观的可靠性图表？欢迎任何建议。

这就是我希望我的图表出现的方式

Answer 1

最简单的方法是

    trace("attribute.default",edit=TRUE)

或其他任何功能。 这样，您可以访问源代码并进行编辑。这些更改仅影响当前的R会话。