为什么 vglm 中的变量顺序对 Hosmer Lemeshow 检验有影响？

问题描述

我正在研究具有序数响应的广义线性模型，但我遇到了与 Hosmer Lemeshow 测试似乎有些偏差的问题。我想知道是否有人可以解释为什么我为同一个模型得到不同的值。

首先，我在 R 中的 faraway 包中加载了债务数据集，并清理了数据以仅包含我想要在我的模型中使用的变量而不会丢失值。然后我在所需的变量上运行 vglm 函数。问题是，我第二次运行时将它们按不同的顺序排列，当我在两个模型上运行 Generalized Hosmer Lemeshow 测试时，我得到了不同的值。请注意，我没有在两个模型中使用不同的变量，当我找到两个模型的摘要时，我得到了相同的参数估计。

library(faraway)
data(debt)

debt <- data.frame(cbind(debt$incomegp,debt$agegp,debt$bankacc,debt$ccarduse))
colnames(debt) <- c("incomegp","agegp","bankacc","ccarduse")
debt <-na.omit(debt)

#Bank Account at the beginning of the predictor list
po_model1 <- vglm(as.ordered(ccarduse)~ bankacc + agegp + incomegp,data=debt,family=cumulative(parallel=TRUE))

#Bank Account at the end of the predictor list
po_model2 <- vglm(as.ordered(ccarduse)~ agegp + incomegp + bankacc,family=cumulative(parallel=TRUE))

library(generalhoslem)
logitgof(debt$ccarduse,fitted(po_model1),g=10,ord=T) #Gives a test statistic of 24.573
logitgof(debt$ccarduse,fitted(po_model2),ord=T) #Gives a test statistic of 19.229

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）