我正在研究一个具有序数响应的广义线性模型，并且我遇到了一些与 Hosmer Lemeshow 测试似乎有些不同的东西。我想知道是否有人可以解释为什么我会为同一模型获得不同的值。
首先，我在 R 中的 faraway 包中加载了债务数据集，并清理了数据以仅包含我想要在我的模型中的变量而没有缺失值。然后我对所需的变量运行 vglm 函数。问题是，我在第二次运行时以不同的顺序排列它们，当我对两个模型运行广义 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 ,data=debt, 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), g=10, ord=T) #Gives a test statistic of 19.229