marginalModelPlots 错误 "need finite 'xlim' values"

Question

我正在尝试运行对二元逻辑回归模型进行一些诊断。具体来说，边际模型图。不幸的是，我一直收到 "need finite 'xlim' values" 错误。下面的代码重现了这个问题。我的模型包括数字变量和分类变量（在模型中转换为虚拟变量）。无论如何，我知道当所有值都是 NA 时会发生此错误，但我的任何数据都不是这种情况，我不确定发生了什么。

set.seed(020275)
df <- data.frame(y=sample(c(0,1), 10, replace=TRUE), 
             cat=sample(c("Red", "Blue", "Green"), 10, replace=TRUE), 
             loc=sample(c("North", "South", "East", "West"), 10,  replace=TRUE), 
             count=runif(10, 0, 10),
             stringsAsFactors = FALSE)

glmModel <- glm(y ~ cat + loc + count, family=binomial(), data=df)
glmModel

library(car)
marginalModelPlots(glmModel)

我收到以下错误：

Error in plot.window(...) : need finite 'xlim' values
In addition: Warning messages:
1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion
2: In min(x) : no non-missing arguments to min; returning Inf
3: In max(x) : no non-missing arguments to max; returning -Inf

正在寻找 ideas/suggestions/guidance 如何处理这个问题。

Answer 1

看来 character 数据类型向量（上例中的 cat 和 loc）与 marginalModelPlots 不兼容，至少对于我当前使用的汽车包版本（ 2.1-1).我发现我可以使用 terms 参数将图限制为变量的子集，同时还包括线性预测图（如下所示）。

marginalModelPlots(glmModel, terms= ~ count)

marginalModelPlots 错误 "need finite 'xlim' values"

marginalModelPlots error "need finite 'xlim' values"

r

r-car