I am trying to plot a calibration plot. There are 329443 observations, and I enter the following:
val.prob.ci.2(Bare$KFREPred, Bare$RRT_in_2YRS, smooth=&