二項ロジスティック回帰モデルでいくつかの診断を実行しようとしています。具体的には、限界モデルがプロットされます。残念ながら、「有限の 'xlim' 値が必要です」というエラーが表示され続けます。以下のコードは問題を再現します。私のモデルには、数値変数とカテゴリ変数の両方が含まれています (モデル内でダミー変数に変換されます)。とにかく、すべての値が 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
これに対処する方法についてのアイデア/提案/ガイダンスを探しています。