Как создать бункеры на диаграмме надежности

Question

Я создал модель логистики / логита с биномиальной переменной ответа, используя

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

а затем я использовал функцию прогнозирования, чтобы создать новый фрейм данных

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.