r - 在 r 中使用 LogisticDx 包中的 gof() 函数时出错

Question

您好，当尝试使用 LogisticDx 包中的 gof() 函数时，我不断收到以下错误消息。

因子错误(G，标签 = dx1[，格式(max(P)，数字 = 3)，by = G]$V1): 无效的“标签”；长度6应该是1或5

我无法找出导致此错误的原因。您可以找到导致错误的代码以及下面有效的代码。 如果您需要更多信息，请告诉我

H

R 版本和加载的库

R版本3.2.3 (2015-12-10) 平台:x86_64-pc-linux-gnu(64位) 运行环境:Ubuntu 15.04

附加基础包: [1] 统计图形 grDevices utils 数据集方法库

其他附加包: [1] ROCR_1.0-7 gplots_2.17.0 LogisticDx_0.2 xtable_1.8-2 pander_0.5.2
[6]plyr_1.8.3Amelia_1.7.4mice_2.25Rcpp_0.11.6knitr_1.11

产生错误的代码

PD <- c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
        1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
        1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
        1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
        0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
        0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
        0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
        0, 0, 0)
E <- c(4, 4, 3, 1, 0, 5, 1, 5, 5, 5, 1, 5, 4, 2, 1, 1, 1, 1, 5,
       0, 5, 5, 5, 5, 0, 0, 4, 5, 2, 4, 5, 4, 5, 3, 0, 5, 3, 3,
       5, 2, 4, 0, 0, 5, 1, 0, 3, 2, 5, 1, 2, 4, 0, 2, 5, 5, 4,
       3, 5, 2, 0, 2, 0, 1, 0, 5, 2, 0, 3, 0, 0, 0, 1, 0, 0, 2,
       3, 0, 0, 1, 0, 0, 0, 0, 0, 2, 0, 1, 3, 0, 2, 0, 3, 0, 0,
       0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 2, 0, 0, 2, 0, 0,
       0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 2, 1, 0, 0)

test.data <- data.frame(PD,E)
factor(test.data$PD)

g.error <- gof(glm(PD ~ E,family=binomial,data=test.data),plot=FALSE)

有效的代码

data(uis)
g.works <- gof(glm(RACE ~ NDRGTX,family=binomial,data=uis),plot=FALSE)

Answer 1

此处的问题与正在执行的 Hosmer–Lemeshow 测试中的组数有关。看来 gof 的默认参数是 g=10，这似乎不适用于 129 个观测值(如果再添加一行，它会起作用)。 如果您将组数设置为 9，它应该可以工作。

gof(glm(PD ~ E,family=binomial,data=test.data), plot=FALSE, g=9)

也许其他人可以回答为什么会出现这种情况？