R中的临床预测规则（校准、鉴别和验证）

这是我的glm模型，我能够创建，但我现在想做Hosmer Lemeshow GOF测试，但得到这个错误，我不理解-

Call:
glm(formula = BC.result ~ Diabetic + Low.diastolic + Pulse, family = 
"binomial", data = Data2)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.0805  -0.4559  -0.3144  -0.2437   2.8259  

 Coefficients:
              Estimate Std. Error z value Pr(>|z|)    
(Intercept)   -5.86311    0.98715  -5.939 2.86e-09 ***
Diabetic       1.21963    0.37395   3.262 0.001108 ** 
Low.diastolic  1.27095    0.35074   3.624 0.000291 ***
Pulse          0.02361    0.00780   3.027 0.002470 ** 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

Null deviance: 276.40  on 485  degrees of freedom
Residual deviance: 249.71  on 482  degrees of freedom
(14 observations deleted due to missingness)
AIC: 257.71
> hl<- hoslem.test(model3$BC.result, fitted(model3, g=10))
Error in model.frame.default(formula = cbind(y0 = 1 - y, y1 = y) ~ cutyhat) 
: 
variable lengths differ (found for 'cutyhat')

呼叫：
glm（公式=BC.result~糖尿病+低舒张压+脉搏，家族=
“二项式”，数据=数据2）
偏差残差：
最小1季度中值3季度最大值
-1.0805  -0.4559  -0.3144  -0.2437   2.8259  
系数：
估计标准误差z值Pr（>z）
（截距）-5.863110.98715-5.9392.86e-09***
糖尿病1.21963 0.37395 3.262 0.001108**
低舒张压1.27095 0.35074 3.624 0.000291***
脉冲0.02361 0.00780 3.027 0.002470**
---
签名。代码：0'***'0.001'***'0.01'*'0.05'.'0.1''1
（二项式族的离散参数取为1）
零偏差：485自由度上的276.40
剩余偏差：482自由度上的249.71
（由于缺失，删除了14项观察结果）
AIC:257.71
>老实说，如果您不熟悉编程，使用一行代码编写一个包会容易得多。尝试
启动
库。否则，学习R的基础知识将非常有用。您的问题似乎是，请尽量缩小范围，并用R语言编写a。我们很想提供帮助，但要求他人为您编写代码不是有效的方法；你能展示一下到目前为止你都做了些什么吗？你在谷歌上搜索过“hosmer lemeshow R”了吗？还是“c指数R”？这两个都给我带来了有用的结果…是的，我已经尝试了所有这些。将编辑包含我收到的所有错误消息的问题。感谢您的数据中缺少值，因此
fitted（。
比您的数据框短14个观察值。尝试
na.omit（）
？再一次，a将是有帮助的。。。