Ola,

Estou tentando construir um calibration plot como o da figura em anexo. Basicamente, ele une a curva de calibração com um gráfico de barras com a quantidade de observações por "bin". Alguem sabe uma maneira fácil de fazer isso?

Para o calibration plot eu ja construi com dois "packages" diferentes mas em nenhum consegui colocar as observações. Os códigos que usei foram:

require(gbm)

calibrate.plot(Database$Event, Database$Predicted)

#ou#

require(caret)

Plot<-calibration(as.factor(Event!=1)~Predicted, data=Database)

plot(Plot)

Desde já agradeço a atenção

Ary Serpa Neto, MD, MSc, PhD
Hospital Israelita Albert Einstein, São Paulo, Brazil (http://www.einstein.br): Department of Critical Care Medicine & Jewish Institute of Education and Research Albert Einstein (http://www.einstein.br/pesquisa)
The PROtective VEntilation NETwork (www.provenet.eu)
Associate Editor: BMC Anesthesiology and BMC Pulmonary Medicine
http://publicationslist.org/a.serpa_neto
http://lattes.cnpq.br/0254664157557375