如何手动计算Hosmer-Lemeshow good of fit test(拟合优度检验)

想采用Hosmer-Lemeshow good of fit test(拟合优度检验)对两个模型进行对比

问题相关代码，请勿粘贴截图

pred<-predict(cox, newdata=test)
hoslem.test(as.numeric(test$complication),pred,g=10)

运行结果及报错内容

    Hosmer and Lemeshow goodness of fit (GOF) test

data: as.numeric(test$complication), pred
X-squared = -313.64, df = 8, p-value = 1

我的解答思路和尝试过的方法

以上得到的卡方值是负数，不知道为什么。我想理解Hosmer-Lemeshow good of fit test的原理，不采用代码，用spss进行Hosmer-Lemeshow good of fit test的卡方值和P值的计算，不知道可不可行。

我想要达到的结果

目前有实际观测数和模型预测数，想用spss计算两个模型的Hosmer-Lemeshow good of fit test的卡方值和P值。

你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答


本次提问扣除的有问必答次数，已经为您补发到账户，我们后续会持续优化，扩大我们的服务范围，为您带来更好地服务。