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

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

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))

чтобы создать график, который выглядит так чтобы создать график, который выглядит так

or

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

от

library(verification)

построить диаграмму надежности. Мне интересно, как я могу разделить бункеры на основе определенных значений. Например, модель, которую я оцениваю, основана на климатологии 19% (0,19), и я хочу, чтобы был интервал (1/3)*климатология, климатология и повышение на (2/3) климатологии. для исходящих бункеров. Как я могу это сделать?

Кроме того, я видел ячейки, представленные в виде кругов, размер которых пропорционален проценту данных, находящихся в этой ячейке. Кто-нибудь знает, как сделать более эстетичную диаграмму надежности в R? Любые рекомендации приветствуются.

Вот как я хочу, чтобы мои диаграммы выглядели Вот как я хочу, чтобы мои диаграммы выглядели


r bin reliability
person Taylor Adams    schedule 14.06.2017    source источник
comment
При обращении за помощью в построении графика обязательно включите воспроизводимый пример с примерными входными данными. Мы понятия не имеем, что вы видите. И расплывчатые вопросы о том, как сделать их более эстетичными, как бы не по теме. Задайте более конкретный вопрос.   -  person MrFlick    schedule 14.06.2017
comment
@MrFlick Я добавил дополнительную информацию к вопросу в ответ на ваш комментарий. Мой вопрос относительно значений бункера является конкретным и не должен вызывать никаких проблем с двусмысленностью. Любая помощь, которую я могу получить по этой теме, приветствуется. Последний вопрос, касающийся внешности, — это просто бонус, если кто-то знает ответ.   -  person Taylor Adams    schedule 14.06.2017

Ответы (1)


arrow_upward
0
arrow_downward

Самым простым может быть использование

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

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

person fcq    schedule 03.07.2019