Hosmer Lemeshow统计错误

Question

我想计算我的多变量逻辑模型的Hosmer Lemeshow统计量。我尝试了hosmerlem（下面的命令）和hosmer.test函数。

  

hosmerlem（y =我的二元因变量，yhat = fit（我的glm））

（逻辑模型的因变量是二进制的，我包含了大约10个预测变量。）

但在这两种情况下，我都会收到此错误： 在quantile.default中出错（yhat，probs = seq（0,1,1 / g））：   如果'na.rm'为FALSE

，则不允许缺少值和NaN

有人可以帮忙吗？ 提前致谢！

Answer 1

要显示 hoslem.test 的R代码，请输入：

ResourceSelection :: hoslem.test

注意，它包含了分位数函数：

qq <- unique(quantile(yhat, probs = seq(0, 1, 1/g)))

分位数函数会抛出错误，因为你的yhat向量包含缺失值。