R语言回归中的Hosmer-Lemeshow以及calibration curve校正曲线

Hosmer-Lemeshow test

详细的Hosmer-Lemeshow test使用方法的介绍请参考：看这里

calibration curve是评价模型的一个重要指标，不信的话请看文章1，文章2，总之很重要的！！！

最近需要画校正曲线，但网上的相关资料实在是太少了，花了点时间研究了一下，跟大家分享一下如何使用Hosmer-Lemeshow test来绘制calibration curve，（p.s. 都是个人理解，可能有错，欢迎讨论）:

首先，按照方法介绍中的样例，我们生成一下数据：

set.seed(43657)
n <- 100
x <- rnorm(n)
xb <- x
pr <- exp(xb)/(1+exp(xb))
y <- 1*(runif(n) < pr)

这里的x是变量，y相当于标签（0，1），这里做的是二分类的。然后我们放到glm函数中fi